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

Journal of Physics A Mathematical and General

Non-linear critical relaxation in the kinetic Ising model

Hiroshi Ikeda

Published Mar 1, 1978• 5 citations
Publisher

Abstract

The authors show, by calculating various linear and non-linear relaxation times exactly, that the scaling law Delta psi (nl) = Delta psi (1) - beta psi is valid for the one-dimensional kinetic Ising model. In the two-dimensional case, the critical exponent of the non-linear relaxation time of the order parameter is estimated from an eleventh-order high-temperature series expansion. The result, Delta M (nl) =2.00+or-0.04, is in agreement with the scaling prediction.

DOI
10.1088/0305-4470/11/3/018

Related papers

Linear and nonlinear critical slowing down in the kinetic Ising model: High-temperature series

Zoltàn Ràcz, Malcolm F. Collins

The difference between the critical exponents (${\ensuremath{\Delta}}^{(\mathrm{l})}$ and ${\ensuremath{\Delta}}^{(\mathrm{nl})}$) of the linear and nonlinear relaxation times of the order parameter (${\ensuremath{\tau}}^{(\mathrm{l})}$ and ${\ensuremath{\tau}}^{(\mathrm{nl})}$) is investigated in the two-dimensional one-spin-flip kinetic Ising model. We have calculated the high-temperature series for ${\ensuremath{\tau}}^{(\mathrm{nl})}$ up to ninth order and made use of the known series for ${\ensuremath{\tau}}^{(\mathrm{l})}$ up to twelfth order. The series are analyzed by the ratio and the Pad\'e-approximant methods. The correlation between the critical-point and critical-exponent estimates in the Pad\'e approximants allows an improvement in the determination of ${\ensuremath{\Delta}}^{(\mathrm{l})}$. The result, ${\ensuremath{\Delta}}^{(\mathrm{l})}=2.125\ifmmode\pm\else\textpm\fi{}0.01$, is higher than previous estimates of 2.0 \ifmmode\pm\else\textpm\fi{} 0.05. The estimate of ${\ensuremath{\Delta}}^{(\mathrm{nl})}$ is less precise but the result ${\ensuremath{\Delta}}^{(\mathrm{nl})}=1.95\ifmmode\pm\else\textpm\fi{}0.15$ leads to the conclusion that ${\ensuremath{\Delta}}^{(\mathrm{l})}\ensuremath{\ne}{\ensuremath{\Delta}}^{(\mathrm{nl})}$ and that the difference between the two is of order $\ensuremath{\beta}$ (critical index of the order parameter). This is in accord with the scaling prediction.

Physical review. B, Solid state•Apr 1, 1976•62 citations

Critical dynamics of the two-dimensional kinetic Ising model: High-temperature series analysis of the autorelaxation time

Jos Rogiers, Joseph O. Indekeu

For the two-dimensional nearest-neighbor kinetic Ising model without conservation laws, we derive on the square lattice the high-temperature series expansion of the autorelaxation time, to the eleventh order. With use of ratio methods and Pad\'e approximants, we obtain for the dynamical critical exponent ${\mathrm{\ensuremath{\Delta}}}_{\mathit{A}}$ the value 2.09\ifmmode\pm\else\textpm\fi{}0.03, which implies the value 2.34\ifmmode\pm\else\textpm\fi{}0.03 for the linear relaxation exponent z. This value is about 10% larger than several other recent estimates. This discrepancy is possibly due to the relatively small number of nonzero coefficients in the series expansion.

Physical review. B, Condensed matter•Apr 1, 1990•22 citations

Non-linear critical slowing down for the time-dependent Ising model

Noam White

Reports the results of a calculation of a high-temperature expansion for the non-linear relaxation time for a square Ising model with zero applied field. The author has also attempted to estimate the corresponding critical exponent.

Journal of Physics C Solid State Physics•Apr 14, 1976•7 citations

Monte Carlo simulation of energy relaxation in the two-dimensional Ising model

Masatoshi Mori

The energy relaxation in the two-dimensional Ising model, with temperature greater than ${\mathit{T}}_{\mathit{c}}$ and in a spin system having a 5120\ifmmode\times\else\texttimes\fi{}5120 lattice, is studied using supercomputer-conducted Monte Carlo simulation. The critical exponents of the linear and nonlinear relaxation times of the internal energy are respectively estimated to be ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$=2.064\ifmmode\pm\else\textpm\fi{}0.026 and ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathrm{nl})}$=1.132\ifmmode\pm\else\textpm\fi{}0.080, being consistent with the Racz scaling law ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$-${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathrm{nl})}$=1-\ensuremath{\alpha}, where \ensuremath{\alpha} is the critical exponent of the specific heat. Based on the author's previous result in which the critical exponent of the linear relaxation for magnetization was determined as ${\mathrm{\ensuremath{\Delta}}}_{\mathit{M}}^{(\mathit{l})}$=2.076\ifmmode\pm\else\textpm\fi{}0.005, it is strongly suggested that ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$=${\mathrm{\ensuremath{\Delta}}}_{\mathit{M}}^{(\mathit{l})}$.

Physical review. B, Condensed matter•May 1, 1993•5 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

Dynamical Critical Exponents of the Four Dimensional Ising Model

Nevzat AKTEKİN

The modified definition of the correlation length for the four dimensional Ising model \xi \propto \epsilon^{-\nu} log^{1/6}\epsilon [1] predicts that the relaxation time for a quantity near the critical temperature T_c should diverge as \tau \propto \epsilon^{-\Delta} log^{\Delta/3}\epsilon where \epsilon, \nu, and \Delta are, respectively, the reduced temperature ( \epsilon=( T_c-T )/T_c, T_c = 6.68 ), the correlation length critical exponent, and a dynamical critical exponent. The \tau-data for the nonlinear relaxation of the order parameter computed on the Creutz cellular automaton [2,3], by using the finite-size lattices of linear dimension L = 12 and 14 as approximations to the infinite lattice for T.

TURKISH JOURNAL OF PHYSICS•Jan 1, 1997•0 citations
The authors show, by calculating various linear and non-linear relaxation times exactly, that the scaling law Delta psi (nl) = Delta psi (1) - beta psi is valid for the one-dimensional kinetic Ising model. In the two-dimensional case, the critical exponent of the non-linear relaxation time of the order parameter is estimated from an eleventh-order high-temperature series expansion. The result, Delta M (nl) =2.00+or-0.04, is in agreement with the scaling prediction.
The difference between the critical exponents (${\ensuremath{\Delta}}^{(\mathrm{l})}$ and ${\ensuremath{\Delta}}^{(\mathrm{nl})}$) of the linear and nonlinear relaxation times of the order parameter (${\ensuremath{\tau}}^{(\mathrm{l})}$ and ${\ensuremath{\tau}}^{(\mathrm{nl})}$) is investigated in the two-dimensional one-spin-flip kinetic Ising model. We have calculated the high-temperature series for ${\ensuremath{\tau}}^{(\mathrm{nl})}$ up to ninth order and made use of the known series for ${\ensuremath{\tau}}^{(\mathrm{l})}$ up to twelfth order. The series are analyzed by the ratio and the Pad\'e-approximant methods. The correlation between the critical-point and critical-exponent estimates in the Pad\'e approximants allows an improvement in the determination of ${\ensuremath{\Delta}}^{(\mathrm{l})}$. The result, ${\ensuremath{\Delta}}^{(\mathrm{l})}=2.125\ifmmode\pm\else\textpm\fi{}0.01$, is higher than previous estimates of 2.0 \ifmmode\pm\else\textpm\fi{} 0.05. The estimate of ${\ensuremath{\Delta}}^{(\mathrm{nl})}$ is less precise but the result ${\ensuremath{\Delta}}^{(\mathrm{nl})}=1.95\ifmmode\pm\else\textpm\fi{}0.15$ leads to the conclusion that ${\ensuremath{\Delta}}^{(\mathrm{l})}\ensuremath{\ne}{\ensuremath{\Delta}}^{(\mathrm{nl})}$ and that the difference between the two is of order $\ensuremath{\beta}$ (critical index of the order parameter). This is in accord with the scaling prediction.
For the two-dimensional nearest-neighbor kinetic Ising model without conservation laws, we derive on the square lattice the high-temperature series expansion of the autorelaxation time, to the eleventh order. With use of ratio methods and Pad\'e approximants, we obtain for the dynamical critical exponent ${\mathrm{\ensuremath{\Delta}}}_{\mathit{A}}$ the value 2.09\ifmmode\pm\else\textpm\fi{}0.03, which implies the value 2.34\ifmmode\pm\else\textpm\fi{}0.03 for the linear relaxation exponent z. This value is about 10% larger than several other recent estimates. This discrepancy is possibly due to the relatively small number of nonzero coefficients in the series expansion.
The energy relaxation in the two-dimensional Ising model, with temperature greater than Tc{\mathit{T}}_{\mathit{c}}Tc​ and in a spin system having a 5120\ifmmode\times\else\texttimes\fi{}5120 lattice, is studied using supercomputer-conducted Monte Carlo simulation. The critical exponents of the linear and nonlinear relaxation times of the internal energy are respectively estimated to be ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$=2.064\ifmmode\pm\else\textpm\fi{}0.026 and ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathrm{nl})}$=1.132\ifmmode\pm\else\textpm\fi{}0.080, being consistent with the Racz scaling law ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$-${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathrm{nl})}$=1-\ensuremath{\alpha}, where \ensuremath{\alpha} is the critical exponent of the specific heat. Based on the author's previous result in which the critical exponent of the linear relaxation for magnetization was determined as ${\mathrm{\ensuremath{\Delta}}}_{\mathit{M}}^{(\mathit{l})}$=2.076\ifmmode\pm\else\textpm\fi{}0.005, it is strongly suggested that ${\mathrm{\ensuremath{\Delta}}}_{\mathit{E}}^{(\mathit{l})}$=${\mathrm{\ensuremath{\Delta}}}_{\mathit{M}}^{(\mathit{l})}$.
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.
The modified definition of the correlation length for the four dimensional Ising model \xi \propto \epsilon^{-\nu} log^{1/6}\epsilon [1] predicts that the relaxation time for a quantity near the critical temperature T_c should diverge as \tau \propto \epsilon^{-\Delta} log^{\Delta/3}\epsilon where \epsilon, \nu, and \Delta are, respectively, the reduced temperature ( \epsilon=( T_c-T )/T_c, T_c = 6.68 ), the correlation length critical exponent, and a dynamical critical exponent. The \tau-data for the nonlinear relaxation of the order parameter computed on the Creutz cellular automaton [2,3], by using the finite-size lattices of linear dimension L = 12 and 14 as approximations to the infinite lattice for T.