A new free weekly publication from APS

Physics - A new free weekly publication from APS
Please visit physics.aps.org
 

SELECTED ARTICLES (0)

Physical Review E

(Statistical, Nonlinear, and Soft Matter Physics)

February 2008

Volume 77, Number 2 , Articles (02xxxx)

Click on icons for further information Rapid Communication icon   Physics Viewpoint icon   Focus icon   Free to Read icon

Rapid Communication icon  Rapid Communication
     Close

     Close
Physics Viewpoint icon  Physics Viewpoint
Selected for a Viewpoint in Physics. Please visit http://physics.aps.org/.
     Close
Focus icon  Featured in Phys. Rev. Focus
Please visit http://focus.aps.org.
     Close
Free to Read icon  Free to Read
Marks articles whose full text is available without a subscription. See the FREE TO READ FAQ.
     Close

Part 1 - Statistical, Soft Matter, and Biological Physics


Part 2 - Nonlinear and Plasma Physics, Fluid Dynamics, and Related Topics

 
 
  • RAPID COMMUNICATIONS
  • ARTICLES
  • BRIEF REPORTS
  • COMMENTS

  • RAPID COMMUNICATIONS
  • ARTICLES
  • BRIEF REPORTS
  • COMMENTS
  • ERRATA
  • Part 1 - Statistical, Soft Matter, and Biological Physics

    RAPID COMMUNICATIONS

    Statistical physics

    Rapid

    Published 6 February 2008 (4 pages)
    020101(R)  Full Text: PDF (97 kB)  | Buy Article
    +
    Show Abstract
    We obtain an effective transport description for the superdiffusive motion of random walkers in stratified flow by projection of the process on the direction of stratification. The effective dimensionally reduced motion is shown to describe a correlated random walk characterized by the Lagrangian velocity correlation. We analyze the projected motion through exact analytical solutions for the distribution density for an arbitrary correlated Gaussian noise and derive an evolution equation for the one-point and conditional two-point displacement densities. The latter gives an explicit effective equation for superdiffusive transport in stratified random flow and demonstrates that the displacement density has a Gaussian scaling form for all times.
    Rapid

    Published 6 February 2008 (4 pages)
    020102(R)  Full Text: PDF (109 kB)  | Buy Article
    +
    Show Abstract
    We present a Brownian inchworm model of a self-propelled elastic dimer in the absence of an external potential. Nonequilibrium noise together with a stretch-dependent damping form the propulsion mechanism. Our model connects three key nonequilibrium features—position-velocity correlations, a nonzero mean internal force, and a drift velocity. Our analytical results, including striking current reversals, compare very well with numerical simulations. The model unifies the propulsion mechanisms of DNA helicases, polar rods on a vibrated surface, crawling keratocytes and Myosin VI. We suggest experimental realizations and tests of the model.
    Rapid

    Published 7 February 2008 (4 pages)
    020103(R)  Full Text: PDF (228 kB)  | Buy Article
    +
    Show Abstract
    We introduce a model for diffusion of two classes of particles (A and B) with priority: where both species are present in the same site the motion of A's takes precedence over that of B's. This describes realistic situations in wireless and communication networks. In regular lattices the diffusion of the two species is normal, but the B particles are significantly slower due to the presence of the A particles. From the fraction of sites where the B particles can move freely, which we compute analytically, we derive the diffusion coefficients of the two species. In heterogeneous networks the fraction of sites where B's are free decreases exponentially with the degree of the sites. This, coupled with accumulation of particles in high-degree nodes, leads to trapping of the low priority particles in scale-free networks.
    Rapid

    Published 25 February 2008 (4 pages)
    020104(R)  Full Text: PDF (111 kB)  | Buy Article
    +
    Show Abstract
    Stochastic resonance (SR) is theoretically investigated for a double quantum dot system represented by two discrete levels in respective wells. The system is driven by a periodic signal and a white noise source with variable amplitude, and thus displays an improved output signal-to-noise ratio, a characteristic signature of SR.

    Granular materials

    Rapid

    Published 20 February 2008 (4 pages)
    020301(R)  Full Text: PDF (292 kB)  | Buy Article
    +
    Show Abstract
    Two neutrally buoyant intruder particles in a granular bed fluidized by vertical, sinusoidal vibration are known to interact with each other over a range of about five intruder diameters. Using molecular dynamics simulations, we investigate in detail the spatial and temporal nature of this interaction. We show that the force of attraction between intruders can be calculated from the local density and kinetic energy using a simple equation of state. Moreover, the interaction can be changed from attractive to repulsive by reducing the coefficient of restitution between the intruders and host particles, one of the key results of this work.

    Colloidal dispersions, suspensions, and aggregates

    Rapid

    Published 7 February 2008 (4 pages)
    020401(R)  Full Text: PDF (131 kB)  | Buy Article
    +
    Show Abstract
    The general form of the electrostatic potential around an arbitrarily charged colloid at a flat interface between a dielectric and a screening phase (such as air and water, respectively) is analyzed in terms of a multipole expansion. The leading term is isotropic in the interfacial plane and varies with d−3 where d is the in-plane distance from the colloid. The effective interaction potential between two arbitrarily charged colloids is likewise isotropic and proportional to d−3, thus generalizing the dipole-dipole repulsion first found for point charges at water interfaces. Anisotropic attractive interaction terms can arise only for higher powers dn with n>=4. The relevance of these findings for recent experiments is discussed.
    Rapid

    Published 20 February 2008 (4 pages)
    020402(R)  Full Text: PDF (232 kB)  | Buy Article
    +
    Show Abstract
    The aging of a charged colloidal system has been studied by small-angle x-ray scattering, in the exchanged momentum range Q=0.03–5  nm−1, and by dynamic light scattering, at different clay concentrations (Cw=0.6–2.8  %). The static structure factor S(Q) has been determined as a function of both aging time and concentration. This is the direct experimental evidence of the existence and evolution with aging time of two different arrested states in a single system simply obtained only by changing its volume fraction: an inhomogeneous state is reached at low concentrations, while a homogeneous one is found at high concentrations.

    Structured and complex fluids

    Rapid

    Published 5 February 2008 (4 pages)
    020501(R)  Full Text: PDF (793 kB)  | Buy Article
    +
    Show Abstract
    We report on a first-principles molecular-dynamics study of structural and dynamical heterogeneity in supercooled liquid silicon. We find that highly tetrahedral configurations are intermittently formed and that spatially heterogeneous dynamics is concurrently induced in the deeply supercooled state (1000  K). This heterogeneity is responsible for the anomalous structural relaxation characterized by the stretched-exponential function. The temporal structural fluctuation is found to give rise to the 1/f dependence in the corresponding power spectral density. In a moderately supercooled state (1600  K), the structural and dynamical heterogeneity is quite weak, in contrast to the deeply supercooled state. The applicability of the Stillinger-Weber potential to the deeply supercooled state is also discussed.

    Liquid crystals

    Rapid

    Published 12 February 2008 (4 pages)
    020701(R)  Full Text: PDF (190 kB)  | Buy Article
    +
    Show Abstract
    A bent-core mesogen containing a tetrathiafulvalene group is studied. The mesophase is an optically isotropic spongelike smectic phase without optical activity. The material responds to uv light in a peculiar way, showing photoconductivity and, in some circumstances, changing its mesomorphic properties. Under an electric field, large chiral domains are segregated from an initial racemic phase within a few minutes.
    Rapid

    Published 20 February 2008 (4 pages)
    020702(R)  Full Text: PDF (202 kB)  | Buy Article
    +
    Show Abstract
    A two-dimensional model within the Q-bar-bar-tensor description of liquid crystals is used to describe the inhomogeneous order reconstruction in a nematic cell driven by tony modulation in the anchoring conditions. Homogeneous and inhomogeneous reconstruction are contrasted: the former is defectless, the latter is defect mediated. While the transition thresholds are comparable in both cases and in good agreement with experimental data, the biaxial wall breaking is considerably slower in the inhomogeneous transition than in the homogeneous one. The shape of the signal given by the electric current flowing through the cell allows us to distinguish the actual path followed by the transition.
    Rapid

    Published 25 February 2008 (4 pages)
    020703(R)  Full Text: PDF (606 kB)  | Buy Article
    +
    Show Abstract
    We measure the dependence of the interparticle force F on the distance R between two colloidal particles with hyperbolic hedgehog defects in a nematic liquid crystal using optical tweezers. The particle-defect pair can be regarded as an elastic “dipole” in the electrostatic analogy. In a parallel configuration, where the dipole vectors are parallel with each other, F is attractive and proportional to R−4. However, F becomes repulsive at small R due to the existence of a defect between the particles. In an antiparallel configuration, where the particles directly face each other, F is repulsive over the whole range of R and proportional to R−3.6. In another antiparallel configuration, where two hyperbolic hedgehog defects directly face each other, F is proportional to R−3.6 and F at small R turns out to be attractive upon tilting the dipoles. Furthermore, we yield the force between particles connected by a stringlike defect called a bubblegum defect.

    Polymers

    Rapid

    Published 13 February 2008 (4 pages)
    020801(R)  Full Text: PDF (498 kB)  | Buy Article
    +
    Show Abstract
    A hexagonal phase was found during the synthesis of interpenetrating polymer networks composed of polystyrene (PS) and poly(methyl methacrylate) (PMMA). By using confocal microscopy, it was found that the regularity of this hexagonal phase further increases upon de-cross-linking of the PS networks in the matrix phase by irradiation with shorter uv wavelengths. We conclude that the cooperation between the cross-link-induced suppression of phase separation and the elastic repulsion between the dispersed PMMA-rich domains is responsible for the emergence of this hexagonal phase.
    Rapid

    Published 15 February 2008 (4 pages)
    020802(R)  Full Text: PDF (242 kB)  | Buy Article
    +
    Show Abstract
    Using Langevin simulations, we find that simple “generic” bead-and-spring homopolymer chains in a sufficiently bad solvent spontaneously develop helical order during the process of collapsing from an initially stretched conformation. The helix formation is initiated by the unstable modes of the straight chain, which drive the system towards a long-lived metastable transient state. The effect is most pronounced if hydrodynamic interactions are screened.

    Biological physics

    Rapid

    Published 1 February 2008 (4 pages)
    020901(R)  Full Text: PDF (240 kB)  | Buy Article
    +
    Show Abstract
    We predict that patterns with correlated surface density of atoms have statistically higher promiscuity (ability to bind stronger to an arbitrary pattern) as compared with noncorrelated patterns with the same average surface density. We suggest that this constitutes a generic design principle for highly connected proteins (hubs) in protein interaction networks. We develop an analytical theory for this effect. We show that our key predictions are generic and independent, qualitatively, on the specific form of the interatomic interaction potential, provided it has a finite range.

    ARTICLES

    Statistical physics

    Published 5 February 2008 (10 pages)
    021101  Full Text: PDF (1171 kB)  | Buy Article
    +
    Show Abstract
    Stochastic coherence (SC) and self-induced stochastic resonance (SISR) are two distinct mechanisms of noise-induced coherent motion. For interacting SC and SISR oscillators, we find that whether or not phase synchronization is achieved depends sensitively on the coupling strength and noise intensities. Specifically, in the case of weak coupling, individual oscillators are insensitive to each other, whereas in the case of strong coupling, one fixed oscillator with optimal coherence can be entrained to the other, adjustable oscillator (i.e., its noise intensity is tunable), achieving phase-locking synchronization, as long as the tunable noise intensity is not beyond a threshold; such synchronization is lost otherwise. For an array lattice of SISR oscillators, except for coupling-enhanced coherence similar to that found in the case of coupled SC oscillators, there is an optimal network topology degree (i.e., number of coupled nodes), such that coherence and synchronization are optimally achieved, implying that the system-size resonance found in an ensemble of noise-driven bistable systems can occur in coupled SISR oscillators.

    Published 6 February 2008 (12 pages)
    021102  Full Text: PDF (900 kB)  | Buy Article
    +
    Show Abstract
    Motivated by the existence of remarkably ordered cluster arrays of bacteria colonies growing in Petri dishes and related problems, we study the spontaneous emergence of clustering and patterns in a simple nonequilibrium system: the individual-based interacting Brownian bug model. We map this discrete model into a continuous Langevin equation which is the starting point for our extensive numerical analyses. For the two-dimensional case we report on the spontaneous generation of localized clusters of activity as well as a melting-freezing transition from a disordered or isotropic phase to an ordered one characterized by hexagonal patterns. We study in detail the analogies and differences with the well-established Kosterlitz-Thouless-Halperin-Nelson-Young theory of equilibrium melting, as well as with another competing theory. For that, we study translational and orientational correlations and perform a careful defect analysis. We find a nonstandard one-stage, defect-mediated transition whose nature is only partially elucidated.

    Published 6 February 2008 (6 pages)
    021103  Full Text: PDF (132 kB)  | Buy Article
    +
    Show Abstract
    We demonstrate control over heat flow in an N-terminal molecular junction. Using simple model Hamiltonians we show that the heat current through two terminals can be tuned by the temperature and coupling parameters of external gating reservoirs. We discuss two models: A fully harmonic system and a model incorporating anharmonic interactions. For both models the control reservoirs induce thermal fluctuations of the transition elements between molecular vibrational states. We find that a fully harmonic model does not show any controllability, while for an anharmonic system the conduction properties of the junction strongly depend on the parameters of the gates. Realizations of the model system within nanodevices and macromolecules are discussed.

    Published 6 February 2008 (8 pages)
    021104  Full Text: PDF (493 kB)  | Buy Article
    +
    Show Abstract
    We study the phase structure of a phantom tethered surface model shedding light on the internal degrees of freedom (IDOF), which correspond to the three-dimensional rodlike structure of the lipid molecules. The so-called tilt order is assumed as IDOF on the surface model. The model is defined by combining the conventional spherical surface model and the XY model, which describes not only the interaction between lipids but also the interaction between the lipids and the surface. The interaction strength between IDOF and the surface varies depending on the interaction strength between the variables of IDOF. We know that the model without IDOF undergoes a first-order transition of surface fluctuations and a first order collapsing transition. We observe in this paper that the order of the surface fluctuation transition changes from first order to second order and to higher order with increasing strength of the interaction between IDOF variables. On the contrary, the order of collapsing transition remains first order and is not influenced by the presence of IDOF.

    Published 6 February 2008 (9 pages)
    021105  Full Text: PDF (183 kB)  | Buy Article
    +
    Show Abstract
    This paper presents a statistical theory of stochastic evaporation and degradation processes in complex polymerlike ring and chain aggregates with multiple degrading bonds between the primary particles (monomers). The exact kinetic solution fully describing fragmentation processes is obtained for such aggregates with arbitrary number of primary particles (monomers) and bonds between them. The effects of additional interaction of multiple bonds with each other is shown to have a drastic impact on the predicted kinetic processes and time-dependent particle size distributions during aggregate degradation. Structural effects associated with different distributions of multiple bonds and bonding configurations in the aggregates are also investigated and shown to have a significant impact on typical fragmentation time and accumulation of fragmenting aggregates in intermediate modes. The developed theory and its results will be important for degradation of multistranded polymers, polymer networks, self-assembling structures, surface nanoclusters and nanotechnology, and formation and evolution of aerosol aggregates resulting from transport and industry emissions.

    Published 7 February 2008 (16 pages)
    021106  Full Text: PDF (382 kB)  | Buy Article
    +
    Show Abstract
    We investigate disordered one- and two-dimensional Heisenberg spin lattices across the transition from integrability to quantum chaos from both statistical many-body and quantum-information perspectives. Special emphasis is devoted to quantitatively exploring the interplay between eigenvector statistics, delocalization, and entanglement in the presence of nontrivial symmetries. The implication of the basis dependence of state delocalization indicators (such as the number of principal components) is addressed, and a measure of relative delocalization is proposed in order to robustly characterize the onset of chaos in the presence of disorder. Both standard multipartite and generalized entanglement are investigated in a wide parameter regime by using a family of spin- and fermion-purity measures, their dependence on delocalization and on energy spectrum statistics being examined. A distinctive correlation between entanglement, delocalization, and integrability is uncovered, which may be generic to systems described by the two-body random ensemble and may point to a new diagnostic tool for quantum chaos. Analytical estimates for typical entanglement of random pure states restricted to a proper subspace of the full Hilbert space are also established and compared with random matrix theory predictions.

    Published 7 February 2008 (9 pages)
    021107  Full Text: PDF (156 kB)  | Buy Article
    +
    Show Abstract
    The time correlation functions for a Gaussian wave-packet preparation of the dissipative harmonic oscillator evolving from three initial conditions for the heat bath are calculated and compared with each other for Ohmic heat baths. The three initial distributions for the bath are the factorized, partially factorized, and unfactorized distributions. Explicit analytical formulas are derived and then used to study the effect of the three initial distributions on the subsequent dynamics. We find that the transient behavior does not depend sensitively on the initial condition as long as the initial Gaussian wave function of the system is centered at the equilibrium point. Differences become noticeable as the center of the wave packet is significantly shifted from the equilibrium point. These observations justify to some extent the prevalent use of factorized initial conditions for studying real time quantum dynamics in dissipative systems. The total energy in the system is also calculated for the three initial states and its relation to features in the decay is pointed out.

    Published 8 February 2008 (5 pages)
    021108  Full Text: PDF (173 kB)  | Buy Article
    +
    Show Abstract
    We study many interacting Brownian particles under a tilted periodic potential. We numerically measure the linear response coefficient of the density field by applying a slowly varying potential transversal to the tilted direction. In equilibrium cases, the linear response coefficient is related to the intensity of density fluctuations in a universal manner, which is called a fluctuation-response relation. We then report numerical evidence that this relation holds even in nonequilibrium cases. This result suggests that Einstein's formula on density fluctuations can be extended to driven diffusive systems when the slowly varying potential is applied in a direction transversal to the driving force.

    Published 11 February 2008 (24 pages)
    021109  Full Text: PDF (334 kB)  | Buy Article
    +
    Show Abstract
    Many of the core concepts and (especially field-theoretic) tools of statistical mechanics have developed within the context of thermodynamic equilibrium, where state variables are all taken to be charges, meaning that their values are inherently preserved under reversal of the direction of time. A principle concern of nonequilibrium statistical mechanics is to understand the emergence and stability of currents, quantities whose values change sign under time reversal. Whereas the correspondence between classical charge-valued state variables and their underlying statistical or quantum ensembles is quite well understood, the study of currents away from equilibrium has been more fragmentary, with classical descriptions relying on the asymmetric auxiliary-field formalism of Martin, Siggia, and Rose (and often restricted to the Markovian assumption of Doi and Peliti), while quantum descriptions employ a symmetric two-field formalism introduced by Schwinger and further clarified by Keldysh. In this paper we demonstrate that for quantum ensembles in which superposition is not violated by very strong conditions of decoherence, there is a large natural generalization of the principles and tools of equilibrium, which not only admits but requires the introduction of current-valued state variables. For these systems, not only do Martin-Siggia-Rose (MSR) and Schwinger-Keldysh (SK) field methods both exist, in some cases they provide inequivalent classical and quantum descriptions of identical ensembles. With these systems for examples, we can both study the correspondence between classical and quantum descriptions of currents, and also clarify the nature of the mapping between the structurally homologous but interpretationally different MSR and SK formalisms.

    Published 11 February 2008 (4 pages)
    021110  Full Text: PDF (76 kB)  | Buy Article
    +
    Show Abstract
    The crucial role of ambient correlations in determining thermodynamic behavior is established. A class of entangled states of two macroscopic systems is constructed such that each component is in a state of thermal equilibrium at a given temperature, and when the two are allowed to interact heat can flow from the colder to the hotter system. A dilute gas model exhibiting this behavior is presented. This reversal of the thermodynamic arrow is a consequence of the entanglement between the two systems, a condition that is opposite to molecular chaos and shown to be unlikely in a low-entropy environment. By contrast, the second law is established by proving Clausius' inequality in a low-entropy environment. These general results strongly support the expectation, first expressed by Boltzmann and subsequently elaborated by others, that the second law is an emergent phenomenon which requires a low-entropy cosmological environment, one that can effectively function as an ideal information sink.

    Published 11 February 2008 (9 pages)
    021111  Full Text: PDF (162 kB)  | Buy Article
    +
    Show Abstract
    We have introduced a set of coupled fractional reaction-diffusion equations to model a multispecies system undergoing anomalous subdiffusion with linear reaction dynamics. The model equations are derived from a mesoscopic continuous time random walk formulation of anomalously diffusing species with linear mean field reaction kinetics. The effect of reactions is manifest in reaction modified spatiotemporal diffusion operators as well as in additive mean field reaction terms. One consequence of the nonseparability of reaction and subdiffusion terms is that the governing evolution equation for the concentration of one particular species may include both reactive and diffusive contributions from other species. The general solution is derived for the multispecies system and some particular special cases involving both irreversible and reversible reaction dynamics are analyzed in detail. We have carried out Monte Carlo simulations corresponding to these special cases and we find excellent agreement with theory.

    Published 13 February 2008 (9 pages)
    021112  Full Text: PDF (172 kB)  | Buy Article
    +
    Show Abstract
    The dynamics of a system formed by a finite number N of globally coupled bistable oscillators and driven by external forces is studied focusing on a global variable defined as the arithmetic mean of each oscillator variable. Several models based on truncation schemes of a hierarchy of stochastic equations for a set of fluctuating cumulant variables are presented. This hierarchy is derived using Itô stochastic calculus, and the noise terms in it are treated using an asymptotic approximation valid for large N. In addition, a simplified one-variable model based on an effective potential is also considered. These models are tested in the framework of the phenomenon of stochastic resonance. In turn, they are used to explain in simple terms the very large gains recently observed in these finite systems.

    Published 14 February 2008 (6 pages)
    021113  Full Text: PDF (380 kB)  | Buy Article
    +
    Show Abstract
    We investigated the nonequilibrium phase transition of the conserved lattice gas model in one dimension using two update methods: i.e., the sequential update and the parallel update. We measured the critical indices of theta, beta, nu||, and nu[perpendicular] and found that, for a parallel update, the exponents were delicately influenced by the hopping rule of active particles. When the hopping rule was designed to be symmetric, the results were found to be consistent with those of the sequential update. The exponents we obtained were precisely the same as the corresponding results of a recently presented lattice model of two species of particles with a conserved field in one dimension, in contrast with the authors' claim. We also found that one of the scaling relations known for absorbing phase transition is violated.

    Published 14 February 2008 (5 pages)
    021114  Full Text: PDF (251 kB)  | Buy Article
    +
    Show Abstract
    Stability of directed transport and molecular separation in a symmetric channel is analyzed. The original mechanism is based on harmonic spatial oscillations of the channel, under which the system exhibits multiple regimes of a directed transport. The particles may be forced to move with different velocities and directions as the amplitude and/or frequency of the oscillations are adjusted to a proper resonance. The advantage of this mechanism in contrast to the ratchet systems is that the average particle velocity is larger than the velocity of the growing of the width of the particle spatial distribution. We have studied the stability of the directed transport with regard to random impacts to the channel parameters and oscillation frequency. Here we present the results of the simulations which show that the ability of the combined longitudinally and transversally vibrating randomized dynamic channel to perform directed molecular transport remains resilient to quite intensive random channel structure fluctuations (50–60  %) and relatively strong random impacts to its oscillations (15–20  %).

    Published 15 February 2008 (9 pages)
    021115  Full Text: PDF (232 kB)  | Buy Article
    +
    Show Abstract
    A two-dimensional square-lattice model for the formation of secondary structures in proteins, the hydrogen-bonding model, is extended to include the effects of solvent quality. This is achieved by allowing configuration-dependent nearest-neighbor interactions. The phase diagram is presented and found to have a much richer variety of phases than either the pure hydrogen-bonding self-avoiding walk model or the standard Theta-point model.

    Published 15 February 2008 (26 pages)
    021116  Full Text: PDF (1011 kB)  | Buy Article
    +
    Show Abstract
    By introducing three well-defined dimensionless numbers, we establish the link between the scale dilatation method able to estimate master (i.e., unique) singular behaviors of the one-component fluid subclass and the universal crossover functions recently estimated [Garrabos and Bervillier, Phys. Rev. E 74, 021113 (2006)] from the bounded results of the massive renormalization scheme applied to the Phid<sup>4</sup>(n) model of scalar order parameter (n=1) and three dimensions (d=3), representative of the Ising-like universality class. The master (i.e., rescaled) crossover functions are then able to fit the singular behaviors of any one-component fluid without adjustable parameter, using only one critical energy scale factor, one critical length scale factor, and two dimensionless asymptotic scale factors, which characterize the fluid critical interaction cell at its liquid-gas critical point. An additional adjustable parameter accounts for quantum effects in light fluids at the critical temperature. The effective extension of the thermal field range along the critical isochore where the master crossover functions seems to be valid corresponds to a correlation length greater than three times the effective range of the microscopic short-range molecular interaction.

    Published 15 February 2008 (5 pages)
    021117  Full Text: PDF (284 kB)  | Buy Article
    +
    Show Abstract
    Continuous-time quantum walks with long-range steps Rgamma (R being the distance between sites) on a discrete line behave in similar ways for all gamma>=2. This is in contrast to classical random walks, which for gamma>3 belong to a different universality class than for gamma<=3. We show that the average probabilities to be at the initial site after time t as well as the mean square displacements are of the same functional form for quantum walks with gamma=2, 4, and with nearest neighbor steps. We interpolate this result to arbitrary gamma>=2.

    Published 15 February 2008 (4 pages)
    021118  Full Text: PDF (210 kB)  | Buy Article
    +
    Show Abstract
    A correlated random walk approach to diffusion is applied to the disordered nonoverlapping Lorentz gas. By invoking the Lu-Torquato theory for chord-length distributions in random media [J. Chem. Phys. 98, 6472 (1993)], an analytic expression for the diffusion constant in arbitrary number of dimensions d is obtained. The result corresponds to an Enskog-like correction to the Boltzmann prediction, being exact in the dilute limit, and better or nearly exact in comparison to renormalized kinetic theory predictions for all allowed densities in d=2,3. Extensive numerical simulations were also performed to elucidate the role of the approximations involved.

    Published 21 February 2008 (5 pages)
    021119  Full Text: PDF (200 kB)  | Buy Article
    +
    Show Abstract
    Cooperative behaviors near the disorder-induced critical point in a random-field Ising model are numerically investigated by analyzing time-dependent magnetization in ordering processes from a special initial condition. We find that the intensity of fluctuations of time-dependent magnetization, chi(t), attains a maximum value at a time t=tau in a normal phase and that chi(tau) and tau exhibit divergences near the disorder-induced critical point. Furthermore, spin configurations around the time tau are characterized by a length scale, which also exhibits a divergence near the critical point. We estimate the critical exponents that characterize these power-law divergences by using a finite-size scaling method.

    Published 25 February 2008 (6 pages)
    021120  Full Text: PDF (90 kB)  | Buy Article
    +
    Show Abstract
    Grand canonical and canonical ensembles become equivalent in the thermodynamic limit, but when the system size is finite the results obtained in the two ensembles deviate from each other. In many important cases, the canonical ensemble provides an appropriate physical description but it is often much easier to perform the calculations in the corresponding grand canonical ensemble. We present a method to compute averages in the canonical ensemble based on calculations of the expectation values in the grand canonical ensemble. The number of particles, which is fixed in the canonical ensemble, is not necessarily the same as the average number of particles in the grand canonical ensemble.

    Published 25 February 2008 (7 pages)
    021121  Full Text: PDF (725 kB)  | Buy Article
    +
    Show Abstract
    In this paper, site percolation on random Phi3 planar graphs is studied by Monte Carlo numerical techniques. The method consists in randomly removing a fraction q=1−p of vertices from graphs generated by Monte Carlo simulations, where p is the occupation probability. The resulting graphs are made of clusters of occupied sites. By measuring several properties of their distribution, it is shown that percolation occurs for an occupation probability above a percolation threshold pc=0.7360(5). Moreover, critical exponents are compatible with those analytically known for bond percolation.

    Published 25 February 2008 (7 pages)
    021122  Full Text: PDF (345 kB)  | Buy Article
    +
    Show Abstract
    We present a numerical method for the Monte Carlo simulation of uncoupled continuous-time random walks with a Lévy alpha-stable distribution of jumps in space and a Mittag-Leffler distribution of waiting times, and apply it to the stochastic solution of the Cauchy problem for a partial differential equation with fractional derivatives both in space and in time. The one-parameter Mittag-Leffler function is the natural survival probability leading to time-fractional diffusion equations. Transformation methods for Mittag-Leffler random variables were found later than the well-known transformation method by Chambers, Mallows, and Stuck for Lévy alpha-stable random variables and so far have not received as much attention; nor have they been used together with the latter in spite of their mathematical relationship due to the geometric stability of the Mittag-Leffler distribution. Combining the two methods, we obtain an accurate approximation of space- and time-fractional diffusion processes almost as easy and fast to compute as for standard diffusion processes.

    Published 26 February 2008 (4 pages)
    021123  Full Text: PDF (91 kB)  | Buy Article
    +
    Show Abstract
    We show that steady-state work fluctuations in periodically modulated systems display universal features, which are not described by the standard fluctuation theorems. Modulated systems often have coexisting stable periodic states. We find that work fluctuations sharply increase near a kinetic phase transition where the state populations are close to each other. We also show that the work variance displays scaling with the distance to a bifurcation point where a stable state disappears and find the critical exponent for a saddle-node bifurcation.

    Published 26 February 2008 (6 pages)
    021124  Full Text: PDF (185 kB)  | Buy Article
    +
    Show Abstract
    An algorithm based on backward induction is devised in order to compute the optimal sequence of games to be played in Parrondo games. The algorithm can be used to find the optimal sequence for any finite number of turns or in the steady state, showing that ABABB… is the sequence with the highest steady state average gain. The algorithm can also be generalized to find the optimal adaptive strategy in a multiplayer version of the games, where a finite number of players may choose, at every turn, the game the whole ensemble should play.

    Published 26 February 2008 (8 pages)
    021125  Full Text: PDF (139 kB)  | Buy Article
    +
    Show Abstract
    We study the ground state energy EG(n) of N classical O(n) vector spins with the Hamiltonian [script H]=−[summation]i>jJijS-vectori·S-vectorj where the coupling constants {Jij} are arbitrary. We prove that EG(n) is independent of n for all n>nmax(N)=[left floor](sqrt(8N + 1)−1)/2[right floor]. We show that this bound is the best possible. We also derive an upper bound for EG(m) in terms of EG(n), for m<n. We obtain an upper bound on the frustration in the system, as measured by F(n)[equivalent][[summation]i>j|Jij|+EG(n)]/[summation]i>j|Jij|. We describe a procedure for constructing a set of Jij's such that an arbitrary given state, {S-vectori}, is the ground state. We show that the problem of finding the ground state for the special case n=N is equivalent to finding the ground state of a corresponding soft-spin problem.

    Published 27 February 2008 (19 pages)
    021126  Full Text: PDF (438 kB)  | Buy Article
    +
    Show Abstract
    We consider the three-dimensional randomly diluted Ising model and study the critical behavior of the static and dynamic spin-spin correlation functions (static and dynamic structure factors) at the paramagnetic-ferromagnetic transition in the high-temperature phase. We consider a purely relaxational dynamics without conservation laws, the so-called model A. We present Monte Carlo simulations and perturbative field-theoretical calculations. While the critical behavior of the static structure factor is quite similar to that occurring in pure Ising systems, the dynamic structure factor shows a substantially different critical behavior. In particular, the dynamic correlation function shows a large-time decay rate which is momentum independent. This effect is not related to the presence of the Griffiths tail, which is expected to be irrelevant in the critical limit, but rather to the breaking of translational invariance, which occurs for any sample and which, at the critical point, is not recovered even after the disorder average.

    Published 27 February 2008 (8 pages)
    021127  Full Text: PDF (553 kB)  | Buy Article
    +
    Show Abstract
    In this paper we propose a generalization of the one-dimensional outflow dynamics (KD). The rule is introduced as a simplification of Galam dynamics (GD) proposed in an earlier paper. We simulate three types of outflow dynamics, GD, Stauffer et al. dynamics, and KD, both on the square and triangular lattices and show whether the outflow dynamics is sensitive to the lattice structure. Moreover, we took into account several types of initial configuration—random, “stripes,” and “circle.” We investigate the dependence between the mean relaxation time and the initial density p of up-spins for each type of initial conditions, as well as dependence between the mean relaxation time and the size of the system. As a result, we show differences and similarities between three types of the outflow dynamics.

    Published 27 February 2008 (6 pages)
    021128  Full Text: PDF (187 kB)  | Buy Article
    +
    Show Abstract
    We calculate analytically the work distribution of a quantum harmonic oscillator with arbitrary time-dependent angular frequency. We provide detailed expressions for the work probability density for adiabatic and nonadiabatic processes, in the limits of low and high temperature. We further verify the validity of the quantum Jarzynski equality.

    Published 27 February 2008 (7 pages)
    021129  Full Text: PDF (201 kB)  | Buy Article
    +
    Show Abstract
    The transport of a dimer, consisting of two Brownian particles bounded by a harmonic potential, moving on a periodic substrate is investigated both numerically and analytically. The mobility and diffusion of the dimer center of mass present distinct properties when compared with those of a monomer under the same transport conditions. Both the average current and the diffusion coefficient are found to be complicated nonmonotonic functions of the driving force. The influence of dimer equilibrium length, coupling strength, and damping constant on the dimer transport properties are also examined in detail.

    Published 27 February 2008 (4 pages)
    021130  Full Text: PDF (159 kB)  | Buy Article
    See Also: Publisher's Note
    +
    Show Abstract
    It is shown analytically that the flashing annihilation term of a Verhulst kinetic leads to the power-law distribution in the stationary state. For the frequency of switching slower than twice the free growth rate this provides the quasideterministic source of a Lévy noise at the macroscopic level.

    Published 28 February 2008 (9 pages)
    021131  Full Text: PDF (266 kB)  | Buy Article
    +
    Show Abstract
    We discuss the application of wavelet transforms to a critical interface model which is known to provide a good description of Barkhausen noise in soft ferromagnets. The two-dimensional version of the model (one-dimensional interface) is considered, mainly in the adiabatic limit of very slow driving. On length scales shorter than a crossover length (which grows with the strength of the surface tension), the effective interface roughness exponent zeta is ~=1.20, close to the expected value for the universality class of the quenched Edwards-Wilkinson model. We find that the waiting times between avalanches are fully uncorrelated, as the wavelet transform of their autocorrelations scales as white noise. Similarly, detrended size-size correlations give a white-noise wavelet transform. Consideration of finite driving rates, still deep within the intermittent regime, shows the wavelet transform of correlations scaling as 1/f1.5 for intermediate frequencies. This behavior is ascribed to intra-avalanche correlations.

    Published 28 February 2008 (16 pages)
    021132  Full Text: PDF (400 kB)  | Buy Article
    +
    Show Abstract
    We study the wetting transition and the directed polymer delocalization transition on diamond hierarchical lattices. These two phase transitions with frozen disorder correspond to the critical points of quadratic renormalizations of the partition function. (These exact renormalizations on diamond lattices can also be considered as approximate Migdal-Kadanoff renormalizations for hypercubic lattices.) In terms of the rescaled partition function z=Z/Ztyp, we find that the critical point corresponds to a fixed point distribution with a power-law tail Pc(z)~Phi(ln z)/z1+µ as z-->+[infinity] [up to some subleading logarithmic correction Phi(ln z)], so that all moments zn with n>µ diverge. For the wetting transition, the first moment diverges [overline z]=+[infinity] (case 0<µ<1), and the critical temperature is strictly below the annealed temperature Tc<Tann. For the directed polymer case, the second moment diverges [overline z[sup 2]]=+[infinity] (case 1<µ<2), and the critical temperature is strictly below the exactly known transition temperature T2 of the second moment. We then consider the correlation length exponent nu: the linearized renormalization around the fixed point distribution coincides with the transfer matrix describing a directed polymer on the Cayley tree, but the random weights determined by the fixed point distribution Pc(z) are broadly distributed. This induces some changes in the traveling wave solutions with respect to the usual case of more narrow distributions.

    Equilibrium and linear transport properties of fluids

    Published 7 February 2008 (4 pages)
    021201  Full Text: PDF (480 kB)  | Buy Article
    +
    Show Abstract
    Alder and Wainwright discovered the slow power decay ~td/2 (d is dimension) of the velocity autocorrelation function in moderately dense hard-sphere fluids using the event-driven molecular dynamics simulations. In the two-dimensional (2D) case, the diffusion coefficient derived using the time correlation expression in linear response theory shows logarithmic divergence, which is called the “2D long-time-tail problem.” We reexamined this problem to perform a large-scale, long-time simulation with 1×106 hard disks using a modern efficient algorithm and found that the decay of the long tail in moderately dense fluids is slightly faster than the power decay (~1/t). We also compared our numerical data with the prediction of the self-consistent mode-coupling theory in the long-time limit [~1/(tsqrt(ln t))].

    Published 8 February 2008 (4 pages)
    021202  Full Text: PDF (234 kB)  | Buy Article
    +
    Show Abstract
    We propose to reconsider the diffusion of atoms in the Knudsen regime in terms of a complex dynamical reflection process. By means of molecular dynamics simulations, we emphasize the asymptotic nature of the cosine law of reflection at the atomic scale, and carefully analyze the resulting strong correlations in the reflection events. A dynamical interpretation of the accommodation coefficient associated with the slip at the wall interface is also proposed. Finally, we show that the first two moments of the stochastic process of reflection depend nonuniformly on the incident angle.

    Published 29 February 2008 (3 pages)
    021203  Full Text: PDF (73 kB)  | Buy Article
    +
    Show Abstract
    The relationship between capillary pressure and saturation in a porous medium often exhibits a power-law dependence. The physical basis for this relation has been substantiated by assuming that capillary pressure is directly related to the pore radius. When the pore space of a medium exhibits fractal structure this approach results in a power-law relation with an exponent of 3−Dv, where Dv is the pore volume fractal dimension. However, larger values of the exponent than are realistically allowed by this result have long been known to occur. Using a thermodynamic formulation for equilibrium capillary pressure we show that the standard result is a special case of the more general exponent (3−Dv)/(3−Ds) where Ds is the surface fractal dimension of the pores. The analysis reduces to the standard result when Ds=2, indicating a Euclidean relationship between a pore's surface area and the volume it encloses, and allows for a larger value for the exponent than the standard result when Ds>2.

    Granular materials

    Published 5 February 2008 (10 pages)
    021301  Full Text: PDF (2147 kB)  | Buy Article
    +
    Show Abstract
    Pattern forming processes are abundant in nature. Here, we report on a particular pattern forming process. Upon withdrawal of fluid from a particle-fluid dispersion in a Hele-Shaw cell, the particles are shown to be left behind in intriguing mazelike patterns. The particles, initially being uniformly spread out in a disc, are slowly pulled inwards and together by capillary and pressure forces. Invading air forms branching fingers, whereas the particles are compiled into comparably narrow branches. These branches are connected in a treelike structure, taking the form of a maze. The characteristic length scale within the structure is found to decrease with the volume fraction of the particles and increase with the plate separation in the Hele-Shaw cell. We present a simulator designed to simulate this phenomenon, which reproduces qualitatively and quantitatively the experiments, as well as a theory that can predict the observed wavelengths.

    Published 7 February 2008 (13 pages)
    021302  Full Text: PDF (814 kB)  | Buy Article
    +
    Show Abstract
    We report on acoustic wave propagation in a regular array of nominally identical beads under isotropic static stress. The weak polydispersity of the beads makes the contact lattice random. Time-frequency analysis of the acoustic signal is performed and allows measurement of the full lattice dispersion relation. Comparison with the theoretical prediction for a perfect triangular lattice gives an indication of the level of randomness in the contact lattice. The results extend, in a consistent way, a previous study restricted to long wavelength propagation [B. Gilles and C. Coste, Phys. Rev. Lett. 90, 174302 (2003)]: The contact lattice is ordered by increasing the stress, and the smaller the wavelength, the higher the stress required to get regular lattice behavior. Measurements involving ballistic propagation of the coherent wave, whatever its frequency, evidence reversible lattice behavior under compression and/or decompression. Nevertheless, correlations of short wavelength incoherent waves are a sensitive probe of disorder, and allow us to exhibit a small irreversible evolution of the lattice.

    Published 12 February 2008 (6 pages)
    021303  Full Text: PDF (836 kB)  | Buy Article
    +
    Show Abstract
    Experiments and computational simulations are carried out to study the behavior of a granular column in a silo whose walls are able to vibrate horizontally. The column is brought to a steady fluidized state and it behaves similar to a hydrostatic system. We study the dynamics of the granular discharge through openings at the bottom of the silo in order to search for a Torricelli-like behavior. We show that the flow rate scales with the wall induced shear rate, and at high rates, the granular bed indeed discharges similar to a viscous fluid.

    Published 26 February 2008 (7 pages)
    021304  Full Text: PDF (341 kB)  | Buy Article
    +
    Show Abstract
    We generalize a recent study of random space-filling bearings to a more realistic situation, where the spacing offset varies randomly during the space-filling procedure, and show that it reproduces well the size distributions observed in recent studies of real fault gouges. In particular, we show that the fractal dimensions of random polydisperse bearings sweep predominantly the low range of values in the spectrum of fractal dimensions observed along real faults, which strengthen the evidence that polydisperse bearings may explain the occurrence of seismic gaps in nature. In addition, the influence of different distributions on the offset is studied and we find that a uniform distribution is the best choice for reproducing the size distribution of fault gouges.

    Published 27 February 2008 (10 pages)
    021305  Full Text: PDF (616 kB)  | Buy Article
    +
    Show Abstract
    It has been recently shown [I. Fouxon et al., Phys. Rev. E 75, 050301(R) (2007); I. Fouxon et al., Phys. Fluids 19, 093303 (2007)] that, in the framework of ideal granular hydrodynamics (IGHD), an initially smooth hydrodynamic flow of a granular gas can produce an infinite gas density in a finite time. Exact solutions that exhibit this property have been derived. Close to the singularity, the granular gas pressure is finite and almost constant. We report molecular dynamics (MD) simulations of a freely cooling gas of nearly elastically colliding hard disks, aimed at identifying the “attempted” density blowup regime. The initial conditions of the simulated flow mimic those of one particular solution of the IGHD equations that exhibits the density blowup. We measure the hydrodynamic fields in the MD simulations and compare them with predictions from the ideal theory. We find a remarkable quantitative agreement between the two over an extended time interval, proving the existence of the attempted blowup regime. As the attempted singularity is approached, the hydrodynamic fields, as observed in the MD simulations, deviate from the predictions of the ideal solution. To investigate the mechanism of breakdown of the ideal theory near the singularity, we extend the hydrodynamic theory by accounting separately for the gradient-dependent transport and for finite density corrections.

    Published 28 February 2008 (9 pages)
    021306  Full Text: PDF (242 kB)  | Buy Article
    +
    Show Abstract
    Granular media are reversible and elastic if the stress increments are small enough. An elastic stress-strain relation, employed previously to determine static stress distributions, in this paper is compared to experiments by Kuwano and Jardine [Geotechnique 52, 727 (2002)] on incremental stress-strain relations, and shown to yield satisfactory agreement. In addition, the yield condition is given a firmer footing.

    Published 28 February 2008 (19 pages)
    021307  Full Text: PDF (721 kB)  | Buy Article
    See Also: Erratum
    +
    Show Abstract
    We employ hydrodynamic equations to investigate nonstationary channel flows of freely cooling dilute gases of hard and smooth spheres with nearly elastic particle collisions. This work focuses on the regime where the sound travel time through the channel is much shorter than the characteristic cooling time of the gas. As a result, the gas pressure rapidly becomes almost homogeneous, while the typical Mach number of the flow drops well below unity. Eliminating the acoustic modes and employing Lagrangian coordinates, we reduce the hydrodynamic equations to a single nonlinear and nonlocal equation of a reaction-diffusion type. This equation describes a broad class of channel flows and, in particular, can follow the development of the clustering instability from a weakly perturbed homogeneous cooling state to strongly nonlinear states. If the heat diffusion is neglected, the reduced equation becomes exactly soluble, and the solution develops a finite-time density blowup. The blowup has the same local features at singularity as those exhibited by the recently found family of exact solutions of the full set of ideal hydrodynamic equations [I. Fouxon et al., Phys. Rev. E 75, 050301(R) (2007); Phys. Fluids 19, 093303 (2007)]. The heat diffusion, however, always becomes important near the attempted singularity. It arrests the density blowup and brings about previously unknown inhomogeneous cooling states (ICSs) of the gas, where the pressure continues to decay with time, while the density profile becomes time-independent. The ICSs represent exact solutions of the full set of granular hydrodynamic equations. Both the density profile of an ICS and the characteristic relaxation time toward it are determined by a single dimensionless parameter [script L] that describes the relative role of the inelastic energy loss and heat diffusion. At [script L]>>1 the intermediate cooling dynamics proceeds as a competition between “holes”: low-density regions of the gas. This competition resembles Ostwald ripening (only one hole survives at the end), and we report a particular regime where the “hole ripening” statistics exhibits a simple dynamic scaling behavior.

    Published 29 February 2008 (14 pages)
    021308  Full Text: PDF (1882 kB)  | Buy Article
    +
    Show Abstract
    Direct measurements of the acceleration of spheres and disks impacting granular media reveal simple power law scalings along with complex dynamics which bear the signatures of both fluid and solid behavior. The penetration depth scales linearly with impact velocity while the collision duration is constant for sufficiently large impact velocity. Both quantities exhibit power law dependence on sphere diameter and density, and gravitational acceleration. The acceleration during impact is characterized by two jumps: a rapid, velocity-dependent increase upon initial contact and a similarly sharp depth-dependent decrease as the impacting object comes to rest. Examination of the measured forces on the sphere in the vicinity of these features leads to an experimentally based granular force model for collision. We discuss our findings in the context of recently proposed phenomenological models that capture qualitative dynamical features of impact but fail both quantitatively and in their inability to capture significant acceleration fluctuations that occur during penetration and which depend on the impacted material.

    Published 29 February 2008 (8 pages)
    021309  Full Text: PDF (498 kB)  | Buy Article
    +
    Show Abstract
    We study the distribution of volume fluctuations in experiments and numerical simulations concerning equal-sized sphere packings prepared with different techniques. We show that the distribution of the local volumes (Voronoï cells) and also the distributions of the global volumes (whole samples) follow remarkably well a shifted and rescaled Gamma distribution that we name a k-Gamma distribution. Such agreement is robust over a broad range of packing fractions and it is observed for several distinct systems. This distribution is characterized by the average packing fraction and a shape parameter “k” which is very sensitive to changes in the structural organization. A statistical mechanics approach predicts such k-Gamma distribution at statistical equilibrium and it links the parameter k with the number of elementary cells which are exchanging volume during the system preparation. The thermodynamical equivalent of k and its relation with the “granular temperature” are also discussed.

    Colloidal dispersions, suspensions, and aggregates

    Published 5 February 2008 (7 pages)
    021401  Full Text: PDF (186 kB)  | Buy Article
    See Also: Publisher's Note
    +
    Show Abstract
    Nonequilibrium steady states in an open system connecting two reservoirs of platelike colloidal particles are investigated by means of a recently proposed phenomenological dynamic density functional theory [M. Bier and R. van Roij, Phys. Rev. E 76, 021405 (2007)]. The platelike colloidal particles are approximated within the Zwanzig model of restricted orientations, which exhibits an isotropic-nematic bulk phase transition. Inhomogeneities of the local chemical potential generate a diffusion current which relaxes to a nonvanishing value if the two reservoirs coupled to the system sustain different chemical potentials. The relaxation process of initial states towards the steady state turns out to comprise two regimes: a smoothening of initial steplike structures followed by an ultimate relaxation of the slowest diffusive mode. The position of a nonequilibrium interface and the particle current of steady states depend nontrivially on the structure of the reservoirs due to the coupling between translational and orientational degrees of freedom of the fluid.

    Published 12 February 2008 (9 pages)
    021402  Full Text: PDF (568 kB)  | Buy Article
    +
    Show Abstract
    Upon drying, colloidal suspensions undergo a phase transformation from a “liquid” to a “gel” state. With further solvent evaporation, tensile stresses develop in the gel, which ultimately leads to fractures. These generally manifest themselves in regular cracking patterns which reflect the physical conditions of the drying process. Here we show experimentally and theoretically how, in the case of a drying droplet of magnetic colloid (ferrofluid), an externally applied magnetic field modifies the stress in the gel and therefore the crack patterns. We find that the analysis of the shape of the cracks allows one to estimate the value of the gel Young's modulus just before the crack nucleation.

    Published 28 February 2008 (5 pages)
    021403  Full Text: PDF (195 kB)  | Buy Article
    +
    Show Abstract
    We study the optical absorbance from a magnetically arranged colloidal structure, and investigate the possibility of creating a magnetically controlled optical sensor using this system. The colloids form chains when exposed to an external magnetic field, which tend to collapse and form a more random particle arrangement when the field is removed. We show that a small magnetic field is able to change the sensor's reflection coefficient by more than 30%, and investigate in detail the relaxation mechanism when the field is turned off.

    Structured and complex fluids

    Published 4 February 2008 (4 pages)
    021501  Full Text: PDF (531 kB)  | Buy Article
    +
    Show Abstract
    We study the aging behavior of aqueous laponite suspension, a model soft glassy material, in creep. We observe that the viscoelastic behavior is time dependent and is strongly influenced by the deformation field; the effect is known to arise due to aging and rejuvenation. We show that irrespective of the strength of the deformation field (shear stress) and age, when the imposed time scale is normalized with a dominating relaxation mode of the system, universal aging behavior is obtained, demonstrating time-stress superposition, a phenomenon that may be generic in a variety of soft materials.

    Published 13 February 2008 (14 pages)
    021502  Full Text: PDF (305 kB)  | Buy Article
    +
    Show Abstract
    The shear-transformation-zone (STZ) theory of plastic deformation in glass-forming materials is reformulated in light of recent progress in understanding the roles played by the effective disorder temperature and entropy flow in nonequilibrium situations. A distinction between fast and slow internal-state variables reduces the theory to just two coupled equations of motion, one describing the plastic response to applied stresses and the other the dynamics of the effective temperature. The analysis leading to these equations contains, as a by-product, a fundamental reinterpretation of the dynamic yield stress in amorphous materials. In order to put all these concepts together in a realistic context, I conclude with a reexamination of the experimentally observed rheological behavior of a bulk metallic glass. That reexamination serves as a test of the STZ dynamics, confirming that system parameters obtained from steady-state properties such as the viscosity can be used to predict transient behaviors.

    Published 14 February 2008 (7 pages)
    021503  Full Text: PDF (702 kB)  | Buy Article
    +
    Show Abstract
    Glassy dynamics was investigated for thin films of atactic polystyrene by complex electric capacitance measurements using dielectric relaxation spectroscopy. During the isothermal aging process the real part of the electric capacitance increased with time, whereas the imaginary part decreased with time. It follows that the aging time dependences of real and imaginary parts of the electric capacitance were primarily associated with change in volume (film thickness) and dielectric permittivity, respectively. Further, dielectric permittivity showed memory and rejuvenation effects in a similar manner to those observed for poly(methyl methacrylate) thin films. On the other hand, volume did not show a strong rejuvenation effect.

    Published 21 February 2008 (17 pages)
    021504  Full Text: PDF (1513 kB)  | Buy Article
    +
    Show Abstract
    We study numerically phase separation in a binary fluid subject to an applied shear flow in two dimensions, with full hydrodynamics. To do so, we introduce a mixed finite-differencing and spectral simulation technique, with a transformation to render trivial the implementation of Lees-Edwards sheared periodic boundary conditions. For systems with inertia, we reproduce the nonequilibrium steady states reported in a recent lattice Boltzmann study. The domain coarsening that would occur in zero shear is arrested by the applied shear flow, which restores a finite-domain-size set by the inverse shear rate. For inertialess systems, in contrast, we find no evidence of nonequilibrium steady states free of finite-size effects: Coarsening persists indefinitely until the typical domain size attains the system size, as in zero shear. We present an analytical argument that supports this observation and that furthermore provides a possible explanation for a hitherto puzzling property of the nonequilibrium steady states with inertia.

    Published 29 February 2008 (7 pages)
    021505  Full Text: PDF (289 kB)  | Buy Article
    +
    Show Abstract
    The influence of a stationary shear flow on the crystallization in a glassy system is studied by means of molecular dynamics simulations and subsequent cluster analysis. The results reveal two opposite effects of the shear flow on the processes of topological ordering in the system. Shear promotes the formation of separated crystallites and suppresses the appearance of the large clusters. The shear-induced ordering proceeds in two stages, where the first stage is related mainly to the growth of crystallites and the second stage is due to an adjustment of the created clusters and a progressive alignment of their lattice directions. The influence of strain and shear rate on the crystallization is also investigated. In particular, we find two plausible phenomenological relations between the shear rate and the characteristic time scale needed for ordering of the amorphous system under shear.

    Films, interfaces, and crystal growth

    Published 14 February 2008 (8 pages)
    021601  Full Text: PDF (167 kB)  | Buy Article
    +
    Show Abstract
    As an alloy surface evolves under capillary forces, differing mobilities of the individual components can lead to kinetic alloy decomposition at the surface. In this paper, we address the relaxation of nanoscale sinusoidal ripples on alloy surfaces by considering the effects of both surface and bulk diffusion. In the absence of bulk diffusion, we derive exact analytical expressions for relaxation rates and identify two natural time scales that govern the relaxation dynamics. Bulk diffusion is shown to reduce kinetic surface segregation and enhance relaxation rates, owing to intermixing near the surface. Our results provide a quantitative framework for the interpretation of relaxation experiments on alloy surfaces.

    Published 27 February 2008 (5 pages)
    021602  Full Text: PDF (967 kB)  | Buy Article
    +
    Show Abstract
    We investigate crystalline order on a two-dimensional paraboloid of revolution by assembling a single layer of millimeter-sized soap bubbles on the surface of a rotating liquid, thus extending the classic work of Bragg and Nye on planar soap bubble rafts. Topological constraints require crystalline configurations to contain a certain minimum number of topological defects such as disclinations or grain boundary scars whose structure is analyzed as a function of the aspect ratio of the paraboloid. We find the defect structure to agree with theoretical predictions and propose a mechanism for scar nucleation in the presence of large Gaussian curvature.

    Published 28 February 2008 (6 pages)
    021603  Full Text: PDF (253 kB)  | Buy Article
    +
    Show Abstract
    We report measurements of the x-ray scattering intensity as mesoporous silica glasses are filled with neon. The intensity of the first peak in the liquidlike diffraction pattern increases nonlinearly with mass adsorbed. We outline a simple model assuming that the major coherent contribution to the first peak in the scattering function S(Q) is due to interference from nearest-neighbor scatterers. This allows us to demonstrate an approach for surface area determination which does not rely on thermodynamic models—and is therefore complementary to existing methods. We also suggest that the overestimation of surface area by the traditional Brunauer-Emmett-Teller method may be resolved by using the capillary, and not the bulk, condensation pressure as the reference pressure p0. Furthermore, the alternative analysis offers an insight into the atomic structure of monatomic sorption, which may be of use for further studies on materials with different surface properties.

    Liquid crystals

    Published 5 February 2008 (6 pages)
    021701  Full Text: PDF (289 kB)  | Buy Article
    +
    Show Abstract
    Using Monte Carlo simulations of the Lebwohl-Lasher model we study the director ordering in a nematic cell where the top and bottom surfaces are patterned with a lattice of ±1 point topological defects of lattice spacing a. As expected on general physical grounds we find that the nematic order depends on the ratio of the height of the cell H to a. For thick cells (H/a>~0.9) we find that the system is very well ordered and the frustration induced by the lattice of defects is relieved in a novel way by a network of half-integer defect lines which emerge from the point defects and hug the top and bottom surfaces of the cell. When H/a<~0.9 the system has zero nematic order parameter and the half-integer defect lines thread through the cell joining point defects on the top and bottom surfaces. We present a simple physical argument in terms of the length of the defect lines to explain these results. To facilitate eventual comparison with experimental systems we also simulate optical textures in the presence of crossed polarizers.

    Published 8 February 2008 (10 pages)
    021702  Full Text: PDF (1146 kB)  | Buy Article
    +
    Show Abstract
    The elastic response of a smectic-C elastomer to three deformations, namely imposed lambdaxx, lambdaxz, and lambdazz, has been modeled using a nonlinear theory of a nematic elastomer with embedded smectic layers, and with the director tilt (in the x direction) at a fixed angle with respect to the smectic layer normal (z direction). The main focus is the elastic response after any soft mode of the sample. It is found that the elastomer contracts in the x direction under lambdaxz shear. On stretching parallel to the layer normal it is found that there is a soft mode that acts to rotate the director toward the z direction. The deformation of the system after this soft mode can be reduced to shear and elongation in the plane of the layers. We make predictions of the mechanical response of the elastomer, in particular the length of the soft plateau and the asymptotic modulus for the elastomer when stretched parallel to the layer normal. Finally, since monodomain Sm-C elastomers are made by the deformation-induced alignment of polydomains, we describe these important systems. Qualitative behavior of the model is then compared to existing experimental literature on the mechanical alignment of polydomains

    Published 8 February 2008 (6 pages)
    021703  Full Text: PDF (331 kB)  | Buy Article
    +
    Show Abstract
    We report a generalized scenario for the formation of modulated smectic phases of bent-core molecules based on locally ferroelectric layering and spontaneous splay of the polarization. Twelve phases are proposed, distinguished by neighboring splay stripes with either syn- or antiorder of the polarization and undulation slope, in addition to layer continuity versus layer discontinuity at the intervening defects. We outline the experimental techniques necessary to differentiate among the phases and interpret previous results in the present context, using high resolution x-ray scattering diffraction and block and undulation models of the layer organization to distinguish among the three 2D lattice types which emerge.

    Published 8 February 2008 (8 pages)
    021704  Full Text: PDF (253 kB)  | Buy Article
    +
    Show Abstract
    Over the last 12 years, the possible existence of a tetrahedratic mesophase, involving a third-rank orientational order parameter and no positional order, has been addressed theoretically and predicted in some cases; no experimental realizations of a purely tetrahedratic phase are known at the time being, but various pieces of evidence suggest that interactions of tetrahedral symmetry do play a significant role in the macroscopic properties of mesophases resulting from banana-shaped (bent-core) mesogens. We address a very simple tetrahedratic mesogenic lattice model, involving continuous interactions; we consider particles possessing Td symmetry, whose centers of mass are associated with a three-dimensional simple-cubic lattice; the pair potential is taken to be isotropic in orientation space and restricted to nearest-neighboring sites; we let the two orthonormal triads {ualpha,  alpha=1,2,3} and {vgamma,  gamma=1,2,3} define the orientations of a pair of interacting particles; we let the unit vectors ualpha be combined to yield four unit vectors {ej,  j=1,2,3,4}, arranged in a tetrahedral fashion; we let the unit vectors vgamma be similarly combined to yield the four unit vectors {fk,  k=1,2,3,4}; and finally we let hjk=(ej·fk). The interaction model studied here is defined by the simplest nontrivial (cubic) polynomial in the scalar products hjk, consistent with the assumed symmetry and favoring orientational order; it is, so to speak, the tetrahedratic counterpart of the Lebwohl-Lasher model for uniaxial nematics. The model was investigated by molecular field (MF) theory and Monte Carlo simulations; MF theory predicts a low-temperature, tetrahedrically ordered phase, undergoing a second-order transition to the isotropic phase at higher temperature; on the other hand, available theoretical treatments point to the transition being driven first order by thermal fluctuations. Simulations showed evidence of a first-order transition.

    Published 13 February 2008 (11 pages)
    021705  Full Text: PDF (255 kB)  | Buy Article
    +
    Show Abstract
    For many years it has been commonly accepted that electroconvection (EC) as primary instability in nematic liquid crystals for the “classical” planar geometry requires a positive anisotropy of the electric conductivity, sigmaa, and a slightly negative dielectric anisotropy, epsilona. This firm belief was supported by many experimental and theoretical studies. Recent experiments, which have surprisingly revealed EC patterns at negative conduction anisotropy as well, have motivated the theoretical studies in this paper. It will be demonstrated that extending the common hydrodynamic description of nematics by the usually neglected flexoelectric effect allows for a simple explanation of EC in the “nonstandard” case sigmaa<0.

    Published 28 February 2008 (12 pages)
    021706  Full Text: PDF (600 kB)  | Buy Article
    +
    Show Abstract
    We present a high-resolution x-ray study of the effects of disorder induced by random cross-linking side-chain smectic elastomers. The influence of variation of the concentration and stiffness of the cross-link units on the disruption of the one-dimensional translational order is reported in detail. Precise analysis of the line shape of the quasi-Bragg peaks associated with the smectic layering indicates a transition from algebraic decaying ordering to disorder. The smectic line shapes can be described by the Caillé correlation function convoluted with a finite-size factor represented by a stretched Gaussian (compressed exponential). The transition to disorder is signaled by a change in the exponent of the stretched Gaussian from 1 (simple Gaussian describing finite-size domains) via 0.5 (Lorentzian describing exponentially decaying short-range correlations) to <0.5 (stretched exponential correlations). For a flexible cross linker the changeover occurs for concentration between 0.15 and 0.20, for a stiff cross linker below about 0.10. Broadening of the higher harmonics of the x-ray peak indicates strong nonuniform strain within finite-size domains and local deformations induced by randomly distributed dislocations.

    Published 29 February 2008 (7 pages)
    021707  Full Text: PDF (929 kB)  | Buy Article
    +
    Show Abstract
    Liquid-crystal anchoring at a polymer surface arises from interactions at several different length scales. At the molecular level, a liquid-crystal molecule may tend to align with the substrate polymer chain, while at the nanometer length scale grooves can exist that arise from the periodic repeat structure of a polymer chain or from nanometer-scale undulations due to surface stresses. On a still longer scale there is the secondary effect of grooves or surface inhomogeneities. We have performed a total of more than 900 ns of atomistic molecular dynamics simulations in order to study the relative importance of the molecular-level interaction and the topography of the polymer surface in liquid-crystal anchoring. Substrates were constructed in which grooves were induced along a direction perpendicular to the constituent molecular chains. In the results presented for the case of 32 5CB molecules on a poly(vinyl alcohol) substrate, the liquid-crystal director orientation appeared to be determined principally by the substrate chain orientation. Only for the deepest grooves did the director align along the grooves and perpendicular to the substrate molecular chain direction.

    Polymers