Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Mathematical Physics

Zeros of the Partition Function for the Heisenberg, Ferroelectric, and General Ising Models

Masuo Suzuki, Michael Ellis Fisher

Published Feb 1, 1971• 164 citations
Publisher

Abstract

The Lee-Yang theorem for the zeros of the partition function of a ferromagnetic Ising model with real pair spin interactions is extended to general Ising models with complex many-spin interactions (satisfying appropriate ``ferromagnetic'' and spin inversion symmetry conditions). When many-spin interactions are present, all zeros lie on the imaginary Hz-axis for sufficiently low (but fixed) T, but, in general, some leave the imaginary axis as T → ∞. The extended Ising theorem is used to prove the same result for a Heisenberg system of arbitrary spin with the real anisotropic pair interaction Hamiltonian Hij=−(JijzSizSjz+JijxSixSjx+JijySiySjy)in an arbitrary transverse field (Hx, Hy) under the ``ferromagnetic'' condition Jijz≥|Jijx| and Jijz≥|Jijy|. The analyticity of the limiting free energy of such a Heisenberg ferromagnet and the absence of a phase transition are thereby established for all (real) nonzero magnetic fields Hz. The Ising theorem is also applied to hydrogen-bonded ferroelectric models to prove, in particular, that the zeros for the KDP model lie on the imaginary electric field axis for all T below the transition temperature Tc.

DOI
10.1063/1.1665583

Related papers

Theorems on the Ising Model with General Spin and Phase Transition

Masuo Suzuki

The theorem of Lee and Yang has been extended to the ferromagnetic Ising model with arbitrarily mixed spin values of Sj = ½, 1, and 32, including the case of equal spin values as a special one. Namely, it has been proved that the zeros of the partition function for the above Ising model with higher spin values lie on the unit circle in the fugacity plane (or complex magnetic-field plane). Expressions for general correlation functions in Ising ferromagnets with higher spin values have been derived in terms of the above generalized theorem. By the use of these expressions, the relations among the critical indices are discussed and the same results are obtained as those predicted by the scaling-law approach.

Journal of Mathematical Physics•Dec 1, 1968•33 citations

Distribution of Zeros of the Partition Function for the Finite Two-Dimensional Heisenberg Model

Vinod Mubayi, Chanchal K. Majumdar, Kanwar Krishan

The exact excitation spectrum of the 3 \ifmmode\times\else\texttimes\fi{} 3 square and the 3 \ifmmode\times\else\texttimes\fi{} 4 rectangular lattices interacting via the nearest-neighbor anisotropic Heisenberg interaction $H=J\ensuremath{\Sigma}{(\mathrm{ij})}^{}[S_{i}^{}{}_{}{}^{z}S_{j}^{}{}_{}{}^{z}+\ensuremath{\gamma}(S_{i}^{}{}_{}{}^{x}S_{j}^{}{}_{}{}^{x}+S_{i}^{}{}_{}{}^{y}S_{j}^{}{}_{}{}^{y})]$ has been obtained for various values of the anisotropy constant $\ensuremath{\gamma}$. For $0\ensuremath{\le}\ensuremath{\gamma}\ensuremath{\le}1$, the zeros of the partition function in the complex $\ensuremath{\mu}=\mathrm{exp}(\ensuremath{-}\frac{m{H}_{e}}{{k}_{B}T})$ plane obey the generalized Lee-Yang theorem for ferromagnetic coupling. The usual ferromagnetic transition seems to persist up to $\ensuremath{\gamma}\ensuremath{\simeq}0.6$. For $0.6\ensuremath{\lesssim}\ensuremath{\gamma}\ensuremath{\le}1$, nothing definite about a ferromagnetic transition can be inferred from the behavior of the zeros. When $\ensuremath{\gamma}>1$, the ferromagnetic zeros do not obey the generalized Lee-Yang theorem at sufficiently low temperatures. For antiferromagnetic coupling, the zeros lie on the negative real axis for all values of $\ensuremath{\gamma}$ and of the temperature of the studied lattices.

Physical review. B, Solid state•Oct 1, 1973•2 citations

Yang-Lee Zeros of the Antiferromagnetic Ising Model

Seung‐Yeon Kim

There exists the famous circle theorem on the Yang-Lee zeros of the ferromagnetic Ising model. However, the Yang-Lee zeros of the antiferromagnetic Ising model are much less well understood than those of the ferromagnetic model. The precise distribution of the Yang-Lee zeros of the antiferromagnetic Ising model only with nearest-neighbor interaction J on LxL square lattices is determined as a function of temperature a=e(2betaJ) (J<0), and its relation to the phase transitions is investigated. In the thermodynamic limit (L-->infinity), the distribution of the Yang-Lee zeros of the antiferromagnetic Ising model cuts the positive real axis in the complex x=e(-2betaH) plane, resulting in the critical magnetic field +/-H(c)(a), where H(c)>0 below the critical temperature a(c)=square root of 2-1. The results suggest that the value of the scaling exponent y(h) is 1 along the critical line for a<a(c).

Physical Review Letters•Sep 23, 2004•68 citations

On properties of the Ising model for complex energy/temperature and magnetic field

V. Matveev, Robert E. Shrock

We study some properties of the Ising model in the plane of the complex (energy/temperature)-dependent variable $u=e^{-4K}$, where $K=J/(k_BT)$, 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/(k_BT)$, these results provide new insights into features of our earlier study of this case. We also consider complex $h=H/(k_BT)$ 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., $μ=e^{iθ}$, 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 ${\cal B}_u$ contains a real line segment extending through part of the physical ferromagnetic interval $0 \le u \le 1$, with a right-hand endpoint $u_{rhe}$ at the temperature for which the Yang-Lee edge singularity occurs at $μ=e^{\pm iθ}$. Conformal field theory arguments are used to relate the singularities at $u_{rhe}$ and the Yang-Lee edge.

arXiv:0711.4639•Journal of Physics A Mathematical and Theoretical•Mar 17, 2008•23 citations

Distribution of Zeros of the Partition Function of the Ising Model

Shigetoshi Katsura, Yoshihiko Abe, Masami Yamamoto

To study the relations between Yang-Lee's theory of the phase transitions and the distribution of the zeros of the partition functions in the antiferromagnetic case, the partition functions of finite (4×4 and 4×6) Ising model with the nearest and the next nearest neighbor interactions are calculated. The distributions of the zeros in the fugacity z plane of the partition function for various combinations of the values of J / k T - J '/ k T or tanh ( J /2 k T )-tanh ( J '/2 k T ) are investigated. The patterns of the distribution of the zeros show several types corresponding to ferromagnetic state, antiferromagnetic state and superantiferromagnetic state. In the antiferromagnetic ordered state the locus of the zeros is found to be nearly two concentric circles and to cross the positive real axis at two points z c and 1/ z c when the next neighbor ferromagnetic interaction is introduced. This means that our system shows the phase transitions at finite critical fields H c and - H c . The locus of zeros of Fisher's model of super-exchange antiferromagnet is also shown to cross the positive real axis at z c and 1/ z c below the Néel temperature.

Journal of the Physical Society of Japan•Feb 1, 1971•23 citations

Theorems on the Partition Functions of the Heisenberg Ferromagnets

Taro Asano

It is shown that all the zeros of the partition function lie on the unit circle in the complex `fugacity' plane and the first inequality of Griffiths, Kelly and Sherman for the spin correlation functions also holds for the anisotropic Heisenberg ferromagnets as in the Ising model. The isotropic Heisenberg ferromagnets and the Ising ferromagnets are involved as special cases and thus the methods of this article provide alternative proof of the corresponding theorems in the Ising model.

Journal of the Physical Society of Japan•Aug 1, 1970•97 citations
The Lee-Yang theorem for the zeros of the partition function of a ferromagnetic Ising model with real pair spin interactions is extended to general Ising models with complex many-spin interactions (satisfying appropriate ferromagnetic'' and spin inversion symmetry conditions). When many-spin interactions are present, all zeros lie on the imaginary Hz-axis for sufficiently low (but fixed) T, but, in general, some leave the imaginary axis as T → ∞. The extended Ising theorem is used to prove the same result for a Heisenberg system of arbitrary spin with the real anisotropic pair interaction Hamiltonian Hij=−(JijzSizSjz+JijxSixSjx+JijySiySjy)in an arbitrary transverse field (Hx, Hy) under the ferromagnetic'' condition Jijz≥|Jijx| and Jijz≥|Jijy|. The analyticity of the limiting free energy of such a Heisenberg ferromagnet and the absence of a phase transition are thereby established for all (real) nonzero magnetic fields Hz. The Ising theorem is also applied to hydrogen-bonded ferroelectric models to prove, in particular, that the zeros for the KDP model lie on the imaginary electric field axis for all T below the transition temperature Tc.
The theorem of Lee and Yang has been extended to the ferromagnetic Ising model with arbitrarily mixed spin values of Sj = ½, 1, and 32, including the case of equal spin values as a special one. Namely, it has been proved that the zeros of the partition function for the above Ising model with higher spin values lie on the unit circle in the fugacity plane (or complex magnetic-field plane). Expressions for general correlation functions in Ising ferromagnets with higher spin values have been derived in terms of the above generalized theorem. By the use of these expressions, the relations among the critical indices are discussed and the same results are obtained as those predicted by the scaling-law approach.
The exact excitation spectrum of the 3 \ifmmode\times\else\texttimes\fi{} 3 square and the 3 \ifmmode\times\else\texttimes\fi{} 4 rectangular lattices interacting via the nearest-neighbor anisotropic Heisenberg interaction $H=J\ensuremath{\Sigma}{(\mathrm{ij})}^{}[S_{i}^{}{}_{}{}^{z}S_{j}^{}{}_{}{}^{z}+\ensuremath{\gamma}(S_{i}^{}{}_{}{}^{x}S_{j}^{}{}_{}{}^{x}+S_{i}^{}{}_{}{}^{y}S_{j}^{}{}_{}{}^{y})]$ has been obtained for various values of the anisotropy constant $\ensuremath{\gamma}$. For $0\ensuremath{\le}\ensuremath{\gamma}\ensuremath{\le}1$, the zeros of the partition function in the complex $\ensuremath{\mu}=\mathrm{exp}(\ensuremath{-}\frac{m{H}_{e}}{{k}_{B}T})$ plane obey the generalized Lee-Yang theorem for ferromagnetic coupling. The usual ferromagnetic transition seems to persist up to $\ensuremath{\gamma}\ensuremath{\simeq}0.6$. For $0.6\ensuremath{\lesssim}\ensuremath{\gamma}\ensuremath{\le}1$, nothing definite about a ferromagnetic transition can be inferred from the behavior of the zeros. When $\ensuremath{\gamma}>1$, the ferromagnetic zeros do not obey the generalized Lee-Yang theorem at sufficiently low temperatures. For antiferromagnetic coupling, the zeros lie on the negative real axis for all values of $\ensuremath{\gamma}$ and of the temperature of the studied lattices.
There exists the famous circle theorem on the Yang-Lee zeros of the ferromagnetic Ising model. However, the Yang-Lee zeros of the antiferromagnetic Ising model are much less well understood than those of the ferromagnetic model. The precise distribution of the Yang-Lee zeros of the antiferromagnetic Ising model only with nearest-neighbor interaction J on LxL square lattices is determined as a function of temperature a=e(2betaJ) (J<0), and its relation to the phase transitions is investigated. In the thermodynamic limit (L-->infinity), the distribution of the Yang-Lee zeros of the antiferromagnetic Ising model cuts the positive real axis in the complex x=e(-2betaH) plane, resulting in the critical magnetic field +/-H(c)(a), where H(c)>0 below the critical temperature a(c)=square root of 2-1. The results suggest that the value of the scaling exponent y(h) is 1 along the critical line for a<a(c).
We study some properties of the Ising model in the plane of the complex (energy/temperature)-dependent variable u=e−4Ku=e^{-4K}u=e−4K, where K=J/(kBT)K=J/(k_BT)K=J/(kB​T), for nonzero external magnetic field, HHH. Exact results are given for the phase diagram in the uuu plane for the model in one dimension and on infinite-length quasi-one-dimensional strips. In the case of real h=H/(kBT)h=H/(k_BT)h=H/(kB​T), these results provide new insights into features of our earlier study of this case. We also consider complex h=H/(kBT)h=H/(k_BT)h=H/(kB​T) and μ=e−2hμ=e^{-2h}μ=e−2h. Calculations of complex-uuu zeros of the partition function on sections of the square lattice are presented. For the case of imaginary hhh, i.e., μ=eiθμ=e^{iθ}μ=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 uuu plane. We find that in this case, the phase boundary Bu{\cal B}_uBu​ contains a real line segment extending through part of the physical ferromagnetic interval 0≤u≤10 \le u \le 10≤u≤1, with a right-hand endpoint urheu_{rhe}urhe​ at the temperature for which the Yang-Lee edge singularity occurs at μ=e±iθμ=e^{\pm iθ}μ=e±iθ. Conformal field theory arguments are used to relate the singularities at urheu_{rhe}urhe​ and the Yang-Lee edge.
To study the relations between Yang-Lee's theory of the phase transitions and the distribution of the zeros of the partition functions in the antiferromagnetic case, the partition functions of finite (4×4 and 4×6) Ising model with the nearest and the next nearest neighbor interactions are calculated. The distributions of the zeros in the fugacity z plane of the partition function for various combinations of the values of J / k T - J '/ k T or tanh ( J /2 k T )-tanh ( J '/2 k T ) are investigated. The patterns of the distribution of the zeros show several types corresponding to ferromagnetic state, antiferromagnetic state and superantiferromagnetic state. In the antiferromagnetic ordered state the locus of the zeros is found to be nearly two concentric circles and to cross the positive real axis at two points z c and 1/ z c when the next neighbor ferromagnetic interaction is introduced. This means that our system shows the phase transitions at finite critical fields H c and - H c . The locus of zeros of Fisher's model of super-exchange antiferromagnet is also shown to cross the positive real axis at z c and 1/ z c below the Néel temperature.