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. B, Solid state

Calculation of dynamic critical properties from a cluster-reaction theory

Kurt Binder, Dietrich Stauffer, Heiner Muller-Krumbhaar

Published Nov 1, 1974• 43 citations
Publisher

Abstract

Critical fluctuations in kinetic Ising models are interpreted in terms of cluster reactions. The basic assumption, that clusters with $l$ spins grow at a rate $\ensuremath{\propto}{l}^{r}$, is tested by Monte Carlo computations in the single-spin-flip case. The dynamic susceptibilities associated with order parameter and energy are then calculated also for nonzero magnetic field, and are shown to fulfill dynamic scaling. The exponent $(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}$ of the relaxation times can be different from the susceptibility exponent $\ensuremath{\gamma}$.

DOI
10.1103/physrevb.10.3853

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

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

Finite-Size Scaling Approach to the Kinetic Ising Model

H. Takano

The kinetic Ising model on a square lattice is studied by means of the finite-size scaling with the help of the Monte Carlo method. The results obtained are consistent with the universality hypothesis that the value of the dynamic critical expotent is universal for certain types of transition probabilities. The present results also suggest that the dynamic critical exponents for the relaxation times of the magnetization and energy are the same. Our estimation of the dynamic critical exponent z gives that z \cong2.2 ±0.1.

Progress of Theoretical Physics•Aug 1, 1982•30 citations

Critical Ising Spin Dynamics near the Percolation Threshold

R. Rammal, A. D. Benoit

A low-temperature approach for the critical Ising spin dynamics on percolation clusters is proposed. The average energy barrier to overturning a cluster of $s$ spins is shown to scale as $\mathrm{ln}s$. This basic result is supported by exact numerical calculations of the energy-barrier distribution for two-and three-dimensional percolation clusters, and is argued to be a consequence of scale invariance. The resulting temperature as well as size dependence of the relaxation-time scales agrees very well with the results of our Monte Carlo simulations.

Physical Review Letters•Aug 12, 1985•49 citations

Critical slowing down in the two-dimensional one-spin-flip Ising model

T. Schneider, E. Stoll

We extend a previous Monte Carlo investigation to include the relaxation times associated with the first time moment of the relaxation function and of the relaxation function squared. Moreover, we introduce, in addition to the order-parameter and energy-relaxation functions, that of the order parameter squared. The asymptotic temperature dependence of the associated relaxation times is found to support the validity of rather extended dynamic scaling hypotheses, predicting that all the associated exponents of slowing down are equal and symmetric with respect to ${T}_{c}$. We also derive various exact lower bounds and inequalities for these relaxation times and exponents. Our numerical data are found to be consistent with these exact results.

Physical review. B, Solid state•Aug 1, 1974•11 citations

Nonequilibrium phase transition in the kinetic Ising model: Violation of Fluctuation-Dissipation relation

Muktish Acharyya

The nonequilibrium dynamic phase transition, in the kinetic Ising model in presence of an oscillating magnetic field, has been studied by Monte Carlo simulation. The fluctuation of dynamic order parameter has been studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined `susceptibility' has also been studied near the dynamic transition point. Similarly, fluctuation of energy and appropriately defined `specific-heat' have been studied as a function of temperature near the dynamic transition point. In both the cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (as power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).

arXiv:cond-mat/9701180•Jan 24, 1997•0 citations
Critical fluctuations in kinetic Ising models are interpreted in terms of cluster reactions. The basic assumption, that clusters with lll spins grow at a rate \ensuremath∝lr\ensuremath{\propto}{l}^{r}\ensuremath∝lr, is tested by Monte Carlo computations in the single-spin-flip case. The dynamic susceptibilities associated with order parameter and energy are then calculated also for nonzero magnetic field, and are shown to fulfill dynamic scaling. The exponent (2\ensuremath−r)\ensuremathβ\ensuremathδ(2\ensuremath{-}r)\ensuremath{\beta}\ensuremath{\delta}(2\ensuremath−r)\ensuremathβ\ensuremathδ of the relaxation times can be different from the susceptibility exponent \ensuremathγ\ensuremath{\gamma}\ensuremathγ.
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}$.
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 kinetic Ising model on a square lattice is studied by means of the finite-size scaling with the help of the Monte Carlo method. The results obtained are consistent with the universality hypothesis that the value of the dynamic critical expotent is universal for certain types of transition probabilities. The present results also suggest that the dynamic critical exponents for the relaxation times of the magnetization and energy are the same. Our estimation of the dynamic critical exponent z gives that z \cong2.2 ±0.1.
A low-temperature approach for the critical Ising spin dynamics on percolation clusters is proposed. The average energy barrier to overturning a cluster of sss spins is shown to scale as lns\mathrm{ln}slns. This basic result is supported by exact numerical calculations of the energy-barrier distribution for two-and three-dimensional percolation clusters, and is argued to be a consequence of scale invariance. The resulting temperature as well as size dependence of the relaxation-time scales agrees very well with the results of our Monte Carlo simulations.
We extend a previous Monte Carlo investigation to include the relaxation times associated with the first time moment of the relaxation function and of the relaxation function squared. Moreover, we introduce, in addition to the order-parameter and energy-relaxation functions, that of the order parameter squared. The asymptotic temperature dependence of the associated relaxation times is found to support the validity of rather extended dynamic scaling hypotheses, predicting that all the associated exponents of slowing down are equal and symmetric with respect to Tc{T}_{c}Tc​. We also derive various exact lower bounds and inequalities for these relaxation times and exponents. Our numerical data are found to be consistent with these exact results.
The nonequilibrium dynamic phase transition, in the kinetic Ising model in presence of an oscillating magnetic field, has been studied by Monte Carlo simulation. The fluctuation of dynamic order parameter has been studied as a function of temperature near the dynamic transition point. The temperature variation of appropriately defined `susceptibility' has also been studied near the dynamic transition point. Similarly, fluctuation of energy and appropriately defined `specific-heat' have been studied as a function of temperature near the dynamic transition point. In both the cases, the fluctuations (of dynamic order parameter and energy) and the corresponding responses diverge (as power law fashion) near the dynamic transition point with similar critical behavior (with identical exponent values).