A careful Monte Carlo investigation of the phase transition very close to the critical point (T -> Tc, H -> 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents beta = 0.3126(4) =~ 5/16, delta^{-1} = 0.1997(4) =~ 1/5 and gamma_{3D} = 1.253(4) =~ 5/4. Our results indicate that, within experimental error, they are given by simple fractions corresponding to the linear interpolations between the respective two-dimensional (Onsager) and four-dimensional (mean field) critical exponents. An analysis of our inverse susceptibility data chi^{-1}(T) vs. /T - Tc/ shows that these data lead to a value of gamma_{3D} compatible with gamma' = gamma and Tc = 4.51152(12), while gamma values obtained recently by high and low temperature series expansions and renormalization group methods are not.
A careful Monte Carlo investigation of the phase transition very close to the critical point (T -> Tc, H -> 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents beta = 0.3126(4) =~ 5/16, delta^{-1} = 0.1997(4) =~ 1/5 and gamma_{3D} = 1.253(4) =~ 5/4. Our results indicate that, within experimental error, they are given by simple fractions corresponding to the linear interpolations between the respective two-dimensional (Onsager) and four-dimensional (mean field) critical exponents. An analysis of our inverse susceptibility data chi^{-1}(T) vs. /T - Tc/ shows that these data lead to a value of gamma_{3D} compatible with gamma' = gamma and Tc = 4.51152(12), while gamma values obtained recently by high and low temperature series expansions and renormalization group methods are not.
Monte Carlo investigations of magnetization versus field, Mc(H), at the critical temperature provide direct accurate results on the critical exponent δ−1 for one, two, three and four-dimensional lattices: δ−1 1D =0, δ−1 2D =0.0666(2)≈1/15, δ−1 3D =0.1997(4)≈1/5, δ−1 4D =0.332(5)≈1/3. This type of Monte Carlo data on δ, which is not easily found in studies of Ising lattices in the current literature, as far as we know, defines extremely well the numerical value of this exponent within very stringent limits.
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.
High temperature series expansion for the critical exponents of the Ising model are reanalysed using a modified ratio method. The analysis shows that a minor modification of the ratio method yields for all lattices a value, for the exponent γ in three dimensions, close to 1.245, therefore lower than the quoted value 1.250, and much closer to the renormalization group (R.G.) value 1.241. The exponent is analysed in two ways : in one method γ is estimated directly while in the other one Tc, is calculated first. With these new values of Tc, the exponent α is recalculated and found to be very close to the R.G. value 0.110. The value of ν is not modified (0.638) and is therefore still a problem for hyperscaling, and is in disagreement with the R.G. value (0.630).
We have tested the theoretical values of critical exponents, predicted for the three--dimensional Heisenberg model, based on the published Monte Carlo (MC) simulation data for the susceptibility. Two different sets of the critical exponents have been considered - one provided by the usual (perturbative) renormalization group (RG) theory, and another predicted by grouping of Feynman diagrams in phi^4 model (our theory). The test consists of two steps. First we determine the critical coupling by fitting the MC data to the theoretical expression, including both confluent and analytical corrections to scaling, the values of critical exponents being taken from theory. Then we use the obtained value of critical coupling to test the agreement between theory and MC data at criticality. As a result, we have found that predictions of our theory (gamma=19/14, eta=1/10, omega=3/5) are consistent, whereas those of the perturbative RG theory (gamma=1.3895, eta=0.0355, omega=0.782) are inconsistent with the MC data. The seemable agreement between the RG prediction for eta and MC results at criticality, reported in literature, appears due to slightly overestimated value of the critical coupling. Estimation of critical exponents of 3D Ising model from complex zeroth of the partition function is discussed. A refined analysis yields the best estimate 1/nu=1.518. We conclude that the recent MC data can be completely explained within our theory (providing 1/nu=1.5 and omega=0.5) rather than within the conventional RG theory.
We have performed large-scale numerical simulations for the ± J Ising model on the simple cubic and the square lattices with asymmetric weight of ferromagnetic and antiferromagnetic bonds. The Monte Carlo renormalization group method is used to estimate T c , γ/ν, β/ν and ν at various points along the phase boundary between ferromagnetic and paramagnetic phases including the tricritical point. The obtained values of T c and the resulting phase diagrams of the simple cubic and the square lattices are much more accurate than those estimated so far. In the simple cubic lattice, we find evidence for the possibility that universality of the ferromagnetic critical exponents holds along the critical line except at the tricritical point, while weak universality holds even at the tricritical point. Weak universality seems to hold also in the square lattice.
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.