We investigate the dynamical behavior of the Ising model under a zero-temperature quench with the initial fraction of up spins $0\ensuremath{\le}x\ensuremath{\le}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\ensuremath{\le}0.5$, persistence for the up (minority) spins shows the behavior ${P}_{\text{min}}(t)\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\gamma}}exp[\ensuremath{-}{(t/\ensuremath{\tau})}^{\ensuremath{\delta}}]$ with time
t, while for the down (majority) spins,
Pmaj(t) approaches a finite value. We find that the timescale $\ensuremath{\tau}$ diverges as ${({x}_{c}\ensuremath{-}x)}^{\ensuremath{-}\ensuremath{\lambda}}$, where
xc=0.5 and $\ensuremath{\lambda}\ensuremath{\simeq}2.31$. The exponent $\ensuremath{\gamma}$ varies as ${\ensuremath{\theta}}_{2d}+{c}_{0}{({x}_{c}\ensuremath{-}x)}^{\ensuremath{\beta}}$ where ${\ensuremath{\theta}}_{2d}\ensuremath{\simeq}0.215$ is very close to the persistence exponent in two dimensions; $\ensuremath{\beta}\ensuremath{\simeq}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\left[(\frac{x\ensuremath{-}{x}_{c}}{{x}_{c}}){L}^{1/\ensuremath{\nu}}\right]$, with $\ensuremath{\nu}\ensuremath{\approx}1.47$. This result suggests that $\ensuremath{\tau}\ensuremath{\sim}{L}^{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{z}}$, where $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{z}=\frac{\ensuremath{\lambda}}{\ensuremath{\nu}}=1.57\ifmmode\pm\else\textpm\fi{}0.11$ is an exponent not explored earlier.
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.
We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of the Ising model, employing the Glauber spin-flip mechanism, in space dimensions
d=2 and
3. Results for the persistence probability and the domain growth are discussed for quenches to various temperatures (
Tf) below the critical one (
Tc), from different initial temperatures
Ti≥Tc. In long time limit, for
Ti>Tc, the persistence probability exhibits power-law decay with exponents
θ≃0.22 and
≃0.18 in
d=2 and
3, respectively. For finite
Ti, the early time behavior is a different power-law whose life-time diverges and exponent decreases as
Ti→Tc. The crossover length between the two steps diverges as the equilibrium correlation length.
Ti=Tc is expected to provide a {\it{new universality class}} for which we obtain
θ≃0.035 in
d=2 and
≃0.10 in
d=3. The time dependence of the average domain size
ℓ, however, is observed to be rather insensitive to the choice of
Ti.
We present very accurate numerical estimates of the time and size dependence of the zero-temperature local persistence in the 2d ferromagnetic Ising model. We show that the effective exponent decays algebraically to an asymptotic value θ that depends upon the initial condition. More precisely, we find that θ takes one universal value 0.199(2) for initial conditions with short-range spatial correlations as in a paramagnetic state, and the value 0.033(1) for initial conditions with the long-range spatial correlations of the critical Ising state. We checked universality by working with a square and a triangular lattice, and by imposing free and periodic boundary conditions. We found that the effective exponent suffers from stronger finite size effects in the former case.