We study some properties of the Ising model in the plane of the complex (energy/temperature)-dependent variable
u=e−4K, where
K=J/(kBT), for nonzero external magnetic field,
H. Exact results are given for the phase diagram in the
u plane for the model in one dimension and on infinite-length quasi-one-dimensional strips. In the case of real
h=H/(kBT), these results provide new insights into features of our earlier study of this case. We also consider complex
h=H/(kBT) and
μ=e−2h. Calculations of complex-
u zeros of the partition function on sections of the square lattice are presented. For the case of imaginary
h, i.e.,
μ=eiθ, we use exact results for the quasi-1D strips together with these partition function zeros for the model in 2D to infer some properties of the resultant phase diagram in the
u plane. We find that in this case, the phase boundary
Bu contains a real line segment extending through part of the physical ferromagnetic interval
0≤u≤1, with a right-hand endpoint
urhe at the temperature for which the Yang-Lee edge singularity occurs at
μ=e±iθ. Conformal field theory arguments are used to relate the singularities at
urhe and the Yang-Lee edge.
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.
Although the physical properties of the 2D and 1D Ising models are quite different, we point out an interesting connection between their complex-temperature phase diagrams. We carry out an exact determination of the complex-temperature phase diagram for the 1D Ising model for arbitrary spin s and show that in the u_s=e^{-K/s^2} plane (i) it consists of N_{c,1D}=4s^2 infinite regions separated by an equal number of boundary curves where the free energy is non-analytic; (ii) these curves extend from the origin to complex infinity, and in both limits are oriented along the angles \\theta_n = (1+2n)\\pi/(4s^2), for n=0,..., 4s^2-1; (iii) of these curves, there are N_{c,NE,1D}=N_{c,NW,1D}=[s^2] in the first and second (NE and NW) quadrants; and (iv) there is a boundary curve (line) along the negative real u_s axis if and only if s is half-integral. We note a close relation between these results and the number of arcs of zeros protruding into the FM phase in our recent calculation of partition function zeros for the 2D spin s Ising model.