The two-dimensional Ising model with a slit is studied. The slit free energy is defined, in which the bulk term, edge terms, and corner terms other than that of the slit are canceled. The bond propagation algorithm is used to calculate the free energy, internal energy, and heat capacity numerically. At the critical point, the slit free energy has a logarithmic correction L^{-1}lnL, in addition to the correction lnL, where L is the typical length of the system. The fitted coefficients of these terms agree with the conformal field theory very well. The integral constant and extrapolation length, which are pending in the conformal field theory, are obtained. In the slit internal energy, the leading term is proportional to L. There is a correction term, ln^{2}L, that has a geometry-independent coefficient, the exact value of which is conjectured to be -sqrt[2]/(2π^{2}). The leading term in the slit heat capacity is proportional to L^{2}. There is also a ln^{2}L term; however, its coefficient depends on the geometry. In the infinitely long strip limit, its exact value is conjectured to be 5/π^{2}. Near the critical point, for different size lattices, the rescaled slit free energies, internal energies, and heat capacities collapse onto the corresponding universal curves. According to the scaling and conformal field theory, the leading term of the slit free energy should be -c/8ln|t| for L|t|≫1 and |t|≪1, where t is the reduced temperature. The amplitude c/8 is universal, and c is the central charge. Our results validate these predictions qualitatively. For L|t|≫1 and |t|≪1, the slit internal energy and heat capacity scale as t^{-1} and t^{-2}, respectively.
The critical 2-dimensional Ising model is studied with four types boundary conditions: free, fixed ferromagnetic, fixed antiferromagnetic, fixed double antiferromagnetic. Using Bond Propagation algorithms with surface fields, we obtained the free energy, internal energy and specific heat numerically on square lattices with square shape and various combinations of the four types boundary conditions. The numerical data are analyzed with finite size scaling. The bulk, edge and corner terms are extracted very accurately. The exact results are conjectured for the corner logarithmic term in the free energy, the edge and corner logarithmic terms in the internal energy and specific heat. The corner logarithmic terms in the free energy agree with the conformal field theory very well.
Using a bond-propagation algorithm, we study the finite-size behavior of the critical two-dimensional Ising model on a finite triangular lattice with free boundaries in five shapes: triangular, rhomboid, trapezoid, hexagonal, and rectangular. The critical free energy, internal energy, and specific heat are calculated. The accuracy of the free energy reaches 10(-26). Based on accurate data on several finite systems with linear size up to N=2000, we extract the bulk, surface, and corner parts of the free energy, internal energy, and specific heat accurately. We confirm the conformal field theory prediction that the corner free energy is universal and find logarithmic corrections in higher-order terms in the critical free energy for the rhomboid, trapezoid, and hexagonal systems, which are absent for the triangular and rectangular systems. The logarithmic edge corrections due to edges parallel or perpendicular to the bond directions in the internal energy are found to be identical, while the logarithmic edge corrections due to corresponding edges in the specific heat are different. The corner internal energy and corner specific heat for angles π/3, π/2, and 2π/3 are obtained, as well as higher-order corrections. Comparing with the corner internal energy and corner specific heat we previously found on a rectangle of the square lattice [Phys. Rev. E 86, 041149 (2012)], we conclude that the corner internal energy and corner specific heat for the rectangular shape are not universal.
The critical two-dimensional Ising model is studied with four types boundary conditions: free, fixed ferromagnetic, fixed antiferromagnetic, and fixed double antiferromagnetic. Using bond propagation algorithms with surface fields, we obtain the free energy, internal energy, and specific heat numerically on square lattices with a square shape and various combinations of the four types of boundary conditions. The calculations are carried out on the square lattices with size N×N and 30<N<1000. The numerical data are analyzed with finite-size scaling. The bulk, edge, and corner terms are extracted very accurately. The exact results are conjectured for the corner logarithmic term in the free energy, the edge logarithmic term in the internal energy, and the corner logarithmic term in the specific heat. The corner logarithmic terms in the free energy agree with the conformal field theory very well.
Using the bond-propagation algorithm, we study the Ising model on a rectangle of size $M\ifmmode\times\else\texttimes\fi{}N$ with free boundaries. For five aspect ratios, $\ensuremath{\rho}=M/N=1$, 2, 4, 8, and 16, the critical free energy, internal energy and specific heat are calculated. The largest size reached is $M\ifmmode\times\else\texttimes\fi{}N=64\ifmmode\times\else\texttimes\fi{}{10}^{6}$. The accuracy of the free energy reaches ${10}^{\ensuremath{-}26}$. Based on these accurate data, we determine exact expansions of the critical free energy, internal energy, and specific heat. With these expansions, we extract the bulk, surface, and corner parts of the free energy, internal energy, and specific heat. The fitted bulk free energy density is given by ${f}_{\ensuremath{\infty}}=0.92969539834161021499(1)$, compared with Onsager's exact result ${f}_{\ensuremath{\infty}}=0.929695398341610214985\ensuremath{\cdots}.$ We confirm the conformal field theory (CFT) prediction of the corner free energy, in which the central charge of the Ising model is found to be $c=0.5\ifmmode\pm\else\textpm\fi{}1\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}10}$, compared with the CFT result
c=0.5. We find that not only the corner free energy but also the corner internal energy and specific heat are geometry independent, i.e., independent of aspect ratio. The implication of this finding for finite-size scaling is discussed. In the second-order correction of the free energy, we confirm the geometry dependence predicted by CFT and determine a geometry-independent constant beyond CFT. High-order corrections are also obtained.
For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with N = l × l spins and the asymptotic Onsager solution ( N → ∞). We derived empirical expressions for critical parameters as functions of N and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when N increased the heat capacity in the critical point increased as ln N . We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis Monte Carlo method.
We test the Gluing Method on the square-lattice two-dimensional Ising model by reconstructing global free-energy generators from complementary local relations. The emphasis is kept explicit throughout: the inputs used as \emph{gluing material} are separated from quantities obtained \emph{after gluing} and from data reserved for validation. At zero magnetic field, Kramers--Wannier duality is absorbed into the coordinate
τ=logsinh2K, so that the reduced free energy decomposes as
Φ(τ)=τ/2+G(∣τ∣). The gluing materials are then only a critical logarithmic local expansion and a low-order high-temperature outer expansion; the low-temperature sector is generated by duality rather than fitted independently. Positive flat weights preserve the prescribed critical asymptotics to every finite algebraic order. With one crossover parameter fixed from three scalar anchors, the resulting bulk free energy agrees with the exact Onsager benchmark at a maximum relative error of
1.915×10−5 on the tested interval. Internal energy and specific heat are derived from the same glued generator. Zero-field finite-size corrections are organized separately as one-variable scaling functions multiplying
L−2,L−4,…. At nonzero field in the thermodynamic limit, the gluing materials are the established high-temperature, critical, and low-temperature local expansions of the universal Ising magnetic scaling free energy
Φmag(η),
η=m/∣h∣8/15. The low-temperature input contains the
∣h∣ cusp associated with spontaneous symmetry breaking. Mapping the glued scaling free energy back to the square lattice through nonlinear scaling fields yields a single generator for the magnetic free-energy increment, magnetization, and susceptibility. The construction reproduces high-temperature linear response, the critical isotherm
M∝∣H∣1/15, and the two one-sided limits
M(H→0±)=±M0 without separately fitting magnetization branches. Finally, finite size and finite field are combined in the scaling variables
xt=Lm and
xh=L15/8h. The two protected one-dimensional axes, a central mixed-response expansion, a low-temperature spectral-rounding relation, and a massive outer description are glued into a two-variable scaling surface
Yglue(xt,xh). The resulting lattice free energy closes in the compact form
ΦLglue(K,H)=Φlatglue(K,H)+L2Yglue(Lm,L15/8h)−Ybulk(Lm,L15/8h)+O(L−4) where
m=m(r,H) and
h=h(r,H) are the nonlinear scaling fields and
Ybulk is the thermodynamic singular scaling surface expressed in the same variables. Thus all finite-size/finite-field information is carried by the excess
Yglue−Ybulk, while the thermodynamic lattice free energy remains the previously constructed
Φlatglue. On eleven transfer-matrix points withheld from the construction, the blind free-energy test gives RMS absolute error
3.90×10−4 and maximum absolute error
9.55×10−4 in
Y. Differentiation of the same free energy produces finite-size magnetization and susceptibility; magnetization is especially stable on the low-temperature side, while the second derivative is more sensitive in crossover regions.