Using previous results from boundary conformal field theory and integrability, a phase diagram is derived for the two-dimensional Ising model at its bulk tri-critical point as a function of boundary magnetic field and boundary spin-coupling constant. A boundary tri-critical point separates phases where the spins on the boundary are ordered or disordered. In the latter range of coupling constant, there is a non-zero critical field where the magnetization is singular. In the former range, as the temperature is lowered, the boundary undergoes a first-order transition while the bulk simultaneously undergoes a second-order transition.
Previous series expansion studies of the magnetization and susceptibility of the Ising model are extended to the higher derivatives of the free energy with respect to the magnetic field. It is found that just above the critical temperature the 2 n th derivative is given by where in two dimensions γ = 1¾, Δ = 1fraction seven-eighths and in three dimensions γ = 1¼, Δ = 1fraction nine-sixteens. In two dimensions γ is exact and Δ is correct to 0.6%; the three-dimensional indices are correct to 0.2%. Just below the critical temperature where in three dimensions 0.303 <or= β <or= 0.318 and 1·55 <or= Δprime <or= 1.65. If the even derivatives are to have the same index above and below T c , then β = fraction five-sixteens and Δprime = 1fraction nine-sixteens. The latter index lies very close to the lower limit so that, although we cannot rule out Δprime = Δ, Δprime = 1fraction five-eighths is not unlikely. In two dimensions our results are consistent with linearity of the index and, if this is assumed, then β = fraction one-eighths and Δprime = 1fraction seven-eighths exactly and there is symmetry about T c . The amplitudes C 2 n + and C n - are estimated for n = 1 to 6.
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.
First, a new method proposed by Baker and Rushbrooke is used to study the simple ferromagnetic Ising model at and below the Curie temperature. Of course, the properties of the Ising model are already well known, so that the main aim here is to assess the potential and reliability of the new method, since it has wide applicability to other models which have not been otherwise studied. Between 8 and 16 coefficients of exact high-temperature expansions for fixed values of the magnetization are derived for various two- and three-dimensional lattices. A Pad\'e-approximant analysis of these expansions at the critical isotherm and magnetic phase boundary enables us to estimate the critical exponents $\ensuremath{\beta}$, ${\ensuremath{\gamma}}^{\ensuremath{'}}$, and $\ensuremath{\delta}$, and plot the spontaneous magnetization. The results are in good agreement with previous calculations. Secondly, an analysis of the exact series expansions provides no support for the conjecture that the phase boundary is a line of essential singularities. However, the same expansions strongly suggest the existence of a "spinodal" curve, whose properties are in reasonable agreement with the predictions of various heuristic arguments (based essentially upon analyticity at the phase boundary and one-phase homogeneity in the critical region). Finally, structure and a mild extension of the proven analyticity of the free energy are used to show the $\ensuremath{\Delta}\ensuremath{\le}{\ensuremath{\Delta}}^{\ensuremath{'}}$, $\ensuremath{\gamma}\ensuremath{\le}{\ensuremath{\gamma}}^{\ensuremath{'}}$.
We consider a semi-infinite two-dimensional Ising model with its spins on the boundary row having a different interaction energy E′1 from the ferromagnetic bulk. We find that the boundary specific heat has two divergent terms: one of which diverges linearly at the bulk critical temperature Tc, and the other, logaritmically. The linearly divergent term is independent of E′1, and the coefficient of the logaritmically divergent term is a decreasing function of E′1. There is a boundary latent heat at Tc, which is identical to McCoy and Wu's result. The boundary spins, which can be either ferromagnetic or antiferromagnetic, are aligned for temperature lower than Tc. The boundary spontaneous magnetization approaches zero in the form of A(E′1)(1−T/Tc)1/2, and the boundary zero field magnetic susceptibility diverges at Tc in the form −B (E′1)ln|1−T/Tc|, where A(E′1) and B(E′) are increasing functions of E′1.