Root-of-unity completeness and Temperley-Lieb structure of the XXZ Bethe ansatz, and analytic limiting curves of six-vertex partition-function zeros

The 2019 exact computation of the torus partition function of the isotropic six-vertex model by Jacobsen, Jiang and Zhang - Bethe ansatz equations solved through Gröbner bases, primary decomposition over cyclotomic extensions, and Galois theory - came with an explicit program of open questions: extend the algebro-geometric method to the quantum-deformed (XXZ) chain, settle the completeness and physicality of Bethe solutions at roots of unity, relate the algebraic decomposition of the solution space to the representation theory of the affine Temperley-Lieb (aTL) algebra, characterize the limiting curves of partition-function zeros analytically, generalize to other topologies, and prove the momentum-sector ideals primary. Follow-up work by the same and adjacent groups resolved the generic-q extension (including the aTL correspondence and the twist quantum number) and the cylinder geometry at the isotropic point, and produced a numerically validated constrained-Q-system classification of all physical solutions at roots of unity together with an exact root-of-unity torus partition function. The open core that survives, and that this problem states, is the mathematical completion of that program: (1) a rigorous proof or refutation of completeness of the primitive-plus-descendant classification of the periodic XXZ Bethe ansatz at roots of unity, which currently rests on a conjectural, numerically tested counting algorithm; (2) the ideal-theoretic geometry of the rationalized root-of-unity equations - the positive-dimensional solution varieties, their primary decomposition, and the fate of the momentum-sector and Galois bookkeeping; (3) the derivation, from the Bethe-solution algebra, of the aTL gluings and Jordan cells expected at roots of unity; and (4) an analytic theory of the condensation curves of partition-function zeros, including their dependence on the deformation parameter, which are presently traced only numerically.

Background

The six-vertex model is an exactly solved lattice model of statistical mechanics in which ice-type rules allow exactly six local vertex configurations on a square lattice; through the transfer-matrix formalism it is equivalent to the spin-1/2 XXZ quantum spin chain, whose anisotropy $\Delta$ is customarily encoded in a deformation parameter $q$ via $\Delta=(q+q^{-1})/2$. The isotropic point $q=1$ ($\Delta=1$) is the Heisenberg XXX chain. Eigenvalues of the transfer matrix are parametrized by sets of Bethe roots solving the Bethe ansatz equations (BAE), a coupled system of algebraic equations. A long-standing issue is the completeness problem: whether the physical solutions of the BAE account for the full Hilbert space, and which solutions must be discarded as unphysical. For the periodic XXX chain the physical solutions divide into regular solutions and singular physical solutions, and the counting was verified numerically by Hao, Nepomechie and Sommese. Marboe and Volin reformulated the rational BAE as a Q-system: polynomial QQ-relations arranged on a Young diagram whose solutions are automatically physical. This turned the solution space into an object of computational algebraic geometry: the BAE generate an ideal in a polynomial ring whose Gröbner basis, quotient ring, and companion matrices encode all solutions and the transfer-matrix eigenvalues algebraically, without solving for individual roots. The algebraic-geometry and Bethe-ansatz program of Jiang and Zhang exploits this. For the torus geometry, Jacobsen, Jiang and Zhang (2019) decomposed the XXX solution space into momentum-sector ideals $I_{M,K,\ell}$ - with $M$ the chain length, $K$ the magnon number, and $\ell$ labeling the total momentum $2\pi\ell/M$, plus a singular sector $\ell=\infty$ - over the algebraic extension $\mathbb{F}_M=\mathbb{Q}(i,\xi_M)$ with $\xi_M=e^{2\pi i/M}$, organized the sectors into Galois orbits so that only one representative per orbit needs computing, and expressed the torus partition function $Z_{M,N}$ as a sum of traces of powers of companion matrices. Closed forms valid for any $N$ were obtained for $M\le 6$ and exact results up to $M=14$. The zeros of these polynomials in the spectral parameter condense onto curves in the partial thermodynamic limit (one lattice direction sent to infinity), a phenomenon governed by the Beraha-Kahane-Weiss (BKW) theorem - the statement that a finite sum $\sum_i c_i(x)\lambda_i(x)^N$ has zeros accumulating where the dominant terms tie in modulus - which underlies Lee-Yang-type zero analyses but yields the accumulation curves only numerically in general. Away from the isotropic point the BAE are trigonometric. After a rationalizing change of variables the same machinery applies, but two regimes must be distinguished. For generic $q$ (not a root of unity) completeness is a mild deformation of the XXX case. When $q$ is a root of unity, $\Delta=\cos(\pi s/p)$, the solution space degenerates: Fabricius-McCoy (FM) strings - bound-state strings of Bethe roots with zero string deviation - can be appended to any solution without changing the energy, producing apparent continua of solutions, exponential degeneracies of the transfer matrix, and descendant-tower structures attributed to an enhanced symmetry (conjecturally a loop algebra or the Onsager algebra, the infinite-dimensional Lie algebra underlying Onsager's exact solution of the two-dimensional Ising model); correspondingly, the Gröbner bases develop vanishing leading coefficients and the quotient rings cease to be zero-dimensional. On the representation-theory side, the six-vertex model is governed by the affine Temperley-Lieb (aTL) algebra, a diagrammatic algebra generated by overlapping arcs acting on link patterns; its standard modules $\mathcal{W}_{j,z}$ are labeled by the number $j$ of through-lines (defect lines) and a twist parameter $z$ recording how the defects spiral around the periodic direction. These modules are irreducible for generic $q$ but become reducible-but-indecomposable at roots of unity, with characteristic gluings and Jordan cells (nontrivial nilpotent action) that feed logarithmic conformal field theories in the continuum limit. In 2022, Boehm, Jacobsen, Jiang and Zhang established the precise correspondence between the aTL module decomposition and the algebraic decomposition of Bethe solutions at generic $q$ with a twist, computed torus partition functions of the loop and Potts models by both geometric-algebra and algebraic-geometry methods, and showed that the resonance condition where $\mathcal{W}_{j,z}$ becomes reducible appears as a degeneration of the Gröbner bases - while leaving the root-of-unity fine structure explicitly open. What remains is the mathematical completion of this program at roots of unity and for the zero curves, as stated below.

Problem Statement

The common objective is to complete, mathematically, the algebro-geometric theory of the solution space of the Bethe ansatz equations (BAE) of the periodic spin-1/2 XXZ chain (equivalently, the q-deformed six-vertex model) at roots of unity, together with the analytic theory of the zeros of the partition function that this solution space produces - bringing the root-of-unity regime to the same status that the 2019 program gave the isotropic case and follow-up work gave generic $q$. The interlocking deliverables that remain open are:

1. Completeness at roots of unity, rigorously. Prove (or refute by an explicit finite counterexample) that the physical solutions of the periodic XXZ BAE at $q=e^{i\pi s/p}$ with $s$ and $p$ coprime consist exactly of the primitive solutions (of the three types: without roots at infinity, with roots at $+\infty$, or with roots at $-\infty$) together with their descendant towers generated by Fabricius-McCoy strings and infinity pairs, for every chain length, magnon number, and root of unity - the constrained-Q-system classification of Hou, Jiang and Miao that is currently validated by extensive numerics and exact finite-size computations but rests on a conjectural counting algorithm.

2. Ideal-theoretic structure of the degenerating solution varieties. At roots of unity the Gröbner bases of the rationalized BAE/Q-system ideals develop vanishing leading coefficients and the quotient rings become positive-dimensional affine varieties. Characterize these ideals: their primary decomposition, dimensions and components, how the continuum of FM-string centers parametrizes them, and what replaces the momentum-sector ideal decomposition and the Galois-orbit bookkeeping of the isotropic case. As an auxiliary mathematical refinement in the isotropic case, prove (or disprove) that each momentum-sector ideal $I_{M,K,\ell}$, including the singular sector, is primary over $\mathbb{F}_M=\mathbb{Q}(i,\xi_M)$ - that is, admits no further nontrivial decomposition over this algebraic extension - and complete the proofs of the Galois classification of the cyclotomic extensions that the source paper announced as future work.

3. Representation-theoretic dictionary. Derive from the algebraic structure of the Bethe-solution space the affine Temperley-Lieb structures at roots of unity: the correspondence between primitive states and the standard modules $\mathcal{W}_{j,z}$, and the way descendant towers realize the gluings of standard modules and the Jordan cells (reducible-but-indecomposable structure) known for the periodic aTL at roots of unity; identify which features of the conjectured enhanced symmetry (loop algebra or Onsager algebra) are visible in the primary decomposition.

4. Analytic limiting curves of partition-function zeros. Replace the numerical Beraha-Kahane-Weiss-type tracing of the condensation curves of the six-vertex torus (and cylinder) partition-function zeros in the partial thermodynamic limit by analytic characterizations: closed-form or implicit algebraic equations for the curves, the scaling of the zero density near curves and branch points (T-points, i.e. points where three or more accumulation curves meet), and the effect of q-deformation on the curves; only isolated curve components have so far been obtained analytically, in the full thermodynamic limit.

A resolution of this problem requires all four deliverables: deliverables 1-3 jointly determine how the root-of-unity solution space decomposes, rigorously and representation-theoretically, while deliverable 4 is the analytic counterpart for the partition-function output of that solution space; the isotropic primary-decomposition and Galois clause is part of deliverable 2 rather than an optional extension.

References

1. Jesper Lykke Jacobsen, Yunfeng Jiang, Yang Zhang, Torus partition function of the six-vertex model from algebraic geometry, JHEP 03 (2019) 152, arXiv:1812.00447, doi:10.1007/JHEP03(2019)152, https://arxiv.org/abs/1812.00447
2. Yunfeng Jiang, Yang Zhang, Algebraic geometry and Bethe ansatz. Part I. The quotient ring for BAE, JHEP 03 (2018) 087, arXiv:1710.04693, doi:10.1007/JHEP03(2018)087, https://arxiv.org/abs/1710.04693
3. Christian Marboe, Dmytro Volin, Fast analytic solver of rational Bethe equations, J. Phys. A 50 (2017) 204002, arXiv:1608.06504, doi:10.1088/1751-8121/aa6b88, https://arxiv.org/abs/1608.06504
4. Wenrui Hao, Rafael I. Nepomechie, Andrew J. Sommese, Completeness of solutions of Bethe's equations, Phys. Rev. E 88 (2013) 052113, arXiv:1308.4645, doi:10.1103/PhysRevE.88.052113, https://arxiv.org/abs/1308.4645
5. A. M. Gainutdinov, Wenrui Hao, Rafael I. Nepomechie, Andrew J. Sommese, Counting solutions of the Bethe equations of the quantum group invariant open XXZ chain at roots of unity, J. Phys. A 48 (2015) 494003, arXiv:1505.02104, doi:10.1088/1751-8113/48/49/494003, https://arxiv.org/abs/1505.02104
6. Zoltan Bajnok, Etienne Granet, Jesper Lykke Jacobsen, Rafael I. Nepomechie, On Generalized Q-systems, JHEP 03 (2020) 177, arXiv:1910.07805, doi:10.1007/JHEP03(2020)177, https://arxiv.org/abs/1910.07805
7. Rafael I. Nepomechie, Q-systems with boundary parameters, J. Phys. A 53 (2020) 294001, arXiv:1912.12702, doi:10.1088/1751-8121/ab9386, https://arxiv.org/abs/1912.12702
8. Zoltan Bajnok, Jesper Lykke Jacobsen, Yunfeng Jiang, Rafael I. Nepomechie, Yang Zhang, Cylinder partition function of the 6-vertex model from algebraic geometry, JHEP 06 (2020) 169, arXiv:2002.09019, doi:10.1007/JHEP06(2020)169, https://arxiv.org/abs/2002.09019
9. Janko Boehm, Jesper Lykke Jacobsen, Yunfeng Jiang, Yang Zhang, Geometric algebra and algebraic geometry of loop and Potts models, JHEP 05 (2022) 068, arXiv:2202.02986, doi:10.1007/JHEP05(2022)068, https://arxiv.org/abs/2202.02986
10. Yuan Miao, Jules Lamers, Vincent Pasquier, On the Q operator and the spectrum of the XXZ model at root of unity, SciPost Phys. 11 (2021) 067, arXiv:2012.10224, doi:10.21468/SciPostPhys.11.3.067, https://arxiv.org/abs/2012.10224
11. Jue Hou, Yunfeng Jiang, Yuan Miao, Rational Q-systems at root of unity I. Closed chains, SciPost Phys. 16 (2024) 129, arXiv:2310.14966, doi:10.21468/SciPostPhys.16.5.129, https://arxiv.org/abs/2310.14966
12. Yunfeng Jiang, Rui Wen, Yang Zhang, Exact quench dynamics from algebraic geometry, Phys. Rev. E 108 (2023) 024128, arXiv:2109.10568, doi:10.1103/PhysRevE.108.024128, https://arxiv.org/abs/2109.10568
13. Arthur Hutsalyuk, Yunfeng Jiang, Balazs Pozsgay, Hefeng Xu, Yang Zhang, Exact Spin Correlators of Integrable Quantum Circuits from Algebraic Geometry, SciPost Phys. 19 (2025) 003, arXiv:2405.16070, doi:10.21468/SciPostPhys.19.1.003, https://arxiv.org/abs/2405.16070
14. Yunfeng Jiang, Yi-Chao Liu, Yuan Miao, Zi-Xi Tan, Rational Q-systems for integrable spin chains without U(1) symmetry, SciPost Phys. 20 (2026) 175, arXiv:2512.01551, doi:10.21468/SciPostPhys.20.6.175, https://arxiv.org/abs/2512.01551
15. Yunfeng Jiang, Chang Liu, Yu Wu, Yang Zhang, Gate Parameter Lee-Yang Zeros and Dynamical Phases in Quantum Circuits, arXiv:2605.29838 (2026), https://arxiv.org/abs/2605.29838
16. Jue Hou, Yunfeng Jiang, Rui-Dong Zhu, Spin-s rational Q-system, SciPost Phys. 16 (2024) 113, arXiv:2303.07640, doi:10.21468/SciPostPhys.16.4.113, https://arxiv.org/abs/2303.07640
17. Dmitry Chernyak, Sebastien Leurent, Dmytro Volin, Completeness of Wronskian Bethe Equations for Rational gl(m|n) Spin Chains, Commun. Math. Phys. 391 (2022) 969-1045, arXiv:2004.02865, doi:10.1007/s00220-021-04275-9, https://arxiv.org/abs/2004.02865
18. Cyril Closset, Osama Khlaif, Twisted indices, Bethe ideals and 3d N=2 infrared dualities, JHEP 05 (2023) 148, arXiv:2301.10753, doi:10.1007/JHEP05(2023)148, https://arxiv.org/abs/2301.10753