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

Progress of Theoretical Physics Supplement

Transfer Matrix and Finite-Size Scaling for the Ising Model on Two- and Three-Dimensional Lattices

Takehiko Oguchi, Y‐h. Taguchi

Published Jan 1, 1986• 8 citations
Publisher

Abstract

The transfer matrices for the ferromagnetic Ising model on the two- and three-dimensional lattices with finite size are formulated using the screw method introduced by Kramers and Wannier, and the largest and next largest eigenvalues of matrices are obtained by a computer. Then the specific heats and magnetic susceptibilities are calculated for each system with various sizes. By use of the finite-size scaling and Roomany-Wyld approximations, the critical exponents of the infinite-size systems are evaluated. They are in quite good agreement with the reliable ones.

DOI
10.1143/ptps.87.23

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

Approximation method for spin-½ Ising models

S. Gartenhaus

An approximation for calculating the thermodynamic properties of spin-\textonehalf{} Ising models is described. The method is based on the formulation of Kramers and Wannier for the square lattice, and involves expressing the free energy in terms of the limiting maximum eigenvalue of a matrix sequence which represents the contribution from a single spin of the lattice. The application of this formulation to various two- and three-dimensional lattices shows that many of these are characterized by an identically structured matrix sequence with exactly two nonvanishing elements in each row and column. To explore the potential of using this similarity and its simplicity for analyzing Ising models we consider the two-dimensional square lattice in an external field and approximate the doubly infinite lattice by a semi-infinite one consisting of $n (\ensuremath{\le}10)$ infinitely long rows. Calculations are carried out and compared with (1) Onsager's formula for the free energy, (2) the Onsanger-Yang formula for the spontaneous magnetization, and (3) values calculated from series. Except for the region around the critical point where the correlation length exceeds $n$ lattice spacings, substantial agreement---to better than 12 places of decimals in some cases---is obtained. As an additional measure of the potential of this method, a value for the critical exponent $\ensuremath{\delta}$ is computed by the extrapolation of the magnetization along the critical isotherm. The result is found to be in excellent agreement with the known value for $\ensuremath{\delta}$.

Physical review. B, Condensed matter•Feb 1, 1983•10 citations

Finite Size Effects and Pseudocritical Behavior. II: Potts Spin Model

Shingo Ohya, Hiroshi Katsumori

Using the transfer matrices constructed by Kramers and Wannier's successive screw method for the 2-d (l × 蜴) and 3-d (l ×m ×蜴) q state Potts model (2 蠆q 蠆7), we examine the extent to which the presently available supercomputers can provide us with finite size effects and pseudocritical properties of the systems. We estimate the bulk coupling parameters Kc from the specific heat anomaly in two and three dimensions and also from the correlation length anomaly in two dimensions. In two dimensions, our results Kc are in good agreement with those given by the Potts duality formula. In three dimensions, however, we need extended precision computations to obtain the numerical results with sufficient accuracy. We regard the behavior of the exponent µ, which characterizes the relation between the specific heat peak Cv0(l) and the lattice size l, as an indicator to discriminate the order of the phase transition by investigating Cv0(l) for several small values of l.

Progress of Theoretical Physics•May 1, 1996•1 citation

Calculation of the Critical Temperature for the Anisotropic Two-Layer Ising Model Using the Transfer Matrix Method

Mehrdad Ghaemi, Mohammad Ghannadi, Behrouz Mirza

A new finite-size scaling approach based on the transfer matrix method is developed to calculate the critical temperature of an anisotropic two-layer Ising ferromagnet on strips of r -wide sites of square lattices. The reduced internal energy per site has been accurately calculated for the ferromagnetic case, with the nearest-neighbor couplings K x and K y (where K x and K y are the nearest-neighbor interactions within each layer in the x and y directions, respectively) and with interlayer coupling K z, using different size-limited lattices. The calculated energies for different lattice sizes intersect at various points when plotted versus the reduced temperature. It is found that the location of the intersection point versus the lattice size can be fitted to a power series in terms of the lattice sizes. The power series is used to obtain the critical temperature of the unlimited two-layer lattice. The results obtained are in good agreement with the accurate values reported by others.

The Journal of Physical Chemistry B•Dec 18, 2002•12 citations

Finite-size scaling in the three-dimensional Ising model

Michael N. Barber, Robert B. Pearson, D. Toussaint +1 more

We have performed extensive studies of finite-size scaling in the three-dimensional Ising model, using a specially designed Ising processor. We show results for the susceptibility, the energy and specific heat, and the connected part of the fourth moment of the magnetization. Corrections to finite-size scaling are found to be small only for lattices of linear dimension greater than about L=24. This effect severely limits the accuracy to which critical parameters can be estimated. We estimate ${K}_{c}$=0.221 650\ifmmode\pm\else\textpm\fi{}0.000 005 for the critical coupling, while exponent estimates are consistent with those obtained by other methods.

Physical review. B, Condensed matter•Aug 1, 1985•134 citations

Finite-size-scaling estimation of critical exponents in the two-dimensional axial next-nearest-neighbor Ising model

Marcelo D. Grynberg, Horacio Ceva

The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.

Physical review. B, Condensed matter•Nov 1, 1988•6 citations
The transfer matrices for the ferromagnetic Ising model on the two- and three-dimensional lattices with finite size are formulated using the screw method introduced by Kramers and Wannier, and the largest and next largest eigenvalues of matrices are obtained by a computer. Then the specific heats and magnetic susceptibilities are calculated for each system with various sizes. By use of the finite-size scaling and Roomany-Wyld approximations, the critical exponents of the infinite-size systems are evaluated. They are in quite good agreement with the reliable ones.
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.
An approximation for calculating the thermodynamic properties of spin-\textonehalf{} Ising models is described. The method is based on the formulation of Kramers and Wannier for the square lattice, and involves expressing the free energy in terms of the limiting maximum eigenvalue of a matrix sequence which represents the contribution from a single spin of the lattice. The application of this formulation to various two- and three-dimensional lattices shows that many of these are characterized by an identically structured matrix sequence with exactly two nonvanishing elements in each row and column. To explore the potential of using this similarity and its simplicity for analyzing Ising models we consider the two-dimensional square lattice in an external field and approximate the doubly infinite lattice by a semi-infinite one consisting of $n (\ensuremath{\le}10)$ infinitely long rows. Calculations are carried out and compared with (1) Onsager's formula for the free energy, (2) the Onsanger-Yang formula for the spontaneous magnetization, and (3) values calculated from series. Except for the region around the critical point where the correlation length exceeds nnn lattice spacings, substantial agreement---to better than 12 places of decimals in some cases---is obtained. As an additional measure of the potential of this method, a value for the critical exponent $\ensuremath{\delta}$ is computed by the extrapolation of the magnetization along the critical isotherm. The result is found to be in excellent agreement with the known value for $\ensuremath{\delta}$.
Using the transfer matrices constructed by Kramers and Wannier's successive screw method for the 2-d (l × 蜴) and 3-d (l ×m ×蜴) q state Potts model (2 蠆q 蠆7), we examine the extent to which the presently available supercomputers can provide us with finite size effects and pseudocritical properties of the systems. We estimate the bulk coupling parameters Kc from the specific heat anomaly in two and three dimensions and also from the correlation length anomaly in two dimensions. In two dimensions, our results Kc are in good agreement with those given by the Potts duality formula. In three dimensions, however, we need extended precision computations to obtain the numerical results with sufficient accuracy. We regard the behavior of the exponent µ, which characterizes the relation between the specific heat peak Cv0(l) and the lattice size l, as an indicator to discriminate the order of the phase transition by investigating Cv0(l) for several small values of l.
A new finite-size scaling approach based on the transfer matrix method is developed to calculate the critical temperature of an anisotropic two-layer Ising ferromagnet on strips of r -wide sites of square lattices. The reduced internal energy per site has been accurately calculated for the ferromagnetic case, with the nearest-neighbor couplings K x and K y (where K x and K y are the nearest-neighbor interactions within each layer in the x and y directions, respectively) and with interlayer coupling K z, using different size-limited lattices. The calculated energies for different lattice sizes intersect at various points when plotted versus the reduced temperature. It is found that the location of the intersection point versus the lattice size can be fitted to a power series in terms of the lattice sizes. The power series is used to obtain the critical temperature of the unlimited two-layer lattice. The results obtained are in good agreement with the accurate values reported by others.
We have performed extensive studies of finite-size scaling in the three-dimensional Ising model, using a specially designed Ising processor. We show results for the susceptibility, the energy and specific heat, and the connected part of the fourth moment of the magnetization. Corrections to finite-size scaling are found to be small only for lattices of linear dimension greater than about L=24. This effect severely limits the accuracy to which critical parameters can be estimated. We estimate Kc{K}_{c}Kc​=0.221 650\ifmmode\pm\else\textpm\fi{}0.000 005 for the critical coupling, while exponent estimates are consistent with those obtained by other methods.
The magnetic critical exponent of the two-dimensional axial next-nearest-neighbor Ising model is estimated by means of the finite-size scaling of the order parameter. For a finite system without ordering field, this magnitude is computed from the two leading eigenvectors of the transfer matrix. The ferromagnetic-paramagnetic transition temperatures estimated with this method exhibit a small lattice-size dependence. Moreover, they are in good agreement with those obtained using phenomenological renormalization, and also with analytical approximations. The correlation-length critical exponent is calculated by means of phenomenological renormalization. The results found with these numerical techniques indicate that the ferromagnetic-paramagnetic transition of the model seems to be continuous, in contradiction with a previous result obtained using the persistence-length criterion, which indicated that the transition is always first order. We found that the phase transition is first order only when the ratio of the next-nearest-neighbor and nearest-neighbor interactions equals one-half.