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

Spectral gap of Lee-Yang Hamiltonians

Chaithanya Rayudu, Jun Takahashi

quant-phcond-mat.stat-mechcs.CCmath-phSubmitted Jul 12, 2026• Announced Jul 14, 2026• 0 citations
arXivHTMLPDF

Abstract

The Lee-Yang theorem and its quantum extensions state that, for a broad class of Hamiltonians on any graph, the partition function's zeros in the complex magnetic field plane lie only on the imaginary axis. For these Hamiltonians, we prove that under a uniform Z-field of any strength h, the ground state has a spectral gap of at least h/4, independent of the system size and of the coupling strengths. The proof uses the zero-freeness of the partition function as given by Asano and Suzuki-Fisher to show exponential decay of the imaginary-time correlations for any product of Z-operators. Our result gives a polynomial-time quantum algorithm for computing the ground state energy of any Lee-Yang Hamiltonian.

AI Summary

The paper asks whether the Lee–Yang zero-freeness property can imply a quantitative ground-state spectral gap and efficient quantum ground-energy estimation. It considers two-local qubit Hamiltonians on arbitrary graphs with ferromagnetic $ZZ$ couplings satisfying a bondwise dominance condition over $XX$, $YY$, $XY$, and $YX$ interactions. Without cross couplings, this condition reduces to $w^z_{ij}\geq\max(|w^x_{ij}|,|w^y_{ij}|)$. For every Hamiltonian in this class and every uniform field $h>0$, it proves that $H_h=H-h\sum_i Z_i$ has a unique ground state and gap $E_1-E_0\geq h/4$, independently of system size and coupling strengths. The proof retains separate field variables at every imaginary-time Trotter layer. Lee–Yang zero-freeness makes the log-partition function analytic, while a Taylor-support lemma forces contributions connecting distant layers to have high degree. Analytic bounds then yield decay of Trotterized correlations. An adaptive Trotter step controls the continuum limit, establishing connected imaginary-time correlation decay for every $Z_S=\prod_{i\in S}Z_i$ at rate $h/2$, with polynomial prefactors. A spectral representation and the completeness of diagonal Pauli products imply the gap and uniqueness. The gap enables polynomial-time adiabatic ground-energy estimation to inverse-polynomial additive accuracy, placing EPR-Hamiltonian energy estimation and bipartite quantum Max-Cut in $\mathsf{BQP}$. The theorem requires a nonzero field and gives no gap bound at $h=0$; energy estimation there uses a small-field approximation. No classical hardness or quantum advantage is established, and containment of the sign-problem-free EPR case in $\mathsf{BPP}$ remains possible.

No sign of AI-generated text

Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.

Comments
23 pages

Related papers

Efficient quantum algorithm for Heisenberg spin systems

Sergey Bravyi, David Gosset, Yinchen Liu +1 more

We consider a broad family of Heisenberg-type quantum spin systems that were studied by Suzuki and Fisher in 1971. This family includes, for example, the Heisenberg antiferromagnet on any bipartite graph. The ground and Gibbs states of these models are Lee-Yang tensors with radius 1: they are associated with multilinear polynomials that possess an extraordinary zero-freeness property inside the unit polydisk in the complex plane. For each Hamiltonian in this family we show that the spectral gap between the first-excited and ground-state energies is lower bounded by $2μ$, where $μ$ is the magnetic field along the $Z$ direction. Using this result we obtain an efficient quantum adiabatic algorithm for the ground energy of any model in this family. The proof is based on a new inequality that relates the spectral gap of a positive semidefinite operator to its Lee-Yang radius -- a quantitative strengthening of prior work of the authors that may find applications elsewhere.

quant-pharXiv:2607.14401•Jul 17, 2026•0 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

The Geometry of Lee-Yang Tensors

Jake Hofgard

We provide a geometric description of Lee-Yang tensors, building upon several recent results of Bravyi, Gosset, Liu, and Wong. By applying tools from complex Perron-Frobenius theory, we generalize one of their key lemmas to bound the spectral ratio of non-normal operators with Lee-Yang radius greater than one. As a direct consequence, we provide a local-to-global principle for bounding spectral gaps of Suzuki-Fisher Hamiltonians, a family of $2$-local Hamiltonians, using only the Lee-Yang radius of each local operator to bound the gap from below. On the other hand, we develop a simple perturbative technique to show that the Gibbs states of $k$-local Hamiltonians, $k \geq 3$, cannot have Lee-Yang radius at least one for all inverse temperatures $β\geq 0$. We apply this same method to conclude that Suzuki-Fisher Hamiltonians are precisely those Hamiltonians whose Gibbs states are Lee-Yang at all inverse temperatures $β\geq 0$.

quant-pharXiv:2610.09485•Oct 8, 2026•0 citations

Lee-Yang tensors and Hamiltonian complexity

Benjamin Wong, Sergey Bravyi, David Gosset +1 more

A complex tensor with $n$ binary indices can be identified with a multilinear polynomial in $n$ complex variables. We say it is a Lee-Yang tensor with radius $r$ if the polynomial is nonzero whenever all variables lie in the open disk of radius $r$. In this work we study quantum states and observables which are Lee-Yang tensors when expressed in the computational basis. We first review their basic properties, including closure under tensor contraction and certain quantum operations. We show that quantum states with Lee-Yang radius $r > 1$ can be prepared by quasipolynomial-sized circuits. We also show that every Hermitian operator with Lee-Yang radius $r > 1$ has a unique principal eigenvector. These results suggest that $r = 1$ is a key threshold for quantum states and observables. Finally, we consider a family of two-local Hamiltonians where every interaction term energetically favors a deformed EPR state $|00\rangle + s|11\rangle$ for some $0 \leq s \leq 1$. We numerically investigate this model and find that on all graphs considered the Lee-Yang radius of the ground state is at least $r = 1/\sqrt{s}$ while the spectral gap between the two smallest eigenvalues is at least $1-s^2$. We conjecture that these lower bounds hold more generally; in particular, this would provide an efficient quantum adiabatic algorithm for the quantum Max-Cut problem on uniformly weighted bipartite graphs.

quant-pharXiv:2602.03605•Feb 4, 2026•0 citations

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

Masuo Suzuki, Michael Ellis Fisher

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.

Journal of Mathematical Physics•Feb 1, 1971•164 citations

Yang-Lee Edge Singularity and ϕ 3 Field Theory

Michael Ellis Fisher

The edge of the gap in the distribution of Yang-Lee zeros at $H=i{H}_{0}(T)$ on the imaginary magnetic field axis in ferromagnets above ${T}_{c}$ is essentially a critical point. In terms of the edge exponents $\ensuremath{\delta}$ and $\ensuremath{\eta}$, the density of zeros obeys $\mathcal{G}({H}^{\ensuremath{'}\ensuremath{'}})\ensuremath{\sim}{[{H}^{\ensuremath{'}\ensuremath{'}}\ensuremath{-}{H}_{0}(T)]}^{\ensuremath{\sigma}}$, with $\ensuremath{\sigma}=\frac{1}{\ensuremath{\delta}}=\frac{(d\ensuremath{-}2+\ensuremath{\eta})}{({d}_{2}\ensuremath{-}\ensuremath{\eta})}$. Classical behavior ($\ensuremath{\sigma}=\frac{1}{2}$) occurs for $d>{d}^{\ifmmode\times\else\texttimes\fi{}}=6$. The appropriate field-theoretic renormalization group entails a $w{\ensuremath{\phi}}^{3}$ coupling and, with $\ensuremath{\epsilon}=6\ensuremath{-}d>~0$, yields $\ensuremath{\eta}\ensuremath{\approx}\ensuremath{-}\frac{\ensuremath{\epsilon}}{9}$ for all $n<\ensuremath{\infty}$. This correlates well with refined series estimates for $d=2$ and $d=3$ and with exact results for $d=1$ ($\ensuremath{\eta}=\ensuremath{-}1$).

Physical Review Letters•Jun 19, 1978•485 citations
The Lee-Yang theorem and its quantum extensions state that, for a broad class of Hamiltonians on any graph, the partition function's zeros in the complex magnetic field plane lie only on the imaginary axis. For these Hamiltonians, we prove that under a uniform Z-field of any strength h, the ground state has a spectral gap of at least h/4, independent of the system size and of the coupling strengths. The proof uses the zero-freeness of the partition function as given by Asano and Suzuki-Fisher to show exponential decay of the imaginary-time correlations for any product of Z-operators. Our result gives a polynomial-time quantum algorithm for computing the ground state energy of any Lee-Yang Hamiltonian.
The paper asks whether the Lee–Yang zero-freeness property can imply a quantitative ground-state spectral gap and efficient quantum ground-energy estimation. It considers two-local qubit Hamiltonians on arbitrary graphs with ferromagnetic ZZZZZZ couplings satisfying a bondwise dominance condition over XXXXXX, YYYYYY, XYXYXY, and YXYXYX interactions. Without cross couplings, this condition reduces to wijz≥max⁡(∣wijx∣,∣wijy∣)w^z_{ij}\geq\max(|w^x_{ij}|,|w^y_{ij}|)wijz​≥max(∣wijx​∣,∣wijy​∣).
For every Hamiltonian in this class and every uniform field h>0h>0h>0, it proves that Hh=H−h∑iZiH_h=H-h\sum_i Z_iHh​=H−h∑i​Zi​ has a unique ground state and gap E1−E0≥h/4E_1-E_0\geq h/4E1​−E0​≥h/4, independently of system size and coupling strengths. The proof retains separate field variables at every imaginary-time Trotter layer. Lee–Yang zero-freeness makes the log-partition function analytic, while a Taylor-support lemma forces contributions connecting distant layers to have high degree. Analytic bounds then yield decay of Trotterized correlations. An adaptive Trotter step controls the continuum limit, establishing connected imaginary-time correlation decay for every ZS=∏i∈SZiZ_S=\prod_{i\in S}Z_iZS​=∏i∈S​Zi​ at rate h/2h/2h/2, with polynomial prefactors. A spectral representation and the completeness of diagonal Pauli products imply the gap and uniqueness.
The gap enables polynomial-time adiabatic ground-energy estimation to inverse-polynomial additive accuracy, placing EPR-Hamiltonian energy estimation and bipartite quantum Max-Cut in BQP\mathsf{BQP}BQP. The theorem requires a nonzero field and gives no gap bound at h=0h=0h=0; energy estimation there uses a small-field approximation. No classical hardness or quantum advantage is established, and containment of the sign-problem-free EPR case in BPP\mathsf{BPP}BPP remains possible.
We consider a broad family of Heisenberg-type quantum spin systems that were studied by Suzuki and Fisher in 1971. This family includes, for example, the Heisenberg antiferromagnet on any bipartite graph. The ground and Gibbs states of these models are Lee-Yang tensors with radius 1: they are associated with multilinear polynomials that possess an extraordinary zero-freeness property inside the unit polydisk in the complex plane. For each Hamiltonian in this family we show that the spectral gap between the first-excited and ground-state energies is lower bounded by 2μ2μ2μ, where μμμ is the magnetic field along the ZZZ direction. Using this result we obtain an efficient quantum adiabatic algorithm for the ground energy of any model in this family. The proof is based on a new inequality that relates the spectral gap of a positive semidefinite operator to its Lee-Yang radius -- a quantitative strengthening of prior work of the authors that may find applications elsewhere.
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.
We provide a geometric description of Lee-Yang tensors, building upon several recent results of Bravyi, Gosset, Liu, and Wong. By applying tools from complex Perron-Frobenius theory, we generalize one of their key lemmas to bound the spectral ratio of non-normal operators with Lee-Yang radius greater than one. As a direct consequence, we provide a local-to-global principle for bounding spectral gaps of Suzuki-Fisher Hamiltonians, a family of 222-local Hamiltonians, using only the Lee-Yang radius of each local operator to bound the gap from below. On the other hand, we develop a simple perturbative technique to show that the Gibbs states of kkk-local Hamiltonians, k≥3k \geq 3k≥3, cannot have Lee-Yang radius at least one for all inverse temperatures β≥0β\geq 0β≥0. We apply this same method to conclude that Suzuki-Fisher Hamiltonians are precisely those Hamiltonians whose Gibbs states are Lee-Yang at all inverse temperatures β≥0β\geq 0β≥0.
A complex tensor with nnn binary indices can be identified with a multilinear polynomial in nnn complex variables. We say it is a Lee-Yang tensor with radius rrr if the polynomial is nonzero whenever all variables lie in the open disk of radius rrr. In this work we study quantum states and observables which are Lee-Yang tensors when expressed in the computational basis. We first review their basic properties, including closure under tensor contraction and certain quantum operations. We show that quantum states with Lee-Yang radius r>1r > 1r>1 can be prepared by quasipolynomial-sized circuits. We also show that every Hermitian operator with Lee-Yang radius r>1r > 1r>1 has a unique principal eigenvector. These results suggest that r=1r = 1r=1 is a key threshold for quantum states and observables. Finally, we consider a family of two-local Hamiltonians where every interaction term energetically favors a deformed EPR state ∣00⟩+s∣11⟩|00\rangle + s|11\rangle∣00⟩+s∣11⟩ for some 0≤s≤10 \leq s \leq 10≤s≤1. We numerically investigate this model and find that on all graphs considered the Lee-Yang radius of the ground state is at least r=1/sr = 1/\sqrt{s}r=1/s​ while the spectral gap between the two smallest eigenvalues is at least 1−s21-s^21−s2. We conjecture that these lower bounds hold more generally; in particular, this would provide an efficient quantum adiabatic algorithm for the quantum Max-Cut problem on uniformly weighted bipartite graphs.
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 edge of the gap in the distribution of Yang-Lee zeros at H=iH0(T)H=i{H}_{0}(T)H=iH0​(T) on the imaginary magnetic field axis in ferromagnets above Tc{T}_{c}Tc​ is essentially a critical point. In terms of the edge exponents $\ensuremath{\delta}$ and $\ensuremath{\eta}$, the density of zeros obeys $\mathcal{G}({H}^{\ensuremath{'}\ensuremath{'}})\ensuremath{\sim}{[{H}^{\ensuremath{'}\ensuremath{'}}\ensuremath{-}{H}_{0}(T)]}^{\ensuremath{\sigma}}$, with $\ensuremath{\sigma}=\frac{1}{\ensuremath{\delta}}=\frac{(d\ensuremath{-}2+\ensuremath{\eta})}{({d}_{2}\ensuremath{-}\ensuremath{\eta})}$. Classical behavior ($\ensuremath{\sigma}=\frac{1}{2}$) occurs for $d>{d}^{\ifmmode\times\else\texttimes\fi{}}=6$. The appropriate field-theoretic renormalization group entails a $w{\ensuremath{\phi}}^{3}$ coupling and, with $\ensuremath{\epsilon}=6\ensuremath{-}d>~0$, yields $\ensuremath{\eta}\ensuremath{\approx}\ensuremath{-}\frac{\ensuremath{\epsilon}}{9}$ for all $n<\ensuremath{\infty}$. This correlates well with refined series estimates for d=2d=2d=2 and d=3d=3d=3 and with exact results for d=1d=1d=1 ($\ensuremath{\eta}=\ensuremath{-}1$).