Scaling and persistence in the two-dimensional Ising model

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.

Zero-temperature coarsening in the two-dimensional long-range Ising model

We investigate the nonequilibrium dynamics following a quench to zero temperature of the nonconserved Ising model with power-law decaying long-range interactions ∝1/r^{d+σ} in d=2 spatial dimensions. The zero-temperature coarsening is always of special interest among nonequilibrium processes, because often peculiar behavior is observed. We provide estimates of the nonequilibrium exponents, viz., the growth exponent α, the persistence exponent θ, and the fractal dimension d_{f}. It is found that the growth exponent α≈3/4 is independent of σ and different from α=1/2, as expected for nearest-neighbor models. In the large σ regime of the tunable interactions only the fractal dimension d_{f} of the nearest-neighbor Ising model is recovered, while the other exponents differ significantly. For the persistence exponents θ this is a direct consequence of the different growth exponents α as can be understood from the relation d-d_{f}=θ/α; they just differ by the ratio of the growth exponents ≈3/2. This relation has been proposed for annihilation processes and later numerically tested for the d=2 nearest-neighbor Ising model. We confirm this relation for all σ studied, reinforcing its general validity.

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.

Minority-spin dynamics in the nonhomogeneous Ising model: Diverging time scales and exponents

We investigate the dynamical behavior of the Ising model under a zero-temperature quench with the initial fraction of up spins 0≤x≤1. In one dimension, the known results for persistence probability are verified; it shows algebraic decay for both up and down spins asymptotically with different exponents. It is found that the conventional finite-size scaling is valid here. In two dimensions, however, the persistence probabilities are no longer algebraic; in particular for x≤0.5, persistence for the up (minority) spins shows the behavior P_{min}(t)∼t^{-γ}exp[-(t/τ)^{δ}] with time t, while for the down (majority) spins, P_{maj}(t) approaches a finite value. We find that the timescale τ diverges as (x_{c}-x)^{-λ}, where x_{c}=0.5 and λ≃2.31. The exponent γ varies as θ_{2d}+c_{0}(x_{c}-x)^{β} where θ_{2d}≃0.215 is very close to the persistence exponent in two dimensions; β≃1. The results in two dimensions can be understood qualitatively by studying the exit probability, which for different system size is found to have the form E(x)=f[(x-x_{c}/x_{c})L^{1/ν}], with ν≈1.47. This result suggests that τ∼L^{z[over ̃]}, where z[over ̃]=λ/ν=1.57±0.11 is an exponent not explored earlier.

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.

Local and global persistence exponents of two quenched continuous-lattice spin models

Local and global persistence exponents associated with zero-temperature quenched dynamics of two-dimensional XY model and three-dimensional Heisenberg model have been estimated using numerical simulations. The method of block persistence has been used to find the global and local exponents simultaneously (in a single simulation). Temperature universality of both the exponents for three-dimensional Heisenberg model has been confirmed by simulating the stochastic (with noise) version of the equation of motion. The noise amplitudes added were small enough to retain the dynamics below criticality. In the second part of our work we have studied scaling associated with correlated persistence sites in the three-dimensional Heisenberg model in the later stages of the dynamics. The relevant length scale associated with correlated persistent sites was found to behave in a manner similar to the dynamic length scale associated with the phase ordering dynamics.

Nonequilibrium relaxations within the ground-state manifold in the antiferromagnetic Ising model on a triangular lattice

We present an extensive Monte Carlo simulation study on the nonequilibrium kinetics of triangular antiferromagnetic Ising model within the ground state ensemble which consists of sectors, each of which is characterized by a unique value of the string density p through a dimer covering method. Building upon our recent work [Phys. Rev. E 68, 066127 (2003)] where we considered the nonequilibrium relaxation observed within the dominant sector with p=2/3, we here focus on the nonequilibrium kinetics within the minor sectors with p<2/3. The initial configurations are chosen as those in which the strings are straight and evenly distributed. In the minor sectors, we observe a characteristic spatial anisotropy in both equilibrium and nonequilibrium spatial correlations. We observe emergence of a critical relaxation region (in the spatial and temporal domain) which grows as p deviates from p=2/3. Spatial anisotropy appears in the equilibrium spatial correlation with the characteristic length scale xi(e,V)(p) diverging with vanishing string density as xi(e,V)(p) approximately p(-2) along the vertical direction, while along the horizontal direction the spatial length scale diverges as xi(e,H) approximately p(-1). Analytic forms for the anisotropic equilibrium correlation functions are given. We also find that the spin autocorrelation function A(t) shows a simple scaling behavior A(t)=A(t/tau(A)(p)), where the time scale tau(A)(p) shows a power-law divergence with vanishing p as tau(A)(p) approximately p(-phi) with phi approximately or equal to 4. These features can be understood in terms of random walk nature of the fluctuations of the strings within the typical separation between neighboring strings.

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.

Persistence in q-state Potts model: a mean-field approach

We study the persistence properties of the T=0 coarsening dynamics of one-dimensional q-state Potts model using a modified mean-field approximation (MMFA). In this approximation, the spatial correlations between the interfaces separating spins with different Potts states is ignored, but the correct time dependence of the mean density P(t) of persistent spins is imposed. For this model, it is known that P(t) follows a power-law decay with time, P(t) approximately t(-theta(q)), where theta(q) is the q-dependent persistence exponent. We study the spatial structure of the persistent region within the MMFA. We show that the persistent site pair correlation function P2(r,t) has the scaling form P2(r,t)=P(t)(2)f(r/t(1/2)) for all values of the persistence exponent theta(q). The scaling function has the limiting behavior f(x) approximately x(-2theta) (x<<1) and f(x)-->1 (x>>1). We then show within the independent interval approximation (IIA) that the distribution n(k,t) of separation k between two consecutive persistent spins at time t has the asymptotic scaling form n(k,t)=t(-2phi)g(t,k/t(phi)), where the dynamical exponent has the form phi=max(1/2,theta). The behavior of the scaling function for large and small values of the arguments is found analytically. We find that for small separations k<<t(phi),n(k,t) approximately P(t)k(-tau), where tau=max[2(1-theta),2theta], while for large separations k>>t(phi), g(t,x) decays exponentially with x. The unusual dynamical scaling form and the behavior of the scaling function is supported by numerical simulations.

Persistence in cluster-cluster aggregation

Persistence is considered in one-dimensional diffusion-limited cluster-cluster aggregation when the diffusion coefficient of a cluster depends on its size s as D(s) approximately s(gamma). The probabilities that a site has been either empty or covered by a cluster all the time define the empty and filled site persistences. The cluster persistence gives the probability of a cluster remaining intact. The empty site and cluster persistences are universal whereas the filled site depends on the initial concentration. For gamma>0 the universal persistences decay algebraically with the exponent 2/(2-gamma). For the empty site case the exponent remains the same for gamma<0 but the cluster persistence shows a stretched exponential behavior as it is related to the small s behavior of the cluster size distribution. The scaling of the intervals between persistent regions demonstrates the presence of two length scales: the one related to the distances between clusters and that between the persistent regions.

Uninfected random walkers in one dimension

We consider a system of unbiased diffusing walkers (A(Phi)<-->(Phi)A) in one dimension with random initial conditions. We investigate numerically the relation between the fraction of walkers U(t) which have never encountered another walker up to time t, calling such walkers "uninfected" and the fraction of sites P(t) which have never been visited by a diffusing particle. We extend our study to include the A+B--> Phi diffusion-limited reaction in one dimension, with equal initial densities of A and B particles distributed homogeneously at t=0. We find U(t) approximately [P(t)]gamma, with gamma approximately 1.39, in both models, though there is evidence that a smaller value of gamma is required for t-->infinity.

Fraction of uninfected walkers in the one-dimensional Potts model

The dynamics of the one-dimensional q-state Potts model, in the zero-temperature limit, can be formulated through the motion of random walkers which either annihilate (A+A-->Phi) or coalesce (A+A-->A) with a q-dependent probability. We consider all of the walkers in this model to be mutually infectious. Whenever two walkers meet, they experience mutual contamination. Walkers which avoid an encounter with another random walker up to time t remain uninfected. The fraction of uninfected walkers is known to obey a power-law decay U(t) approximately t(-phi(q)), with a nontrivial exponent phi(q) [C. Monthus, Phys. Rev. E 54, 4844 (1996); S. N. Majumdar and S. J. Cornell, ibid. 57, 3757 (1998)]. We probe the numerical values of phi(q) to a higher degree of accuracy than previous simulations and relate the exponent phi(q) to the persistence exponent theta(q) [B. Derrida, V. Hakim, and V. Pasquier, Phys. Rev. Lett. 75, 751 (1995)], through the relation phi(q)=gamma(q)theta(q) where gamma is an exponent introduced in [S. J. O'Donoghue and A. J. Bray, preceding paper, Phys. Rev. E 65, 051113 (2002)]. Our study is extended to include the coupled diffusion-limited reaction A+A-->B, B+B-->A in one dimension with equal initial densities of A and B particles. We find that the density of walkers decays in this model as rho(t) approximately t(-1/2). The fraction of sites unvisited by either an A or a B particle is found to obey a power law, P(t) approximately t(-theta) with theta approximately 1.33. We discuss these exponents within the context of the q-state Potts model and present numerical evidence that the fraction of walkers which remain uninfected decays as U(t) approximately t(-phi), where phi approximately 1.13 when infection occurs between like particles only, and phi approximately 1.93 when we also include cross-species contamination. We find that the relation between phi and theta in this model can also be characterized by an exponent gamma, where similarly, phi=gamma(theta).