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

Modern Physics Letters B

The Glauber and Metropolis Transition Rates on the Stationary States of the Ising Model

G. L. S. Paula, W. Figueiredo

Published Jun 10, 1997• 0 citations
Publisher

Abstract

We have applied the Glauber and Metropolis prescriptions to investigate the stationary states of the Ising model in one and two dimensions. We have employed the formalism of the master equation to follow the evolution of the system towards the stationary states. Although the Glauber and Metropolis transition rates lead the system to the same equilibrium states for the Ising model in the Monte Carlo simulations, we show that they can predict different results if we disregard the correlations between spins. The critical temperature of the one-dimensional Ising model cannot even be found by using the Metropolis algorithm and the mean field approximation. However, taking into account only correlations between nearest neighbor spins, the resulting stationary states become identical for both Glauber and Metropolis transition rates.

DOI
10.1142/s0217984997000694

Related papers

Critical behavior of the classical spin-1 Ising model: a combined low-temperature series expansion and Metropolis Monte Carlo analysis

Amir Taheridehkordi, Roberto Zivieri

In this paper, we theoretically study the critical properties of the classical spin-1 Ising model using two approaches: 1) the analytical low-temperature series expansion and 2) the numerical Metropolis Monte Carlo technique. Within this analysis, we discuss the critical behavior of one-, two- and three-dimensional systems modeled by the first-neighbor spin-1 Ising model for different types of exchange interactions. The comparison of the results obtained according the Metropolis Monte Carlo simulations allows us to highlight the limits of the widely used mean-field theory approach. We also show, via a simple transformation, that for the special case where the bilinear and bicubic terms are set equal to zero in the Hamiltonian the partition function of the spin-1 Ising model can be reduced to that of the spin-1/2 Ising model with temperature dependent external field and temperature independent exchange interaction times an exponential factor depending on the other terms of the Hamiltonian and confirm this result numerically by using the Metropolis Monte Carlo simulation. Finally, we investigate the dependence of the critical temperature on the strength of long-range interactions included in the Ising Hamiltonian comparing it with that of the first-neighbor spin-1/2 Ising model.

arXiv:2007.08593•Jul 16, 2020•0 citations

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

Nonequilibrium phase transition in an Ising model without detailed balance

Manoj Kumar, Chandan Dasgupta

We study a two-dimensional ferromagnetic Ising model in which spins are updated using modified versions of the Metropolis and Glauber algorithms. These update rules do not obey the detailed balance condition. The steady-state behavior of the model is studied using molecular field theory and Monte Carlo simulations. This model is found to exhibit a nonequilibrium phase transition from a "paramagnetic" state with zero magnetization to a "ferromagnetic" state with nonzero magnetization as the variable that plays the role of temperature in the spin updates is decreased. From detailed Monte Carlo simulations using the modified Metropolis algorithm, we demonstrate explicitly the nonequilibrium nature of the transition and show that it cannot be described as an equilibrium transition with an effective temperature different from the temperature used in the spin updates. The critical exponents that characterize the singular behavior near this continuous phase transition are calculated from finite size scaling of specific heat, magnetization, susceptibility, and correlation length. We find that the values of these exponents are the same (within error bars) as those of the equilibrium Ising model in two dimensions.

Physical review. E•Nov 9, 2020•17 citations

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

Time-Dependent Statistics of the Ising Model in a Magnetic Field

Huey Wen Huang

It is of interest to solve the master equation of a one-dimensional Ising model in a magnetic field (Glauber's equation) not only as a mathematical problem, but also for some practical applications. While an exact treatment of Glauber's equation seems impossible, we consider its approximate solutions. We review the applications of the perturbation-expansion method and applications of the mean-field theory to this problem. We also introduce a local-equilibrium method that improves the mean-field theory by taking the short-range correlation into account. The new method yields an analytically soluble differential equation, which, contrary to the preceding two methods, has an exact steady-state solution.

Physical Review A•Nov 1, 1973•19 citations

Monte Carlo simulation of an antiferromagnetic Ising model at two competing temperatures

B. C. S. Grandi, W. Figueiredo

We consider a two-dimensional antiferromagnet Ising system interacting with a heat bath at temperature T. The dynamics of the system is simulated by two competing stochastic processes: the two-spin-exchange Kawasaki kinetics at temperature T>0 and the one-spin-flip Glauber dynamics at T(G)-->0(-), which mimics the increase of the energy of the system. These two processes have probabilities 1-p and p, respectively. Monte Carlo simulations were employed to determine the phase diagram for the stationary states of the model and the corresponding critical exponents. Contrary to the ferromagnetic case, the phase diagram obtained does not exhibit the phenomenon of self-organization: for any nonzero value of the competing parameter p, and for any value of T, the only stationary phase which remains is the ferromagnetic one. At the phase transition between the antiferromagnetic and paramagnetic phases, at p=0, the values found for the critical exponents agree with those of the corresponding equilibrium Ising model.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•May 1, 1999•1 citation
We have applied the Glauber and Metropolis prescriptions to investigate the stationary states of the Ising model in one and two dimensions. We have employed the formalism of the master equation to follow the evolution of the system towards the stationary states. Although the Glauber and Metropolis transition rates lead the system to the same equilibrium states for the Ising model in the Monte Carlo simulations, we show that they can predict different results if we disregard the correlations between spins. The critical temperature of the one-dimensional Ising model cannot even be found by using the Metropolis algorithm and the mean field approximation. However, taking into account only correlations between nearest neighbor spins, the resulting stationary states become identical for both Glauber and Metropolis transition rates.
In this paper, we theoretically study the critical properties of the classical spin-1 Ising model using two approaches: 1) the analytical low-temperature series expansion and 2) the numerical Metropolis Monte Carlo technique. Within this analysis, we discuss the critical behavior of one-, two- and three-dimensional systems modeled by the first-neighbor spin-1 Ising model for different types of exchange interactions. The comparison of the results obtained according the Metropolis Monte Carlo simulations allows us to highlight the limits of the widely used mean-field theory approach. We also show, via a simple transformation, that for the special case where the bilinear and bicubic terms are set equal to zero in the Hamiltonian the partition function of the spin-1 Ising model can be reduced to that of the spin-1/2 Ising model with temperature dependent external field and temperature independent exchange interaction times an exponential factor depending on the other terms of the Hamiltonian and confirm this result numerically by using the Metropolis Monte Carlo simulation. Finally, we investigate the dependence of the critical temperature on the strength of long-range interactions included in the Ising Hamiltonian comparing it with that of the first-neighbor spin-1/2 Ising model.
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.
We study a two-dimensional ferromagnetic Ising model in which spins are updated using modified versions of the Metropolis and Glauber algorithms. These update rules do not obey the detailed balance condition. The steady-state behavior of the model is studied using molecular field theory and Monte Carlo simulations. This model is found to exhibit a nonequilibrium phase transition from a "paramagnetic" state with zero magnetization to a "ferromagnetic" state with nonzero magnetization as the variable that plays the role of temperature in the spin updates is decreased. From detailed Monte Carlo simulations using the modified Metropolis algorithm, we demonstrate explicitly the nonequilibrium nature of the transition and show that it cannot be described as an equilibrium transition with an effective temperature different from the temperature used in the spin updates. The critical exponents that characterize the singular behavior near this continuous phase transition are calculated from finite size scaling of specific heat, magnetization, susceptibility, and correlation length. We find that the values of these exponents are the same (within error bars) as those of the equilibrium Ising model in two dimensions.
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.
It is of interest to solve the master equation of a one-dimensional Ising model in a magnetic field (Glauber's equation) not only as a mathematical problem, but also for some practical applications. While an exact treatment of Glauber's equation seems impossible, we consider its approximate solutions. We review the applications of the perturbation-expansion method and applications of the mean-field theory to this problem. We also introduce a local-equilibrium method that improves the mean-field theory by taking the short-range correlation into account. The new method yields an analytically soluble differential equation, which, contrary to the preceding two methods, has an exact steady-state solution.
We consider a two-dimensional antiferromagnet Ising system interacting with a heat bath at temperature T. The dynamics of the system is simulated by two competing stochastic processes: the two-spin-exchange Kawasaki kinetics at temperature T>0 and the one-spin-flip Glauber dynamics at T(G)-->0(-), which mimics the increase of the energy of the system. These two processes have probabilities 1-p and p, respectively. Monte Carlo simulations were employed to determine the phase diagram for the stationary states of the model and the corresponding critical exponents. Contrary to the ferromagnetic case, the phase diagram obtained does not exhibit the phenomenon of self-organization: for any nonzero value of the competing parameter p, and for any value of T, the only stationary phase which remains is the ferromagnetic one. At the phase transition between the antiferromagnetic and paramagnetic phases, at p=0, the values found for the critical exponents agree with those of the corresponding equilibrium Ising model.