An extended chain Ising model and its Glauber dynamics

It was first proposed that an extended chain Ising (ECI) model contains the Ising chain model, single spin double-well potentials and a pure phonon heat bath of a specific energy exchange with the spins. The extension method is easy to apply to high dimensional cases. Then the single spin-flip probability (rate) of the ECI model is deduced based on the Boltzmann principle and general statistical principles of independent events and the model is simplified to an extended chain Glauber-Ising (ECGI) model. Moreover, the relaxation dynamics of the ECGI model were simulated by the Monte Carlo method and a comparison with the predictions of the special chain Glauber-Ising (SCGI) model was presented. It was found that the results of the two models are consistent with each other when the Ising chain length is large enough and temperature is relative low, which is the most valuable case of the model applications. These show that the ECI model will provide a firm physical base for the widely used single spin-flip rate proposed by Glauber and a possible route to obtain the single spin-flip rate of other form and even the multi-spin-flip rate.     

Metastability and the Ising model

We discuss a recent theorem which establishes a precise connection between (i) the approximate degeneracy of the zero eigenvalue for the generator of the Glauber dynamics of the Ising model in a small nonzero field and below the critical temperature, (ii) the existence of a partition of the configuration space into a normal region and a metastable region. This enables us to demonstrate that the recent approach to metastability of Davies and Martin may be viewed as a simple (although in some ways fairly crude) approximation to the conventional approach. We also obtain what appear to be the first results concerning the stability of metastable states under small perturbations.     

Investigation of the behavior of an Ising metamagnet under the influence of a field using Glauber dynamics

In this study, the kinetics of the Ising metamagnet where the interlayer interactions are ferromagnetic has been investigated under the mean field approximation. In describing the kinetics of the system, Glauber stochastic dynamics in the presence of an external field which performs time-dependent oscillations, has been utilized. Obtained results could be identified by two distinct types: the asymmetric solutions oscillating in the vicinity of finite values where the lattice magnetization has different values and the symmetric solutions being zero where the sublattice magnetizations are equal to each other. On the other hand, it has been observed that in the case where the system's initial state has a homogenous magnetization it exhibits two different periodical behaviors in the course of time.     

Zero-temperature dynamics for the ferromagnetic Ising model on random graphs

We consider Glauber dynamics at zero temperature for the ferromagnetic Ising model on the usual random graph model on N vertices, with on average γ edges incident to each vertex, in the limit as N→∞. Based on numerical simulations, Svenson (Phys. Rev. E 64 (2001) 036122) reported that the dynamics fails to reach a global energy minimum for a range of values of γ. The present paper provides a mathematically rigorous proof that this failure to find the global minimum in fact happens for all γ>0. A lower bound on the residual energy is also given.     

Anomalous dynamics in the Ising chain

We analyze a 1D Ising system with anomalous distributions of nearest neighbor interactions and show that the single-spin-flip dynamics exhibit breakdown of dynamic scaling. The results are obtained by a real-space numerical method applied to the exact equations of motion and they may be explained by domain wall motion arguments reformulated in terms of extreme value statistics.     

Damage spreading on the 3–12 lattice with competing Glauber and Kawasaki dynamics

The damage spreading of the Ising model on the 3–12 lattice with competing Glauber and Kawasaki dynamics is studied. The difference between the two kinds of nearest-neighboring spin interactions (interaction between two 12-gons, or interaction between a 12-gon and a triangle) are considered in the Hamiltonian. It is shown that the ratio of the interaction strengthF between the two kinds of interactions plays an important role in determining the critical temperature T_d of phase transition from frozen to chaotic. Two methods are used to introduce the bond dilution on the Ising model on the 3–12 lattice: regular and random. The maximum of the average damage spreading 〈D〉_max can approach values lower than 0.5 in both cases and the reason can be attributed to the ’survivors’ among the spins. We have also, for the first time, presented the phase diagram of the mixed G-K dynamics in the 3–12 lattice which shows what happens when going from pure Glauber to pure Kawasaki     

Some new results on the two-dimensional kinetic Ising model in the phase coexistence region

We consider a Glauber dynamics reversible with respect to the two-dimensional Ising model in a finite square of sideL with open boundary conditions, in the absence of an external field and at large inverse temperature . We prove that the gap in the spectrum of the generator restricted to the invariant subspace of functions which are even under global spin flip is much larger than the true gap. As a consequence we are able to show that there exists a new time scalet _even, much smaller than the global relaxation timet _rel, such that, with large probability, any initial configuration first relaxes to one of the two phases in a time scale of ordert _even and only after a time scale of the order oft _rel does it reach the final equilibrium by jumping, via a large deviation, to the opposite phase. It also follows that, with large probability, the time spent by the system during the first jump from one phase to the opposite one is much shorter than the relaxation time.     

Metastability and exponential approach to equilibrium for low-temperature stochastic Ising models

We present a proof of the exponential convergence to equilibrium of single-spin-flip stochastic dynamics for the two-dimensional Ising ferromagnet in the low-temperature case with not too small external magnetic fieldh uniformly in the volume and in the boundary conditions.     

Effect of disorder on relaxation in the one-dimensional Glauber model

We study the long-time relaxation of magnetization in a disordered linear chain of Ising spins from an initially aligned state. The coupling constants are ferromagnetic and nearest-neighbor only, taking valuesJ ₀ andJ ₁ with probabilitiesp and 1–p, respectively. The time evolution of the system is governed by the Glauber master equation. It is shown that for large timest, the magnetizationM(t) varies as [exp(–₀ t](t), where ₀ is a function of the stronger bond strengthJ ₀ only, and (t) decreases slower than an exponential. For very long times, we find that ln (t) varies as –t ^1/3. For low enough temperatures, there is an intermediate time regime when ln (t) varies as –t ^1/2. The results can be extended to more general probability distributions of ferromagnetic coupling constants, assuming thatM(t) can only increase if any bond in the chain is strengthened. If the coupling constants have a continuous distribution in which the probability density varies as a power law near some maximum valueJ ₀, we find that ln (t) varies as –t ^1/3(lnt)^2/3 for large times.     

Phase transition in a two-dimensional Ising ferromagnet based on the generalized zero-temperature Glauber dynamics

At zero temperature, based on the Ising model, the phase transition in a two-dimensional square lattice is studied using the generalized zero-temperature Glauber dynamics. Using Monte Carlo （MC） renormalization group methods, the static critical exponents and the dynamic exponent are studied; the type of phase transition is found to be of the first order.     

Droplet dynamics for asymmetric Ising model

Nucleation from a metastable state is studied for an anisotropic Ising model at very low temperatures. It turns out that the critical nucleus as well as configurations on a typical path to it differ from the Wulff shape of an equilibrium droplet.     

Metastability in Glauber Dynamics in the Low-Temperature Limit: Beyond Exponential Asymptotics

We consider Glauber dynamics of classical spin systems of Ising type in the limit when the temperature tends to zero in finite volume. We show that information on the structure of the most profound minima and the connecting saddle points of the Hamiltonian can be translated into sharp estimates on the distribution of the times of metastable transitions between such minima as well as the low lying spectrum of the generator. In contrast with earlier results on such problems, where only the asymptotics of the exponential rates is obtained, we compute the precise pre-factors up to multiplicative errors that tend to 1 as T0. As an example we treat the nearest neighbor Ising model on the 2 and 3 dimensional square lattice. Our results improve considerably earlier estimates obtained by Neves–Schonmann,⁽¹⁾ Ben Arous–Cerf,⁽²⁾ and Alonso–Cerf.⁽³⁾ Our results employ the methods introduced by Bovier, Eckhoff, Gayrard, and Klein in refs. 4 and 5.     

Ferroelectricity in diatomic Ising chain as investigated by elastic Ising model

An elastic Ising model for a one-dimensional diatomic spin chain is proposed to explain the ferroelectricity induced by the collinear magnetic order with a low-excited energy state. A statistical theory based on this model is developed to calculate the electrical and magnetic properties of Ca₃CoMnO₆, a typical quasi-one-dimensional diatomic spin chain system. The calculated ferroelectric polarization and dielectric susceptibility show a good agreement with recently reported data on Ca₃Co_2-xMn_xO₆ (x ≈0.96) (Phys. Rev. Lett. 100 047601 (2008)), although the predicted magnetic susceptibility does not coincide well with experiment. We also address the rationality and deficiency of this model by including a first-order correction which improves the consistency between the model and experiment.     

The Spectral Gap of the 2-D Stochastic Ising Model with Mixed Boundary Conditions

We establish upper bounds for the spectral gap of the stochastic Ising model at low temperatures in an l×l box with boundary conditions which are not purely plus or minus; specifically, we assume the magnitude of the sum of the boundary spins over each interval of length l in the boundary is bounded by l, where <1. We show that for any such boundary condition, when the temperature is sufficiently low (depending on ), the spectral gap decreases exponentially in l.     

The Spectral Gap of the 2-D Stochastic Ising Model with Nearly Single-Spin Boundary Conditions

We establish upper bounds for the spectral gap of the stochastic Ising model at low temperature in an N×N box, with boundary conditions which are plus except for small regions at the corners which are either free or minus. The spectral gap decreases exponentially in the size of the corner regions, when these regions are of size at least of order logN. This means that removing as few as O(logN) plus spins from the corners produces a spectral gap far smaller than the order N ^–2 gap believed to hold under the all-plus boundary condition. Our results are valid at all subcritical temperatures.     

Critical dynamics of the three-dimensional Ising model: A Monte Carlo study

Extensive Monte Carlo simulations have been performed to analyze the dynamical behavior of the three-dimensional Ising model with local dynamics. We have studied the equilibrium correlation functions and the power spectral densities of odd and even observables. The exponential relaxation times have been calculated in the asymptotic one-exponential time region. We find that the critical exponentz=2.09 ±0.02 characterizes the algebraic divergence with lattice size for all observables. The influence of scaling corrections has been analyzed. We have determined integrated relaxation times as well. Their dynamical exponentz _int agrees withz for correlations of the magnetization and its absolute value, but it is different for energy correlations. We have applied a scaling method to analyze the behavior of the correlation functions. This method verifies excellent scaling behavior and yields a dynamical exponentz _scal which perfectly agrees withz.     

Critical dynamics of finite Ising model

The Glauber kinetics of Ising spins is considered as a queueing process and simulated event by event as first proposed by Bortz, Kalos, and Lebowitz. The advantage of this algorithm compared to the standard single-flip Monte Carlo method is discussed for the situation of slowing down of dynamics. This process is used to generate fluctuations of magnetization and energy in the critical regimeT=T_c of two-dimensional Ising models. The analysis of these fluctuations leads to numerical determination of the critical exponents for dynamics: for the size dependence of correlation time atT _c , and for frequency dependence of the power spectrumS()~ ^–µ . From the finite-size scaling hypothesis, scaling relations are settled which are confirmed by this numerical experiment.     

Ferroelectricity in a diatomic Ising chain as investigated by the elastic Ising model

An elastic Ising model for a one-dimensional diatomic spin chain is proposed to explain the ferroelectricity induced by the collinear magnetic order with a low-excited energy state. A statistical theory based on this model is developed to calculate the electrical and magnetic properties of Ca₃CoMnO₆, a typical quasi-one-dimensional diatomic spin chain system. The calculated ferroelectric polarization and dielectric susceptibility show a good agreement with recently reported data on Ca₃Co_2-xMn_xO₆ (x ≈0.96) (Phys. Rev. Lett. 100 047601 (2008)), although the predicted magnetic susceptibility does not coincide well with experiment. We also address the rationality and deficiency of this model by including a first-order correction which improves the consistency between the model and experiment.     

Some New Results on the Kinetic Ising Model in a Pure Phase

We consider a general class of Glauber dynamics reversible with respect to the standard Ising model in ^d with zero external field and inverse temperature strictly larger than the critical value _c in dimension 2 or the so called slab threshold in dimension d 3. We first prove that the inverse spectral gap in a large cube of side N with plus boundary conditions is, apart from logarithmic corrections, larger than N in d = 2 while the logarithmic Sobolev constant is instead larger than N ² in any dimension. Such a result substantially improves over all the previous existing bounds and agrees with a similar computations obtained in the framework of a one dimensional toy model based on mean curvature motion. The proof, based on a suggestion made by H. T. Yau some years ago, explicitly constructs a subtle test function which forces a large droplet of the minus phase inside the plus phase. The relevant bounds for general d 2 are then obtained via a careful use of the recent –approach to the Wulff construction. Finally we prove that in d = 2 the probability that two independent initial configurations, distributed according to the infinite volume plus phase and evolving under any coupling, agree at the origin at time t is bounded from below by a stretched exponential , again apart from logarithmic corrections. Such a result should be considered as a first step toward a rigorous proof that, as conjectured by Fisher and Huse some years ago, the equilibrium time auto-correlation of the spin at the origin decays as a stretched exponential in d = 2.     

On the ergodic properties of Glauber dynamics

We show that if there is an infinite volume Gibbs measure which satisfies a logarithmic Sobolev inequality with local coefficients of moderate growth, then the corresponding stochastic dynamics decays to equilibrium exponentially fast in the uniform norm.