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

Off-equilibrium computation of the dynamic critical exponent of the three-dimensional Heisenberg model

Antonio Astillero, Juan Jesus Ruiz-Lorenzo

Published Jun 11, 2019• 0 citations
OpenAlex

Abstract

Working in the out of equilibrium regime and using state-of-the-art techniques we have computed the dynamic critical exponent of the three dimensional Heisenberg model. We have run very large lattices ($L\le 250$) in CPUs and GPUs obtaining $z=2.041(16)$ from the growth of the correlation length and $z=2.034(22)$ for the decay of the energy. We compare our values with that previously computed at equilibrium with relatively small lattices ($L\le 24$), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, $\eta$ and $\nu$, in this out of equilibrium regime.

Related papers

Computation of the dynamic critical exponent of the three-dimensional Heisenberg model

Antonio Astillero, Juan Jesus Ruiz-Lorenzo

Working in and out of equilibrium and using state-of-the-art techniques we have computed the dynamic critical exponent of the three-dimensional Heisenberg model. By computing the integrated autocorrelation time at equilibrium, for lattice sizes L≤64, we have obtained z=2.033(5). In the out-of-equilibrium regime we have run very large lattices (L≤250) obtaining z=2.04(2) from the growth of the correlation length. We compare our values with that previously computed at equilibrium with relatively small lattices (L≤24), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, η and ν, in the out-of-equilibrium regime.

arXiv:1906.04518•Physical review. E•Dec 16, 2019•8 citations

Equilibrium and out-of-equilibrium critical dynamics of the three-dimensional Heisenberg model with random cubic anisotropy

Antonio Astillero, Juan Jesus Ruiz-Lorenzo

We study the critical dynamics of the three-dimensional Heisenberg model with random cubic anisotropy in the out-of-equilibrium and equilibrium regimes. Analytical approaches based on field theory predict that the universality class of this model is that of the three-dimensional site-diluted Ising model. We have been able to estimate the dynamic critical exponent by working in the equilibrium regime and by computing the integrated autocorrelation times obtaining z=2.50(5) (without taking into account scaling corrections) and z=2.29(11) (by fixing the scaling corrections to that predicted by field theory). In the out-of-equilibrium regime we have focused on the study of the dynamic correlation length, allowing us to compute the dynamic critical exponent obtaining z=2.38(2), which is compatible with the equilibrium ones. Finally, both estimates are also compatible with the most accurate prediction z=2.35(2), from numerical simulations of the 3D site-diluted Ising model, in agreement with the predictions based on field theory.

arXiv:2410.14275•Physical review. E•Jun 5, 2025•1 citation

Critical exponents and equation of state of the three-dimensional Heisenberg universality class

Massimo Campostrini B, Martin Hasenbusch B, Andrea Pelissetto C +2 more

We improve the theoretical estimates of the critical exponents for the three-dimensional Heisenberg universality class. We find $\ensuremath{\gamma}=1.3960(9),$ $\ensuremath{\nu}=0.7112(5),$ $\ensuremath{\eta}=0.0375(5),$ $\ensuremath{\alpha}=\ensuremath{-}0.1336(15),$ $\ensuremath{\beta}=0.3689(3),$ and $\ensuremath{\delta}=4.783(3).$ We consider an improved lattice ${\ensuremath{\varphi}}^{4}$ Hamiltonian with suppressed leading scaling corrections. Our results are obtained by combining Monte Carlo simulations based on finite-size scaling methods and high-temperature expansions. The critical exponents are computed from high-temperature expansions specialized to the ${\ensuremath{\varphi}}^{4}$ improved model. By the same technique we determine the coefficients of the small-magnetization expansion of the equation of state. This expansion is extended analytically by means of approximate parametric representations, obtaining the equation of state in the whole critical region. We also determine a number of universal amplitude ratios.

arXiv:cond-mat/0110336•Physical review. B, Condensed matter•Apr 2, 2002•482 citations

Critical exponents of the classical three-dimensional Heisenberg model: A single-cluster Monte Carlo study

Christian Holm, Wolfhard Janke E

We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to ${100}^{3}$ spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.

arXiv:hep-lat/9301002•Physical review. B, Condensed matter•Jul 1, 1993•232 citations

Short-time critical and coarsening dynamics of the classical three-dimensional Heisenberg model

Henrique A. Fernandes, Roberto da Silva, J. R. Drugowich de Felício

The dynamic critical behaviour of the three-dimensional classical Heisenberg model was studied by using short-time Monte Carlo simulations. The dynamic critical exponents z , θ and θ g , as well as the static critical exponents ν and β, were estimated by following the time evolution of the z -component of the total magnetization and its moments of higher order. The anomalous dimension of the magnetization x 0 was also calculated from a relation obtained through finite size scaling in the non-equilibrium stage. Afterwards, the local persistence exponent θ l was obtained by following the power law behaviour of the local persistence probability when the system was quenched from high temperatures to the null one ( T = 0). The time evolution of P l ( t ) was also analysed for T >0 and the expected exponential stretch was observed.

Journal of Statistical Mechanics Theory and Experiment•Oct 5, 2006•17 citations

Critical Exponents from Other Expansions

It is useful to compare the results obtained so far with other approaches to the critical exponents. One is a similar field-theoretic approach based on perturbation expansions of φ 4-theories. But instead of working in D = 4−ε dimensions and continuing the results to ε = 1 to obtain critical exponents in the physical dimension D = 3, perturbation expansions can be derived directly in three dimensions, and studied in the regime of small masses. Fundamentally, this approach is less satisfactory than the previous one since there is no small parameter analogous to ε which permits us to prove the existence of a scaling limit at least for small ε. Another disadvantage is that it is impossible to calculate the expansion coefficients analytically, except for one- and two loop-diagrams. In one respect the three-dimensional approach has, however, an advantage: the power series have long been available up to six loops for ω, η, and ν, and there exist recent results up to seven loops for the critical exponents η and ν. Although numerical expansions are not as esthetic as the exact expansions derived in this book, their increased order leads to a higher accuracy in the critical exponents, as we shall see below in Section 20.4. Another approach is based on high-temperature expansions for lattice models of the O(N)symmetric classical Heisenberg model with the energy of Eq. (1.52). The critical exponents will be derived from such expansions in Section 20.8. Finally we take advantage of the fact that the O(N)-symmetric classical Heisenberg model is exactly solvable in the limit of large N, where it reduces to the spherical model [recall the remark after Eq. (1.69)]. It is therefore possible to expand the φ 4-theory around this limit in powers of 1/N. The results will be compared with those from the perturbation expansions in

WORLD SCIENTIFIC eBooks•Jul 1, 2001•0 citations
Working in the out of equilibrium regime and using state-of-the-art techniques we have computed the dynamic critical exponent of the three dimensional Heisenberg model. We have run very large lattices (L≤250L\le 250L≤250) in CPUs and GPUs obtaining z=2.041(16)z=2.041(16)z=2.041(16) from the growth of the correlation length and z=2.034(22)z=2.034(22)z=2.034(22) for the decay of the energy. We compare our values with that previously computed at equilibrium with relatively small lattices (L≤24L\le 24L≤24), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, η\etaη and ν\nuν, in this out of equilibrium regime.
Working in and out of equilibrium and using state-of-the-art techniques we have computed the dynamic critical exponent of the three-dimensional Heisenberg model. By computing the integrated autocorrelation time at equilibrium, for lattice sizes L≤64, we have obtained z=2.033(5). In the out-of-equilibrium regime we have run very large lattices (L≤250) obtaining z=2.04(2) from the growth of the correlation length. We compare our values with that previously computed at equilibrium with relatively small lattices (L≤24), with that provided by means a three-loops calculation using perturbation theory and with experiments. Finally we have checked previous estimates of the static critical exponents, η and ν, in the out-of-equilibrium regime.
We study the critical dynamics of the three-dimensional Heisenberg model with random cubic anisotropy in the out-of-equilibrium and equilibrium regimes. Analytical approaches based on field theory predict that the universality class of this model is that of the three-dimensional site-diluted Ising model. We have been able to estimate the dynamic critical exponent by working in the equilibrium regime and by computing the integrated autocorrelation times obtaining z=2.50(5) (without taking into account scaling corrections) and z=2.29(11) (by fixing the scaling corrections to that predicted by field theory). In the out-of-equilibrium regime we have focused on the study of the dynamic correlation length, allowing us to compute the dynamic critical exponent obtaining z=2.38(2), which is compatible with the equilibrium ones. Finally, both estimates are also compatible with the most accurate prediction z=2.35(2), from numerical simulations of the 3D site-diluted Ising model, in agreement with the predictions based on field theory.
We improve the theoretical estimates of the critical exponents for the three-dimensional Heisenberg universality class. We find $\ensuremath{\gamma}=1.3960(9),$ $\ensuremath{\nu}=0.7112(5),$ $\ensuremath{\eta}=0.0375(5),$ $\ensuremath{\alpha}=\ensuremath{-}0.1336(15),$ $\ensuremath{\beta}=0.3689(3),$ and $\ensuremath{\delta}=4.783(3).$ We consider an improved lattice ${\ensuremath{\varphi}}^{4}$ Hamiltonian with suppressed leading scaling corrections. Our results are obtained by combining Monte Carlo simulations based on finite-size scaling methods and high-temperature expansions. The critical exponents are computed from high-temperature expansions specialized to the ${\ensuremath{\varphi}}^{4}$ improved model. By the same technique we determine the coefficients of the small-magnetization expansion of the equation of state. This expansion is extended analytically by means of approximate parametric representations, obtaining the equation of state in the whole critical region. We also determine a number of universal amplitude ratios.
We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to 1003{100}^{3}1003 spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.
The dynamic critical behaviour of the three-dimensional classical Heisenberg model was studied by using short-time Monte Carlo simulations. The dynamic critical exponents z , θ and θ g , as well as the static critical exponents ν and β, were estimated by following the time evolution of the z -component of the total magnetization and its moments of higher order. The anomalous dimension of the magnetization x 0 was also calculated from a relation obtained through finite size scaling in the non-equilibrium stage. Afterwards, the local persistence exponent θ l was obtained by following the power law behaviour of the local persistence probability when the system was quenched from high temperatures to the null one ( T = 0). The time evolution of P l ( t ) was also analysed for T >0 and the expected exponential stretch was observed.
It is useful to compare the results obtained so far with other approaches to the critical exponents. One is a similar field-theoretic approach based on perturbation expansions of φ 4-theories. But instead of working in D = 4−ε dimensions and continuing the results to ε = 1 to obtain critical exponents in the physical dimension D = 3, perturbation expansions can be derived directly in three dimensions, and studied in the regime of small masses. Fundamentally, this approach is less satisfactory than the previous one since there is no small parameter analogous to ε which permits us to prove the existence of a scaling limit at least for small ε. Another disadvantage is that it is impossible to calculate the expansion coefficients analytically, except for one- and two loop-diagrams. In one respect the three-dimensional approach has, however, an advantage: the power series have long been available up to six loops for ω, η, and ν, and there exist recent results up to seven loops for the critical exponents η and ν. Although numerical expansions are not as esthetic as the exact expansions derived in this book, their increased order leads to a higher accuracy in the critical exponents, as we shall see below in Section 20.4. Another approach is based on high-temperature expansions for lattice models of the O(N)symmetric classical Heisenberg model with the energy of Eq. (1.52). The critical exponents will be derived from such expansions in Section 20.8. Finally we take advantage of the fact that the O(N)-symmetric classical Heisenberg model is exactly solvable in the limit of large N, where it reduces to the spherical model [recall the remark after Eq. (1.69)]. It is therefore possible to expand the φ 4-theory around this limit in powers of 1/N. The results will be compared with those from the perturbation expansions in