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/9710098 · Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Monte Carlo method to calculate the central charge and critical exponents

Paul J. M. Bastiaansen, H. J. F. Knops

Published Apr 1, 1998• 16 citations
arXivHTMLPDF

Abstract

We present a finite size scaling technique to calculate the central charge and some critical exponents of two-dimensional critical models with a Monte Carlo simulation. We use systems with dimensions $L\ifmmode\times\else\texttimes\fi{}M,$ and focus on the scaling behavior in $M/L.$ The finite size scaling relation that we use is the universal expression for the stress tensor on the torus. The stress tensor is the operator that governs the anisotropy of the system, and stems for the theory of conformal invariance. We show that a lattice representation of the stress tensor can easily be constructed, such that its expectation value on the torus can be calculated in a Monte Carlo simulation. In doing so, we observe that the stress tensor turns out to be remarkably insensitive to critical slowing down. We show that the typical simulation time scales with the linear system dimension $L$ roughly as ${L}^{4},$ and that this scaling holds for the ordinary Metropolis algorithm as well as for more sophisticated cluster algorithms, such that it is fruitless to invoke the latter. We test the method for the Ising model (with central charge $c=$$\frac{1}{2}$), the Ashkin-Teller model $(c=1),$ and the F model (also $c=1$).

DOI
10.1103/physreve.57.3784

Related papers

Determination of the conformal-field-theory central charge by the Wang-Landau algorithm

P. Belov, Anton Nazarov, А. О. Сорокин

We present a simple method to estimate the central charge of the conformal field theory corresponding to a critical point of a two-dimensional lattice model from Monte Carlo simulations. The main idea is to use the Wang-Landau flat-histogram algorithm, which allows us to obtain the free energy of a lattice model on a torus as a function of torus radii. The central charge is calculated with good precision from a free-energy scaling at the critical point. We apply the method to the Ising, tricritical Ising (Blume-Capel), Potts, and site-diluted Ising models, and we also discuss an estimation of the conformal weights.

arXiv:1612.06374•Physical review. E•Jun 14, 2017•5 citations

Single-cluster Monte Carlo study of the Ising model on two-dimensional random lattices

Wolfhard Janke E, Mohammad W. Katoot, Ramon Villanova

We use the single-cluster Monte Carlo update algorithm to simulate the Ising model on two-dimensional Poissonian random lattices with up to 80 000 sites which are linked together according to the Voronoi-Delaunay prescription. In one set of simulations we use reweighting techniques and finite-size scaling analysis to investigate the critical properties of the model in the very vicinity of the phase transition. In the other set of simulations we study the approach to criticality in the disordered phase, making use of improved estimators for measurements. From both sets of simulations we obtain clear evidence that the critical exponents agree with the exactly known exponents for regular lattices, i.e., that (lattice) universality holds for the two-dimensional Ising model.

arXiv:hep-lat/9402008•Physical review. B, Condensed matter•Apr 1, 1994•60 citations

Transfer matrix and Monte Carlo tests of critical exponents in Ising model

J. Kaupužs

The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction ~1/L^0.25 of a very small magnitude also has been detected, as it can be expected from our recently developed GFD (grouping of Feynman diagrams) theory. Monte Carlo simulations of the squared magnetization of 3D Ising model have been performed by Wolff's algorithm in the range of the reduced temperatures t =< 0.000086 and system sizes L =< 410. The effective critical exponent beta_eff(t) tends to increase above the currently accepted numerical values. The critical coupling K_c=0.22165386(51) has been extracted from the Binder cumulant data within 96 =< L =< 384. The critical exponent 1/nu, estimated from the finite-size scaling of the derivatives of the Binder cumulant, tends to decrease slightly below the RG value 1.587 for the largest system sizes. The finite-size scaling of accurately simulated maximal values of the specific heat C_v in 3D Ising model confirms a logarithmic rather than power-like critical singularity of C_v.

arXiv:cond-mat/0405197•May 10, 2004•1 citation

Monte Carlo method for critical systems in infinite volume: The planar Ising model

Victor Herdeiro, Benjamin Doyon

In this paper we propose a Monte Carlo method for generating finite-domain marginals of critical distributions of statistical models in infinite volume. The algorithm corrects the problem of the long-range effects of boundaries associated to generating critical distributions on finite lattices. It uses the advantage of scale invariance combined with ideas of the renormalization group in order to construct a type of "holographic" boundary condition that encodes the presence of an infinite volume beyond it. We check the quality of the distribution obtained in the case of the planar Ising model by comparing various observables with their infinite-plane prediction. We accurately reproduce planar two-, three-, and four-point of spin and energy operators. We also define a lattice stress-energy tensor, and numerically obtain the associated conformal Ward identities and the Ising central charge.

arXiv:1605.05350•Physical review. E•Oct 31, 2016•5 citations

Conformal anomaly and critical exponents of theXYIsing model

M. P. Nightingale, Enzo Granato, J. Michael Kosterlitz

We use extensive Monte Carlo transfer-matrix calculations on infinite strips of widths L up to 30 lattice spacing and a finite-size scaling analysis to obtain critical exponents and conformal anomaly number c for the two-dimensional XY Ising model. This model is expected to describe the critical behavior of a class of systems with simultaneous U(1) and ${\mathit{Z}}_{2}$ symmetries of which the fully frustrated XY model is a special case. The effective values obtained for c show a significant decrease with L at different points along the line where the transition to the ordered phase takes place in a single transition. Extrapolations based on power-law corrections give values consistent with c=3/2 although larger values cannot be ruled out. Critical exponents are obtained more accurately and are consistent with previous Monte Carlo simulations suggesting critical behavior and with recent calculations for the frustrated XY model.

arXiv:cond-mat/9505048•Physical review. B, Condensed matter•Sep 1, 1995•55 citations

A NEW APPROACH TO DYNAMIC FINITE-SIZE SCALING

MEHMET DİLAVER, SEMRA GÜNDÜÇ, MERAL AYDIN +1 more

In this work we have considered the Taylor series expansion of the dynamic scaling relation of the magnetization with respect to small initial magnetization values in order to study the dynamic scaling behavior of two- and three-dimensional Ising models. We have used the literature values of the critical exponents and of the new dynamic exponent x0 to observe the dynamic finite-size scaling behavior of the time evolution of the magnetization during early stages of the Monte Carlo simulation. For the three-dimensional Ising model we have also presented that this method opens the possibility of calculating z and x0 separately. Our results show good agreement with the literature values. Measurements done on lattices with different sizes seem to give very good scaling.

arXiv:cond-mat/0302571•International Journal of Modern Physics C•Sep 1, 2003•1 citation
We present a finite size scaling technique to calculate the central charge and some critical exponents of two-dimensional critical models with a Monte Carlo simulation. We use systems with dimensions L\ifmmode×\else\texttimes\fiM,L\ifmmode\times\else\texttimes\fi{}M,L\ifmmode×\else\texttimes\fiM, and focus on the scaling behavior in M/L.M/L.M/L. The finite size scaling relation that we use is the universal expression for the stress tensor on the torus. The stress tensor is the operator that governs the anisotropy of the system, and stems for the theory of conformal invariance. We show that a lattice representation of the stress tensor can easily be constructed, such that its expectation value on the torus can be calculated in a Monte Carlo simulation. In doing so, we observe that the stress tensor turns out to be remarkably insensitive to critical slowing down. We show that the typical simulation time scales with the linear system dimension LLL roughly as L4,{L}^{4},L4, and that this scaling holds for the ordinary Metropolis algorithm as well as for more sophisticated cluster algorithms, such that it is fruitless to invoke the latter. We test the method for the Ising model (with central charge c=$$\frac{1}{2}), the Ashkin-Teller model (c=1),(c=1),(c=1), and the F model (also c=1c=1c=1).
We present a simple method to estimate the central charge of the conformal field theory corresponding to a critical point of a two-dimensional lattice model from Monte Carlo simulations. The main idea is to use the Wang-Landau flat-histogram algorithm, which allows us to obtain the free energy of a lattice model on a torus as a function of torus radii. The central charge is calculated with good precision from a free-energy scaling at the critical point. We apply the method to the Ising, tricritical Ising (Blume-Capel), Potts, and site-diluted Ising models, and we also discuss an estimation of the conformal weights.
We use the single-cluster Monte Carlo update algorithm to simulate the Ising model on two-dimensional Poissonian random lattices with up to 80 000 sites which are linked together according to the Voronoi-Delaunay prescription. In one set of simulations we use reweighting techniques and finite-size scaling analysis to investigate the critical properties of the model in the very vicinity of the phase transition. In the other set of simulations we study the approach to criticality in the disordered phase, making use of improved estimators for measurements. From both sets of simulations we obtain clear evidence that the critical exponents agree with the exactly known exponents for regular lattices, i.e., that (lattice) universality holds for the two-dimensional Ising model.
The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction ~1/L^0.25 of a very small magnitude also has been detected, as it can be expected from our recently developed GFD (grouping of Feynman diagrams) theory. Monte Carlo simulations of the squared magnetization of 3D Ising model have been performed by Wolff's algorithm in the range of the reduced temperatures t =< 0.000086 and system sizes L =< 410. The effective critical exponent beta_eff(t) tends to increase above the currently accepted numerical values. The critical coupling K_c=0.22165386(51) has been extracted from the Binder cumulant data within 96 =< L =< 384. The critical exponent 1/nu, estimated from the finite-size scaling of the derivatives of the Binder cumulant, tends to decrease slightly below the RG value 1.587 for the largest system sizes. The finite-size scaling of accurately simulated maximal values of the specific heat C_v in 3D Ising model confirms a logarithmic rather than power-like critical singularity of C_v.
In this paper we propose a Monte Carlo method for generating finite-domain marginals of critical distributions of statistical models in infinite volume. The algorithm corrects the problem of the long-range effects of boundaries associated to generating critical distributions on finite lattices. It uses the advantage of scale invariance combined with ideas of the renormalization group in order to construct a type of "holographic" boundary condition that encodes the presence of an infinite volume beyond it. We check the quality of the distribution obtained in the case of the planar Ising model by comparing various observables with their infinite-plane prediction. We accurately reproduce planar two-, three-, and four-point of spin and energy operators. We also define a lattice stress-energy tensor, and numerically obtain the associated conformal Ward identities and the Ising central charge.
We use extensive Monte Carlo transfer-matrix calculations on infinite strips of widths L up to 30 lattice spacing and a finite-size scaling analysis to obtain critical exponents and conformal anomaly number c for the two-dimensional XY Ising model. This model is expected to describe the critical behavior of a class of systems with simultaneous U(1) and Z2{\mathit{Z}}_{2}Z2​ symmetries of which the fully frustrated XY model is a special case. The effective values obtained for c show a significant decrease with L at different points along the line where the transition to the ordered phase takes place in a single transition. Extrapolations based on power-law corrections give values consistent with c=3/2 although larger values cannot be ruled out. Critical exponents are obtained more accurately and are consistent with previous Monte Carlo simulations suggesting critical behavior and with recent calculations for the frustrated XY model.
In this work we have considered the Taylor series expansion of the dynamic scaling relation of the magnetization with respect to small initial magnetization values in order to study the dynamic scaling behavior of two- and three-dimensional Ising models. We have used the literature values of the critical exponents and of the new dynamic exponent x0 to observe the dynamic finite-size scaling behavior of the time evolution of the magnetization during early stages of the Monte Carlo simulation. For the three-dimensional Ising model we have also presented that this method opens the possibility of calculating z and x0 separately. Our results show good agreement with the literature values. Measurements done on lattices with different sizes seem to give very good scaling.