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

Computational Methods in Applied Mathematics

Transfer Matrix Study of Finite-size Corrections in the 2D Ising Model

J. Kaupužs

Published Jan 1, 2005• 2 citations
Publisher

Abstract

Abstract Effective exact transfer matrix algorithms have been developed to compute the two-point correlation function G(r) of the 2D Ising model on a square finite size lattice. Systems including up to 800 spins have been considered and corrections to the finite-size scaling at the critical point have been analysed.As a new result, we have found that the correlation function has a nontrivial amplitude correction of a very small magnitude.

DOI
10.2478/cmam-2005-0003

Related papers

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

Transfer matrix and Monte Carlo tests of critical exponents in lattice models

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 of square geometry with periodic boundaries oriented either along <10> or along <11> direction 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 proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.

arXiv:cond-mat/0201221•Jan 14, 2002•0 citations

FINITE-SIZE CORRECTIONS TO CORRELATION FUNCTION AND SUSCEPTIBILITY IN 2D ISING MODEL

J. Kaupužs

Transfer matrix calculations of the critical two-point correlation function in 2D Ising model on a finite-size [Formula: see text] lattice with periodic boundaries along 〈11〉 direction are extended to L = 21. A refined analysis of the correlation function in 〈10〉 crystallographic direction at the distance r = L indicates the existence of a nontrivial finite-size correction of a very small amplitude with correction-to-scaling exponent ω < 2 in agreement with our foregoing study for L ≤ 20. Here we provide an additional evidence and show that amplitude a of the multiplicative correction term 1 + aL-ωis about -3.5·10-8if ω = 1/4 (the expected value). We calculate also the susceptibility for L ≤ 18 in order to compare our numerical estimates for the constant background contribution with the known very precise value and to look for possible nontrivial corrections to scaling. The numerical analysis reveals a perfect agreement for the background term, as well as shows that the nontrivial correction term, detected by our analysis in the correlation function, likely cancels in the susceptibility.

International Journal of Modern Physics C•Aug 1, 2006•5 citations

Finite-size analysis of the hard-square lattice gas

Wenan Guo, Henk W. J. Blöte

We investigate the hard-square lattice-gas model by means of transfer-matrix calculations and a finite-size-scaling analysis. Using a minimal set of assumptions we find that the spectrum of correction-to-scaling exponents is consistent with that of the exactly solved Ising model, and that the critical exponents and correlation-length amplitudes closely follow the relation predicted by conformal invariance. Assuming that these spectra are exactly identical, and conformal invariance, we determine the critical point, the conformal anomaly, and the temperature and magnetic exponents with numerical margins of 10(-11) or less. These results are in a perfect agreement with the exactly known Ising universal parameters in two dimensions. In order to obtain this degree of precision, we included system sizes as large as feasible, and used extended-precision floating-point arithmetic. The latter resource provided a substantial improvement of the analysis, despite the fact that it restricted the transfer-matrix calculations to finite sizes of at most 34 lattice units.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 30, 2002•52 citations

Finite-size scaling corrections in two-dimensional Ising and Potts ferromagnets

S L A de Queiroz

Finite-size corrections to scaling of critical correlation lengths and free energies of Ising and three-state Potts ferromagnets are analysed by numerical methods, on strips of width N sites of square, triangular and honeycomb lattices. Strong evidence is given that the amplitudes of the `analytical' correction terms, N -2 , are identically zero for triangular and honeycomb Ising systems. For Potts spins, our results are broadly consistent with this lattice-dependent pattern of cancellations, though for correlation lengths non-vanishing (albeit rather small) amplitudes cannot be entirely ruled out.

arXiv:cond-mat/9912090•Journal of Physics A Mathematical and General•Jan 24, 2000•27 citations

Finite-size scaling study of the two-dimensional Ising spin glass

David A. Huse, I. Morgenstern

We have calculated finite-size correlation lengths for strips of the square-lattice Ising spin glass using the transfer-matrix method. In order to minimize sampling errors we study strip lengths up to 5\ifmmode\times\else\texttimes\fi{}${10}^{6}$ lattice spacings. A phenomenological renormalization-group analysis indicates that there are strong corrections to simple power-law scaling near the zero-temperature critical point, as is to be expected near a lower critical dimension. We examine models with Gaussian, exp(-${J}^{2}$/2), and exponential, exp(-\ensuremath{\Vert}J\ensuremath{\Vert}), distributions of couplings; the Gaussian distribution shows stronger finite-size corrections. The correlation-length exponent is estimated to be \ensuremath{\nu}=4.2\ifmmode\pm\else\textpm\fi{}0.5, although we do not want to rule out the possibility that \ensuremath{\nu} is significantly larger than this.

Physical review. B, Condensed matter•Sep 1, 1985•40 citations
Abstract Effective exact transfer matrix algorithms have been developed to compute the two-point correlation function G(r) of the 2D Ising model on a square finite size lattice. Systems including up to 800 spins have been considered and corrections to the finite-size scaling at the critical point have been analysed.As a new result, we have found that the correlation function has a nontrivial amplitude correction of a very small magnitude.
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.
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 of square geometry with periodic boundaries oriented either along <10> or along <11> direction 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 proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.
Transfer matrix calculations of the critical two-point correlation function in 2D Ising model on a finite-size [Formula: see text] lattice with periodic boundaries along 〈11〉 direction are extended to L = 21. A refined analysis of the correlation function in 〈10〉 crystallographic direction at the distance r = L indicates the existence of a nontrivial finite-size correction of a very small amplitude with correction-to-scaling exponent ω < 2 in agreement with our foregoing study for L ≤ 20. Here we provide an additional evidence and show that amplitude a of the multiplicative correction term 1 + aL-ωis about -3.5·10-8if ω = 1/4 (the expected value). We calculate also the susceptibility for L ≤ 18 in order to compare our numerical estimates for the constant background contribution with the known very precise value and to look for possible nontrivial corrections to scaling. The numerical analysis reveals a perfect agreement for the background term, as well as shows that the nontrivial correction term, detected by our analysis in the correlation function, likely cancels in the susceptibility.
We investigate the hard-square lattice-gas model by means of transfer-matrix calculations and a finite-size-scaling analysis. Using a minimal set of assumptions we find that the spectrum of correction-to-scaling exponents is consistent with that of the exactly solved Ising model, and that the critical exponents and correlation-length amplitudes closely follow the relation predicted by conformal invariance. Assuming that these spectra are exactly identical, and conformal invariance, we determine the critical point, the conformal anomaly, and the temperature and magnetic exponents with numerical margins of 10(-11) or less. These results are in a perfect agreement with the exactly known Ising universal parameters in two dimensions. In order to obtain this degree of precision, we included system sizes as large as feasible, and used extended-precision floating-point arithmetic. The latter resource provided a substantial improvement of the analysis, despite the fact that it restricted the transfer-matrix calculations to finite sizes of at most 34 lattice units.
Finite-size corrections to scaling of critical correlation lengths and free energies of Ising and three-state Potts ferromagnets are analysed by numerical methods, on strips of width N sites of square, triangular and honeycomb lattices. Strong evidence is given that the amplitudes of the `analytical' correction terms, N -2 , are identically zero for triangular and honeycomb Ising systems. For Potts spins, our results are broadly consistent with this lattice-dependent pattern of cancellations, though for correlation lengths non-vanishing (albeit rather small) amplitudes cannot be entirely ruled out.
We have calculated finite-size correlation lengths for strips of the square-lattice Ising spin glass using the transfer-matrix method. In order to minimize sampling errors we study strip lengths up to 5\ifmmode\times\else\texttimes\fi{}106{10}^{6}106 lattice spacings. A phenomenological renormalization-group analysis indicates that there are strong corrections to simple power-law scaling near the zero-temperature critical point, as is to be expected near a lower critical dimension. We examine models with Gaussian, exp(-J2{J}^{2}J2/2), and exponential, exp(-\ensuremath{\Vert}J\ensuremath{\Vert}), distributions of couplings; the Gaussian distribution shows stronger finite-size corrections. The correlation-length exponent is estimated to be \ensuremath{\nu}=4.2\ifmmode\pm\else\textpm\fi{}0.5, although we do not want to rule out the possibility that \ensuremath{\nu} is significantly larger than this.