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

Le Journal de Physique Colloques

MONTE CARLO STUDIES OF DYNAMIC CRITICAL PHENOMENA

D. P. Landau, Simon Y. Tang, Stephan Wansleben

Published Dec 1, 1988• 10 citations
Publisher

Abstract

We review the current state of knowledge of dynamic critical phenomena in simple magnetic models and present high statistics Monte Carlo data for critical slowing down in Ising and Potts models in two and three dimensions. By fitting the relaxation functions to multi-exponential decays at the critical temperature and applying finite size scaling we are able to extract very accurate estimates for the dynamical critical exponent z. Our results are compared with previous numerical estimates and with the predictions of renormalization group theory.

DOI
10.1051/jphyscol:19888701

Related papers

Monte Carlo studies of phase transitions and critical phenomena

Ibragimkhan K. Kamilov, Akai Kurbanovich Murtazaev, Kh. K. Aliev

The current state of Monte Carlo work on phase transitions and critical phenomena is reviewed. Both classical and quantum Monte Carlo results are discussed with emphasis on statistical lattice-model studies involving the highly accurate calculation of critical exponents. It is shown that finite-size scaling is an effective method whereby not only simple lattice models but also those of complex practical magnetic materials can be treated. For models allowing crossover phenomena, it is shown that both finite-size scaling theory and conventional power-law functions must be used in carrying out a Monte Carlo analysis.

Physics-Uspekhi•Jul 31, 1999•65 citations

Monte Carlo investigation of critical dynamics in the three-dimensional Ising model

Stephan Wansleben, D. P. Landau

We report the results of a Monte Carlo investigation of the (equilibrium) time-displaced correlation functions for the magnetization and energy of a simple cubic Ising model as a function of time, temperature, and lattice size. The simulations were carried out on a CDC CYBER 205 supercomputer employing a high-speed, vectorized multispin coding program and using a total of 5\ifmmode\times\else\texttimes\fi{}${10}^{12}$ Monte Carlo spin-flip trials. We used L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices with periodic boundary conditions and L as large as 96. The short-time and long-time behaviors of the correlation functions are analyzed by fits to a sum of exponential decays, and the critical exponent z for the largest relaxation time is extracted using a finite-size-scaling analysis. Our estimate z=2.04\ifmmode\pm\else\textpm\fi{}0.03 resolves an intriguing contradiction in the literature; it satisfies the theoretical lower bound and is in agreement with the prediction obtained by \ensuremath{\epsilon} expansion. We also consider various small systematic errors that typically occur in the analysis of relaxation functions and show how they can lead to spurious results if sufficient care is not exercised.

Physical review. B, Condensed matter•Mar 1, 1991•129 citations

Critical dynamics of the four-state Potts model

S. Jain

The critical relaxation of the four-state Potts model on L*L lattices with L<or=64 is studied by the Monte Carlo method. The decay of the magnetisation from the nonequilibrium initial state exhibits a crossover from a power-law relaxation for t<<t c (L) to an exponential one for t>>t c (L), where t c (L) is some crossover time. This crossover is discussed within the context of dynamic finite-size scaling theory. As a result, the dynamic critical exponent is estimated to be Z=2.94+or-0.29.

Journal of Physics A Mathematical and General•Jun 21, 1987•1 citation

DYNAMIC CRITICAL EXPONENTS OF THE ISING MODEL WITH MULTISPIN INTERACTIONS

César Simões, J. R. Drugowich de Felício

We revisit the short-time dynamics of 2D Ising model with three spin interactions in one direction and estimate the critical exponents z, θ, β and ν. Taking properly into account the symmetry of the Hamiltonian, we obtain results completely different from those obtained by Wang et al.10 For the dynamic exponent z our result coincides with that of the 4-state Potts model in two dimensions. In addition, results for the static exponents ν and β agree with previous estimates obtained from finite size scaling combined with conformal invariance. Finally, for the new dynamic exponent θ we find a negative and close to zero value, a result also expected for the 4-state Potts model according to Okano et al.

arXiv:cond-mat/0103165•Modern Physics Letters B•Jun 30, 2001•15 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

Critical dynamics of disordered two-dimensional Ising systems: a Monte Carlo study

Vladimir V. Prudnikov, O. N. Markov

The critical relaxation of the magnetization in a two-dimensional Ising model with quenched random non-magnetic impurities has been studied by numerical simulation. A squared lattice of size 400 2 with spin concentrations p=1.0; 0.95, 0,9, 0.85, 0.8, 0.75 and 0.7 was considered. The dynamical critical exponent z was determined by the Monte Carlo method combined with the dynamical renormalization group method. The following values were received for z(p): z(1)=2.24+or-0.07, z(0.95)=2.24+or-0.06, z(0.9)=2.24+or-0.06, z(0.85)=2.38+or-0.05, z(0.8)=2.51+or-0.06, z(0.75)=2.66+or-0.07 and z(0.7)=2.88+or-0.06. A singular dynamic scaling behaviour with z=A' mod ln(p-p c ) mod +B' was found for systems with spin concentrations p<or=0.85 near the percolation threshold p c , where A'=0.56+or-0.07 and B'=1.62+or-0.07 are temperature-independent constants.

Journal of Physics A Mathematical and General•Mar 23, 1995•5 citations
We review the current state of knowledge of dynamic critical phenomena in simple magnetic models and present high statistics Monte Carlo data for critical slowing down in Ising and Potts models in two and three dimensions. By fitting the relaxation functions to multi-exponential decays at the critical temperature and applying finite size scaling we are able to extract very accurate estimates for the dynamical critical exponent z. Our results are compared with previous numerical estimates and with the predictions of renormalization group theory.
The current state of Monte Carlo work on phase transitions and critical phenomena is reviewed. Both classical and quantum Monte Carlo results are discussed with emphasis on statistical lattice-model studies involving the highly accurate calculation of critical exponents. It is shown that finite-size scaling is an effective method whereby not only simple lattice models but also those of complex practical magnetic materials can be treated. For models allowing crossover phenomena, it is shown that both finite-size scaling theory and conventional power-law functions must be used in carrying out a Monte Carlo analysis.
We report the results of a Monte Carlo investigation of the (equilibrium) time-displaced correlation functions for the magnetization and energy of a simple cubic Ising model as a function of time, temperature, and lattice size. The simulations were carried out on a CDC CYBER 205 supercomputer employing a high-speed, vectorized multispin coding program and using a total of 5\ifmmode\times\else\texttimes\fi{}1012{10}^{12}1012 Monte Carlo spin-flip trials. We used L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices with periodic boundary conditions and L as large as 96. The short-time and long-time behaviors of the correlation functions are analyzed by fits to a sum of exponential decays, and the critical exponent z for the largest relaxation time is extracted using a finite-size-scaling analysis. Our estimate z=2.04\ifmmode\pm\else\textpm\fi{}0.03 resolves an intriguing contradiction in the literature; it satisfies the theoretical lower bound and is in agreement with the prediction obtained by \ensuremath{\epsilon} expansion. We also consider various small systematic errors that typically occur in the analysis of relaxation functions and show how they can lead to spurious results if sufficient care is not exercised.
The critical relaxation of the four-state Potts model on L*L lattices with L<or=64 is studied by the Monte Carlo method. The decay of the magnetisation from the nonequilibrium initial state exhibits a crossover from a power-law relaxation for t<<t c (L) to an exponential one for t>>t c (L), where t c (L) is some crossover time. This crossover is discussed within the context of dynamic finite-size scaling theory. As a result, the dynamic critical exponent is estimated to be Z=2.94+or-0.29.
We revisit the short-time dynamics of 2D Ising model with three spin interactions in one direction and estimate the critical exponents z, θ, β and ν. Taking properly into account the symmetry of the Hamiltonian, we obtain results completely different from those obtained by Wang et al.10 For the dynamic exponent z our result coincides with that of the 4-state Potts model in two dimensions. In addition, results for the static exponents ν and β agree with previous estimates obtained from finite size scaling combined with conformal invariance. Finally, for the new dynamic exponent θ we find a negative and close to zero value, a result also expected for the 4-state Potts model according to Okano et al.
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 critical relaxation of the magnetization in a two-dimensional Ising model with quenched random non-magnetic impurities has been studied by numerical simulation. A squared lattice of size 400 2 with spin concentrations p=1.0; 0.95, 0,9, 0.85, 0.8, 0.75 and 0.7 was considered. The dynamical critical exponent z was determined by the Monte Carlo method combined with the dynamical renormalization group method. The following values were received for z(p): z(1)=2.24+or-0.07, z(0.95)=2.24+or-0.06, z(0.9)=2.24+or-0.06, z(0.85)=2.38+or-0.05, z(0.8)=2.51+or-0.06, z(0.75)=2.66+or-0.07 and z(0.7)=2.88+or-0.06. A singular dynamic scaling behaviour with z=A' mod ln(p-p c ) mod +B' was found for systems with spin concentrations p<or=0.85 near the percolation threshold p c , where A'=0.56+or-0.07 and B'=1.62+or-0.07 are temperature-independent constants.