arXiv:cond-mat/0405197

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

J. Kaupužs

Published May 10, 2004• 1 citation

Abstract

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.

DOI: 10.48550/arxiv.cond-mat/0405197

Related papers

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/0201221Jan 14, 20020 citations

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

J. Kaupužs

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.

Computational Methods in Applied MathematicsJan 1, 20052 citations

Monte Carlo test of critical exponents and amplitudes in 3D Ising and phi^4 lattice models

J. Kaupužs

We have tested the leading correction-to-scaling exponent omega in O(n)-symmetric models on a three-dimensional lattice by analysing the recent Monte Carlo (MC) data. We have found that the effective critical exponent, estimated at finite sizes of the system L and L/2, decreases remarkably within the range of the simulated L values. This shows the incorrectness of some claims that omega has a very accurate value 0.845(10) at n=1. A selfconsistent infinite volume extrapolation yields row estimates omega=0.547, omega=0.573, and omega=0.625 at n=1, 2, and 3, respectively, in approximate agreement with the corresponding exact values 1/2, 5/9, and 3/5 predicted by our recently developed GFD (grouping of Feynman diagrams) theory. We have fitted the MC data for the susceptibility of 3D Ising model at criticality showing that the effective critical exponent eta tends to increase well above the usually accepted values around 0.036. We have fitted the data within [L;8L], including several terms in the asymptotic expansion with fixed exponents, to obtain the effective amplitudes depending on L. This method clearly demonstrates that the critical exponents of GFD theory are correct (the amplitudes converge to certain asymptotic values at L tending to infinity), whereas those of the perturbative renormalization group (RG) theory are incorrect (the amplitudes diverge). A modification of the standard Ising model by introducing suitable "improved" action (Hamiltonian) does not solve the problem in favour of the perturbative RG theory.

arXiv:cond-mat/0106271Jun 14, 20010 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 CAug 1, 20065 citations

Energy fluctuations and the singularity of specific heat in a 3D Ising model

J. Kaupužs

We study the energy fluctuations in 3D Ising model near the phase transition point. Specific heat is a relevant quantity which is directly related to the mean squared amplitude of the energy fluctuations in the system. We have made extensive Monte Carlo simulations in 3D Ising model to clarify the character of the singularity of the specific heat Cv based on the finite-size scaling of its maximal values Cvmax depending on the linear size of the lattice L. An original iterative method has been used which automatically finds the pseudocritical temperature corresponding to the maximum of Cv. The simulations made up to L ≤ 128 with application of the Wolff's cluster algorithm allowed us to verify the possible power-like as well as logarithmic singularity of the specific heat predicted by different theoretical treatments. The most challenging and interesting result we have obtained is that the finite-size scaling of Cvmax in 3D Ising model is well described by a logarithmic rather than power-like ansatz, just like in 2D case. Another modification of our iterative method has been considered to estimate the critical coupling of 3D Ising model from the Binder cumulant data within L ε [96; 384]. Furthermore, the critical exponent β has been evaluated from the simulated magnetization data within the range of reduced temperatures t ≥ 0.000086 and system sizes L ≤ 410.

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIEMay 25, 20042 citations

Corrections to scaling in the 3D Ising model: a comparison between MC and MCRG results

J. Kaupužs, Roderick Melnik

Corrections to scaling in the 3D Ising model are studied based on Monte Carlo (MC) simulation results for very large lattices with linear lattice sizes up to L=3456. Our estimated values of the correction-to-scaling exponent omega tend to decrease below the usually accepted value about 0.83 when the smallest lattice sizes are discarded from the fits. This behavior apparently confirms some of the known estimates of the Monte Carlo renormalization group (MCRG) method, i.e., omega about 0.7 and omega = 0.75(5). We discuss the possibilities that omega is either really smaller than usually expected or these values of omega describe some transient behavior. We propose refining MCRG simulations and analysis to resolve this issue. In distinction from omega, our actual MC estimations of the critical exponents eta and nu provide stable values eta=0.03632(13) and nu=0.63017(31), which well agree with those of the conformal bootstrap method, i.e., eta=0.0362978(20) and nu=0.6299709(40).

arXiv:2206.04393Jun 9, 20220 citations