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.
In a previous paper of this series, we studied $\mathfrak{P}(\overline{\ensuremath{\sigma}})$, the probability that the average boundary spin for a half-plane of Ising spins is $\overline{\ensuremath{\sigma}}$. We extend our study of such functions by generalizing the previous derivation to any magnetic system and using them to relate the magnetic behavior at the critical isotherm to the spin-spin correlations at the critical temperature. We also clarify the meaning of the secondary maxima previously found in $\mathfrak{P}(\overline{\ensuremath{\sigma}})$ by a more accurate calculation and by examining the magnetization on interior rows of the half-plane. These considerations show that $\mathfrak{P}(\overline{\ensuremath{\sigma}})$ is more properly to be considered as the sum of two separate spin probability functions. The various ways in which the thermodynamic limit may be taken and the effects that these different limiting procedures have on $\mathfrak{P}(\overline{\ensuremath{\sigma}})$ are discussed in detail. The influence of boundary conditions on spin probability functions is studied by reversing the sign of one column of horizontal bonds. We compute the additional free energy resulting from the misfit bonds and use this to show that $\mathfrak{P}(\overline{\ensuremath{\sigma}})$ is no longer bimodal but is rectangular for $|\overline{\ensuremath{\sigma}}|$ less than the boundary spontaneous magnetization. Finally, we present in graphical form numerical integrations of the boundary magnetization and 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.
We study in detail the asymptotic behavior, for large separations, of the correlation between two spins in the same row of the two-dimensional Ising model with the particular value of a purely imaginary magnetic field $H=\frac{1}{2}i\ensuremath{\pi}\mathrm{kT}$. We derive the partition function and the magnetization per spin first proposed by Lee and Yang and find that the correlation function approaches its limiting value exponentially. In contrast to the case of no magnetic field, the approach is shown to be oscillatory rather than monotonic. The first three terms of the asymptotic behavior are explicitly given. The correlation function is obtained by expressing it as a block Toeplitz determinant whose behavior for large separation is studied by solving a pair of coupled Wiener-Hopf sum equations.
We compute exactly the spin-spin correlation functions $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉$ for the two-dimensional Ising model on a square lattice in zero magnetic field for
T>Tc and
T<Tc. We then analyze the correlation functions in the scaling limit $T\ensuremath{\rightarrow}{T}_{c},{M}^{2}+{N}^{2}\ensuremath{\rightarrow}\ensuremath{\infty}$ such that $(T\ensuremath{-}{T}_{c})$ is fixed. In this scaling limit $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉={R}^{\ensuremath{-}\frac{1}{4}}{F}_{\ifmmode\pm\else\textpm\fi{}}(t)+{R}^{\ensuremath{-}\frac{5}{4}}{F}_{1\ifmmode\pm\else\textpm\fi{}}(t)+o({R}^{\ensuremath{-}\frac{5}{4}})$, where
t is the scaling variable $\frac{R}{\ensuremath{\xi}}$ and ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ and ${F}_{1\ifmmode\pm\else\textpm\fi{}}(t)$ are the scaling functions ($\ensuremath{\xi}$ is the correlation length). We derive exact expressions for these scaling functions, in terms of a Painlev\'e function of the third kind and analyze both the small- and large-
t behavior. A table of values for ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ (good to ten significant digits) is also given. As an application we computer the coefficients ${C}_{0\ifmmode\pm\else\textpm\fi{}}$ and ${C}_{1\ifmmode\pm\else\textpm\fi{}}$ in the expansion ${k}_{B}T\ensuremath{\chi}(T)={C}_{0\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{7}{4}}+{C}_{1\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{3}{4}}+O(1)$ of the zero-field susceptibility $\ensuremath{\chi}(T)$ as $T\ensuremath{\rightarrow}{T}_{c}^{\ifmmode\pm\else\textpm\fi{}}$.
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 partition functions of the two- and three-dimensional finite Ising models in the presence of a magnetic field have been calculated with use of a high-speed computer. Periodic boundary conditions were imposed. All the numbers of configurations for the 3×3×3, 5×5, and 4×6 lattices are tabulated. The distribution of zeros of the partition functions in the complex magnetic-field plane has been obtained both for the ferro- and antiferromagnetic cases. The zeros in the ferromagnetic cases are distributed on the unit circle according to the Yang-Lee theorem, and the dependence on temperature of the distribution functions has been studied. In the antiferromagnetic case most of zeros are on the negative real axis, but a certain number of complex roots appear in the left half of the plane. The magnetization and the magnetic susceptibility have been calculated as functions of temperature.