The Ising model at magnetic field H=0 is one of the most important exactly solved models in statistical mechanics, solved on the square lattice in 1944 by Onsager, and later by others on the triangular lattice; its magnetic susceptibility at H=0 continues to be an unsolved aspect of the model, however. The susceptibility can either be viewed as a sum over all correlation functions of the integrable model at H=0, or else as a second derivative of the free energy of the non-integrable Ising model in a field. Therefore, its analytic properties, though derived from the non-integrable model, can be studied through series expansion of the integrable correlation functions. We begin this process by analyzing the first four terms in the form factor expansion of the diagonal correlation functions, and after summing over the diagonal form factor expansion, the first four terms of the diagonal susceptibility expansion. We have been able to reduce the form factor and susceptibility expansion terms, given as multi-dimensional integrals, to closed-form functions in all cases. Under the limit of H going to infinity with the interaction energy E going to minus infinity, the isotropic Ising model becomes a hard particle lattice gas model, with exclusion of nearest neighbor lattice sites. In this limit, the triangular lattice Ising model becomes the exactly solved hard hexagon model, while the square lattice becomes the non-integrable hard squares model. We study in detail these models for finite lattice sizes in order to understand the differences between integrable and non-integrable models. We consider the partition functions zeros and their density for different boundary conditions, and find notable differences in the density which is attributed to an extra factorization in the transfer matrices of hard hexagons which is absent in hard squares. We also study the special point at fugacity z=-1 of hard squares where all eigenvalues of the transfer matrix are equimodular and where the grand partition function's value depends on boundary conditions. | 300 pages
The Ising model at magnetic field H=0 is one of the most important exactly solved models in statistical mechanics, solved on the square lattice in 1944 by Onsager, and later by others on the triangular lattice; its magnetic susceptibility at H=0 continues to be an unsolved aspect of the model, however. The susceptibility can either be viewed as a sum over all correlation functions of the integrable model at H=0, or else as a second derivative of the free energy of the non-integrable Ising model in a field. Therefore, its analytic properties, though derived from the non-integrable model, can be studied through series expansion of the integrable correlation functions. We begin this process by analyzing the first four terms in the form factor expansion of the diagonal correlation functions, and after summing over the diagonal form factor expansion, the first four terms of the diagonal susceptibility expansion. We have been able to reduce the form factor and susceptibility expansion terms, given as multi-dimensional integrals, to closed-form functions in all cases. Under the limit of H going to infinity with the interaction energy E going to minus infinity, the isotropic Ising model becomes a hard particle lattice gas model, with exclusion of nearest neighbor lattice sites. In this limit, the triangular lattice Ising model becomes the exactly solved hard hexagon model, while the square lattice becomes the non-integrable hard squares model. We study in detail these models for finite lattice sizes in order to understand the differences between integrable and non-integrable models. We consider the partition functions zeros and their density for different boundary conditions, and find notable differences in the density which is attributed to an extra factorization in the transfer matrices of hard hexagons which is absent in hard squares. We also study the special point at fugacity z=-1 of hard squares where all eigenvalues of the transfer matrix are equimodular and where the grand partition function's value depends on boundary conditions. | 300 pages
We use the recently derived form factor expansions of the diagonal two-point correlation function of the square Ising model to study the susceptibility for a magnetic field applied only to one diagonal of the square lattice, for the isotropic Ising model. We exactly evaluate the one and two particle contributions χ (1) d and χ (2) d of the corresponding susceptibility, and obtain linear differential equations for the three and four particle contributions, as well as the five particle contribution χ (5) d , but only modulo a given prime. We use these exact linear differential equations to show that not only the Russian-doll structure but also the direct sum structure on the linear differential operators for the n -particle contributions χ ( n ) d are quite directly inherited from the direct sum structure on the form factors f ( n ) . We show that the n th particle contributions χ ( n ) d have their singularities at roots of unity. These singularities become dense on the unit circle |sinh 2 E v / kT sinh 2 E h / kT | = 1 as n → ∞ .
The partition functions of ferromagnetic Ising models of square lattices in a finite magnetic field is deduced using topological considerations within a heuristic graph-theoretical approach. These equations are derived separately for low and high temperature regimes while the exact solution of Onsager is obtained therefrom when the magnetic field is zero. The derived partition function equations here are almost similar to those given by Onsager, thus indicating a straight-forward protocol, even when the magnetic field is present. The spontaneous magnetization derived here using the Helmholtz free energy is identical with that arising from the exact solution. The partition functions lead to the known series expansions of the magnetization and zero-field susceptibility.
We consider the rectangular Ising model on a half-plane of infinite extent and study some of the consequences connected with the presence of the boundary. Only the spins on the boundary row are allowed to interact with a magnetic field H. The method of Pfaffians is employed to obtain exact expressions for the partition function. It is found that the free energy is the sum of two terms, one of which is independent of H and proportional to the total number of lattice sites, while the other depends on H and is proportional to the number of lattice sites on the boundary. This separation makes it possible to define various thermodynamic quantities associated with the boundary. In particular, the boundary magnetization is shown to be discontinuous, in the ferromagnetic case, at zero magnetic field for temperatures below the bulk critical temperature Tc. This discontinuity, which is the spontaneous boundary magnetization, goes to zero as ${(1\ensuremath{-}\frac{T}{{T}_{c}})}^{\frac{1}{2}}$ as $T\ensuremath{\rightarrow}{T}_{c}\ensuremath{-}$. For T=Tc, the discontinuity is of course absent, and the boundary magnetization behaves as $\ensuremath{-}\mathfrak{H}\mathrm{ln}\mathfrak{H}$ for small H. The boundary susceptibility at zero magnetic field in the ferromagnetic case exhibits a logarithmic singularity at T=Tc, both above and below transition. An interesting feature is that the ferromagnetic boundary magnetization, although discontinuous for T<Tc, may be analytically continued beyond the point H=0. We interpret this as a hystersis phenomenon which we study in detail by computing the probability distribution function for the average boundary spin. The correlation function for two spins, both on the boundary row, is also obtained exactly and its asymptotic behavior is given. Finally, we derive an expression for the magnetization in any row and explicitly evaluate it for the second row, i.e., the row next to the boundary.
Abstract This chapter summarizes the many results of the exactly solved Ising model in two dimensions at H = 0 including: the partition function, free energy, spontaneous magnetization, the Toeplitz determinant representations of the correlation functions, the Painlevé VI differential equation for the diagonal two spin correlations, the form factor expansion of the magnetic susceptibility, the scaling limit, the boundary properties of the half plane lattice, and the random layered lattice. For H ≠ 0 the Lee-Yang circle theorem, the free energy and magnetization for H/kBT = iπ/2 expansions for small H, the results of Zamolodchikov at T = Tc and a discussion of extended analyticity are given.
For various Ising models two approaches are discussed, one is that of simulating lattices, also called gauging on exact equations, the other is that of calculating analytical expressions for the boundary free energy of Ising lattices. The first approach allows to conjecture a solution for some Ising models, that have sofar not been solved, once some exact partial result for the problem is known. The second approach aims at furnishing such a partial result in the form of a condition for the critical temperature. An example of such a result was recently given for the 2D Ising square lattice with nearest and next-nearest-neighbor interactions. The critical line that separates the ordered (ferromagnetic) phase from the disordered (paramagnetic) phase showed good agreement in the moderate and strong nearest neighbor coupling limit with several results obtained by Monte Carlo, transfer matrix and series expansion results. We extend the discussion of the critical line, finding an excellent fit, now also in other points, like the Padé point, as well as cusp behavior at the Onsager point where the lattice decouples into two 2D square lattices with only nearest-neighbor interaction. Combination of this result with a geometrical argument in the simulation approach leads to a critical exponent 2−2≈0.5858, comparable to the exponent 4/7≈0.5714 found from renormalization arguments.
An example of a two-dimensional crystal in which each lattice point interacts with six neighbors is the so-called \triangular Ising lattice. In this work the transfer matrix method is used to nd the thermodynamic properties of the triangular Ising model with a limited number of rows, in the absence or presence of a magnetic eld. In the absence of a magnetic eld the results have been used to nd an analytical expression for the partition function of a model with any number of rows, from which the thermodynamic properties of the triangular planar lattice model are obtained as a polynomial in terms of x = e 2j tanh 2j=2 cosh2j, where J = jkT is the interaction energy of the two nearest neighbor spins. The degree of the polynomial depends on the number of rows in the model. The results are compared with the exact values. The calculated free energy is found to be very accurate, but the calculated internal energy and heat capacity show some small deviations from the exact values. In the presence of a magnetic eld the solution to the matrices is too complicated for a general analytical expression to be proposed for a model with any number of rows. However, our exact solution for the model with n = 4, 5, and 7 reveals an important point: that such results are independent of n when the eld-spin interaction energy is almost equal to or larger than that of the spin-spin interaction.