© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:cond-mat/0603038 · Journal of Statistical Mechanics Theory and Experiment

Dynamic critical exponents for Swendsen–Wang and Wolff algorithms obtained by a nonequilibrium relaxation method

Jianqing Du, Bo Zheng, Jian‐Sheng Wang

Published May 11, 2006• 14 citations

Abstract

Using a nonequilibrium relaxation method, we calculate the dynamic critical exponent z of the two-dimensional Ising model for the Swendsen–Wang and Wolff algorithms. We examine dynamic relaxation processes following a quench from a disordered or an ordered initial state to the critical temperature T c , and measure the exponential relaxation time of the system energy. For the Swendsen–Wang algorithm with an ordered or a disordered initial state, and for the Wolff algorithm with an ordered initial state, the exponential relaxation time fits well to a logarithmic size dependence up to a lattice size L = 8192. For the Wolff algorithm with a disordered initial state, we obtain an effective dynamic exponent z exp = 1.19(2) up to L = 2048. For comparison, we also compute the effective dynamic exponents through the integrated correlation times. In addition, an exact result of the Swendsen–Wang dynamic spectrum of a one-dimensional Ising chain is derived.

DOI: 10.1088/1742-5468/2006/05/p05004

Related papers

Critical dynamics of cluster algorithms in the random-bond Ising model

Ulvi Kanbur, Zeynep Demir Vatansever

In the present work, we present an extensive Monte Carlo simulation study on the dynamical properties of the two-dimensional random-bond Ising model. The correlation time τ of the Swendsen-Wang and Wolff cluster algorithms is calculated at the critical point. The dynamic critical exponent z of both algorithms is also measured by using the numerical data for several lattice sizes up to L=512. It is found for both algorithms that the autocorrelation time decreases considerably and the critical slowing-down effect reduces upon the introduction of bond disorder. Additionally, simulations with the Metropolis algorithm are performed, and the critical slowing-down effect is observed to be more pronounced in the presence of disorder, confirming the previous findings in the literature. Moreover, the existence of the non-self-averaging property of the model is demonstrated by calculating the scaled form of the standard deviation of autocorrelation times. Finally, the critical exponent ratio of the magnetic susceptibility is estimated by using the average cluster size of the Wolff algorithm.

Physical review. E•Feb 28, 2024•2 citations

Nonequilibrium behaviors of the three-dimensional Heisenberg model in the Swendsen-Wang algorithm

Yoshihiko Nonomura, Yusuke Tomita

Recently, it was shown [Y. Nonomura, J. Phys. Soc. Jpn. 83, 113001 (2014)JUPSAU0031-901510.7566/JPSJ.83.113001] that the nonequilibrium critical relaxation of the two-dimensional (2D) Ising model from a perfectly ordered state in the Wolff algorithm is described by stretched-exponential decay, and a universal scaling scheme was found to connect nonequilibrium and equilibrium behaviors. In the present study we extend these findings to vector spin models, and the 3D Heisenberg model could be a typical example. To evaluate the critical temperature and critical exponents precisely using the above scaling scheme, we calculate nonequilibrium ordering from the perfectly disordered state in the Swendsen-Wang algorithm, and we find that the critical ordering process is described by stretched-exponential growth with a comparable exponent to that of the 3D XY model. The critical exponents evaluated in the present study are consistent with those in previous studies.

arXiv:1508.05218•Physical review. E•Jan 4, 2016•8 citations

Empirical relations between static and dynamic exponents for Ising model cluster algorithms

Paul David Coddington, CLIVE F. BAILLIE

We have measured the autocorrelations for the Swendsen-Wang and the Wolff cluster update algorithms for the Ising model in two, three, and four dimensions. The data for the Wolff algorithm suggest that the autocorrelations are linearly related to the specific heat, in which case the dynamic critical exponent is ${\mathit{z}}_{\mathrm{i}\mathrm{n}\mathrm{t},}$${\mathit{E}}^{\mathrm{W}}$=\ensuremath{\alpha}/\ensuremath{\nu}. For the Swendsen-Wang algorithm, scaling the autocorrelations by the average maximum cluster size gives either a constant or a logarithm, which implies that ${\mathit{z}}_{\mathrm{i}\mathrm{n}\mathrm{t},}$${\mathit{E}}^{\mathrm{SW}}$=\ensuremath{\beta}/\ensuremath{\nu} for the Ising model.

Physical Review Letters•Feb 17, 1992•45 citations

Dynamic exponents for potts model cluster algorithms

Paul David Coddington, CLIVE F. BAILLIE

We have studied the Swendsen-Wang and Wolff cluster update algorithms for the Ising model in 2, 3 and 4 dimensions. The data indicate simple relations between the specific heat and the Wolff autocorrelations, and between the magnetization and the Swendsen-Wang autocorrelations. This implies that the dynamic critical exponents are related to the static exponents of the Ising model. We also investigate the possibility of similar relationships for the Q-state Potts model.

Nuclear Physics B - Proceedings Supplements•Jan 1, 1992•1 citation

Finite-size scaling and critical exponents in critical relaxation

Zhibing Li, Lothar Schülke, Bo Zheng

We simulate the critical relaxation process of the two-dimensional Ising model with the initial state both completely disordered or completely ordered. Results of a method to measure both the dynamic and static critical exponents are reported, based on the finite-size scaling for the dynamics at the early time. From the time-dependent Binder cumulant, the dynamical exponent z is extracted independently, while the static exponents \ensuremath{\beta}/\ensuremath{\nu} and \ensuremath{\nu} are obtained from the time evolution of the magnetization and its higher moments. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9508148•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Mar 1, 1996•77 citations

Critical nonequilibrium relaxation in the Swendsen-Wang algorithm in the Berezinsky-Kosterlitz-Thouless and weak first-order phase transitions

Yoshihiko Nonomura, Yusuke Tomita

Recently we showed that the critical nonequilibrium relaxation in the Swendsen-Wang algorithm is widely described by the stretched-exponential relaxation of physical quantities in the Ising or Heisenberg models. Here we make a similar analysis in the Berezinsky-Kosterlitz-Thouless phase transition in the two-dimensional (2D) XY model and in the first-order phase transition in the 2D q=5 Potts model and find that these phase transitions are described by the simple exponential relaxation and power-law relaxation of physical quantities, respectively. We compare the relaxation behaviors of these phase transitions with those of the second-order phase transition in the three- and four-dimensional XY models and in the 2D q-state Potts models for 2≤q≤4 and show that the species of phase transitions can be clearly characterized by the present analysis. We also compare the size dependence of relaxation behaviors of the first-order phase transition in the 2D q=5 and 6 Potts models and propose a quantitative criterion on "weakness" of the first-order phase transition.

arXiv:1509.08352•Physical Review E•Dec 10, 2015•7 citations

Using a nonequilibrium relaxation method, we calculate the dynamic critical exponent z of the two-dimensional Ising model for the Swendsen–Wang and Wolff algorithms. We examine dynamic relaxation processes following a quench from a disordered or an ordered initial state to the critical temperature T c , and measure the exponential relaxation time of the system energy. For the Swendsen–Wang algorithm with an ordered or a disordered initial state, and for the Wolff algorithm with an ordered initial state, the exponential relaxation time fits well to a logarithmic size dependence up to a lattice size L = 8192. For the Wolff algorithm with a disordered initial state, we obtain an effective dynamic exponent z exp = 1.19(2) up to L = 2048. For comparison, we also compute the effective dynamic exponents through the integrated correlation times. In addition, an exact result of the Swendsen–Wang dynamic spectrum of a one-dimensional Ising chain is derived.

In the present work, we present an extensive Monte Carlo simulation study on the dynamical properties of the two-dimensional random-bond Ising model. The correlation time τ of the Swendsen-Wang and Wolff cluster algorithms is calculated at the critical point. The dynamic critical exponent z of both algorithms is also measured by using the numerical data for several lattice sizes up to L=512. It is found for both algorithms that the autocorrelation time decreases considerably and the critical slowing-down effect reduces upon the introduction of bond disorder. Additionally, simulations with the Metropolis algorithm are performed, and the critical slowing-down effect is observed to be more pronounced in the presence of disorder, confirming the previous findings in the literature. Moreover, the existence of the non-self-averaging property of the model is demonstrated by calculating the scaled form of the standard deviation of autocorrelation times. Finally, the critical exponent ratio of the magnetic susceptibility is estimated by using the average cluster size of the Wolff algorithm.

Recently, it was shown [Y. Nonomura, J. Phys. Soc. Jpn. 83, 113001 (2014)JUPSAU0031-901510.7566/JPSJ.83.113001] that the nonequilibrium critical relaxation of the two-dimensional (2D) Ising model from a perfectly ordered state in the Wolff algorithm is described by stretched-exponential decay, and a universal scaling scheme was found to connect nonequilibrium and equilibrium behaviors. In the present study we extend these findings to vector spin models, and the 3D Heisenberg model could be a typical example. To evaluate the critical temperature and critical exponents precisely using the above scaling scheme, we calculate nonequilibrium ordering from the perfectly disordered state in the Swendsen-Wang algorithm, and we find that the critical ordering process is described by stretched-exponential growth with a comparable exponent to that of the 3D XY model. The critical exponents evaluated in the present study are consistent with those in previous studies.

We have measured the autocorrelations for the Swendsen-Wang and the Wolff cluster update algorithms for the Ising model in two, three, and four dimensions. The data for the Wolff algorithm suggest that the autocorrelations are linearly related to the specific heat, in which case the dynamic critical exponent is

{\mathit{z}}_{\mathrm{i}\mathrm{n}\mathrm{t},}

{\mathit{E}}^{\mathrm{W}}

=\ensuremath{\alpha}/\ensuremath{\nu}. For the Swendsen-Wang algorithm, scaling the autocorrelations by the average maximum cluster size gives either a constant or a logarithm, which implies that

{\mathit{z}}_{\mathrm{i}\mathrm{n}\mathrm{t},}

{\mathit{E}}^{\mathrm{SW}}

=\ensuremath{\beta}/\ensuremath{\nu} for the Ising model.

We simulate the critical relaxation process of the two-dimensional Ising model with the initial state both completely disordered or completely ordered. Results of a method to measure both the dynamic and static critical exponents are reported, based on the finite-size scaling for the dynamics at the early time. From the time-dependent Binder cumulant, the dynamical exponent z is extracted independently, while the static exponents \ensuremath{\beta}/\ensuremath{\nu} and \ensuremath{\nu} are obtained from the time evolution of the magnetization and its higher moments. \textcopyright{} 1996 The American Physical Society.

Recently we showed that the critical nonequilibrium relaxation in the Swendsen-Wang algorithm is widely described by the stretched-exponential relaxation of physical quantities in the Ising or Heisenberg models. Here we make a similar analysis in the Berezinsky-Kosterlitz-Thouless phase transition in the two-dimensional (2D) XY model and in the first-order phase transition in the 2D q=5 Potts model and find that these phase transitions are described by the simple exponential relaxation and power-law relaxation of physical quantities, respectively. We compare the relaxation behaviors of these phase transitions with those of the second-order phase transition in the three- and four-dimensional XY models and in the 2D q-state Potts models for 2≤q≤4 and show that the species of phase transitions can be clearly characterized by the present analysis. We also compare the size dependence of relaxation behaviors of the first-order phase transition in the 2D q=5 and 6 Potts models and propose a quantitative criterion on "weakness" of the first-order phase transition.