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

arXiv:2503.14602 · Physical review. B./Physical review. B

Yang-Lee zeros of two-dimensional nearest-neighbor antiferromagnetic Ising models: A numerical linked cluster expansion study

Mahmoud Abdelshafy, Muhammad Sedik

Published Aug 8, 2025• 2 citations
arXivHTMLPDF

Abstract

We study Yang-Lee zeros as they approach the thermodynamic limit of the two-dimensional (2D) nearest-neighbor antiferromagnetic Ising model on square and triangular lattices. We employ the numerical linked cluster expansion (NLCE) equipped with exact enumeration (EE) of the partition function to estimate the Laplacian of the free energy, which is proportional to the zeros density. Using a modified NLCE, where the expansion can be carried directly on the Yang-Lee zeros of the involved clusters, we estimate the density of Yang-Lee zeros in the thermodynamic limit. NLCE gives significantly more zeros than EE in the complex field plane providing more insights on how the root curves look in the thermodynamic limit. For the square lattice at $T\ensuremath{\ll}{T}_{c}$, the results suggest that two vertical lines at $\ifmmode\pm\else\textpm\fi{}{h}_{c}(T)$ in the complex field plane (i.e., two concentric circles in the complex fugacity plane) are the thermodynamic root curves. A similar picture is found for the triangular lattice for phase transitions at large values of magnetic field while further study is needed for phase transitions at smaller values of magnetic field. The convergence of the NLCE and (EE) calculations of the partition function to the thermodynamic limit is studied in both lattices and the temperature-field phase diagram is obtained from Yang-Lee zeros using both methods. This NLCE-based approach will facilitate the study of different types of phase transitions using Yang-Lee zeros in future research.

DOI
10.1103/x38c-w4z2

Related papers

Yang-Lee zeros of certain antiferromagnetic models

Muhammad Sedik, Junaid Majeed Bhat, Abhishek Dhar +1 more

We revisit the somewhat less studied problem of Yang-Lee zeros of the Ising antiferromagnet. For this purpose, we study two models, the nearest-neighbor model on a square lattice and the more tractable mean-field model corresponding to infinite-ranged coupling between all sites. In the high-temperature limit, we show that the logarithm of the Yang-Lee zeros can be written as a series in half odd integer powers of the inverse temperature, k, with the leading term ∼k^{1/2}. This result is true in any dimension and for arbitrary lattices. We also show that the coefficients of the expansion satisfy simple identities (akin to sum rules) for the nearest-neighbor case. These identities are verified numerically by computing the exact partition function for a two-dimensional square lattice of size 16×16. For the mean-field model, we write down the partition function (termed the mean-field polynomials) for the ferromagnetic (FM) and antiferromagnetic (AFM) cases and derive from them the mean-field equations. We analytically show that at high temperatures the zeros of the AFM mean-field polynomial scale as ∼k^{1/2} as well. Using a simple numerical method, we find the roots lie on certain curves (the root curves), in the thermodynamic limit for the mean-field polynomials for the AFM case as well as for the FM one. Our results show a new root curve that was not found earlier. Our results also clearly illustrate the phase transition expected for the FM and AFM cases, in the language of Yang-Lee zeros. Moreover, for the AFM case, we observe that the root curves separate two distinct phases of zero and nonzero complex staggered magnetization, and thus depict a complex phase boundary.

arXiv:2309.14562•Physical review. E•Jul 11, 2024•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

Complex-temperature properties of the two-dimensional Ising model for nonzero magnetic field

V. Matveev, Robert E. Shrock

We study the complex-temperature phase diagram of the square-lattice Ising model for nonzero external magnetic field H, i.e., for 0\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}\ensuremath{\infty}, where \ensuremath{\mu}=${\mathit{e}}^{\mathrm{\ensuremath{-}}2\mathrm{\ensuremath{\beta}}\mathit{H}}$. We also carry out a similar analysis for -\ensuremath{\infty}\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}0. The results for the interval -1\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}1 provide a way of continuously connecting the two known exact solutions of this model, viz., at \ensuremath{\mu}=1 (Onsager, Yang) and \ensuremath{\mu}=-1 (Lee and Yang). Our methods include numerical calculations of complex-temperature zeros of the partition function and an analysis of low-temperature series expansions. For real nonzero H, the inner branch of a lima\ifmmode \mbox{\c{c}}\else \c{c}\fi{}on bounding the FM phase breaks and forms two complex-conjugate arcs. We study the singularities and associated exponents of thermodynamic functions at the endpoints of these arcs. For \ensuremath{\mu}0, there are two line segments of singularities on the negative and positive u axis, and we carry out a similar study of the behavior at the inner endpoints of these arcs, which constitute the nearest singularities to the origin in this case. Finally, we also determine the exact complex-temperature phase diagrams at \ensuremath{\mu}=-1 on the honeycomb and triangular lattices and discuss the relation between these and the corresponding zero-field phase diagrams. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9507120•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1996•53 citations

Density of Yang-Lee zeros in the thermodynamic limit from tensor network methods

Artur García-Sáez, Tzu-Chieh Wei

The distribution of Yang-Lee zeros in the ferromagnetic Ising model in both two and three dimensions is studied on the complex field plane directly in the thermodynamic limit via the tensor network methods. The partition function is represented as a contraction of a tensor network and is efficiently evaluated with an iterative tensor renormalization scheme. The free-energy density and the magnetization are computed on the complex field plane. Via the discontinuity of the magnetization, the density of the Yang-Lee zeros is obtained to lie on the unit circle, consistent with the Lee-Yang circle theorem. Distinct features are observed at different temperatures---below, above, and at the critical temperature. Application of the tensor network approach is also made to the $q$-state Potts models in both two and three dimensions and a previous debate on whether, in the thermodynamic limit, the Yang-Lee zeros lie on a unit circle for $q>2$ is resolved: they clearly do not lie on a unit circle except at the zero temperature. For the Potts models $(q=3,4,5,6)$ investigated in two dimensions, as the temperature is lowered the radius of the zeros at a fixed angle from the real axis shrinks exponentially towards unity with the inverse temperature.

arXiv:1508.07751•Physical Review B•Sep 17, 2015•15 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

Yang–Lee zeros of the Yang–Lee model

Giuseppe Mussardo, Riccarda Bonsignori, Andrea Trombettoni

Abstract To understand the distribution of the Yang–Lee zeros in quantum integrable field theories we analyse the simplest of these systems given by the 2D Yang–Lee model. The grand-canonical partition function of this quantum field theory, as a function of the fugacity z and the inverse temperature β , can be computed in terms of the thermodynamics Bethe Ansatz based on its exact S -matrix. We extract the Yang–Lee zeros in the complex plane by using a sequence of polynomials of increasing order N in z which converges to the grand-canonical partition function. We show that these zeros are distributed along curves which are approximate circles as it is also the case of the zeros for purely free theories. There is though an important difference between the interactive theory and the free theories, for the radius of the zeros in the interactive theory goes continuously to zero in the high-temperature limit β → 0 while in the free theories it remains close to 1 even for small values of β , jumping to 0 only at β = 0 .

arXiv:1708.06444•Journal of Physics A Mathematical and Theoretical•Sep 29, 2017•7 citations
We study Yang-Lee zeros as they approach the thermodynamic limit of the two-dimensional (2D) nearest-neighbor antiferromagnetic Ising model on square and triangular lattices. We employ the numerical linked cluster expansion (NLCE) equipped with exact enumeration (EE) of the partition function to estimate the Laplacian of the free energy, which is proportional to the zeros density. Using a modified NLCE, where the expansion can be carried directly on the Yang-Lee zeros of the involved clusters, we estimate the density of Yang-Lee zeros in the thermodynamic limit. NLCE gives significantly more zeros than EE in the complex field plane providing more insights on how the root curves look in the thermodynamic limit. For the square lattice at T\ensuremath≪TcT\ensuremath{\ll}{T}_{c}T\ensuremath≪Tc​, the results suggest that two vertical lines at \ifmmode±\else\textpm\fihc(T)\ifmmode\pm\else\textpm\fi{}{h}_{c}(T)\ifmmode±\else\textpm\fihc​(T) in the complex field plane (i.e., two concentric circles in the complex fugacity plane) are the thermodynamic root curves. A similar picture is found for the triangular lattice for phase transitions at large values of magnetic field while further study is needed for phase transitions at smaller values of magnetic field. The convergence of the NLCE and (EE) calculations of the partition function to the thermodynamic limit is studied in both lattices and the temperature-field phase diagram is obtained from Yang-Lee zeros using both methods. This NLCE-based approach will facilitate the study of different types of phase transitions using Yang-Lee zeros in future research.
We revisit the somewhat less studied problem of Yang-Lee zeros of the Ising antiferromagnet. For this purpose, we study two models, the nearest-neighbor model on a square lattice and the more tractable mean-field model corresponding to infinite-ranged coupling between all sites. In the high-temperature limit, we show that the logarithm of the Yang-Lee zeros can be written as a series in half odd integer powers of the inverse temperature, k, with the leading term ∼k^{1/2}. This result is true in any dimension and for arbitrary lattices. We also show that the coefficients of the expansion satisfy simple identities (akin to sum rules) for the nearest-neighbor case. These identities are verified numerically by computing the exact partition function for a two-dimensional square lattice of size 16×16. For the mean-field model, we write down the partition function (termed the mean-field polynomials) for the ferromagnetic (FM) and antiferromagnetic (AFM) cases and derive from them the mean-field equations. We analytically show that at high temperatures the zeros of the AFM mean-field polynomial scale as ∼k^{1/2} as well. Using a simple numerical method, we find the roots lie on certain curves (the root curves), in the thermodynamic limit for the mean-field polynomials for the AFM case as well as for the FM one. Our results show a new root curve that was not found earlier. Our results also clearly illustrate the phase transition expected for the FM and AFM cases, in the language of Yang-Lee zeros. Moreover, for the AFM case, we observe that the root curves separate two distinct phases of zero and nonzero complex staggered magnetization, and thus depict a complex phase boundary.
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 the complex-temperature phase diagram of the square-lattice Ising model for nonzero external magnetic field H, i.e., for 0\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}\ensuremath{\infty}, where \ensuremath{\mu}=${\mathit{e}}^{\mathrm{\ensuremath{-}}2\mathrm{\ensuremath{\beta}}\mathit{H}}$. We also carry out a similar analysis for -\ensuremath{\infty}\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}0. The results for the interval -1\ensuremath{\le}\ensuremath{\mu}\ensuremath{\le}1 provide a way of continuously connecting the two known exact solutions of this model, viz., at \ensuremath{\mu}=1 (Onsager, Yang) and \ensuremath{\mu}=-1 (Lee and Yang). Our methods include numerical calculations of complex-temperature zeros of the partition function and an analysis of low-temperature series expansions. For real nonzero H, the inner branch of a lima\ifmmode \mbox{\c{c}}\else \c{c}\fi{}on bounding the FM phase breaks and forms two complex-conjugate arcs. We study the singularities and associated exponents of thermodynamic functions at the endpoints of these arcs. For \ensuremath{\mu}0, there are two line segments of singularities on the negative and positive u axis, and we carry out a similar study of the behavior at the inner endpoints of these arcs, which constitute the nearest singularities to the origin in this case. Finally, we also determine the exact complex-temperature phase diagrams at \ensuremath{\mu}=-1 on the honeycomb and triangular lattices and discuss the relation between these and the corresponding zero-field phase diagrams. \textcopyright{} 1996 The American Physical Society.
The distribution of Yang-Lee zeros in the ferromagnetic Ising model in both two and three dimensions is studied on the complex field plane directly in the thermodynamic limit via the tensor network methods. The partition function is represented as a contraction of a tensor network and is efficiently evaluated with an iterative tensor renormalization scheme. The free-energy density and the magnetization are computed on the complex field plane. Via the discontinuity of the magnetization, the density of the Yang-Lee zeros is obtained to lie on the unit circle, consistent with the Lee-Yang circle theorem. Distinct features are observed at different temperatures---below, above, and at the critical temperature. Application of the tensor network approach is also made to the qqq-state Potts models in both two and three dimensions and a previous debate on whether, in the thermodynamic limit, the Yang-Lee zeros lie on a unit circle for q>2q>2q>2 is resolved: they clearly do not lie on a unit circle except at the zero temperature. For the Potts models (q=3,4,5,6)(q=3,4,5,6)(q=3,4,5,6) investigated in two dimensions, as the temperature is lowered the radius of the zeros at a fixed angle from the real axis shrinks exponentially towards unity with the inverse temperature.
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.
Abstract To understand the distribution of the Yang–Lee zeros in quantum integrable field theories we analyse the simplest of these systems given by the 2D Yang–Lee model. The grand-canonical partition function of this quantum field theory, as a function of the fugacity z and the inverse temperature β , can be computed in terms of the thermodynamics Bethe Ansatz based on its exact S -matrix. We extract the Yang–Lee zeros in the complex plane by using a sequence of polynomials of increasing order N in z which converges to the grand-canonical partition function. We show that these zeros are distributed along curves which are approximate circles as it is also the case of the zeros for purely free theories. There is though an important difference between the interactive theory and the free theories, for the radius of the zeros in the interactive theory goes continuously to zero in the high-temperature limit β → 0 while in the free theories it remains close to 1 even for small values of β , jumping to 0 only at β = 0 .