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. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Finite-size scaling of the Glauber model of critical dynamics

James H. Luscombe, Marshall Luban, Joseph P. Reynolds

Published Jun 1, 1996• 49 citations
Publisher

Abstract

We obtain the exact critical relaxation time ${\mathrm{\ensuremath{\tau}}}_{\mathit{L}}$(\ensuremath{\xi}), where \ensuremath{\xi} is the bulk correlation length, for the Glauber kinetic Ising model of spins on a one-dimensional lattice of finite length L for both periodic and free boundary conditions (BC's). We show that, independent of the BC's, the dynamic critical exponent has the well-known value z=2, and we comment on a recent claim that z=1 for this model. The ratio ${\mathrm{\ensuremath{\tau}}}_{\mathit{L}}$(\ensuremath{\xi})/${\mathrm{\ensuremath{\tau}}}_{\mathrm{\ensuremath{\infty}}}$(\ensuremath{\xi}), in the double limit L,\ensuremath{\xi}\ensuremath{\rightarrow}\ensuremath{\infty} for fixed x=L/\ensuremath{\xi}, approaches a limiting functional form, ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$(L/\ensuremath{\xi}), the finite-size scaling function. For free BC's we derive the exact scaling function ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$(x)=[1+(\ensuremath{\omega}(x)/x${)}^{2}$${]}^{\mathrm{\ensuremath{-}}1}$, where \ensuremath{\omega}(x) is the smallest root of the transcendental equation \ensuremath{\omega}tan(\ensuremath{\omega}/2)=x. We provide expansions of \ensuremath{\omega}(x) in powers of x and ${\mathit{x}}^{\mathrm{\ensuremath{-}}1}$ for the regimes of small and large x, respectively, and establish their radii of convergence. The scaling function shows anomalous behavior at small x, ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$(x)\ensuremath{\approxeq}x, instead of the usual ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$(x)\ensuremath{\approxeq}${\mathit{x}}^{\mathit{z}}$, as x\ensuremath{\rightarrow}0. This is because, even for finite L, the lifetime of the slowest dynamical mode diverges for T\ensuremath{\rightarrow} 0 K. For periodic BC's, with the exception of one system, ${\mathrm{\ensuremath{\tau}}}_{\mathit{L}}$ is independent of L, and hence ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$=1. The exceptional system, that with an odd number of spins and antiferromagnetic couplings, exhibits frustration at T=0 K, and the scaling function is given by ${\mathit{f}}_{\mathrm{\ensuremath{\tau}}}$(x)=[1+(\ensuremath{\pi}/x${)}^{2}$${]}^{\mathrm{\ensuremath{-}}1}$. \textcopyright{} 1996 The American Physical Society.

DOI
10.1103/physreve.53.5852

Related papers

Universal scaling form of the correlation length in Ising strips with periodic, free, fixed, and mixed boundary conditions

Theodore W. Burkhardt, Ihnsouk Guim

According to finite-size scaling theory, the correlation length in an Ising strip with infinite length and a width of L lattice constants has the form ${\ensuremath{\xi}}_{L}$(T)=Lf(c(T-${T}_{c}$)L) in the scaling regime L\ensuremath{\gg}1,\ensuremath{\Vert}T-${T}_{c}$\ensuremath{\Vert}/${T}_{c}$\ensuremath{\ll}1. Here f(x) is a universal scaling function, which does, however, depend on the boundary conditions at the edges of the strip. Utilizing the correspondence between the two-dimensional Ising model and the one-dimensional quantum Ising model in a transverse field, we calculate the explicit form of the function f(x) for spin-spin and energy-energy correlations in strips with periodic, free, fixed, and mixed free-fixed boundary conditions.

Physical review. B, Condensed matter•Feb 1, 1987•23 citations

Critical dynamics of the open Ising chain

J. Kamphorst Leal da Silva, Adriana G. Moreira, M. S. Soares +1 more

The one-dimensional ferromagnetic Ising model with open boundary conditions within the Glauber dynamics is studied. From the behavior of the order parameter, the explicit scaling form of the relaxation time is obtained. For finite systems, the dynamical critical exponent is z=1, in contrast to the value of the infinite open chain (z=2).

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 1, 1995•44 citations

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

Dynamic exponent for 2D Ising model by power spectra method

Kwan-tai Leung

We consider power spectra for fluctuating quantities at a critical point, P( omega ) varies as omega - phi , where phi is determined by the dynamic exponent z and static exponents specific to the quantity considered. From the magnetization spectra obtained with simulations, z is found to be 2.13+or-0.03 for the two-dimensional Glauber kinetic Ising model, hence in accord with recent Monte Carlo relaxational studies, but at odds with those obtained by series expansion and damage spreading methods. Excellent dynamic finite-size scaling for system sizes L=10 to 128 supports our estimate.

Journal of Physics A Mathematical and General•Dec 7, 1993•6 citations

Finite Size Effects in the Critical Dynamics of the Kinetic Ising Model

M. Cemal Yalabık, J. D. Gunton

The time constant which characterizes the critical slowing down of the kinetic Ising model is studied as a function of system size for n×n lattices with periodic boundary conditions, for n from two to ten. In particular, it is shown that the behavior of this time constant at the temperature at which the specific heat exhibits a maximum is consistent with Suzuki's finite size dynamical scaling theory and yields a dynamical critical exponent of z≃2.

Progress of Theoretical Physics•Dec 1, 1979•24 citations

Spin-spin correlations in finite systems: Scaling hypothesis and corrections to bulk behavior

Surjit Singh, Raj K. Pathria

We study the correlation function G(R,T;L) and the correlation length \ensuremath{\xi}(T;L) in a finite spherical model if size ${L}^{d\mathrm{\ensuremath{-}}d\mathcal{'}}$\ifmmode\times\else\texttimes\fi{}${\ensuremath{\infty}}^{d\mathcal{'}}$ under periodic boundary conditions and emphasize the role of the variable L/\ensuremath{\xi}(L), rather than L/\ensuremath{\xi}(\ensuremath{\infty}), as a natural scaling variable of the system throughout the transition region, including temperatures below ${T}_{c}$(\ensuremath{\infty}). We obtain a variety of finite-size effects, some of which may have validity for all O(n) models with n>2.

Physical review. B, Condensed matter•Jan 1, 1986•16 citations
We obtain the exact critical relaxation time \ensuremathτL{\mathrm{\ensuremath{\tau}}}_{\mathit{L}}\ensuremathτL​(\ensuremath{\xi}), where \ensuremath{\xi} is the bulk correlation length, for the Glauber kinetic Ising model of spins on a one-dimensional lattice of finite length L for both periodic and free boundary conditions (BC's). We show that, independent of the BC's, the dynamic critical exponent has the well-known value z=2, and we comment on a recent claim that z=1 for this model. The ratio \ensuremathτL{\mathrm{\ensuremath{\tau}}}_{\mathit{L}}\ensuremathτL​(\ensuremath{\xi})/\ensuremathτ\ensuremath∞{\mathrm{\ensuremath{\tau}}}_{\mathrm{\ensuremath{\infty}}}\ensuremathτ\ensuremath∞​(\ensuremath{\xi}), in the double limit L,\ensuremath{\xi}\ensuremath{\rightarrow}\ensuremath{\infty} for fixed x=L/\ensuremath{\xi}, approaches a limiting functional form, f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​(L/\ensuremath{\xi}), the finite-size scaling function. For free BC's we derive the exact scaling function f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​(x)=[1+(\ensuremath{\omega}(x)/x{)}^{2}$${]}^{\mathrm{\ensuremath{-}}1}, where \ensuremath{\omega}(x) is the smallest root of the transcendental equation \ensuremath{\omega}tan(\ensuremath{\omega}/2)=x. We provide expansions of \ensuremath{\omega}(x) in powers of x and x\ensuremath−1{\mathit{x}}^{\mathrm{\ensuremath{-}}1}x\ensuremath−1 for the regimes of small and large x, respectively, and establish their radii of convergence. The scaling function shows anomalous behavior at small x, f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​(x)\ensuremath{\approxeq}x, instead of the usual f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​(x)\ensuremath{\approxeq}xz{\mathit{x}}^{\mathit{z}}xz, as x\ensuremath{\rightarrow}0. This is because, even for finite L, the lifetime of the slowest dynamical mode diverges for T\ensuremath{\rightarrow} 0 K. For periodic BC's, with the exception of one system, \ensuremathτL{\mathrm{\ensuremath{\tau}}}_{\mathit{L}}\ensuremathτL​ is independent of L, and hence f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​=1. The exceptional system, that with an odd number of spins and antiferromagnetic couplings, exhibits frustration at T=0 K, and the scaling function is given by f\ensuremathτ{\mathit{f}}_{\mathrm{\ensuremath{\tau}}}f\ensuremathτ​(x)=[1+(\ensuremath{\pi}/x{)}^{2}$${]}^{\mathrm{\ensuremath{-}}1}. \textcopyright{} 1996 The American Physical Society.
According to finite-size scaling theory, the correlation length in an Ising strip with infinite length and a width of L lattice constants has the form ${\ensuremath{\xi}}_{L}$(T)=Lf(c(T-Tc{T}_{c}Tc​)L) in the scaling regime L\ensuremath{\gg}1,\ensuremath{\Vert}T-Tc{T}_{c}Tc​\ensuremath{\Vert}/Tc{T}_{c}Tc​\ensuremath{\ll}1. Here f(x) is a universal scaling function, which does, however, depend on the boundary conditions at the edges of the strip. Utilizing the correspondence between the two-dimensional Ising model and the one-dimensional quantum Ising model in a transverse field, we calculate the explicit form of the function f(x) for spin-spin and energy-energy correlations in strips with periodic, free, fixed, and mixed free-fixed boundary conditions.
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}$.
We consider power spectra for fluctuating quantities at a critical point, P( omega ) varies as omega - phi , where phi is determined by the dynamic exponent z and static exponents specific to the quantity considered. From the magnetization spectra obtained with simulations, z is found to be 2.13+or-0.03 for the two-dimensional Glauber kinetic Ising model, hence in accord with recent Monte Carlo relaxational studies, but at odds with those obtained by series expansion and damage spreading methods. Excellent dynamic finite-size scaling for system sizes L=10 to 128 supports our estimate.
We study the correlation function G(R,T;L) and the correlation length \ensuremath{\xi}(T;L) in a finite spherical model if size ${L}^{d\mathrm{\ensuremath{-}}d\mathcal{'}}$\ifmmode\times\else\texttimes\fi{}${\ensuremath{\infty}}^{d\mathcal{'}}$ under periodic boundary conditions and emphasize the role of the variable L/\ensuremath{\xi}(L), rather than L/\ensuremath{\xi}(\ensuremath{\infty}), as a natural scaling variable of the system throughout the transition region, including temperatures below Tc{T}_{c}Tc​(\ensuremath{\infty}). We obtain a variety of finite-size effects, some of which may have validity for all O(n) models with n>2.