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:2606.07415

Fermion sign problem and the structure of Lee-Yang zeros. II. Finite temperature results for a model system without interactions

Ran-Chen He, Jia-Xi Zeng, Shu Yang, Cong Wang, Qi-Jun Ye, Xin-Zheng Li

Published Jun 5, 2026• 0 citations
arXivHTMLPDF

Abstract

Beyond the analysis of the Lee-Yang (LY) zero of $ξ$ at $0$ K presented by our previous work [He et. al. Phys. Rev. E 113, 24115 (2026)], it is important but intricate to understand how these zeros evolve with temperature ($T$). Here, we use an analytically solvable noninteracting one-dimensional particle-on-a-ring model to address this. We determine the trajectories of these zeros and analyze how their evolution with $T$ reshapes the analytic structure of the partition function. In particular, the zero originating from $ξ=-1$ at $T=0$ remains close to $-1$ at low $T$, where it governs the sign factor and strongly constrains continuation along the real $ξ$ axis. This explains why both direct extrapolation and implicit schemes such as contour-based fitting can fail in the low-$T$ regime, even at high fitting order, while becoming reasonable again once the relevant zeros move away at higher $T$s. Furthermore, based on the polynomial structure of the partition function, we propose a new fitting strategy for low-$T$ fermionic properties. The key is to first obtain reliable high-$T$ fermionic properties by continuing sign-problem-free data in $ξ\in[0,1]$ to $ξ=-1$, and then extend this information toward lower $T$ through $T$-fitting of the $ξ$-independent remainder $ϕ(β)=Z_{\text{F}}$. These results provide a solvable benchmark for diagnosing the validity of analytic continuation and suggest a possible route toward treating more realistic interacting fermionic systems.

DOI
10.48550/arxiv.2606.07415

Related papers

Fermion sign problem and the structure of Lee-Yang zeros: The form of the partition function for indistinguishable particles and its zeros at 0 K

Ran-Chen He, Jia-Xi Zeng, Shu Yang +3 more

To simulate indistinguishable particles, recent studies of path-integral molecular dynamics formulated their partition function Z as a recurrence relation involving a variable ξ, with ξ=1(-1) for bosons (fermions). Inspired by Lee-Yang phase transition theory, we extend ξ into the complex plane and reformulate Z as a polynomial in ξ. By analyzing the distribution of the partition function zeros, we gain insights into the analytical properties of indistinguishable particles, particularly regarding the fermion sign problem (FSP). We found that at 0 K, the partition function zeros for N particles are located at ξ=-1, -1/2, -1/3, ..., -1/(N-1). This distribution disrupts the analytic continuation of thermodynamic quantities, expressed as functions of ξ and typically performed along ξ=1→-1, whenever the paths intersect these zeros. Moreover, we highlight the zero at ξ=-1, which induces an extra term in the free energy of the fermionic systems compared to ones at other ξ=e^{iθ} values. If a path connects this zero to a bosonic system with identical potential energies, it brings a transition resembling a phase transition. These findings provide a fresh perspective on the successes and challenges of emerging FSP studies based on analytic continuation techniques.

arXiv:2507.22779•Physical review. E•Dec 23, 2025•6 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

Lee-Yang theory of the two-dimensional quantum Ising model

Pascal M. Vecsei, José Luis Lado, Christian Flindt

Determining the phase diagram of interacting quantum many-body systems is an important task for a wide range of problems such as the understanding and design of quantum materials. For classical equilibrium systems, the Lee-Yang formalism provides a rigorous foundation of phase transitions, and these ideas have also been extended to the quantum realm. Here, we develop a Lee-Yang theory of quantum phase transitions that can include thermal fluctuations caused by a finite temperature, and it thereby provides a link between the classical Lee-Yang formalism and recent theories of phase transitions at zero temperature. Our methodology exploits analytic properties of the moment generating function of the order parameter in systems of finite size, and it can be implemented in combination with tensor-network calculations. Specifically, the onset of a symmetry-broken phase is signaled by the zeros of the moment generating function approaching the origin in the complex plane of a counting field that couples to the order parameter. Moreover, the zeros can be obtained by measuring or calculating the high cumulants of the order parameter. We determine the phase diagram of the two-dimensional quantum Ising model and thereby demonstrate the potential of our method to predict the critical behavior of two-dimensional quantum systems at finite temperatures.

arXiv:2204.08223•Physical review. B./Physical review. B•Aug 1, 2022•20 citations

Partition function zeros of quantum many-body systems: perturbative results

Muhammad Sedik, B. Sriram Shastry

A systematic method to find the Yang-Lee partition function zeros of quantum many-body systems based on perturbation theory at finite temperatures was recently introduced in arXiv:2504.01880. This method identifies wave-vector and temperature-dependent complex virtual energies obtainable from the thermal electronic Greens function. The collection of virtual energies over all $\vec{k}$ yield the Yang-Lee zeroes. We apply this approach to the one-dimensional Hubbard model for different boundary conditions. We compare the results obtained by this method up to second-order in perturbation theory with the results found by exact diagonalization. We also propose a quantity that could be used for experimental detection of these zeros of a Hubbard ring. An example of the detection method is presented using exact diagonalization of the 8-site Hubbard ring.

arXiv:2606.20808•Jun 18, 2026•0 citations

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

Lee-Yang theorem for fermions

Chaithanya Rayudu, Takahiro Misawa, Andrew Zhao +1 more

Lee-Yang theorems are a powerful tool for studying many-body systems, with applications ranging from analyzing phase transitions to proving the efficiency of certain classical and quantum algorithms. In this work, we prove a Lee-Yang zero-freeness theorem for the partition function of a broad class of interacting fermion models, implying the existence of a provably efficient quantum algorithm for estimating their ground-state energies. This class includes several well-known models such as the attractive Hubbard model, repulsive Hubbard model on bipartite graphs, and the interacting Hofstadter model. Our results also rigorously establish the nonexistence of phase transitions in these models in the presence of a nonzero local external field.

quant-pharXiv:2609.23942•Sep 22, 2026•0 citations
Beyond the analysis of the Lee-Yang (LY) zero of ξξξ at 000 K presented by our previous work [He et. al. Phys. Rev. E 113, 24115 (2026)], it is important but intricate to understand how these zeros evolve with temperature (TTT). Here, we use an analytically solvable noninteracting one-dimensional particle-on-a-ring model to address this. We determine the trajectories of these zeros and analyze how their evolution with TTT reshapes the analytic structure of the partition function. In particular, the zero originating from ξ=−1ξ=-1ξ=−1 at T=0T=0T=0 remains close to −1-1−1 at low TTT, where it governs the sign factor and strongly constrains continuation along the real ξξξ axis. This explains why both direct extrapolation and implicit schemes such as contour-based fitting can fail in the low-TTT regime, even at high fitting order, while becoming reasonable again once the relevant zeros move away at higher TTTs. Furthermore, based on the polynomial structure of the partition function, we propose a new fitting strategy for low-TTT fermionic properties. The key is to first obtain reliable high-TTT fermionic properties by continuing sign-problem-free data in ξ∈[0,1]ξ\in[0,1]ξ∈[0,1] to ξ=−1ξ=-1ξ=−1, and then extend this information toward lower TTT through TTT-fitting of the ξξξ-independent remainder ϕ(β)=ZFϕ(β)=Z_{\text{F}}ϕ(β)=ZF​. These results provide a solvable benchmark for diagnosing the validity of analytic continuation and suggest a possible route toward treating more realistic interacting fermionic systems.
To simulate indistinguishable particles, recent studies of path-integral molecular dynamics formulated their partition function Z as a recurrence relation involving a variable ξ, with ξ=1(-1) for bosons (fermions). Inspired by Lee-Yang phase transition theory, we extend ξ into the complex plane and reformulate Z as a polynomial in ξ. By analyzing the distribution of the partition function zeros, we gain insights into the analytical properties of indistinguishable particles, particularly regarding the fermion sign problem (FSP). We found that at 0 K, the partition function zeros for N particles are located at ξ=-1, -1/2, -1/3, ..., -1/(N-1). This distribution disrupts the analytic continuation of thermodynamic quantities, expressed as functions of ξ and typically performed along ξ=1→-1, whenever the paths intersect these zeros. Moreover, we highlight the zero at ξ=-1, which induces an extra term in the free energy of the fermionic systems compared to ones at other ξ=e^{iθ} values. If a path connects this zero to a bosonic system with identical potential energies, it brings a transition resembling a phase transition. These findings provide a fresh perspective on the successes and challenges of emerging FSP studies based on analytic continuation techniques.
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 .
Determining the phase diagram of interacting quantum many-body systems is an important task for a wide range of problems such as the understanding and design of quantum materials. For classical equilibrium systems, the Lee-Yang formalism provides a rigorous foundation of phase transitions, and these ideas have also been extended to the quantum realm. Here, we develop a Lee-Yang theory of quantum phase transitions that can include thermal fluctuations caused by a finite temperature, and it thereby provides a link between the classical Lee-Yang formalism and recent theories of phase transitions at zero temperature. Our methodology exploits analytic properties of the moment generating function of the order parameter in systems of finite size, and it can be implemented in combination with tensor-network calculations. Specifically, the onset of a symmetry-broken phase is signaled by the zeros of the moment generating function approaching the origin in the complex plane of a counting field that couples to the order parameter. Moreover, the zeros can be obtained by measuring or calculating the high cumulants of the order parameter. We determine the phase diagram of the two-dimensional quantum Ising model and thereby demonstrate the potential of our method to predict the critical behavior of two-dimensional quantum systems at finite temperatures.
A systematic method to find the Yang-Lee partition function zeros of quantum many-body systems based on perturbation theory at finite temperatures was recently introduced in arXiv:2504.01880. This method identifies wave-vector and temperature-dependent complex virtual energies obtainable from the thermal electronic Greens function. The collection of virtual energies over all k⃗\vec{k}k yield the Yang-Lee zeroes. We apply this approach to the one-dimensional Hubbard model for different boundary conditions. We compare the results obtained by this method up to second-order in perturbation theory with the results found by exact diagonalization. We also propose a quantity that could be used for experimental detection of these zeros of a Hubbard ring. An example of the detection method is presented using exact diagonalization of the 8-site Hubbard ring.
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.
Lee-Yang theorems are a powerful tool for studying many-body systems, with applications ranging from analyzing phase transitions to proving the efficiency of certain classical and quantum algorithms. In this work, we prove a Lee-Yang zero-freeness theorem for the partition function of a broad class of interacting fermion models, implying the existence of a provably efficient quantum algorithm for estimating their ground-state energies. This class includes several well-known models such as the attractive Hubbard model, repulsive Hubbard model on bipartite graphs, and the interacting Hofstadter model. Our results also rigorously establish the nonexistence of phase transitions in these models in the presence of a nonzero local external field.