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

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Cluster Monte Carlo simulation of the transverse Ising model

Henk W. J. Blöte, Youjin Deng

Published Dec 10, 2002• 282 citations
Publisher

Abstract

We formulate a cluster Monte Carlo method for the anisotropic limit of Ising models on (d+1)-dimensional lattices, which in effect, are equivalent with d-dimensional quantum transverse Ising models. Using this technique, we investigate the transverse Ising models on the square, triangular, Kagome, honeycomb, and simple-cubic lattices. The Monte Carlo data are analyzed by finite-size scaling. In each case we find, as expected, that the critical behavior fits well in the (d+1)-dimensional Ising universality class. For the transverse Ising model on the square lattice, we determine the Binder cumulant of the classical counterpart for a range of aspect ratios between the system sizes in the third or "classical" direction and that in the other two directions. Matching this universal function with the case of the isotropic Ising model yields the length ratio relating the isotropic Ising model with the anisotropic limit. The efficiency of the present algorithm is reflected by the precision of its results, which improves significantly on earlier analyses.

DOI
10.1103/physreve.66.066110

Related papers

Finite size scaling study of lattice models in the three-dimensional Ising universality class

Martin Hasenbusch B

We study the spin-1/2 Ising model and the Blume-Capel model at various values of the parameter $D$ on the simple cubic lattice. To this end we perform Monte Carlo simulations using a hybrid of the local Metropolis, the single cluster and the wall cluster algorithm. Using finite size scaling we determine the value ${D}^{\ensuremath{\ast}}=0.656(20)$ of the parameter $D$, where leading corrections to scaling vanish. We find $\ensuremath{\omega}=0.832(6)$ for the exponent of leading corrections to scaling. In order to compute accurate estimates of critical exponents, we construct improved observables that have a small amplitude of the leading correction for any model. Analyzing data obtained for $D=0.641$ and 0.655 on lattices of a linear size up to $L=360$ we obtain $\ensuremath{\nu}=0.63002(10)$ and $\ensuremath{\eta}=0.03627(10)$. We compare our results with those obtained from previous Monte Carlo simulations and high-temperature series expansions of lattice models, by using field-theoretic methods and experiments.

arXiv:1004.4486•Physical Review B•Nov 23, 2010•333 citations

Ising model on three-dimensional random lattices: A Monte Carlo study

Wolfhard Janke E, Ramon Villanova

We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.

arXiv:cond-mat/0211446•Physical review. B, Condensed matter•Oct 31, 2002•31 citations

Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model

Alan M. Ferrenberg, Jiahao Xu, David Paul Landau

While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.

arXiv:1806.03558•Physical review. E•Apr 4, 2018•191 citations

Critical behavior of (1+1)-dimensional transverse Ising models for spin S≥12

K. A. Penson, M. Kolb

The phase transitions in the ground state of (1+1)-dimensional Ising models in a transverse field are investigated for spin $S=\frac{1}{2}, 1, \frac{3}{2}, \mathrm{and} 2$. The critical fields and "thermal" and correlation exponents are calculated as a function of $S$ by the finite-size scaling method. The universality of the transition with respect to $S$ is fully confirmed. Furthermore, we elaborate on a recent conjecture by Luck and others relating the correlation exponent $\ensuremath{\eta}$ to the amplitude of a correlation length for a two-dimensional (2D) classical isotropic system. The extension of this relation to quantum chains as extremely anisotropic 2D systems has been considered. We find that while no simple prescription can be given to calculate various correlation exponents from the amplitudes, the ratios of the amplitudes of different quantities still have universal character.

Physical review. B, Condensed matter•Aug 1, 1984•10 citations

Quantum-critical properties of the long-range transverse-field Ising model from quantum Monte Carlo simulations

Jan Alexander Koziol, Anja Langheld, Sebastian C. Kapfer +1 more

The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.

arXiv:2103.09469•Physical review. B./Physical review. B•Jun 23, 2021•47 citations

Ising universality in three dimensions: a Monte Carlo study

Henk W. J. Blöte, Erik Luijten, Jouke R. Heringa

We investigate three Ising models on the simple cubic lattice by means of Monte Carlo methods and finite-size scaling. These models are the spin-1/2 Ising model with nearest-neighbour interactions, a spin-1/2 model with nearest-neighbour and third-neighbour interactions, and a spin-1 model with nearest-neighbour interactions. The results are in accurate agreement with the hypothesis of universality. Analysis of the finite-size scaling behaviour reveals corrections beyond those caused by the leading irrelevant scaling field. We find that the correction-to-scaling amplitudes are strongly dependent on the introduction of further-neighbour interactions or a third spin state. In a spin-1 Ising model, these corrections appear to be very small. This is very helpful for the determination of the universal constants of the Ising model. The renormalization exponents of the Ising model are determined as y t =1.587 (2), y h =2.4815 (15) and y i =-0.82 (6). The universal ratio Q=(m 2 ) 2 /(m 4 ) is equal to 0.6233 (4) for periodic systems with cubic symmetry. The critical point of the nearest-neighbour spin-1/2 model is K c =0.2216546 (10).

arXiv:cond-mat/9509016•Journal of Physics A Mathematical and General•Nov 21, 1995•308 citations
We formulate a cluster Monte Carlo method for the anisotropic limit of Ising models on (d+1)-dimensional lattices, which in effect, are equivalent with d-dimensional quantum transverse Ising models. Using this technique, we investigate the transverse Ising models on the square, triangular, Kagome, honeycomb, and simple-cubic lattices. The Monte Carlo data are analyzed by finite-size scaling. In each case we find, as expected, that the critical behavior fits well in the (d+1)-dimensional Ising universality class. For the transverse Ising model on the square lattice, we determine the Binder cumulant of the classical counterpart for a range of aspect ratios between the system sizes in the third or "classical" direction and that in the other two directions. Matching this universal function with the case of the isotropic Ising model yields the length ratio relating the isotropic Ising model with the anisotropic limit. The efficiency of the present algorithm is reflected by the precision of its results, which improves significantly on earlier analyses.
We study the spin-1/2 Ising model and the Blume-Capel model at various values of the parameter DDD on the simple cubic lattice. To this end we perform Monte Carlo simulations using a hybrid of the local Metropolis, the single cluster and the wall cluster algorithm. Using finite size scaling we determine the value ${D}^{\ensuremath{\ast}}=0.656(20)$ of the parameter DDD, where leading corrections to scaling vanish. We find $\ensuremath{\omega}=0.832(6)$ for the exponent of leading corrections to scaling. In order to compute accurate estimates of critical exponents, we construct improved observables that have a small amplitude of the leading correction for any model. Analyzing data obtained for D=0.641D=0.641D=0.641 and 0.655 on lattices of a linear size up to L=360L=360L=360 we obtain $\ensuremath{\nu}=0.63002(10)$ and $\ensuremath{\eta}=0.03627(10)$. We compare our results with those obtained from previous Monte Carlo simulations and high-temperature series expansions of lattice models, by using field-theoretic methods and experiments.
We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.
While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.
The phase transitions in the ground state of (1+1)-dimensional Ising models in a transverse field are investigated for spin S=12,1,32,and2S=\frac{1}{2}, 1, \frac{3}{2}, \mathrm{and} 2S=21​,1,23​,and2. The critical fields and "thermal" and correlation exponents are calculated as a function of SSS by the finite-size scaling method. The universality of the transition with respect to SSS is fully confirmed. Furthermore, we elaborate on a recent conjecture by Luck and others relating the correlation exponent $\ensuremath{\eta}$ to the amplitude of a correlation length for a two-dimensional (2D) classical isotropic system. The extension of this relation to quantum chains as extremely anisotropic 2D systems has been considered. We find that while no simple prescription can be given to calculate various correlation exponents from the amplitudes, the ratios of the amplitudes of different quantities still have universal character.
The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.
We investigate three Ising models on the simple cubic lattice by means of Monte Carlo methods and finite-size scaling. These models are the spin-1/2 Ising model with nearest-neighbour interactions, a spin-1/2 model with nearest-neighbour and third-neighbour interactions, and a spin-1 model with nearest-neighbour interactions. The results are in accurate agreement with the hypothesis of universality. Analysis of the finite-size scaling behaviour reveals corrections beyond those caused by the leading irrelevant scaling field. We find that the correction-to-scaling amplitudes are strongly dependent on the introduction of further-neighbour interactions or a third spin state. In a spin-1 Ising model, these corrections appear to be very small. This is very helpful for the determination of the universal constants of the Ising model. The renormalization exponents of the Ising model are determined as y t =1.587 (2), y h =2.4815 (15) and y i =-0.82 (6). The universal ratio Q=(m 2 ) 2 /(m 4 ) is equal to 0.6233 (4) for periodic systems with cubic symmetry. The critical point of the nearest-neighbour spin-1/2 model is K c =0.2216546 (10).