Aging in ferromagnetic ordering: full decay and finite-size scaling of autocorrelation

Nonuniversal aging during phase separation with long-range interaction

Issues concerning the kinetics of phase transitions are not well established for the cases where the order parameter remains conserved with time, particularly when the interatomic interactions are long-range in nature. Here we present results on structure, growth, and aging from Monte Carlo simulations of the two-dimensional long-range Ising model. In our computer simulations, random initial configurations, for 50:50 compositions of up and down spins, mimicking high-temperature equilibrium states, have been quenched to temperatures inside the coexistence curve. Our analyses of the simulation data, for such a protocol, show interesting dependence of the aging exponent, λ, on σ, the parameter, within the Hamiltonian, that controls the range of interaction. These nonuniversal values of λ are compared with a theoretical result for lower bounds. For this purpose, we extracted information on relevant aspects of structural properties during the evolution. To estimate λ, as is necessary, we also calculated the average domain size and analyzed its time dependence to obtain the growth exponent α which also is nonuniversal. The trends in the values of λ and α, as well as an anomaly in structure, suggest that a crossover from the long-range to the short-range variety occurs at σ≃1. The location of this boundary and the nonuniversality provide a picture that is surprisingly different from that of the corresponding static critical phenomena. Furthermore, our results suggest an important scaling law combining α and λ.

Aging following a zero-temperature quench in the d=3 Ising model

Aging in phase-ordering kinetics of the d=3 Ising model following a quench from infinite to zero temperature is studied by means of Monte Carlo simulations. In this model the two-time spin-spin autocorrelator C_{ag} is expected to obey dynamical scaling and to follow asymptotically a power-law decay with the autocorrelation exponent λ. Previous work indicated that the lower Fisher-Huse bound of λ≥d/2=1.5 is violated in this model. Using much larger systems than previously studied, the instantaneous exponent for λ we obtain at late times does not disagree with this bound. By conducting systematic fits to the data of C_{ag} using different Ansätze for the leading correction term, we find λ=1.58(14), with most of the error attributed to the systematic uncertainty regarding the Ansätze. This result is in contrast to the recent report that below the roughening transition universality might be violated.

Influence of roughening transition on magnetic ordering

In the literature of magnetic phase transitions, in addition to a critical point, the existence of another special point has been discussed. This is related to the broadening of the interface between two different ordering phases and is referred to as the point of roughening transition. While the equilibrium properties associated with this transition are well understood, the influence of this on nonequilibrium dynamics still needs to be investigated. In this paper we present comprehensive results, from Monte Carlo simulations, on coarsening dynamics in a system, over a wide range of temperature, in space dimension d=3, for which there exists a roughening transition at a nonzero temperature T_{R}. An advanced analysis of the simulation data, on structure, growth, and aging, shows that the onset of unexpected glasslike slow dynamics in this system, that has received attention in recent times, for quenches to zero temperature, actually occurs at this transition point. This implies that the structure and aging depend upon the final temperature, when the latter lies between 0 and T_{R}. This is a very interesting exception to universality in coarsening dynamics. The results also demonstrate an important structure-dynamics connection in the phase-ordering dynamics. We compare the key results with those from d=2, for which there exists no nonzero roughening transition temperature. The absence of the above-mentioned anomalous features in the latter dimension places our conjecture on the role of the roughening transition on a firmer footing.

Phase-ordering kinetics in the Allen-Cahn (Model A) class: Universal aspects elucidated by electrically induced transition in liquid crystals

The two-dimensional (2D) Ising model is the statistical physics textbook example for phase transitions and their kinetics. Quenched through the Curie point with Glauber rates, the late-time description of the ferromagnetic domain coarsening finds its place at the scalar sector of the Allen-Cahn (or Model A) class, which encompasses phase-ordering kinetics endowed with a nonconserved order parameter. Resisting exact results sought for theoreticians since Lifshitz's first account in 1962, the central quantities of 2D Model A-most scaling exponents and correlation functions-remain known up to approximate theories whose disparate outcomes urge experimental assessment. Here we perform such assessment based on a comprehensive study of the coarsening of 2D twisted nematic liquid crystals whose kinetics is induced by a superfast electrical switching from a spatiotemporally chaotic (disordered) state to a two-phase concurrent, equilibrium one. Tracking the dynamics via optical microscopy, we first show the sharp evidence of well-established Model A aspects, such as the dynamic exponent z=2 and the dynamic scaling hypothesis, to then move forward. We confirm the Bray-Humayun theory for Porod's regime describing intradomain length scales of the two-point spatial correlators and show that their nontrivial decay beyond the Porod's scale can be captured in a free-from-parameter fashion by Gaussian theories, namely the Ohta-Jasnow-Kawasaki (OJK) and Mazenko theories. Regarding time-related statistics, we corroborate the aging hypothesis in Model A systems, which includes the collapse of two-time correlators into a master curve whose format is, actually, best accounted for by a solution of the local scaling invariance theory: the same solution that fits the 2D nonconserved Ising model correlator along with the Fisher-Huse conjecture. We also suggest the true value for the local persistence exponent in Model A class, in disfavor of the exact outcome for the diffusion and OJK equations. Finally, we observe a fractal morphology for persistence clusters and extract their universal dimension. Given its accuracy and possibilities, this experimental setup may work as a prototype to address further universality issues in the realm of nonequilibrium systems.

Kinetics of domain growth and aging in a two-dimensional off-lattice system

We have used molecular dynamics simulations for a comprehensive study of phase separation in a two-dimensional single-component off-lattice model where particles interact through the Lennard-Jones potential. Via state-of-the-art methods we have analyzed simulation data on structure, growth, and aging for nonequilibrium evolutions in the model. These data were obtained following quenches of well-equilibrated homogeneous configurations, with density close to the critical value, to various temperatures inside the miscibility gap, having vapor-"liquid" as well as vapor-"solid" coexistence. For the vapor-liquid phase separation we observe that ℓ, the average domain length, grows with time (t) as t^{1/2}, a behavior that has connection with hydrodynamics. At low-enough temperature, a sharp crossover of this time dependence to a much slower, temperature-dependent, growth is identified within the timescale of our simulations, implying "solid"-like final state of the high-density phase. This crossover is, interestingly, accompanied by strong differences in domain morphology and other structural aspects between the two situations. For aging, we have presented results for the order-parameter autocorrelation function. This quantity exhibits data collapse with respect to ℓ/ℓ_{w}, ℓ, and ℓ_{w} being the average domain lengths at times t and t_{w} (≤t), respectively, the latter being the age of a system. Corresponding scaling function follows a power-law decay: ∼(ℓ/ℓ_{w})^{-λ} for t≫t_{w}. The decay exponent λ, for the vapor-liquid case, is accurately estimated via the application of an advanced finite-size scaling method. The obtained value is observed to satisfy a bound.

Aging in the Long-Range Ising Model

The current understanding of aging phenomena is mainly confined to the study of systems with short-ranged interactions. Little is known about the aging of long-ranged systems. Here, the aging in the phase-ordering kinetics of the two-dimensional Ising model with power-law long-range interactions is studied via Monte Carlo simulations. The dynamical scaling of the two-time spin-spin autocorrelator is well described by simple aging for all interaction ranges studied. The autocorrelation exponents are consistent with λ=1.25 in the effectively short-range regime, while for stronger long-range interactions the data are consistent with λ=d/2=1. For very long-ranged interactions, strong finite-size effects are observed. We discuss whether such finite-size effects could be misinterpreted phenomenologically as subaging.

Aging exponents for nonequilibrium dynamics following quenches from critical points

Via Monte Carlo simulations we study nonequilibrium dynamics in the nearest-neighbor Ising model, following quenches to points inside the ordered region of the phase diagram. With the broad objective of quantifying the nonequilibrium universality classes corresponding to spatially correlated and uncorrelated initial configurations, in this paper we present results for the decay of the order-parameter autocorrelation function for quenches from the critical point. This autocorrelation is an important probe for the aging dynamics in far-from-equilibrium systems and typically exhibits power-law scaling. From the state-of-the-art analysis of the simulation results, we quantify the corresponding exponents (λ) for both conserved and nonconserved (order-parameter) dynamics of the model in space dimension d=3. Via structural analysis we demonstrate that the exponents satisfy a bound. We also revisit the d=2 case to obtain more accurate results. It appears that irrespective of the dimension, λ is approximately the same for both conserved and nonconserved dynamics.

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.

Zero-temperature ordering dynamics in a two-dimensional biaxial next-nearest-neighbor Ising model

We investigate the dynamics of a two-dimensional biaxial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H=-J_{0}∑_{i,j=1}^{L}[(S_{i,j}S_{i+1,j}+S_{i,j}S_{i,j+1})-κ(S_{i,j}S_{i+2,j}+S_{i,j}S_{i,j+2})]. For κ<1, the system does not reach the equilibrium ground state and keeps evolving in active states forever. For κ≥1, though, the system reaches a final state, but it does not reach the ground state always and freezes to a striped state with a finite probability like a two-dimensional ferromagnetic Ising model and an axial next-nearest-neighbor Ising (ANNNI) model. The overall dynamical behavior for κ>1 and κ=1 is quite different. The residual energy decays in a power law for both κ>1 and κ=1 from which the dynamical exponent z has been estimated. The persistence probability shows algebraic decay for κ>1 with an exponent θ=0.22±0.002 while the dynamical exponent for ordering z=2.33±0.01. For κ=1, the system belongs to a completely different dynamical class with θ=0.332±0.002 and z=2.47±0.04. We have computed the freezing probability for different values of κ. We have also studied the decay of autocorrelation function with time for a different regime of κ values. The results have been compared with those of the two-dimensional ANNNI model.

Direct test of the Gaussian auxiliary field ansatz in nonconserved order parameter phase ordering dynamics

The assumption that the local order parameter is related to an underlying spatially smooth auxiliary field, u(r[over ⃗],t), is a common feature in theoretical approaches to non-conserved order parameter phase separation dynamics. In particular, the ansatz that u(r[over ⃗],t) is a Gaussian random field leads to predictions for the decay of the autocorrelation function which are consistent with observations, but distinct from predictions using alternative theoretical approaches. In this paper, the auxiliary field is obtained directly from simulations of the time-dependent Ginzburg-Landau equation in two and three dimensions. The results show that u(r[over ⃗],t) is equivalent to the distance to the nearest interface. In two dimensions, the probability distribution, P(u), is well approximated as Gaussian except for small values of u/L(t), where L(t) is the characteristic length-scale of the patterns. The behavior of P(u) in three dimensions is more complicated; the non-Gaussian region for small u/L(t) is much larger than that in two dimensions but the tails of P(u) begin to approach a Gaussian form at intermediate times. However, at later times, the tails of the probability distribution appear to decay faster than a Gaussian distribution.

Pattern, growth, and aging in aggregation kinetics of a Vicsek-like active matter model

Via molecular dynamics simulations, we study kinetics in a Vicsek-like phase-separating active matter model. Quantitative results, for isotropic bicontinuous pattern, are presented on the structure, growth, and aging. These are obtained via the two-point equal-time density-density correlation function, the average domain length, and the two-time density autocorrelation function. Both the correlation functions exhibit basic scaling properties, implying self-similarity in the pattern dynamics, for which the average domain size exhibits a power-law growth in time. The equal-time correlation has a short distance behavior that provides reasonable agreement between the corresponding structure factor tail and the Porod law. The autocorrelation decay is a power-law in the average domain size. Apart from these basic similarities, the overall quantitative behavior of the above-mentioned observables is found to be vastly different from those of the corresponding passive limit of the model which also undergoes phase separation. The functional forms of these have been quantified. An exceptionally rapid growth in the active system occurs due to fast coherent motion of the particles, mean-squared-displacements of which exhibit multiple scaling regimes, including a long time ballistic one.

Pathways to equilibrium orientation fluctuations in finite stripe-forming systems

Small-angle orientation fluctuations in ordered stripe-forming systems free of topological defects can exhibit aging and anisotropic growth of two length scales. In infinitely extended systems, the stripe orientation field develops a dominant modulation length λ_{∥}^{*}(t) in the direction parallel to the stripes, which increases with time t as λ_{∥}^{*}(t)∼t^{1/4}. Simultaneously, the orientation correlation length ξ_{⊥}(t) in the direction perpendicular to the stripes increases as ξ_{⊥}(t)∼t^{1/2} [Riesch et al., Interface Focus 7, 20160146 (2017)2042-889810.1098/rsfs.2016.0146]. Here we show that finite systems of size L_{⊥}×L_{∥} with periodic boundary conditions reach equilibrium when the dominant modulation length λ_{∥}^{*}(t) reaches the system size L_{∥} in the stripe direction. The equilibration time τ_{eq}^{∥} is solely determined by L_{∥}, with τ_{eq}^{∥}∼L_{∥}^{4}. In systems with L_{⊥}<L_{∥}^{2}/2πλ_{p}, where λ_{p} is the undulation penetration length, the initial aging and coarsening dynamics changes at the crossover time τ_{C}^{⊥}∼L_{⊥}^{2} to an aging and coarsening dynamics described by the one-dimensional Mullins-Herring equation, before reaching equilibrium at τ_{∥}^{eq}. Our work reveals the two pathways to equilibrium in stripe phases with periodic boundary conditions, the finite-size scaling behavior of equilibrium orientation fluctuations, and the characteristic exponents associated with the influence of a finite system size.

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.

Evidence of aging and dynamic scaling in the collapse of a polymer

We investigate a newly framed two-time property for the nonequilibrium evolution dynamics during the collapse of a homopolymer via Monte Carlo simulations of a model polymer. Our results show evidence of aging effects, as observed in the slow dynamics of structural and spin glasses, along with the presence of a dynamic scaling of the autocorrelation function ∼x(-λ(c)) (x being the ratio of the cluster sizes at two different times). We estimate the value of λ(c) unambiguously by applying a finite-size scaling analysis to the numerical data. The value thus obtained obeys a bound which we predict via general theoretical arguments. The results presented should be of general validity and may trigger direct experimental verification in single-polymer dynamics.

Fractality in persistence decay and domain growth during ferromagnetic ordering: Dependence upon initial correlation

The dynamics of ordering in the Ising model, following quench to zero temperature, has been studied via Glauber spin-flip Monte Carlo simulations in space dimensions d=2 and 3. One of the primary objectives has been to understand phenomena associated with the persistent spins, viz., time decay in the number of unaffected spins, growth of the corresponding pattern, and its fractal dimensionality for varying correlation length in the initial configurations, prepared at different temperatures, at and above the critical value. It is observed that the fractal dimensionality and the exponent describing the power-law decay of persistence probability are strongly dependent upon the relative values of nonequilibrium domain size and the initial equilibrium correlation length. Via appropriate scaling analyses, these quantities have been estimated for quenches from infinite and critical temperatures. The above-mentioned dependence is observed to be less pronounced in the higher dimension. In addition to these findings for the local persistence, we present results for the global persistence as well. Furthermore, important observations on the standard domain growth problem are reported. For the latter, a controversy in d=3, related to the value of the exponent for the power-law growth of the average domain size with time, has been resolved.

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.