Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Master Equation for Ising Model

Steve P. Heims

Published Apr 19, 1965• 43 citations
Publisher

Abstract

A spin system with unperturbed Hamiltonian ${\mathcal{H}}_{0}=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{B}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}{{\ensuremath{\sigma}}_{x}}^{k}\ensuremath{-}\ensuremath{\gamma}H\ensuremath{\Sigma}{{\ensuremath{\sigma}}_{x}}^{j}$ relaxing via the spin-lattice coupling $G=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{C}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}({{\ensuremath{\sigma}}_{y}}^{k}+{{\ensuremath{\sigma}}_{z}}^{k})$ is studied by means of the general density-matrix theory of magnetic relaxation. By making some assumptions about the magnitude and time constants of the lattice correlation functions $〈{C}_{\mathrm{ij}}(t){C}_{\mathrm{kl}}(0)〉$, a master equation is obtained. It agrees at high temperatures with a master equation previously suggested by Glauber for the one-dimensional nearest-neighbor case. At high temperatures the magnetic moment relaxes with a single relaxation time, and the spin pair-correlation functions satisfy a closed set of equations. At low temperatures, however, the equations for the magnetization and the correlation functions are coupled to higher-order moments.

DOI
10.1103/physrev.138.a587

Related papers

On the dynamics of the Ising model (master equation/spin correlation functions/hierarchy of correlation e

Elliott W. Montroll

A two-dimensional (and to some degree three-di- mensional) version of Glauber's one-dimensional relaxation model is described. The model is constructed to yield the Ising model of cooperative phenomena at equilibrium. A complete hi- erarchy of differential equations for multispin correlation func- tions is constructed. Some remarks are made concerning the solu- tion of them for the initial value problem of determining the relaxation of an initial set of distributions. Investigations of the equilibrium thermodynamic properties of the Ising model have expanded our intuition on the general na- ture of phase transitions. They have also provided an accurate representation of the behavior of certain real systems well char- acterized by the model. Dynamical models of cooperative phe- nomena have proven harder to develop. Glauber (1) introduced an interesting master equation for the dynamics of the Ising spin variables on a one-dimensional chain. This equation is constructed so that the spin-flipping process at a given time is determined by two competing influences: (i) temperature ef- fects which at higher temperatures tend to randomize the sys- tem, and (ii) correlation effects of a given with its neighbors which tend to align a given with its neighbors. Although some literature has been generated for the investigation of two- dimensional and three-dimensional generalizations of the Glauber model (cf. Kawasaki (2) and Matsudaira (3, 4)), a de- tailed analysis of the appropriate hierarchy of rate equations for the dynamics of correlation functions has not yet been given. The purpose of this report is to make a contribution to such a program.

•Jan 1, 2017•0 citations

Relaxation of the Ising spin system coupled to a bosonic bath and the time dependent mean field equation

Máté Tibor Veszeli, Gábor Vattay

The Ising model does not have strictly defined dynamics, only a spectrum. There are different ways to equip it with time dependence, e.g., the Glauber or the Kawasaki dynamics, which are both stochastic, but it means there is a master equation that can also describe their dynamics. These equations can be derived from the Redfield equation, where the spin system is weakly coupled to a bosonic bath. In this paper, we investigate the temperature dependence of the relaxation time of a Glauber-type master equation, especially in the case of the fully connected, uniform Ising model. The finite-size effects were analyzed with a reduced master equation and the thermodynamic limit with a time-dependent mean field equation.

arXiv:2109.02341•PLoS ONE•Feb 28, 2022•1 citation

Glauber dynamics of the kinetic Ising model

Zhao‐Di Yang

In this work we study the Glauber dynamics of the one-dimensional Ising model with nearest-neighbor and next-nearest-neighbor interactions, for which an approximate solution of the magnetization per site is obtained. When the dynamical critical exponent z is investigated following the treatment of Cordery, Sarker, and Tobochnik [Phys. Rev. B 24, 5402 (1981)], our observation shows that its upper-bound value is the same as the known value, thus implying that z is independent of the range of the interaction. We also suggest a high-temperature expansion approximation which is then used to solve the two-dimensional Glauber dynamics governed by a master equation; this solution is compared with that of the decoupling method. The time-delayed correlation function is also calculated.

Physical review. B, Condensed matter•Nov 1, 1992•16 citations

Glauber dynamics in Ising systems: A general approach

Shaon Sahoo, Soumya Kanti Ganguly

In the Glauber model, the only restriction on the choice of the transition rates ($w_j$'s) is that these rates must satisfy the equation of detailed balance at equilibrium for all possible configurations. Though a linear form for $w_j$'s works best for analytical study of dynamics, generally it is not possible to get this form for which the above condition is obeyed. In this paper we discuss the best possible way of writing $w_j$'s in a linear form for a generic system where a spin is coupled to $z$ neighbors with different coupling constants. In our approach, the linear form of a $w_j$ involves $z$ independent parameters whose optimal values are derived by a linear regression. Generically, the Moore-Penrose pseudoinverse matrix involved in the regression is obtained solely from the configuration matrix and takes a simple form of dimension $z \times 2^z$. The fact that the pseudoinverse matrix does not depend on the Hamiltonian parameters or dimensionality of system makes our approach more appealing. In the second part of the paper we use our new approach to study Ising systems in different dimensions; in particular, we analyze the relaxation time and find out the critical temperatures ($T_C$) for two and three dimensional systems. It is easy to study the effect of magnetic field ($H$) within our approach; our analysis shows that $T_C$ decreases quadratically with (weak) magnetic field ($\delta T_C/T_C\propto -H^2$).

•Jan 21, 2014•0 citations

Dynamic Behavior of a Near-Ising System Near the Critical Point

Paul H. E. Meijer, John C. Edwards

Applying the master equation to an Ising model, equations of motion are obtained for the average of a single spin and for the average of a pair of spins. These equations relax to equilibrium values that predict a self-consistency relationship among the internal field, temperature, and exchange energy. The critical temperature obtained is that of the Bethe approximation, and is extended to include a static external field in the magnetization direction. With this, the relaxation of an Ising model near the Curie temperature has been investigated in a self-consistent manner. The results are similar to those of Suzuki and Kubo, but in the Bethe approximation instead of the Bragg-Williams approximation, and hence predict the relaxation of the correlation as well. In particular we studied a nearly Ising spin system that is a system with relatively small isotropic transverse coupling. An oscillating, linear, and external magnetic field applied transverse to the strongly coupled direction of the spins leads to lowering of the critical temperature.

Journal of Applied Physics•Mar 1, 1969•2 citations

II. Exact solutions of disordered Ising spin chains in a magnetic field

Stephen B. Haley

The one-dimensional Ising model with arbitrary nearest neighbor interactions and a magnetic field is formulated in terms of the spin correlation function relations derived in the preceding article I. The local magnetization and the two-spin correlation functions are shown to satisfy a set of coupled algebraic summation equations. The solution of this set is given in terms of a power series in hyberbolic tangents of the magnetic field strength, for chains of arbitrary length. For cases in which the interaction constants are random variables, explicit expressions are derived for the configurational averaged magnetization and the two-spin correlation functions to the lowest order in the magnetic field.

Journal of Mathematical Physics•May 1, 1978•3 citations
A spin system with unperturbed Hamiltonian H0=12\ensuremathΣ\ensuremathΣBjk\ensuremathσxj\ensuremathσxk\ensuremath−\ensuremathγH\ensuremathΣ\ensuremathσxj{\mathcal{H}}_{0}=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{B}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}{{\ensuremath{\sigma}}_{x}}^{k}\ensuremath{-}\ensuremath{\gamma}H\ensuremath{\Sigma}{{\ensuremath{\sigma}}_{x}}^{j}H0​=21​\ensuremathΣ\ensuremathΣBjk​\ensuremathσx​j\ensuremathσx​k\ensuremath−\ensuremathγH\ensuremathΣ\ensuremathσx​j relaxing via the spin-lattice coupling G=12\ensuremathΣ\ensuremathΣCjk\ensuremathσxj(\ensuremathσyk+\ensuremathσzk)G=\frac{1}{2}\ensuremath{\Sigma}\ensuremath{\Sigma}{C}_{\mathrm{jk}}{{\ensuremath{\sigma}}_{x}}^{j}({{\ensuremath{\sigma}}_{y}}^{k}+{{\ensuremath{\sigma}}_{z}}^{k})G=21​\ensuremathΣ\ensuremathΣCjk​\ensuremathσx​j(\ensuremathσy​k+\ensuremathσz​k) is studied by means of the general density-matrix theory of magnetic relaxation. By making some assumptions about the magnitude and time constants of the lattice correlation functions 〈Cij(t)Ckl(0)〉〈{C}_{\mathrm{ij}}(t){C}_{\mathrm{kl}}(0)〉〈Cij​(t)Ckl​(0)〉, a master equation is obtained. It agrees at high temperatures with a master equation previously suggested by Glauber for the one-dimensional nearest-neighbor case. At high temperatures the magnetic moment relaxes with a single relaxation time, and the spin pair-correlation functions satisfy a closed set of equations. At low temperatures, however, the equations for the magnetization and the correlation functions are coupled to higher-order moments.
A two-dimensional (and to some degree three-di- mensional) version of Glauber's one-dimensional relaxation model is described. The model is constructed to yield the Ising model of cooperative phenomena at equilibrium. A complete hi- erarchy of differential equations for multispin correlation func- tions is constructed. Some remarks are made concerning the solu- tion of them for the initial value problem of determining the relaxation of an initial set of distributions. Investigations of the equilibrium thermodynamic properties of the Ising model have expanded our intuition on the general na- ture of phase transitions. They have also provided an accurate representation of the behavior of certain real systems well char- acterized by the model. Dynamical models of cooperative phe- nomena have proven harder to develop. Glauber (1) introduced an interesting master equation for the dynamics of the Ising spin variables on a one-dimensional chain. This equation is constructed so that the spin-flipping process at a given time is determined by two competing influences: (i) temperature ef- fects which at higher temperatures tend to randomize the sys- tem, and (ii) correlation effects of a given with its neighbors which tend to align a given with its neighbors. Although some literature has been generated for the investigation of two- dimensional and three-dimensional generalizations of the Glauber model (cf. Kawasaki (2) and Matsudaira (3, 4)), a de- tailed analysis of the appropriate hierarchy of rate equations for the dynamics of correlation functions has not yet been given. The purpose of this report is to make a contribution to such a program.
The Ising model does not have strictly defined dynamics, only a spectrum. There are different ways to equip it with time dependence, e.g., the Glauber or the Kawasaki dynamics, which are both stochastic, but it means there is a master equation that can also describe their dynamics. These equations can be derived from the Redfield equation, where the spin system is weakly coupled to a bosonic bath. In this paper, we investigate the temperature dependence of the relaxation time of a Glauber-type master equation, especially in the case of the fully connected, uniform Ising model. The finite-size effects were analyzed with a reduced master equation and the thermodynamic limit with a time-dependent mean field equation.
In this work we study the Glauber dynamics of the one-dimensional Ising model with nearest-neighbor and next-nearest-neighbor interactions, for which an approximate solution of the magnetization per site is obtained. When the dynamical critical exponent z is investigated following the treatment of Cordery, Sarker, and Tobochnik [Phys. Rev. B 24, 5402 (1981)], our observation shows that its upper-bound value is the same as the known value, thus implying that z is independent of the range of the interaction. We also suggest a high-temperature expansion approximation which is then used to solve the two-dimensional Glauber dynamics governed by a master equation; this solution is compared with that of the decoupling method. The time-delayed correlation function is also calculated.
In the Glauber model, the only restriction on the choice of the transition rates (wjw_jwj​'s) is that these rates must satisfy the equation of detailed balance at equilibrium for all possible configurations. Though a linear form for wjw_jwj​'s works best for analytical study of dynamics, generally it is not possible to get this form for which the above condition is obeyed. In this paper we discuss the best possible way of writing wjw_jwj​'s in a linear form for a generic system where a spin is coupled to zzz neighbors with different coupling constants. In our approach, the linear form of a wjw_jwj​ involves zzz independent parameters whose optimal values are derived by a linear regression. Generically, the Moore-Penrose pseudoinverse matrix involved in the regression is obtained solely from the configuration matrix and takes a simple form of dimension z×2zz \times 2^zz×2z. The fact that the pseudoinverse matrix does not depend on the Hamiltonian parameters or dimensionality of system makes our approach more appealing. In the second part of the paper we use our new approach to study Ising systems in different dimensions; in particular, we analyze the relaxation time and find out the critical temperatures (TCT_CTC​) for two and three dimensional systems. It is easy to study the effect of magnetic field (HHH) within our approach; our analysis shows that TCT_CTC​ decreases quadratically with (weak) magnetic field (δTC/TC∝−H2\delta T_C/T_C\propto -H^2δTC​/TC​∝−H2).
Applying the master equation to an Ising model, equations of motion are obtained for the average of a single spin and for the average of a pair of spins. These equations relax to equilibrium values that predict a self-consistency relationship among the internal field, temperature, and exchange energy. The critical temperature obtained is that of the Bethe approximation, and is extended to include a static external field in the magnetization direction. With this, the relaxation of an Ising model near the Curie temperature has been investigated in a self-consistent manner. The results are similar to those of Suzuki and Kubo, but in the Bethe approximation instead of the Bragg-Williams approximation, and hence predict the relaxation of the correlation as well. In particular we studied a nearly Ising spin system that is a system with relatively small isotropic transverse coupling. An oscillating, linear, and external magnetic field applied transverse to the strongly coupled direction of the spins leads to lowering of the critical temperature.
The one-dimensional Ising model with arbitrary nearest neighbor interactions and a magnetic field is formulated in terms of the spin correlation function relations derived in the preceding article I. The local magnetization and the two-spin correlation functions are shown to satisfy a set of coupled algebraic summation equations. The solution of this set is given in terms of a power series in hyberbolic tangents of the magnetic field strength, for chains of arbitrary length. For cases in which the interaction constants are random variables, explicit expressions are derived for the configurational averaged magnetization and the two-spin correlation functions to the lowest order in the magnetic field.