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:2412.07328 · Physical review. B./Physical review. B

Low-temperature series expansion of the square-lattice Ising model: A study based on Fisher zeros

De-Zhang Li, Xin Wang, Xiao‐Bao Yang

Published Apr 21, 2025• 3 citations
arXivHTMLPDF

Abstract

Low-temperature expansion of Ising model has long been a topic of significant interest in condensed matter and statistical physics. In this paper we present new results of the coefficients in the low-temperature series of the Ising partition function on the square lattice, in the cases of a zero field and of an imaginary field $i(\ensuremath{\pi}/2){k}_{B}T$. The coefficients in the low-temperature series of the free energy in the thermodynamic limit are represented using the explicit expression of the density function of the Fisher zeros. The asymptotic behavior of the sequence of the coefficients when the order goes to infinity is determined exactly, for both the series of the free energy and of the partition function. Our analytic and numerical results demonstrate that the convergence radius of the sequence is dependent on the accumulation points of the Fisher zeros, which have the smallest modulus. In the zero field case this accumulation point is the physical critical point, while in the imaginary field case it corresponds to a nonphysical singularity. We further discuss the relation between the series coefficients and the energy state degeneracies, using the combinatorial expression of the coefficients and the subgraph expansion.

DOI
10.1103/physrevb.111.155432

Related papers

Exact low-temperature series expansion for the partition function of the two-dimensional zero-field s=1/2 Ising model on the infinite square lattice

Grzegorz Siudem, Agata Fronczak, Piotr Fronczak

In this paper, we provide the exact expression for the coefficients in the low-temperature series expansion of the partition function of the two-dimensional Ising model on the infinite square lattice. This is equivalent to exact determination of the number of spin configurations at a given energy. With these coefficients, we show that the ferromagnetic--to--paramagnetic phase transition in the square lattice Ising model can be explained through equivalence between the model and the perfect gas of energy clusters model, in which the passage through the critical point is related to the complete change in the thermodynamic preferences on the size of clusters. The combinatorial approach reported in this article is very general and can be easily applied to other lattice models.

arXiv:1410.7963•Oct 29, 2014•0 citations

Free-fermion approach to the partition function zeros: Special boundary conditions and product form of solution

De-Zhang Li, Xin Wang

Partition function zeros are powerful tools in understanding critical behavior. In this paper, we present results of the Fisher zeros of two-dimensional Ising models, in the framework of free-fermion eight-vertex model. First, we succeed in finding special boundary conditions for the free-fermion model, under which the partition function of a finite lattice can be expressed in a double product form. Using appropriate mappings, these boundary conditions are transformed into the corresponding versions of the square, triangular, and honeycomb lattice Ising models. Each Ising model is studied in the cases of a zero field and of an imaginary field i ( π / 2 ) k B T . For the square lattice model, we rediscover the famous Brascamp-Kunz (B-K) boundary conditions. For the triangular and honeycomb lattice models, we obtain the B-K-type boundary conditions, and the Fisher zeros are conveniently solved from the product form of partition function. The advantage of B-K-type boundary conditions is that the Fisher zeros of any finite lattice exactly lie on certain loci, and the accumulation points of zeros can be easily determined in the thermodynamic limit. Our finding and method would be very helpful in studying the partition function zeros of vertex and Ising models.

arXiv:2507.21943•Physical Review Research•Nov 13, 2025•1 citation

Fisher zeros of the Ising antiferromagnet in an arbitrary nonzero magnetic field

Seung‐Yeon Kim

From the exact partition functions of the Ising model on L x L square lattices (up to L = 14) in an arbitrary nonzero external magnetic field H, the precise distributions of the Fisher zeros in the complex temperature (a = e(2betaJ)) plane are obtained. Our results indicate that in the limit L-->infinity the Fisher zeros for H not = 0 do cut the positive real axis at the antiferromagnetic critical point ac(H); these values are compared with the results of closed-form approximations for the antiferromagnetic critical line. From the Fisher zeros the thermal scaling exponent yt = 1 is also obtained along the critical line, indicating the logarithmic singularity of the specific heat for the antiferromagnetic Ising model even in a strong uniform magnetic field.

Physical Review E•Jan 10, 2005•15 citations

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

Mahmoud Abdelshafy, Muhammad Sedik

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.

arXiv:2503.14602•Physical review. B./Physical review. B•Aug 8, 2025•2 citations

The susceptibility of the square lattice Ising model: New developments

William P. Orrick, Bernie G. Nickel, Anthony John Guttmann +1 more

We have made substantial advances in elucidating the properties of the susceptibility of the square lattice Ising model. We discuss its analyticity properties, certain closed form expressions for subsets of the coefficients, and give an algorithm of complexity O(N6) to determine its first N coefficients. As a result, we have generated and analyzed series with more than 300 terms in both the high- and low-temperature regime. We quantify the effect of irrelevant variables to the scaling-amplitude functions. In particular, we find and quantify the breakdown of simple scaling, in the absence of irrelevant 1 scaling fields, arising first at order |T −Tc | 9/4, though high-low temperature symmetry is still preserved. At terms of order |T −Tc | 17/4 and beyond, this symmetry is no longer present. The short-distance terms are shown to have the form (T − Tc) p (log |T − Tc|) q with p ≥ q2. Conjectured exact expressions for some correlation functions and series coefficients in terms of elliptic theta functions also foreshadow future developments.

•Nov 15, 2012•69 citations

Low-temperature behavior of the finite-size one-dimensional Ising model and the partition function zeros

Julian Lee

In contrast to an infinite chain, the low-temperature expansion of a one-dimensional free-field Ising model has a strong dependence on boundary conditions. I derive an explicit formula for the leading term of the expansion both under open and periodic boundary conditions and show that they are related to different distributions of the partition function zeros on the complex temperature plane. In particular, when a periodic boundary condition is imposed, the leading coefficient of the expansion grows with increasing size of the chain, due to the zeros approaching the origin.

arXiv:1407.6919•Journal of the Korean Physical Society•Sep 1, 2014•10 citations
Low-temperature expansion of Ising model has long been a topic of significant interest in condensed matter and statistical physics. In this paper we present new results of the coefficients in the low-temperature series of the Ising partition function on the square lattice, in the cases of a zero field and of an imaginary field i(\ensuremathπ/2)kBTi(\ensuremath{\pi}/2){k}_{B}Ti(\ensuremathπ/2)kB​T. The coefficients in the low-temperature series of the free energy in the thermodynamic limit are represented using the explicit expression of the density function of the Fisher zeros. The asymptotic behavior of the sequence of the coefficients when the order goes to infinity is determined exactly, for both the series of the free energy and of the partition function. Our analytic and numerical results demonstrate that the convergence radius of the sequence is dependent on the accumulation points of the Fisher zeros, which have the smallest modulus. In the zero field case this accumulation point is the physical critical point, while in the imaginary field case it corresponds to a nonphysical singularity. We further discuss the relation between the series coefficients and the energy state degeneracies, using the combinatorial expression of the coefficients and the subgraph expansion.
In this paper, we provide the exact expression for the coefficients in the low-temperature series expansion of the partition function of the two-dimensional Ising model on the infinite square lattice. This is equivalent to exact determination of the number of spin configurations at a given energy. With these coefficients, we show that the ferromagnetic--to--paramagnetic phase transition in the square lattice Ising model can be explained through equivalence between the model and the perfect gas of energy clusters model, in which the passage through the critical point is related to the complete change in the thermodynamic preferences on the size of clusters. The combinatorial approach reported in this article is very general and can be easily applied to other lattice models.
Partition function zeros are powerful tools in understanding critical behavior. In this paper, we present results of the Fisher zeros of two-dimensional Ising models, in the framework of free-fermion eight-vertex model. First, we succeed in finding special boundary conditions for the free-fermion model, under which the partition function of a finite lattice can be expressed in a double product form. Using appropriate mappings, these boundary conditions are transformed into the corresponding versions of the square, triangular, and honeycomb lattice Ising models. Each Ising model is studied in the cases of a zero field and of an imaginary field i ( π / 2 ) k B T . For the square lattice model, we rediscover the famous Brascamp-Kunz (B-K) boundary conditions. For the triangular and honeycomb lattice models, we obtain the B-K-type boundary conditions, and the Fisher zeros are conveniently solved from the product form of partition function. The advantage of B-K-type boundary conditions is that the Fisher zeros of any finite lattice exactly lie on certain loci, and the accumulation points of zeros can be easily determined in the thermodynamic limit. Our finding and method would be very helpful in studying the partition function zeros of vertex and Ising models.
From the exact partition functions of the Ising model on L x L square lattices (up to L = 14) in an arbitrary nonzero external magnetic field H, the precise distributions of the Fisher zeros in the complex temperature (a = e(2betaJ)) plane are obtained. Our results indicate that in the limit L-->infinity the Fisher zeros for H not = 0 do cut the positive real axis at the antiferromagnetic critical point ac(H); these values are compared with the results of closed-form approximations for the antiferromagnetic critical line. From the Fisher zeros the thermal scaling exponent yt = 1 is also obtained along the critical line, indicating the logarithmic singularity of the specific heat for the antiferromagnetic Ising model even in a strong uniform magnetic field.
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.
We have made substantial advances in elucidating the properties of the susceptibility of the square lattice Ising model. We discuss its analyticity properties, certain closed form expressions for subsets of the coefficients, and give an algorithm of complexity O(N6) to determine its first N coefficients. As a result, we have generated and analyzed series with more than 300 terms in both the high- and low-temperature regime. We quantify the effect of irrelevant variables to the scaling-amplitude functions. In particular, we find and quantify the breakdown of simple scaling, in the absence of irrelevant 1 scaling fields, arising first at order |T −Tc | 9/4, though high-low temperature symmetry is still preserved. At terms of order |T −Tc | 17/4 and beyond, this symmetry is no longer present. The short-distance terms are shown to have the form (T − Tc) p (log |T − Tc|) q with p ≥ q2. Conjectured exact expressions for some correlation functions and series coefficients in terms of elliptic theta functions also foreshadow future developments.
In contrast to an infinite chain, the low-temperature expansion of a one-dimensional free-field Ising model has a strong dependence on boundary conditions. I derive an explicit formula for the leading term of the expansion both under open and periodic boundary conditions and show that they are related to different distributions of the partition function zeros on the complex temperature plane. In particular, when a periodic boundary condition is imposed, the leading coefficient of the expansion grows with increasing size of the chain, due to the zeros approaching the origin.