Growth of long-range correlations after a quench in conserved-order-parameter systems

Hyperuniformity in phase ordering: the roles of activity, noise, and non-constant mobility

Hyperuniformity emerges generically in the coarsening regime of phase-separating fluids. Numerical studies of active and passive systems have shown that the structure factorS(q) behaves asqςforq → 0, with hyperuniformity exponentς = 4. For passive systems, this result was explained in 1991 by a qualitative scaling analysis of Tomita, exploiting isotropy at scales much larger than the coarsening length. Here we reconsider and extend Tomita's argument to address cases of active phase separation and of non-constant mobility, again findingς = 4. We further show that dynamical noise of varianceDcreates a transientς = 2 regime forq^≪q^∗∼Dt[1-(d+2)ν]/2, crossing over toς = 4 at largerq^. Here,νis the coarsening exponent for the domain sizeℓ, such thatℓ(t)∼tν, andq^∝qℓis the rescaled wavenumber. In diffusive coarseningν=1/3, so the rescaled crossover wavevectorq^∗vanishes at large times whend⩾2. The slowness of this decay suggests a natural explanation for experiments that observe a long-livedς = 2 scaling in phase-separatingactivefluids (where noise is typically large). Conversely, ind = 1, we demonstrate that with noise theς = 2 regime survives ast→∞, withq^∗∼D5/6. (The structure factor is not then determined by the zero-temperature fixed point.) We confirm our analytical predictions by numerical simulations of continuum theories for active and passive phase separation in the deterministic case and of Model B for the stochastic case. We also compare them with related findings for a system near an absorbing-state transition rather than undergoing phase separation. A central role is played throughout by the presence or absence of a conservation law for the centre of mass positionRof the order parameter field.

Initial correlation dependence of aging in phase separating solid binary mixtures and ordering ferromagnets

Following quenches of initial configurations having long range spatial correlations, prepared at the demixing critical point, to points inside the miscibility gap, we study aging phenomena in solid binary mixtures. Results on the decay of the two-time order-parameter autocorrelation functions, obtained from Monte Carlo simulations of the two-dimensional Ising model, with Kawasaki exchange kinetics, are analyzed via state-of-the art methods. The outcome is compared with that obtained for the ordering in uniaxial ferromagnets. For the latter, we have performed Monte Carlo simulations of the same model using the Glauber mechanism. For both types of systems we provide comparative discussion of our results with reference to those concerning quenches with configurations having no spatial correlation. We also discuss the role of structure on the decay of these correlations.

Ferromagnetism-induced phase separation in a two-dimensional spin fluid

We study the liquid-gas phase separation observed in a system of repulsive particles dressed with ferromagnetically aligning spins, a so-called "spin fluid." Microcanonical ensemble numerical simulations of finite-size systems reveal that magnetization sets in and induces a liquid-gas phase separation between a disordered gas and a ferromagnetic dense phase at low enough energies and large enough densities. The dynamics after a quench into the coexistence region show that the order parameter associated with the liquid-vapor phase separation follows an algebraic law with an unusual exponent, as it is forced to synchronize with the growth of the magnetization: this suggests that for finite size systems the magnetization sets in along a Curie line, which is also the gas-side spinodal line, and that the coexistence region ends at a tricritical point. This picture is confirmed at the mean-field level with different approximation schemes, namely, a Bethe lattice resolution and a virial expansion complemented by the introduction of a self-consistent Weiss-like molecular field. However, a detailed finite-size scaling analysis shows that in two dimensions the ferromagnetic phase escapes the Berezinskii-Kosterlitz-Thouless scenario and that the long-range order is not destroyed by the unbinding of topological defects. The Curie line thus becomes a magnetic crossover in the thermodynamic limit. Finally, the effects of the magnetic interaction range and those of the interaction softness are characterized within a mean-field semianalytical low-density approach.

Ballistic aggregation in systems of inelastic particles: Cluster growth, structure, and aging

We study far-from-equilibrium dynamics in models of freely cooling granular gas and ballistically aggregating compact clusters. For both the cases, from event-driven molecular dynamics simulations, we have presented detailed results on structure and dynamics in space dimensions d=1 and 2. Via appropriate analyses it has been confirmed that the ballistic aggregation mechanism applies in d=1 granular gases as well. Aging phenomena for this mechanism, in both the dimensions, have been studied via the two-time density autocorrelation function. This quantity is demonstrated to exhibit scaling property similar to that in the standard phase transition kinetics. The corresponding functional forms have been quantified and the outcomes have been discussed in connection with the structural properties. Our results on aging establish a more complete equivalence between the granular gas and the ballistic aggregation models in d=1.

Dimensionality dependence of aging in kinetics of diffusive phase separation: Behavior of order-parameter autocorrelation

Behavior of two-time autocorrelation during the phase separation in solid binary mixtures is studied via numerical solutions of the Cahn-Hilliard equation as well as Monte Carlo simulations of the Ising model. Results are analyzed via state-of-the-art methods, including the finite-size scaling technique. Full forms of the autocorrelation in space dimensions 2 and 3 are obtained empirically. The long-time behavior is found to be power law, with exponents unexpectedly higher than the ones for the ferromagnetic ordering. Both Cahn-Hilliard and Ising models provide consistent results.

Critical domain-wall dynamics of model B

With Monte Carlo methods, we simulate the critical domain-wall dynamics of model B, taking the two-dimensional Ising model as an example. In the macroscopic short-time regime, a dynamic scaling form is revealed. Due to the existence of the quasirandom walkers, the magnetization shows intrinsic dependence on the lattice size L . An exponent which governs the L dependence of the magnetization is measured to be sigma=0.243(8) .

Autocorrelation exponent of conserved spin systems in the scaling regime following a critical quench

We study the autocorrelation function of a conserved spin system following a quench at the critical temperature. Defining the correlation length L(t) approximately t(1/z), we find that for times t' and t satisfying L(t')<<L(t)<<L(t')(phi) well inside the scaling regime, the spin autocorrelation function behaves like s(t)s(t') approximately L(t')(-(d-2+eta))[L(t')/L(t)](lambda(')(c)). For the O(n) model in the n-->infinity limit, we show that lambda(')(c)=d+2 and phi=z/2. We give a heuristic argument suggesting that this result is, in fact, valid for any dimension d and spin vector dimension n. We present numerical simulations for the conserved Ising model in d=1 and d=2, which are fully consistent with the present theory.

Experimental study of spinodal decomposition in a 1D conserved order parameter system

We study the zigzag instability coarsening of splay-bend walls formed in a nematic liquid crystal under external fields. The vertexes of zigzag can be considered as kinks in a one-dimensional order parameter system and the geometrical constraints associated with the necessary equal length sum of zig and zag segments impose a conserved quantity in this Cahn-Hilliard-type problem. In the late stage of coarsening, the characteristic length of the system L(t) shows a logarithmic increase in time and the dynamical scaling law holds. We then try to extract the nontrivial asymptotic scaling exponent lambda of the two-time correlation function, defined by lim(<phi(0,t)phi(0,t('))> approximately [L(t)/L(t('))](-lambda). The scaling exponents with respective time references, t(')=32 and 64 s, after quench are found to be lambda approximately 2 which is larger than the value with respective time reference t(')=0, predicted by numerical simulation.

Nonequilibrium critical dynamics of the triangular antiferromagnetic Ising model

We investigate the nonequilibrium critical dynamics of antiferromagnetic Ising model on a two-dimensional triangular lattice via dynamic Monte Carlo simulation employing spin-flip kinetics. Macroscopic degeneracy of the ground state originating from geometric frustration fundamentally affects the nonequilibrium dynamics of the system. In particular, the defects and the loose spins (whose flip costs no energy) play key roles in the dynamics. The long-time evolution is characterized by a critical dynamic scaling with a growing length scale xi(t). With random initial configurations, xi(t) exhibits a subdiffusive growth in time, xi(t) approximately t(1/z) with 1/z approximately 0.43, while xi(t) shows a diffusive growth with z=2 for the relaxation within the dominant sector of the ground-state manifold. The nonequilibrium critical dynamics therefore exhibits an interesting initial-state dependence. Persistence and the two-time temporal properties are also discussed.

Universality of the off-equilibrium response function in the kinetic Ising chain

The off-equilibrium response function chi(t,t(w)) and autocorrelation function C(t,t(w)) of an Ising chain with spin-exchange dynamics are studied numerically and compared with the same quantities in the case of spin-flip dynamics. It is found that, even though these quantities are different in the two cases, the parametric plot of chi(t,t(w)) versus C(t,t(w)) is the same. While this result could be expected in higher dimensionality, where chi(C) is related to the equilibrium state, it is far from trivial in the one-dimensional case where this relation does not hold. The origin of the universality of chi(C) is traced back to the optimization of domains position with respect to the perturbing external field. This mechanism is investigated resorting to models with a single domain moving in a random environment.

Ordering dynamics of heisenberg spins with torque: crossover, spin waves, and defects

We study the effect of a torque induced by the local molecular field on the phase ordering dynamics of the Heisenberg model when the total magnetization is conserved. The torque drives the zero-temperature ordering dynamics to a new fixed point, characterized by exponents z=2 and lambda approximately 5. This "torque-driven" fixed point is approached at times such that t>>g(2), where g is the strength of the torque. All physical quantities, like the domain size L(t) and the equal and unequal time correlation functions, obey a crossover scaling form over the entire range of g. An attempt to understand this crossover behavior from the approximate Gaussian closure scheme fails completely, implying that the dynamics at late times cannot be understood from the dynamics of defects alone. We provide convincing arguments that the spin configurations can be decomposed in terms of defects and spin waves which interact with each other even at late times. In the absence of the torque term, the spin waves decay faster, but even so we find that the Gaussian closure scheme is inconsistent. In the latter case the inconsistency may be remedied by including corrections to a simple Gaussian distribution. For completeness we include a discussion of the ordering dynamics at T(c), where the torque is shown to be relevant, with exponents z=4-varepsilon/2 and lambda=d (where varepsilon=6-d). We show to all orders in perturbation theory that lambda=d as a consequence of the conservation law.