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

Invaded Cluster Simulations of the XY Model in Two and Three

Ilija Dukovski, Jonathan L. Machta, Lincoln Chayes

Published Jan 1, 2002• 0 citations
OpenAlex

Abstract

The invaded cluster algorithm is used to study the XY model in two and three dimensions up to sizes 20002 and 1203, respectively. A soft spin O(2) model, in the same universality class as the three-dimensional XY model, is also studied. The static critical properties of the model and the dynamical properties of the algorithm are reported. The results are Kc=0.45412(2) for the three-dimensional XY model and η=0.037(2) for the three-dimensional XY universality class. For the two-dimensional XY model the results are Kc=1.120(1) and η=0.251(5). The invaded cluster algorithm does not show any critical slowing for the magnetization or critical temperature estimator for the two-dimensional or three-dimensional XY models.

Related papers

Invaded cluster simulations of theXYmodel in two and three dimensions

Ilija Dukovski, Jonathan L. Machta, Lincoln Chayes

The invaded cluster algorithm is used to study the XY model in two and three dimensions up to sizes 2000(2) and 120(3), respectively. A soft spin O(2) model, in the same universality class as the three-dimensional XY model, is also studied. The static critical properties of the model and the dynamical properties of the algorithm are reported. The results are K(c)=0.454 12(2) for the three-dimensional XY model and eta=0.037(2) for the three-dimensional XY universality class. For the two-dimensional XY model the results are K(c)=1.120(1) and eta=0.251(5). The invaded cluster algorithm does not show any critical slowing for the magnetization or critical temperature estimator for the two-dimensional or three-dimensional XY models.

arXiv:cond-mat/0105143•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 3, 2002•24 citations

Cluster Monte Carlo and network flow algorithms for critical phenomena

Ilija Dukovski

Phase transitions and critical phenomena in several spin models are studied. These models include a two and three dimensional XY model and Ising models in structured and random field. The methods of study are cluster Monte Carlo simulations and network flow optimization algorithms. The XY model and a related soft spin model belonging to the O(2) universality class are studied using a modified version of the Invaded Cluster algorithm. The static critical properties of the models and the dynamic characteristics of the algorithm are studied. The critical exponent η for the O(2) universality class in two and three dimensions is measured. The Invaded Cluster algorithm for the XY model does not show any critical slowing for global variables. Benzene adsorption in a nanoporous solid is modeled with an Ising model with both internal and external magnetic fields. The first order phase transition and the critical properties of the system are studied. The system is simulated using the two replica cluster algorithm. The critical temperature of the system is found to decrease with increasing difference of the values of the internal fields. It is found that the system undergoes a phase transition for a relatively wide, experimentally obtainable, region of values of its parameters. The random field Ising model in three dimensions is studied at zero temperature. The nature of randomness induced phase transition is studied using a numerical network flow optimization method. It is found that the latent heat at the phase transition vanishes in the limit of infinite system size which can be taken as a proof of a second order phase transition.

•Jan 1, 2001•0 citations

PARALLEL INVADED CLUSTER ALGORITHM FOR THE ISING MODEL

Y. S. CHOI, Jonathan L. Machta, Pablo Tamayo +1 more

A parallel version of the invaded cluster algorithm is described. Results from large scale (up to 40962 and 5123) simulations of the Ising model are reported. No evidence of critical slowing down is found for the three-dimensional Ising model. The magnetic exponent is estimated to be 2.482±0.001(β/ν=0.518±0.001) for the three-dimensional Ising model.

arXiv:cond-mat/9806127•International Journal of Modern Physics C•Feb 1, 1999•16 citations

Invaded cluster algorithm for Potts models

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

The invaded cluster algorithm, a method for simulating phase transitions, is described in detail. Theoretical, albeit nonrigorous, justification of the method is presented and the algorithm is applied to Potts models in two and three dimensions. The algorithm is shown to be useful for both first-order and continuous transitions and evidently provides an efficient way to distinguish between these possibilities. The dynamic properties of the invaded cluster algorithm are studied. Numerical evidence suggests that the algorithm has no critical slowing for Ising models. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9605103•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 1, 1996•80 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, 1999•4 citations

Dynamic and static properties of the invaded cluster algorithm

K. Moriarty, Jonathan L. Machta, L. Chayes

Simulations of the two-dimensional Ising and three-state Potts models at their critical points are performed using the invaded cluster (IC) algorithm. It is argued that observables measured on a sublattice of size l should exhibit a crossover to Swendsen-Wang (SW) behavior for l sufficiently less than the lattice size L, and a scaling form is proposed to describe the crossover phenomenon. It is found that the energy autocorrelation time ${\ensuremath{\tau}}_{\ensuremath{\varepsilon}}(l,L)$ for an $l\ifmmode\times\else\texttimes\fi{}l$ sublattice attains a maximum in the crossover region, and a dynamic exponent ${z}^{\mathrm{IC}}$ for the IC algorithm is defined according to ${\ensuremath{\tau}}_{\ensuremath{\varepsilon},\mathrm{max}}\ensuremath{\sim}{L}^{{z}^{\mathrm{IC}}}.$ Simulation results for the three-state model yield ${z}^{\mathrm{IC}}=0.346\ifmmode\pm\else\textpm\fi{}0.002,$ which is smaller than values of the dynamic exponent found for the SW and Wolff algorithms and also less than the Li-Sokal bound. The results are less conclusive for the Ising model, but it appears that ${z}^{\mathrm{IC}}<0.21$ and possibly that ${\ensuremath{\tau}}_{\ensuremath{\varepsilon},\mathrm{max}}\ensuremath{\sim}\mathrm{ln}L$ so that ${z}^{\mathrm{IC}}=0$---similar to previous results for the SW and Wolff algorithms.

arXiv:cond-mat/9811138•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Feb 1, 1999•8 citations
The invaded cluster algorithm is used to study the XY model in two and three dimensions up to sizes 20002 and 1203, respectively. A soft spin O(2) model, in the same universality class as the three-dimensional XY model, is also studied. The static critical properties of the model and the dynamical properties of the algorithm are reported. The results are Kc=0.45412(2) for the three-dimensional XY model and η=0.037(2) for the three-dimensional XY universality class. For the two-dimensional XY model the results are Kc=1.120(1) and η=0.251(5). The invaded cluster algorithm does not show any critical slowing for the magnetization or critical temperature estimator for the two-dimensional or three-dimensional XY models.
The invaded cluster algorithm is used to study the XY model in two and three dimensions up to sizes 2000(2) and 120(3), respectively. A soft spin O(2) model, in the same universality class as the three-dimensional XY model, is also studied. The static critical properties of the model and the dynamical properties of the algorithm are reported. The results are K(c)=0.454 12(2) for the three-dimensional XY model and eta=0.037(2) for the three-dimensional XY universality class. For the two-dimensional XY model the results are K(c)=1.120(1) and eta=0.251(5). The invaded cluster algorithm does not show any critical slowing for the magnetization or critical temperature estimator for the two-dimensional or three-dimensional XY models.
Phase transitions and critical phenomena in several spin models are studied. These models include a two and three dimensional XY model and Ising models in structured and random field. The methods of study are cluster Monte Carlo simulations and network flow optimization algorithms. The XY model and a related soft spin model belonging to the O(2) universality class are studied using a modified version of the Invaded Cluster algorithm. The static critical properties of the models and the dynamic characteristics of the algorithm are studied. The critical exponent η for the O(2) universality class in two and three dimensions is measured. The Invaded Cluster algorithm for the XY model does not show any critical slowing for global variables. Benzene adsorption in a nanoporous solid is modeled with an Ising model with both internal and external magnetic fields. The first order phase transition and the critical properties of the system are studied. The system is simulated using the two replica cluster algorithm. The critical temperature of the system is found to decrease with increasing difference of the values of the internal fields. It is found that the system undergoes a phase transition for a relatively wide, experimentally obtainable, region of values of its parameters. The random field Ising model in three dimensions is studied at zero temperature. The nature of randomness induced phase transition is studied using a numerical network flow optimization method. It is found that the latent heat at the phase transition vanishes in the limit of infinite system size which can be taken as a proof of a second order phase transition.
The invaded cluster algorithm, a method for simulating phase transitions, is described in detail. Theoretical, albeit nonrigorous, justification of the method is presented and the algorithm is applied to Potts models in two and three dimensions. The algorithm is shown to be useful for both first-order and continuous transitions and evidently provides an efficient way to distinguish between these possibilities. The dynamic properties of the invaded cluster algorithm are studied. Numerical evidence suggests that the algorithm has no critical slowing for Ising models. \textcopyright{} 1996 The American Physical Society.
. 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...
Simulations of the two-dimensional Ising and three-state Potts models at their critical points are performed using the invaded cluster (IC) algorithm. It is argued that observables measured on a sublattice of size l should exhibit a crossover to Swendsen-Wang (SW) behavior for l sufficiently less than the lattice size L, and a scaling form is proposed to describe the crossover phenomenon. It is found that the energy autocorrelation time ${\ensuremath{\tau}}_{\ensuremath{\varepsilon}}(l,L)$ for an $l\ifmmode\times\else\texttimes\fi{}l$ sublattice attains a maximum in the crossover region, and a dynamic exponent zIC{z}^{\mathrm{IC}}zIC for the IC algorithm is defined according to ${\ensuremath{\tau}}_{\ensuremath{\varepsilon},\mathrm{max}}\ensuremath{\sim}{L}^{{z}^{\mathrm{IC}}}.$ Simulation results for the three-state model yield ${z}^{\mathrm{IC}}=0.346\ifmmode\pm\else\textpm\fi{}0.002,$ which is smaller than values of the dynamic exponent found for the SW and Wolff algorithms and also less than the Li-Sokal bound. The results are less conclusive for the Ising model, but it appears that zIC<0.21{z}^{\mathrm{IC}}<0.21zIC<0.21 and possibly that ${\ensuremath{\tau}}_{\ensuremath{\varepsilon},\mathrm{max}}\ensuremath{\sim}\mathrm{ln}L$ so that zIC=0{z}^{\mathrm{IC}}=0zIC=0---similar to previous results for the SW and Wolff algorithms.