Physical Review B

Locally converging algorithms for determining the critical temperature in Ising systems

Eshel Faraggi, Daniel T. Robb

Published Oct 20, 2008• 4 citations

Abstract

We introduce a class of algorithms that converge to criticality automatically, in a way similar to the invaded cluster algorithm. Unlike the invaded cluster algorithm which uses global percolation as a test for criticality, these local algorithms use an average over local observables, specifically the number of satisfied bonds, in a feedback loop which drives the system toward criticality. Two specific algorithms are introduced, the average algorithm and the locally converging Wolff algorithm. We apply these algorithms to study the Ising square lattice and the Ising Bethe lattice. We find reasonable convergence to the critical temperature for both systems under the locally converging Wolff algorithm. We also re-examine the phase diagram of the dilute two-dimensional (2D) Ising model and find results supporting our previously reported conclusions regarding the existence of a local regime of magnetization below the percolations threshold. In addition, the presented algorithms are computationally more efficient than the invaded cluster algorithm, requiring less CPU time and memory.

DOI: 10.1103/physrevb.78.134416

Related papers

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/9507094Physical Review LettersOct 9, 199577 citations

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. EFeb 28, 20242 citations

Single-cluster algorithm for the site-bond-correlated Ising model

Paulo R. A. Campos, Roberto Nicolau Onody

We extend the Wolff algorithm to include correlated spin interactions in diluted magnetic systems. This algorithm is applied to study the site-bond-correlated Ising model on a two-dimensional square lattice. We use a finite-size scaling procedure to obtain the phase diagram in the temperature-concentration space. We also have verified that the autocorrelation time diminishes in the presence of dilution and correlation, showing that the Wolff algorithm performs even better in such situations.

arXiv:cond-mat/9709317Physical review. B, Condensed matterDec 1, 19974 citations

Invaded cluster algorithm for critical properties of periodic and aperiodic planar Ising models

Oliver Redner, Michael Baake

We demonstrate that the invaded cluster algorithm, introduced by Machta et al (1995 Phys. Rev. Lett. 75 2792-5), is a fast and reliable tool for determining the critical temperature and the magnetic critical exponent of periodic and aperiodic ferromagnetic Ising models in two dimensions. The algorithm is shown to reproduce the known values of the critical temperature on various periodic and quasiperiodic graphs with an accuracy of more than three significant digits, but only modest computational effort. On two quasiperiodic graphs which were not investigated in this respect before, the 12-fold symmetric square-triangle tiling and the 10-fold symmetric Tübingen triangle tiling, we determine the critical temperature. Furthermore, a generalization of the algorithm to non-identical coupling strengths is presented and applied to a class of Ising models on the Labyrinth tiling. For generic cases in which the heuristic Harris-Luck criterion predicts deviations from the Onsager universality class, we find a magnetic critical exponent different from the Onsager value. But notable exceptions to the criterion are found which consist not only of the exactly solvable cases, in agreement with a recent exact result, but also of the self-dual ones and maybe more.

arXiv:cond-mat/9910136Journal of Physics A Mathematical and GeneralApr 28, 20006 citations

Comparison of cluster algorithms for the bond-diluted Ising model

Arnold H. Kole, G. T. Barkema, Lars Fritz

Monte Carlo cluster algorithms are popular for their efficiency in studying the Ising model near its critical temperature. We might expect that this efficiency extends to the bond-diluted Ising model. We show, however, that this is not always the case by comparing how the correlation times τ_{w} and τ_{sw} of the Wolff and Swendsen-Wang cluster algorithms scale as a function of the system size L when applied to the two-dimensional bond-diluted Ising model. We demonstrate that the Wolff algorithm suffers from a much longer correlation time than in the pure Ising model, caused by isolated (groups of) spins which are infrequently visited by the algorithm. With a simple argument we prove that these cause the correlation time τ_{w} to be bounded from below by L^{z_{w}} with a dynamical exponent z_{w}=γ/ν≈1.75 for a bond concentration p<1. Furthermore, we numerically show that this lower bound is actually taken for several values of p in the range 0.5<p<1. Moreover, we show that the Swendsen-Wang algorithm does not suffer from the same problem. Consequently, it has a much shorter correlation time, shorter than in the pure Ising model even. Numerically at p=0.6, we find that its dynamical exponent is z_{sw}=0.09(4).

arXiv:2107.08534Physical review. EJan 31, 20224 citations

Invaded cluster algorithm for critical properties of periodic and aperiodic planar Ising models

Oliver Redner, Michael Baake

. We demonstrate that the invaded cluster algorithm, recently introduced by Machta et al, is a fast and reliable tool for determining the critical temperature and the magnetic critical exponent of periodic and aperiodic ferromagnetic Ising models in two dimensions. The algorithm is shown to reproduce the known values of the critical temperature on various periodic and quasiperiodic graphs with an accuracy of more than three significant digits. On two quasiperiodic graphs which were not investigated in this respect before, the twelvefold symmetric square-triangle tiling and the tenfold symmetric Tubingen triangle tiling, we determine the critical temperature. Furthermore, a generalization of the algorithm to non-identical coupling strengths is presented and applied to a class of Ising models on the Labyrinth tiling. For generic cases in which the heuristic Harris-Luck criterion predicts deviations from the Onsager universality class, we find a magnetic critical exponent different from t...

Jan 1, 19994 citations