A kinetic one-dimensional Ising model on a ring evolves according to a generalization of Glauber rates, such that spins at even (odd) lattice sites experience a temperature T(e) (T(o)). Detailed balance is violated so that the spin chain settles into a nonequilibrium stationary state, characterized by multiple interactions of increasing range and spin order. We derive the equations of motion for arbitrary correlation functions and solve them to obtain an exact representation of the steady state. Two nontrivial amplitudes reflect the sublattice symmetries; otherwise, correlations decay exponentially, modulo the periodicity of the ring. In the long-chain limit, they factorize into products of two-point functions, in precise analogy to the equilibrium Ising chain. The exact solution confirms the expectation, based on simulations and renormalization group arguments, that the long-time, long-distance behavior of this two-temperature model is Ising-like, in spite of the apparent complexity of the stationary distribution.
A kinetic one-dimensional Ising model is coupled to two heat baths, such that spins at even (odd) lattice sites experience a temperature T-e (T-o). Spin-flips occur with Glauber-type rates generalized to the case of two temperatures. Driven by the temperature differential, the spin chain settles into a nonequilibrium steady state which corresponds to the stationary solution of a master equation. We construct a perturbation expansion of this master equation in terms of the temperature difference and compute explicitly the first two corrections to the equilibrium Boltzmann distribution. The key result is the emergence of additional spin operators in the steady state, increasing in spatial range and order of spin products. We comment on the violation of detailed balance and entropy production in the steady state.
A kinetic one-dimensional Ising model is coupled to two heat baths, such that spins at even (odd) lattice sites experience a temperature T e ( T o ). Spin-flips occur with Glauber-type rates generalized to the case of two temperatures. Driven by the temperature differential, the spin chain settles into a non-equilibrium steady state which corresponds to the stationary solution of a master equation. We construct a perturbation expansion of this master equation in terms of the temperature difference and compute explicitly the first two corrections to the equilibrium Boltzmann distribution. The key result is the emergence of additional spin operators in the steady state, increasing in spatial range and order of spin products. We comment on the violation of detailed balance and entropy production in the steady state.
The kinetic Glauber–Ising spin chain is one of the very few exactly solvable models of non-equilibrium statistical mechanics. Nevertheless, existing solutions do not yield tractable expressions for two-time correlation and response functions of observables involving products of more than one or two spins. We use a new approach to solve explicitly the full hierarchy of differential equations for the correlation and response functions. From this general solution follow closed expressions for arbitrary multispin two-time correlation and response functions, for the case where the system is quenched from equilibrium at T i > 0 to some arbitrary T ⩾ 0. By way of application, we give the results for two- and four-spin two-time correlation and response functions. From the standard mapping, these also imply new exact results for two-time particle correlation and response functions in one-dimensional diffusion limited annihilation.
Time-dependent statistics of the Ising model proposed by Glauber for the one-dimensional chain are extended to the example of a two-dimensional square lattice. Each spin is assumed to change its state through the interaction with a heat bath. The equations of motion for both the single spin and the spin correlation functions are solved approximately by using a decoupling procedure where the many-body correlation functions are taken as sums of products of pair correlation functions. As a special case, our theory allows the approximate calculation of the equilibrium properties of the system and it turns out that, in this case, our result is an improvement over the Bethe approximation. Both the frequency-dependent magnetic susceptibility and the decay of the magnetic moment to the equilibrium state are calculated above and below the Curie temperature. The fluctuation–dissipation theorem developed by Glauber for the linear chain is shown to hold in the two-dimensional case also.
The nonequilibrium dynamics of classical random Ising spin chains with nonconserved magnetization are studied using an asymptotically exact real space renormalization group (RSRG). We focus on random field Ising model (RFIM) spin chains with and without a uniform applied field, as well as on Ising spin glass chains in an applied field. For the RFIM we consider a universal regime where the random field and the temperature are both much smaller than the exchange coupling. In this regime, the Imry-Ma length that sets the scale of the equilibrium correlations is large and the coarsening of domains from random initial conditions (e.g., a quench from high temperature) occurs over a wide range of length scales. The two types of domain walls that occur diffuse in opposite random potentials, of the form studied by Sinai, and domain walls annihilate when they meet. Using the RSRG we compute many universal asymptotic properties of both the nonequilibrium dynamics and the equilibrium limit. We find that the configurations of the domain walls converge rapidly toward a set of system-specific time-dependent positions that are independent of the initial conditions. Thus the behavior of this nonequilibrium system is pseudodeterministic at long times because of the broad distributions of barriers that occur on the long length scales involved. Specifically, we obtain the time dependence of the energy, the magnetization, and the distribution of domain sizes (found to be statistically independent). The equilibrium limits agree with known exact results. We obtain the exact scaling form of the two-point equal time correlation function <S0(t)Sx(t)> and the two-time autocorrelations <S0(t('))S0(t)>. We also compute the persistence properties of a single spin, of local magnetization, and of domains. The analogous quantities for the +/-J Ising spin glass in an applied field are obtained from the RFIM via a gauge transformation. In addition to these we compute the two-point two-time correlation function <S0(t)Sx(t)><S0(t('))Sx(t('))> which can in principle be measured by experiments on spin-glass-like systems. The thermal fluctuations are studied and found to be dominated by rare events; in particular all moments of truncated equal time correlations are computed. Physical properties which are typically measured in aging experiments are also studied, focusing on the response to a small magnetic field which is applied after waiting for the system to equilibrate for a time t(w). The nonequilibrium fluctuation-dissipation ratio X(t,t(w)) is computed. We find that for (t-t(w)) approximately t(alpha)(w) with alpha;<1, the ratio equal to its equilibrium value X=1, although time translational invariance does not hold in this regime. For t-t(w) approximately t(w) the ratio exhibits an aging regime with a nontrivial X=X(t/t(w)) not equal to 1, but the behavior is markedly different from mean field theory. Finally the distribution of the total magnetization and of the number of domains is computed for large finite size systems. General issues about convergence toward equilibrium and the possibilities of weakly history-dependent evolution in other random systems are discussed.
We consider a one-dimensional lattice of Ising-type variables where the ferromagnetic exchange interaction
J between neighboring sites is frustrated by a long-ranged antiferromagnetic interaction of strength
g between the sites
i and
j, decaying as ${|i\ensuremath{-}j|}^{\ensuremath{-}\ensuremath{\alpha}}$, with $\ensuremath{\alpha}>1$. For $\ensuremath{\alpha}$ smaller than a certain threshold ${\ensuremath{\alpha}}_{0}$, which is larger than 2 and depends on the ratio
J/g, the ground state consists of an ordered sequence of segments with equal length and alternating magnetization. The width of the segments depends on both $\ensuremath{\alpha}$ and the ratio
J/g. Our Monte Carlo study shows that the on-site magnetization vanishes at finite temperatures and finds no indication of any phase transition. Yet, the modulation present in the ground state is recovered at finite temperatures in the two-point correlation function, which oscillates in space with a characteristic spatial period. The latter depends on $\ensuremath{\alpha}$ and
J/g and decreases smoothly from the ground-state value as the temperature is increased. Such an oscillation of the correlation function is exponentially damped over a characteristic spatial scale, the correlation length, which asymptotically diverges roughly as the inverse of the temperature as
T=0 is approached. This suggests that the long-range interaction causes the Ising chain to fall into a universality class consistent with an underlying continuous symmetry. The ${e}^{\ensuremath{\Delta}/T}$ temperature dependence of the correlation length and the uniform ferromagnetic ground state, characteristic of the
g=0 discrete Ising symmetry, are recovered for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{0}$.