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

arXiv:cond-mat/9708063 · Physical Review Letters

Universal Dynamics of Independent Critical Relaxation Modes

M. P. Nightingale, Henk W. J. Blöte

Published Feb 2, 1998• 13 citations
arXivHTMLPDF

Abstract

We obtain the relaxation times of several, progressively rapid, independent modes of three models in a two-dimensional Ising universality class. Their size dependence can be described by one single dynamic exponent and universal amplitude ratios. This analysis is based on variational approximations of the eigenstates of the Markov matrix describing heat-bath, single-spin-flip dynamics. Monte Carlo computation of the corresponding autocorrelations and cross correlations, in which the variational error is systematically reduced, yields eigenvalues and the associated relaxation times with considerably higher statistical accuracy than is the case for traditional correlations.

DOI
10.1103/physrevlett.80.1007

Related papers

Monte Carlo computation of correlation times of independent relaxation modes at criticality

M. P. Nightingale, Henk W. J. Blöte

We investigate aspects of the universality of Glauber critical dynamics in two dimensions. We compute the critical exponent z and numerically corroborate its universality for three different models in the static Ising universality class and for five independent relaxation modes. We also present evidence for universality of amplitude ratios, which shows that, as far as dynamic behavior is concerned, each model in a given universality class is characterized by a single nonuniversal metric factor which determines the overall time scale. This paper also discusses in detail the variational and projection methods that are used to compute relaxation times with high accuracy.

arXiv:cond-mat/0001251•Physical review. B, Condensed matter•Jul 1, 2000•88 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

Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model

Martin Hasenbusch B

We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.

arXiv:1908.01702•Physical review. E•Feb 24, 2020•39 citations

Relaxational dynamics in 3D randomly diluted Ising models

Martin Hasenbusch B, Andrea Pelissetto C, Ettore Vicari

We study the purely relaxational dynamics (model A) at criticality in three-dimensional disordered Ising systems whose static critical behaviour belongs to the randomly diluted Ising universality class. We consider the site-diluted and bond-diluted Ising models, and the ± J Ising model along the paramagnetic–ferromagnetic transition line. We perform Monte Carlo simulations at the critical point using the Metropolis algorithm and study the dynamic behaviour in equilibrium at various values of the disorder parameter. The results provide a robust evidence of the existence of a unique model-A dynamic universality class which describes the relaxational critical dynamics in all considered models. In particular, the analysis of the size dependence of suitably defined autocorrelation times at the critical point provides the estimate z = 2.35(2) for the universal dynamic critical exponent. We also study the off-equilibrium relaxational dynamics following a quench from to T c . In agreement with the field-theory scenario, the analysis of the off-equilibrium dynamic critical behaviour gives an estimate of z that is perfectly consistent with the equilibrium estimate z = 2.35(2).

arXiv:0709.4179•Journal of Statistical Mechanics Theory and Experiment•Nov 23, 2007•23 citations

Dynamic Exponent of the Two-Dimensional Ising Model and Monte Carlo Computation of the Subdominant Eigenvalue of the Stochastic Matrix

M. P. Nightingale, Henk W. J. Blöte

We introduce a novel variance-reducing Monte Carlo algorithm for accurate determination of correlation times. We apply this method to two-dimensional Ising systems with sizes up to $15\ifmmode\times\else\texttimes\fi{}15$, using single-spin flip dynamics, random site selection, and transition probabilities according to the heat-bath method. From a finite-size scaling analysis of these correlation times, the dynamic critical exponent $z$ is determined as $z\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}2.1665(12)$.

arXiv:cond-mat/9601059•Physical Review Letters•Jun 10, 1996•143 citations

Critical nonequilibrium cluster-flip relaxations in Ising models

Yusuke Tomita, Yoshihiko Nonomura

We investigate nonequilibrium relaxations of Ising models at the critical point by using a cluster update. While preceding studies imply that nonequilibrium cluster-flip dynamics at the critical point are universally described by the stretched-exponential function, we find that the dynamics changes from the stretched exponential to the power function as the dimensionality is increased: The two-, three-, four-, and infinite-dimensional Ising models are numerically studied, and the four- and infinite-dimensional Ising models exhibit the power-law relaxation. We also show that the finite-size scaling analysis using the normalized correlation length is markedly effective for the analysis of relaxational processes rather than the direct use of the Monte Carlo step.

arXiv:1808.07761•Physical review. E•Nov 12, 2018•5 citations
We obtain the relaxation times of several, progressively rapid, independent modes of three models in a two-dimensional Ising universality class. Their size dependence can be described by one single dynamic exponent and universal amplitude ratios. This analysis is based on variational approximations of the eigenstates of the Markov matrix describing heat-bath, single-spin-flip dynamics. Monte Carlo computation of the corresponding autocorrelations and cross correlations, in which the variational error is systematically reduced, yields eigenvalues and the associated relaxation times with considerably higher statistical accuracy than is the case for traditional correlations.
We investigate aspects of the universality of Glauber critical dynamics in two dimensions. We compute the critical exponent z and numerically corroborate its universality for three different models in the static Ising universality class and for five independent relaxation modes. We also present evidence for universality of amplitude ratios, which shows that, as far as dynamic behavior is concerned, each model in a given universality class is characterized by a single nonuniversal metric factor which determines the overall time scale. This paper also discusses in detail the variational and projection methods that are used to compute relaxation times with high accuracy.
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.
We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.
We study the purely relaxational dynamics (model A) at criticality in three-dimensional disordered Ising systems whose static critical behaviour belongs to the randomly diluted Ising universality class. We consider the site-diluted and bond-diluted Ising models, and the ± J Ising model along the paramagnetic–ferromagnetic transition line. We perform Monte Carlo simulations at the critical point using the Metropolis algorithm and study the dynamic behaviour in equilibrium at various values of the disorder parameter. The results provide a robust evidence of the existence of a unique model-A dynamic universality class which describes the relaxational critical dynamics in all considered models. In particular, the analysis of the size dependence of suitably defined autocorrelation times at the critical point provides the estimate z = 2.35(2) for the universal dynamic critical exponent. We also study the off-equilibrium relaxational dynamics following a quench from to T c . In agreement with the field-theory scenario, the analysis of the off-equilibrium dynamic critical behaviour gives an estimate of z that is perfectly consistent with the equilibrium estimate z = 2.35(2).
We introduce a novel variance-reducing Monte Carlo algorithm for accurate determination of correlation times. We apply this method to two-dimensional Ising systems with sizes up to $15\ifmmode\times\else\texttimes\fi{}15$, using single-spin flip dynamics, random site selection, and transition probabilities according to the heat-bath method. From a finite-size scaling analysis of these correlation times, the dynamic critical exponent zzz is determined as z=2.1665(12)z\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}2.1665(12)z=2.1665(12).
We investigate nonequilibrium relaxations of Ising models at the critical point by using a cluster update. While preceding studies imply that nonequilibrium cluster-flip dynamics at the critical point are universally described by the stretched-exponential function, we find that the dynamics changes from the stretched exponential to the power function as the dimensionality is increased: The two-, three-, four-, and infinite-dimensional Ising models are numerically studied, and the four- and infinite-dimensional Ising models exhibit the power-law relaxation. We also show that the finite-size scaling analysis using the normalized correlation length is markedly effective for the analysis of relaxational processes rather than the direct use of the Monte Carlo step.