Persistence and dynamics in the ANNNI chain

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.

Zero-temperature coarsening in the Ising model with asymmetric second-neighbor interactions in two dimensions

We consider the zero-temperature coarsening in the Ising model in two dimensions where the spins interact within the Moore neighborhood. The Hamiltonian is given by H=-∑_{〈i,j〉}S_{i}S_{j}-κ∑_{〈i,j^{'}〉}S_{i}S_{j^{'}}, where the two terms are for the first neighbors and second neighbors, respectively, and κ≥0. The freezing phenomenon, already noted in two dimensions for κ=0, is seen to be present for any κ. However, the frozen states show more complicated structure as κ is increased; e.g., local antiferromagnetic motifs can exist for κ>2. Finite-sized systems also show the existence of an isoenergetic active phase for κ>2, which vanishes in the thermodynamic limit. The persistence probability shows universal behavior for κ>0; however, it is clearly different from the κ=0 results when a nonhomogeneous initial condition is considered. Exit probability shows universal behavior for all κ≥0. The results are compared with other models in two dimensions having interactions beyond the first neighbor.

Nonequilibrium critical dynamics of the two-dimensional Ashkin-Teller model at the Baxter line

We investigate the short-time universal behavior of the two-dimensional Ashkin-Teller model at the Baxter line by performing time-dependent Monte Carlo simulations. First, as preparatory results, we obtain the critical parameters by searching the optimal power-law decay of the magnetization. Thus, the dynamic critical exponents θ_{m} and θ_{p}, related to the magnetic and electric order parameters, as well as the persistence exponent θ_{g}, are estimated using heat-bath Monte Carlo simulations. In addition, we estimate the dynamic exponent z and the static critical exponents β and ν for both order parameters. We propose a refined method to estimate the static exponents that considers two different averages: one that combines an internal average using several seeds with another, which is taken over temporal variations in the power laws. Moreover, we also performed the bootstrapping method for a complementary analysis. Our results show that the ratio β/ν exhibits universal behavior along the critical line corroborating the conjecture for both magnetization and polarization.

Low-temperature Glauber dynamics under weak competing interactions

We consider the low but nonzero-temperature regimes of the Glauber dynamics in a chain of Ising spins with first- and second-neighbor interactions J1,J2. For 0<-J2/|J1|<1 it is known that at T=0 the dynamics is both metastable and noncoarsening, while being always ergodic and coarsening in the limit of T→0+. Based on finite-size scaling analyses of relaxation times, here we argue that in that latter situation the asymptotic kinetics of small or weakly frustrated -J2/|J1| ratios is characterized by an almost ballistic dynamic exponent z≃1.03(2) and arbitrarily slow velocities of growth. By contrast, for noncompeting interactions the coarsening length scales are estimated to be almost diffusive.

Effect of the nature of randomness on quenching dynamics of the Ising model on complex networks

Randomness is known to affect the dynamical behavior of many systems to a large extent. In this paper we investigate how the nature of randomness affects the dynamics in a zero-temperature quench of the Ising model on two types of random networks. In both networks, which are embedded in a one-dimensional space, the first-neighbor connections exist and the average degree is 4 per node. In random model A the second-neighbor connections are rewired with a probability p, while in random model B additional connections between neighbors at a Euclidean distance l(l > 1) are introduced with a probability P(l) proportionally l(-α). We find that for both models, the dynamics leads to freezing such that the system gets locked in a disordered state. The point at which the disorder of the nonequilibrium steady state is maximum is located. The behavior of dynamical quantities such as residual energy, order parameter, and persistence are discussed and compared. Overall, the behavior of physical quantities are similar, although subtle differences are observed due to the difference in the nature of randomness.

Zero-temperature dynamics in the two-dimensional axial next-nearest-neighbor Ising model

We investigate the dynamics of a two-dimensional axial next-nearest-neighbor Ising model following a quench to zero temperature. The Hamiltonian is given by H= -J_(0) summation operator(L)_(i,j=1)S_(i,j)S_(i+1,j)-J_(1)summation operator_(i,j=1)(S_{i,j}S_{i,j+1}-kappaS_{i,j}S_{i,j+2}) . For kappa<1 , the system does not reach the equilibrium ground state but slowly evolves to a metastable state. For kappa>1 , the system shows a behavior similar to that of the two-dimensional ferromagnetic Ising model in the sense that it freezes to a striped state with a finite probability. The persistence probability shows algebraic decay here with an exponent theta=0.235+/-0.001 while the dynamical exponent of growth z=2.08+/-0.01 . For kappa=1 , the system belongs to a completely different dynamical class; it always evolves to the true ground state with the persistence and dynamical exponent having unique values. Much of the dynamical phenomena can be understood by studying the dynamics and distribution of the number of domain walls. We also compare the dynamical behavior to that of a Ising model in which both the nearest and next-nearest-neighbor interactions are ferromagnetic.

Global persistence exponent of the double-exchange model

We obtained the global persistence exponent for a continuous spin model on the simple cubic lattice with double-exchange interaction by using two different methods. First, we estimated the exponent theta(g) by following the time evolution of probability P(t) that the order parameter of the model does not change its sign up to time t[P(t) approximately t(-theta(g)]. Afterwards, that exponent was estimated through the scaling collapse of the universal function L(theta(g)(z)P(t) for different lattice sizes. Our results for both approaches are in very good agreement with each other.

Persistence and the random bond Ising model in two dimensions

We study the zero-temperature persistence phenomenon in the random bond +/-J Ising model on a square lattice via extensive numerical simulations. We find strong evidence for "blocking" regardless of the amount disorder present in the system. The fraction of spins which never flips displays interesting nonmonotonic, double-humped behavior as the concentration of ferromagnetic bonds is varied from zero to one. The peak is identified with the onset of the zero-temperature spin glass transition in the model. The residual persistence is found to decay algebraically and the persistence exponent theta(p) approximately = 0.9 over the range 0.1< or =p< or =0.9. Our results are completely consistent with the result of Gandolfi, Newman, and Stein for infinite systems that this model has "mixed" behavior, namely positive fractions of spins that flip finitely and infinitely often, respectively. [Gandolfi, Newman and Stein, Commun. Math. Phys. 214, 373 (2000).].