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/0011276 · Physical Review Letters

Probability-Changing Cluster Algorithm for Potts Models

Yusuke Tomita, Yutaka Okabe

Published Jan 22, 2001• 56 citations
arXivHTMLPDF

Abstract

We propose a new effective cluster algorithm of tuning the critical point automatically, which is an extended version of the Swendsen-Wang algorithm. We change the probability of connecting spins of the same type, p = 1-e(-J/k(B)T), in the process of the Monte Carlo spin update. Since we approach the canonical ensemble asymptotically, we can use the finite-size scaling analysis for physical quantities near the critical point. Simulating the two-dimensional Potts models to demonstrate the validity of the algorithm, we have obtained the critical temperatures and critical exponents which are consistent with the exact values; the comparison has been made with the invaded cluster algorithm.

DOI
10.1103/physrevlett.86.572

Related papers

Comparison of cluster algorithms for two-dimensional Potts models

CLIVE F. BAILLIE, Paul David Coddington

We have measured the dynamical critical exponent z for the Swendsen-Wang and the Wolff cluster update algorithms, as well as a number of variants of these algorithms, for the q=2 and q=3 Potts models in two dimensions. We find that although the autocorrelation times differ considerably between algorithms, the critical exponents are the same. For q=2, we find that although the data are better fitted by a logarithmic increase in the autocorrelation time with lattice size, they are also consistent with a power law with exponent z\ensuremath{\approxeq}0.25, especially if there are non-negligible corrections to scaling.

Physical review. B, Condensed matter•May 1, 1991•74 citations

Probability-Changing Cluster Algorithm: Study of Three-Dimensional Ising Model and Percolation Problem

Yusuke Tomita, Yutaka Okabe

We present a detailed description of the idea and procedure for the newly proposed Monte Carlo algorithm of tuning the critical point automatically, which is called the probability-changing cluster (PCC) algorithm [Y. Tomita and Y. Okabe, Phys. Rev. Lett. {\bf 86} (2001) 572]. Using the PCC algorithm, we investigate the three-dimensional Ising model and the bond percolation problem. We employ a refined finite-size scaling analysis to make estimates of critical point and exponents. With much less efforts, we obtain the results which are consistent with the previous calculations. We argue several directions for the application of the PCC algorithm.

arXiv:cond-mat/0203454•Journal of the Physical Society of Japan•Jun 15, 2002•15 citations

Single-cluster dynamics for the random-cluster model

Youjin Deng, Xiaofeng Qian, Henk W. J. Blöte

We formulate a single-cluster Monte Carlo algorithm for the simulation of the random-cluster model. This algorithm is a generalization of the Wolff single-cluster method for the q-state Potts model to noninteger values q>1. Its results for static quantities are in a satisfactory agreement with those of the existing Swendsen-Wang-Chayes-Machta (SWCM) algorithm, which involves a full-cluster decomposition of random-cluster configurations. We explore the critical dynamics of this algorithm for several two-dimensional Potts and random-cluster models. For integer q, the single-cluster algorithm can be reduced to the Wolff algorithm, for which case we find that the autocorrelation functions decay almost purely exponentially, with dynamic exponents z(exp)=0.07 (1), 0.521 (7), and 1.007 (9) for q=2, 3, and 4, respectively. For noninteger q, the dynamical behavior of the single-cluster algorithm appears to be very dissimilar to that of the SWCM algorithm. For large critical systems, the autocorrelation function displays a range of power-law behavior as a function of time. The dynamic exponents are relatively large. We provide an explanation for this peculiar dynamic behavior.

arXiv:0907.1711•Physical Review E•Sep 30, 2009•5 citations

Monte Carlo simulation of phase transitions in a two-dimensional random-bond Potts model

S. Chen, Alan M. Ferrenberg, D. P. Landau

Using the ``multihit'' Swendsen-Wang cluster flipping method, we performed extensive Monte Carlo simulations to investigate the critical behavior of the two-dimensional (2D) eight-state random-bond Potts model. We applied finite-size-scaling techniques to extract the critical exponents for two different sets of bond strengths, from which we concluded that the transition is second order with critical exponents for both sets falling into same universality class, that of a 2D Ising model. A variation of the Lee-Kosterlitz method for determining the order of a phase transition was also applied. The double-peaked structure of the specific heat, which was found in some of the bond configurations, was also studied by simulation on periodic arrangements of strong and weak bonds.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 1, 1995•105 citations

Replica-exchange cluster algorithm

Wolfhard Janke E, Elmar Bittner

In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature/energy range around the critical point.By combining the replica-exchange algorithm with cluster updates and an adaptive routine to find the range of interest, we introduce a new flexible and powerful method for systematic investigations of critical phenomena.As a result, we gain two further orders of magnitude in the performance for 2D and 3D Ising models in comparison with the recently proposed Wang-Landau recursion for cluster algorithms based on the multibondic algorithm, which is already a great improvement over the standard multicanonical variant.

arXiv:1108.0354•International Symposium on Lattice Field Theory (LATTICE)•Jun 23, 2010•0 citations

Invaded Cluster Algorithm for Equilibrium Critical Points

Jonathan L. Machta, Y. S. Choi, A. Lücke +2 more

A new cluster algorithm based on invasion percolation is described. The algorithm samples the critical point of a spin system without a priori knowledge of the critical temperature and provides an efficient way to determine the critical temperature and other observables in the critical region. The method is illustrated for the two- and three-dimensional Ising models. The algorithm equilibrates spin configurations much faster than the closely related Swendsen-Wang algorithm.

arXiv:cond-mat/9507094•Physical Review Letters•Oct 9, 1995•77 citations
We propose a new effective cluster algorithm of tuning the critical point automatically, which is an extended version of the Swendsen-Wang algorithm. We change the probability of connecting spins of the same type, p = 1-e(-J/k(B)T), in the process of the Monte Carlo spin update. Since we approach the canonical ensemble asymptotically, we can use the finite-size scaling analysis for physical quantities near the critical point. Simulating the two-dimensional Potts models to demonstrate the validity of the algorithm, we have obtained the critical temperatures and critical exponents which are consistent with the exact values; the comparison has been made with the invaded cluster algorithm.
We have measured the dynamical critical exponent z for the Swendsen-Wang and the Wolff cluster update algorithms, as well as a number of variants of these algorithms, for the q=2 and q=3 Potts models in two dimensions. We find that although the autocorrelation times differ considerably between algorithms, the critical exponents are the same. For q=2, we find that although the data are better fitted by a logarithmic increase in the autocorrelation time with lattice size, they are also consistent with a power law with exponent z\ensuremath{\approxeq}0.25, especially if there are non-negligible corrections to scaling.
We present a detailed description of the idea and procedure for the newly proposed Monte Carlo algorithm of tuning the critical point automatically, which is called the probability-changing cluster (PCC) algorithm [Y. Tomita and Y. Okabe, Phys. Rev. Lett. {\bf 86} (2001) 572]. Using the PCC algorithm, we investigate the three-dimensional Ising model and the bond percolation problem. We employ a refined finite-size scaling analysis to make estimates of critical point and exponents. With much less efforts, we obtain the results which are consistent with the previous calculations. We argue several directions for the application of the PCC algorithm.
We formulate a single-cluster Monte Carlo algorithm for the simulation of the random-cluster model. This algorithm is a generalization of the Wolff single-cluster method for the q-state Potts model to noninteger values q>1. Its results for static quantities are in a satisfactory agreement with those of the existing Swendsen-Wang-Chayes-Machta (SWCM) algorithm, which involves a full-cluster decomposition of random-cluster configurations. We explore the critical dynamics of this algorithm for several two-dimensional Potts and random-cluster models. For integer q, the single-cluster algorithm can be reduced to the Wolff algorithm, for which case we find that the autocorrelation functions decay almost purely exponentially, with dynamic exponents z(exp)=0.07 (1), 0.521 (7), and 1.007 (9) for q=2, 3, and 4, respectively. For noninteger q, the dynamical behavior of the single-cluster algorithm appears to be very dissimilar to that of the SWCM algorithm. For large critical systems, the autocorrelation function displays a range of power-law behavior as a function of time. The dynamic exponents are relatively large. We provide an explanation for this peculiar dynamic behavior.
Using the ``multihit'' Swendsen-Wang cluster flipping method, we performed extensive Monte Carlo simulations to investigate the critical behavior of the two-dimensional (2D) eight-state random-bond Potts model. We applied finite-size-scaling techniques to extract the critical exponents for two different sets of bond strengths, from which we concluded that the transition is second order with critical exponents for both sets falling into same universality class, that of a 2D Ising model. A variation of the Lee-Kosterlitz method for determining the order of a phase transition was also applied. The double-peaked structure of the specific heat, which was found in some of the bond configurations, was also studied by simulation on periodic arrangements of strong and weak bonds.
In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature/energy range around the critical point.By combining the replica-exchange algorithm with cluster updates and an adaptive routine to find the range of interest, we introduce a new flexible and powerful method for systematic investigations of critical phenomena.As a result, we gain two further orders of magnitude in the performance for 2D and 3D Ising models in comparison with the recently proposed Wang-Landau recursion for cluster algorithms based on the multibondic algorithm, which is already a great improvement over the standard multicanonical variant.