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

International Journal of Modern Physics C

SWENDSEN-WANG DYNAMICS ON LARGE 2D CRITICAL ISING MODELS

János Kertész, Dietrich Stauffer

Published Dec 1, 1992• 12 citations
Publisher

Abstract

Lattices up to size 6400*6400 are simulated on a 32-processor Intel Hypercube parallel computer. At the critical temperature of the square-lattice Ising model, the magnetization decays roughly exponentially with time if we start with all spins parallel and if the lattice is large enough. This result is consistent with those of Heermann and Burkitt.

DOI
10.1142/s0129183192000877

Related papers

LARGE SCALE SIMULATIONS OF THE KINETIC ISING MODEL

Christian Münkel

We present Monte Carlo simulation results for the dynamical critical exponent z of the two- and three-dimensional kinetic Ising model. The z-values were calculated from the magnetization relaxation from an ordered state into the equilibrium state at Tc for very large systems with up to (169984)2 and (3072)3 spins. To our knowledge, these are the largest Ising-systems simulated todate. We also report the successful simulation of very large lattices on a massively parallel MIMD computer with high speedups of approximately 1000 and an efficiency of about 0.93.

International Journal of Modern Physics C•Dec 1, 1993•4 citations

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

Jianqing Du, Bo Zheng, Jian‐Sheng Wang

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.

arXiv:cond-mat/0603038•Journal of Statistical Mechanics Theory and Experiment•May 11, 2006•14 citations

Critical Behavior Simulation of the Two-Dimensional Ising Model—A Computational Physics Study Based on the Monte Carlo Method

佳钢 李

To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.

Modern Physics•Jan 1, 2026•0 citations

The critical region of the random-bond Ising model

A. L. Talapov, Щур Лев Николаевич

We describe results of the cluster algorithm special purpose processor simulations of the 2D Ising model with impurity bonds. We use large lattices with up to 10 6 spins to define critical temperature regions where both finite-size corrections and corrections to scaling are small. High-accuracy data unambiguously show an increase in the magnetization and magnetic susceptibility effective exponents beta and gamma , caused by impurities. The Mand chi -singularities become sharper, while the specific heat singularity is smoothed.

arXiv:hep-lat/9404011•Journal of Physics Condensed Matter•Oct 3, 1994•25 citations

HIGH SPEED ALGORITHMS AND SIMULATIONS OF THE CRITICAL DYNAMICS OF THE ISING MODEL

Hans-Otto Heuer

We review the development of Monte Carlo algorithms for the simulation of the Ising model with local dynamics. The multi lattice coding algorithm for the pure and for the site-disordered Ising model is described in detail. Application of these algorithms on fast vector computers have led to new progress concerning the critical dynamics of the pure Ising model leading to the critical exponent z = 2:05. The importance of disorder for the critical dynamics of the Ising model has been investigated recently. We review the numerical results and compare with theoretical predictions. A crossover scenario is found as for the static behavior with an effective dynamical exponent z = 1:75 fl. 1 Introduction The Ising model is a standard model in statistical physics for the description of a ferromagnet in an extremely simplified manner. It is assumed that quantum spin operators may be replaced by classical spins S i = \\Sigma 1 which are oriented parallel or antiparallel to a magnetic filed H: H...

WORLD SCIENTIFIC eBooks•May 1, 1996•1 citation

Dependence of Critical Parameters of 2D Ising Model on Lattice Size

B. V. Kryzhanovsky, M. Yu. Malsagov, Iakov Karandashev

For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with N = l × l spins and the asymptotic Onsager solution ( N → ∞). We derived empirical expressions for critical parameters as functions of N and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when N increased the heat capacity in the critical point increased as ln N . We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis Monte Carlo method.

arXiv:2002.00615•Optical Memory and Neural Networks•Jan 1, 2018•7 citations
Lattices up to size 6400*6400 are simulated on a 32-processor Intel Hypercube parallel computer. At the critical temperature of the square-lattice Ising model, the magnetization decays roughly exponentially with time if we start with all spins parallel and if the lattice is large enough. This result is consistent with those of Heermann and Burkitt.
We present Monte Carlo simulation results for the dynamical critical exponent z of the two- and three-dimensional kinetic Ising model. The z-values were calculated from the magnetization relaxation from an ordered state into the equilibrium state at Tc for very large systems with up to (169984)2 and (3072)3 spins. To our knowledge, these are the largest Ising-systems simulated todate. We also report the successful simulation of very large lattices on a massively parallel MIMD computer with high speedups of approximately 1000 and an efficiency of about 0.93.
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.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.
We describe results of the cluster algorithm special purpose processor simulations of the 2D Ising model with impurity bonds. We use large lattices with up to 10 6 spins to define critical temperature regions where both finite-size corrections and corrections to scaling are small. High-accuracy data unambiguously show an increase in the magnetization and magnetic susceptibility effective exponents beta and gamma , caused by impurities. The Mand chi -singularities become sharper, while the specific heat singularity is smoothed.
We review the development of Monte Carlo algorithms for the simulation of the Ising model with local dynamics. The multi lattice coding algorithm for the pure and for the site-disordered Ising model is described in detail. Application of these algorithms on fast vector computers have led to new progress concerning the critical dynamics of the pure Ising model leading to the critical exponent z = 2:05. The importance of disorder for the critical dynamics of the Ising model has been investigated recently. We review the numerical results and compare with theoretical predictions. A crossover scenario is found as for the static behavior with an effective dynamical exponent z = 1:75 fl. 1 Introduction The Ising model is a standard model in statistical physics for the description of a ferromagnet in an extremely simplified manner. It is assumed that quantum spin operators may be replaced by classical spins S i = \\Sigma 1 which are oriented parallel or antiparallel to a magnetic filed H: H...
For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with N = l × l spins and the asymptotic Onsager solution ( N → ∞). We derived empirical expressions for critical parameters as functions of N and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when N increased the heat capacity in the critical point increased as ln N . We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis Monte Carlo method.