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

Progress of Theoretical Physics

Time-Dependent Ising Model with Long Range Interaction

Kyozi Kawasaki, 知司 山田

Published Jan 1, 1968• 18 citations
Publisher

Abstract

A systematic method is developed for studying the time-dependent Ising model proposed by Glauber and extended by others. The graphical technique developed for the equilibrium statistical mechanics has been found to be useful in this nonequilibrium problem as well with some modifications. A hierarchy of equations linking various spin cluster functions is obtained. If we choose as the expansion parameter the reciprocal range of the spin interaction which is taken to be of Kac's type, the hierarchy is terminated, and closed sets of nonlinear kinetic equations satisfied by spin cluster functions follow. The dynamical susceptibility is studied using the hierachy, and the polydispersive nature of the spin relaxation is found to be closely connected to the deviation of the spin pair correlation from the Ornstein-Zernike form via multi-spin cluster functions. A sum rule on the distribution of relaxation times is found, which provides a precise expression for the critical slowing-down at the Curie point. The closed set of kinetic equations involving at most three spin cluster functions is used to evaluate explicitly the approximate dynamical susceptibility above the Curie point, and the characteristic time-scale of the polydispersive relaxation is found to increase indefinitely near the Curie point.

DOI
10.1143/ptp.39.1

Related papers

Theory for the dynamics of clusters near the critical point. I. Relaxation of the Glauber kinetic Ising model

Kurt Binder, Dietrich Stauffer, Heiner Muller-Krumbhaar

A semiphenomenological cluster theory is developed for dynamic critical properties, which is not limited to small deviations from equilibrium. Explicit numerical expressions are derived for linear and nonlinear response functions of the kinetic Ising model, which are compatible with dynamic scaling and recent calculations of dynamic critical exponents. An uniaxial magnet (or Ising system) and its critical fluctuations are approximated by a "droplet" or "cluster" model, the critical part of the free energy being the sum of the contributions of clusters each containing $l$ spins. These clusters are assumed to grow and shrink with a phenomenological rate $\ensuremath{\propto}{l}^{r}$, $r\ensuremath{\approx}1$. These rates are compatible with direct Monte Carlo simulations in two and three dimensions. The resulting cluster reaction and diffusion equation is an approximation to the master equation of the Glauber kinetic Ising model (in which the magnetization is not conserved), and has the same structure as used in nucleation theories. In linear response to space- and time-dependent fields the relaxation times of energy and magnetization are expressed as triple integrals. All relaxation times diverge as ${|1\ensuremath{-}\frac{T}{{T}_{c}}|}^{\ensuremath{-}(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}}$. This treatment is consistent with dynamic scaling and universality, and is more general than the conventional (Van Hove) theory of critical slowing down. Also the (smaller) exponents found for the response to localized variations ("autocorrelation functions") agree with dynamic scaling. The wave-vector dependence of the relaxation times arises from the static correlation function only, and not from the cluster diffusion term. If the equilibrium cluster distribution is assumed to be that of the Fisher droplet model, and for special values of $r$, the (discrete) eigenvalue spectrum and eigenfunctions (generalized Laguerre polynomials) of the cluster reaction equation are found explicitly. Relaxation functions and frequency-dependent susceptibilities are expressed by hypergeometric series, even in the case of nonlinear response. At the critical point, the spectrum becomes continuous and the exponential decay for large times is replaced by power-law behavior. Similarly, it is found that dynamic scaling applies also to the nucleation rate for this model, and the scaling exponent of the nucleation rate is $2\ensuremath{-}\ensuremath{\alpha}+(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}$.

Physical review. B, Solid state•Dec 1, 1975•98 citations

Time-Dependent Spin-Correlation Function in Ferromagnet and Antiferromagnet

Takehiko Oguchi, Ikuo Ono

A statistical theory of the localized model of ferromagnetism is developed by considering the effective density matrices for a single spin, and two nearest-neighbor spins, respectively. The Ising Model and the Heisenberg model, and even the case of the anisotropic exchange interaction, are treated within a unified mathematical scheme. The result is reduced to that of the Bethe approximation in the case of the Ising model, and the constant coupling approximation in the case of the Heisenberg model. The method is generalized to calculate the static spin-correlation function. It is obtained in Ornstein-Zernike type. In the nonequilibrium state, the expectation value of the spin decays by the homogeneous damping as well as the spin-diffusion process. This relaxation equation is different from that by Van Hove because of the existence of the homogeneous damping. The diffusion constant does not vanish even at the Curie point TC, contrary to Van Hove's result, although the damping constant vanishes at TC. The latter fact causes the critical slowing down of the regression of fluctuation. The solution of this relaxation equation gives the time-dependent spin-correlation function. The method can be applied to the antiferromagnetic case, too.

Journal of Applied Physics•Feb 1, 1968•4 citations

Investigation of the Ising-type models within cluster approach

Р. Р. Левицкий, С. И. Сороков

The method for calculation of the correlation functions of the Ising-type systems with short-range interaction and with arbitrary value of spin is developed within cluster approximation. For the Ising model (spin $S^z=\pm1$) the expressions for pair and ternary correlation functions within two-particle approximation in $\set{q}$-space are obtained for the hypercubic Bravais lattices. In the 1D case the exact expressions for them in the site space is obtained as well. On the basis of the Glauber equation within two-particle cluster approximation the longitudinal dynamical susceptibility $χ(\set{q},E)$ is found. In the 1D case and in the absence of external field the expression for $χ(\set{q},E)$ is exact. For the Emery-Blume-Griffiths model ($S^z=-2,0,2$) within two-particle approximation the pair correlation functions are calculated. The four-particle cluster approximation is used for calculation of static susceptibility $χ(\set{q})$ of $KD_2PO_4$ ferroelectrics.

arXiv:cond-mat/9606134•Jun 18, 1996•0 citations

SOME DYNAMICAL PROPERTIES OF THE ISING FERROMAGNET

Noboru Matsudaira

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.

Canadian Journal of Physics•Jun 1, 1967•28 citations

Monte Carlo Investigation of Dynamic Critical Phenomena in the Two-Dimensional Kinetic Ising Model

E. Stoll, Kurt Binder, T. Schneider

Extending the Monte Carlo method to dynamic critical phenomena we investigated the time-dependent correlation functions in the two-dimensional one-spin-flip Ising model and the critical behavior of the associated relaxation times. These relaxation times are the following: ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}}^{\ensuremath{\Delta}T}$, characterizing the approach of the order parameter to equilibrium after a change of temperature $\ensuremath{\Delta}T$ of the system; ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}$ and ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}^{A}$ characterizing the slowing down of the order-parameter correlation and autocorrelation functions, respectively; ${\ensuremath{\tau}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}$ and ${\ensuremath{\tau}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}^{A}$, characterizing the slowing down of the energy correlation and autocorrelation functions; and finally ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}$, characterizing the cross-correlation function. We give estimates for the associated exponents ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}}^{\ensuremath{\Delta}T}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}\ensuremath{\approx}1.90\ifmmode\pm\else\textpm\fi{}0.10$, and ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}^{A}\ensuremath{\approx}1.60\ifmmode\pm\else\textpm\fi{}0.10$, ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}^{A}\ensuremath{\approx}0.95\ifmmode\pm\else\textpm\fi{}0.10$, ${\ensuremath{\Delta}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}^{A}\ensuremath{\approx}0$, which are consistent with the dynamic scaling hypothesis and with exact inequalities. A detailed comparison with recent high-temperature-expansion estimates is performed, and the reliability of the Monte Carlo results is carefully analyzed.

Physical review. B, Solid state•Oct 1, 1973•181 citations

Critical Slowing Down in the Kinetic Ising Model

Hideo Yahata, Masuo Suzuki

A simple dynamical model of interacting Ising spins is discussed near the critical point. It is shown that all the moments are positive and never vanish even at the critical point. Lower and upper bounds of the susceptibility are obtained, which yield the non-vanishing of the real part of the susceptibility. This gives a simple proof for the poly-dispersive nature of our system. The method of high temperature expansions for the dynamical susceptibility is presented and is applied to the two-dimensional Ising model up to the ninth order, the results of which show that the critical index of slowing down (the relaxation time τ (1) ∝( T - T c ) - Δ ; Δ =2±0.05) is different from that of the static susceptibility.

Journal of the Physical Society of Japan•Dec 1, 1969•143 citations
A systematic method is developed for studying the time-dependent Ising model proposed by Glauber and extended by others. The graphical technique developed for the equilibrium statistical mechanics has been found to be useful in this nonequilibrium problem as well with some modifications. A hierarchy of equations linking various spin cluster functions is obtained. If we choose as the expansion parameter the reciprocal range of the spin interaction which is taken to be of Kac's type, the hierarchy is terminated, and closed sets of nonlinear kinetic equations satisfied by spin cluster functions follow. The dynamical susceptibility is studied using the hierachy, and the polydispersive nature of the spin relaxation is found to be closely connected to the deviation of the spin pair correlation from the Ornstein-Zernike form via multi-spin cluster functions. A sum rule on the distribution of relaxation times is found, which provides a precise expression for the critical slowing-down at the Curie point. The closed set of kinetic equations involving at most three spin cluster functions is used to evaluate explicitly the approximate dynamical susceptibility above the Curie point, and the characteristic time-scale of the polydispersive relaxation is found to increase indefinitely near the Curie point.
A semiphenomenological cluster theory is developed for dynamic critical properties, which is not limited to small deviations from equilibrium. Explicit numerical expressions are derived for linear and nonlinear response functions of the kinetic Ising model, which are compatible with dynamic scaling and recent calculations of dynamic critical exponents. An uniaxial magnet (or Ising system) and its critical fluctuations are approximated by a "droplet" or "cluster" model, the critical part of the free energy being the sum of the contributions of clusters each containing lll spins. These clusters are assumed to grow and shrink with a phenomenological rate $\ensuremath{\propto}{l}^{r}$, $r\ensuremath{\approx}1$. These rates are compatible with direct Monte Carlo simulations in two and three dimensions. The resulting cluster reaction and diffusion equation is an approximation to the master equation of the Glauber kinetic Ising model (in which the magnetization is not conserved), and has the same structure as used in nucleation theories. In linear response to space- and time-dependent fields the relaxation times of energy and magnetization are expressed as triple integrals. All relaxation times diverge as ${|1\ensuremath{-}\frac{T}{{T}_{c}}|}^{\ensuremath{-}(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}}$. This treatment is consistent with dynamic scaling and universality, and is more general than the conventional (Van Hove) theory of critical slowing down. Also the (smaller) exponents found for the response to localized variations ("autocorrelation functions") agree with dynamic scaling. The wave-vector dependence of the relaxation times arises from the static correlation function only, and not from the cluster diffusion term. If the equilibrium cluster distribution is assumed to be that of the Fisher droplet model, and for special values of rrr, the (discrete) eigenvalue spectrum and eigenfunctions (generalized Laguerre polynomials) of the cluster reaction equation are found explicitly. Relaxation functions and frequency-dependent susceptibilities are expressed by hypergeometric series, even in the case of nonlinear response. At the critical point, the spectrum becomes continuous and the exponential decay for large times is replaced by power-law behavior. Similarly, it is found that dynamic scaling applies also to the nucleation rate for this model, and the scaling exponent of the nucleation rate is $2\ensuremath{-}\ensuremath{\alpha}+(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}$.
A statistical theory of the localized model of ferromagnetism is developed by considering the effective density matrices for a single spin, and two nearest-neighbor spins, respectively. The Ising Model and the Heisenberg model, and even the case of the anisotropic exchange interaction, are treated within a unified mathematical scheme. The result is reduced to that of the Bethe approximation in the case of the Ising model, and the constant coupling approximation in the case of the Heisenberg model. The method is generalized to calculate the static spin-correlation function. It is obtained in Ornstein-Zernike type. In the nonequilibrium state, the expectation value of the spin decays by the homogeneous damping as well as the spin-diffusion process. This relaxation equation is different from that by Van Hove because of the existence of the homogeneous damping. The diffusion constant does not vanish even at the Curie point TC, contrary to Van Hove's result, although the damping constant vanishes at TC. The latter fact causes the critical slowing down of the regression of fluctuation. The solution of this relaxation equation gives the time-dependent spin-correlation function. The method can be applied to the antiferromagnetic case, too.
The method for calculation of the correlation functions of the Ising-type systems with short-range interaction and with arbitrary value of spin is developed within cluster approximation. For the Ising model (spin Sz=±1S^z=\pm1Sz=±1) the expressions for pair and ternary correlation functions within two-particle approximation in { q }\set{q}{q}-space are obtained for the hypercubic Bravais lattices. In the 1D case the exact expressions for them in the site space is obtained as well. On the basis of the Glauber equation within two-particle cluster approximation the longitudinal dynamical susceptibility χ({ q },E)χ(\set{q},E)χ({q},E) is found. In the 1D case and in the absence of external field the expression for χ({ q },E)χ(\set{q},E)χ({q},E) is exact. For the Emery-Blume-Griffiths model (Sz=−2,0,2S^z=-2,0,2Sz=−2,0,2) within two-particle approximation the pair correlation functions are calculated. The four-particle cluster approximation is used for calculation of static susceptibility χ({ q })χ(\set{q})χ({q}) of KD2PO4KD_2PO_4KD2​PO4​ ferroelectrics.
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.
Extending the Monte Carlo method to dynamic critical phenomena we investigated the time-dependent correlation functions in the two-dimensional one-spin-flip Ising model and the critical behavior of the associated relaxation times. These relaxation times are the following: ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}}^{\ensuremath{\Delta}T}$, characterizing the approach of the order parameter to equilibrium after a change of temperature $\ensuremath{\Delta}T$ of the system; ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}$ and ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}^{A}$ characterizing the slowing down of the order-parameter correlation and autocorrelation functions, respectively; ${\ensuremath{\tau}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}$ and ${\ensuremath{\tau}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}^{A}$, characterizing the slowing down of the energy correlation and autocorrelation functions; and finally ${\ensuremath{\tau}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}$, characterizing the cross-correlation function. We give estimates for the associated exponents ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}}^{\ensuremath{\Delta}T}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}\ensuremath{\approx}{\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}\ensuremath{\approx}1.90\ifmmode\pm\else\textpm\fi{}0.10$, and ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\ensuremath{\mu}}^{A}\ensuremath{\approx}1.60\ifmmode\pm\else\textpm\fi{}0.10$, ${\ensuremath{\Delta}}_{\ensuremath{\delta}\ensuremath{\mu}\ensuremath{\delta}\mathcal{H}}^{A}\ensuremath{\approx}0.95\ifmmode\pm\else\textpm\fi{}0.10$, ${\ensuremath{\Delta}}_{\ensuremath{\delta}\mathcal{H}\ensuremath{\delta}\mathcal{H}}^{A}\ensuremath{\approx}0$, which are consistent with the dynamic scaling hypothesis and with exact inequalities. A detailed comparison with recent high-temperature-expansion estimates is performed, and the reliability of the Monte Carlo results is carefully analyzed.
A simple dynamical model of interacting Ising spins is discussed near the critical point. It is shown that all the moments are positive and never vanish even at the critical point. Lower and upper bounds of the susceptibility are obtained, which yield the non-vanishing of the real part of the susceptibility. This gives a simple proof for the poly-dispersive nature of our system. The method of high temperature expansions for the dynamical susceptibility is presented and is applied to the two-dimensional Ising model up to the ninth order, the results of which show that the critical index of slowing down (the relaxation time τ (1) ∝( T - T c ) - Δ ; Δ =2±0.05) is different from that of the static susceptibility.