Ground state of the supersymmetric XYZ spin chain: unconditional proof of the Bazhanov–Mangazeev, Razumov–Stroganov and Zinn-Justin eigenvector, sum-rule and positivity conjectures

The spin-$1/2$ XYZ Hamiltonian on a periodic chain of odd length $L=2n+1$ is supersymmetric on the coupling line $J_xJ_y+J_xJ_z+J_yJ_z=0$ (eight-vertex crossing parameter $\eta=\pi/3$). The doubly degenerate ground-state eigenvalue $\Theta_n=(a+b)^{2n+1}$ of the transfer matrix is now a theorem, but the ground-state eigenvector itself is known only conjecturally. Bazhanov and Mangazeev, Razumov and Stroganov, Mangazeev and Bazhanov, and Zinn-Justin conjectured that the components of this eigenvector, the sums of its components, and its square norm are given by special polynomials ($s_n$, $\bar{s}_n$, $p_n$, $q_n$) that Rosengren unified in a family $T_n^{(\mathbf{k})}$ of symmetric polynomials satisfying a non-stationary Lamé equation and identified with tau functions of the Painlevé VI equation obtained from Picard's algebraic solution by a lattice of Bäcklund transformations. Only the tau-function side of these identifications is rigorous: the identification with eigenvector components is still empirical, Zinn-Justin's square-norm sum rule was proved only under an unproved conjecture on the inhomogeneous eigenvalue problem, and Hagendorf–Rosengren's finite-size nearest-neighbour correlation formula holds only under an unproved analyticity assumption on Baxter's Q-operator. The problem is to make this conjecture cluster (the eigenvector-component, sum-rule and positivity conjectures) unconditional: existence and characterization of the inhomogeneous eigenvector, closed-form component formulas and the sum rule, and the conjectured positivity of the polynomial coefficients with a combinatorial interpretation.

Background

The square-lattice eight-vertex model of Baxter assigns Boltzmann weights $a,b,c,d$ to the sixteen arrow configurations on a vertex (organized into Baxter's $4\times 4$ R-matrix, the local solution of the Yang--Baxter equation that makes successive transfer matrices commute); its row-to-row transfer matrix generates a commuting family, and a logarithmic derivative of the transfer matrix gives the spin-$1/2$ XYZ Hamiltonian, the anisotropic Heisenberg chain with couplings $J_x,J_y,J_z$. The trigonometric (six-vertex) limit of this machinery is the XXZ chain with anisotropy $\Delta$. On the parameter line $(a^2+ab)(b^2+ab)=(c^2+ab)(d^2+ab)$, equivalently $J_xJ_y+J_xJ_z+J_yJ_z=0$ (the elliptic analogue of $\Delta=-1/2$, crossing parameter $\eta=\pi/3$), the chain carries a non-local $\mathcal{N}=(2,2)$ lattice supersymmetry discovered by Hagendorf and Fendley (the Hamiltonian is the anticommutator of two fermionic nilpotent charges), and the line is called supersymmetric or combinatorial.

For a periodic chain of odd length $L=2n+1$ on this line, Stroganov conjectured in 2001 and Hagendorf and Liénardy proved in 2018 that the transfer matrix has the doubly degenerate eigenvalue $\Theta_n=(a+b)^{2n+1}$, which is the largest eigenvalue for positive weights; the corresponding two-dimensional eigenspace is the ground state of the XYZ Hamiltonian. Supersymmetry fixes the ground-state energy and the degeneracy, but not the components of the eigenvector — the last missing piece of the exact finite-size solution.

The conjecture cluster concerns those components. Baxter's Q-operator is an auxiliary transfer-matrix-like object whose eigenvalues satisfy, together with the transfer-matrix eigenvalues, the functional TQ relation underlying his solution of the eight-vertex model. Bazhanov and Mangazeev (2005) found that the ground-state Q-operator eigenvalue on the supersymmetric line is expressed through polynomials $s_n,\bar{s}_n$ that appear to have positive integer coefficients. Mangazeev and Bazhanov (2010) conjectured closed formulas for components of the ground-state eigenvector in terms of further polynomials $p_n,q_n$ (related to $s_n$ by factorizations such as $s_{2n+1}(y^2)=s_{2n+1}(0)\,p_n(y)p_n(-y)$), and a square-norm sum rule: with the couplings parametrized by the combinatorial parameter $\zeta$ as $J_x=2(1+\zeta)/(\zeta^2+3)$, $J_y=2(1-\zeta)/(\zeta^2+3)$, $J_z=(\zeta^2-1)/(\zeta^2+3)$ (equivalently $\zeta=cd/ab$ in terms of the eight-vertex weights), and with the ground-state vector $\Psi_-$ normalized so that its components are integer-coefficient polynomials in $\zeta$ with no common polynomial factor and a fixed leading component at $\zeta=0$, the square norm is $|\Psi_-|^2=(4/3)^n\,\zeta^{n(n+1)}\,s_n(\zeta^{-2})\,s_{-n-1}(\zeta^{-2})$. Razumov and Stroganov (2010) conjectured companion component and component-sum rules measuring the antiferromagneticity of the chain, in direct analogy with the Razumov–Stroganov conjectures for the six-vertex model at $\Delta=-1/2$, which were proved by Cantini and Sportiello and which identify ground-state components with enumerations of alternating sign matrices — square arrays of $0$ and $\pm1$ with row and column sums equal to $1$ and partial sums in $\{0,1\}$, the central objects of integrable combinatorics. In the trigonometric case the proofs rest on integral representations of eigenvectors via the quantum Knizhnik–Zamolodchikov equation (a system of difference equations whose contour-integral solutions yield the components); no elliptic counterpart of these representations is known, which is the essential new difficulty here.

Rosengren (2013–2015) unified these polynomials in a family $T_n^{(\mathbf{k})}$ of symmetric polynomials indexed by $\mathbf{k}\in\mathbb{Z}^4$ with $|\mathbf{k}|+m=2n$ ($m$ the number of variables): after a change of variables they satisfy a non-stationary Schrödinger equation with elliptic (Lamé) potential; specializations of $T_n^{(\mathbf{k})}$ are tau functions (analytic functions whose logarithmic derivatives encode the invariants of an isomonodromic deformation) of the sixth Painlevé equation — the isomonodromic deformation whose tau functions here arise from Picard's algebraic seed solution by a four-dimensional lattice of Bäcklund transformations (transformations mapping solutions of the equation to new solutions); and in special cases $T_n^{(\mathbf{k})}$ agrees with characters (weight-generating functions of representations) of the affine Lie algebra $C_{2n}^{(1)}$. Rosengren proved rigorously the identifications on the tau-function side (recursions, Painlevé VI structure) and the identification of the Q-operator-eigenvalue polynomials and of three-colour domain-wall partition functions with $T$-specializations. What could not be proved is the identification with eigenvector components of the Hamiltonian, which remains empirical; Rosengren explicitly poses this as the open problem.

Zinn-Justin (2013) initiated the rigorous analysis of the inhomogeneous model on the combinatorial line and derived the square-norm sum rule from a single unproved conjecture (labelled Conjecture 2.1 in later work): the space of solutions of the inhomogeneous transfer-matrix eigenvalue problem with the spin-reversal constraint is spanned by a vector whose components are entire functions of the inhomogeneity parameters and satisfy quasi-periodicity and exchange relations. Brasseur and Hagendorf (2021) derived sum rules and many components conditionally on this conjecture, but showed that even granting it the comparison with the Bazhanov–Mangazeev and Razumov–Stroganov vectors remains non-rigorous because the normalization conventions differ and the exact link between the polynomial families is still lacking. Hagendorf and Rosengren (2024) proved a nearest-neighbour correlation formula in terms of $s_n,\bar{s}_n$, but only under an assumption (their Assumption 3.4) that Baxter's TQ relation admits a solution analytic in the crossing parameter $\eta$ near $\pi/3$, the correlation formula itself having been checked numerically for chains up to length $L=11$; they state that the major obstruction is that the ground-state eigenvectors are not known explicitly, and that integral formulas for eight-vertex eigenvectors would presumably complete the proofs of the ground-state properties conjectured by Bazhanov–Mangazeev, Razumov–Stroganov and Zinn-Justin.

Problem Statement

Establish, without unproved auxiliary assumptions, the conjectured exact structure of the ground state of the periodic supersymmetric XYZ spin chain of odd length $L=2n+1$ (equivalently, the eight-vertex model on its combinatorial line $\eta=\pi/3$), for arbitrary $n$ and all admissible weights (all weight-tuples on the combinatorial line). The required deliverables form one program — an unconditional, explicit finite-size theory of the ground state — and are:

(1) Existence and characterization of the inhomogeneous ground-state eigenvector. Prove Zinn-Justin's conjecture (Conjecture 2.1 in the later literature): the space of solutions of the inhomogeneous transfer-matrix eigenvalue problem with the inhomogeneous generalization of the eigenvalue $\Theta_n$ and the spin-reversal constraint (the vector is an eigenvector of the operator $F$ reversing all spins, $F|\Psi\rangle=(-1)^{n}|\Psi\rangle$) is spanned by a vector whose components are entire functions of the inhomogeneity parameters, are generically non-vanishing with no common factor, and obey the quasi-periodicity relations under shifts of any inhomogeneity by $\pi$ and by $2\pi\tau$ (where $p=e^{\mathrm{i}\pi\tau}$ is the elliptic nome of the weight parametrization), together with the exchange (R-matrix) relations.

(2) Component formulas and sum rules. Deduce unconditionally that, in the homogeneous limit, the components of the ground-state basis vectors and the sums of their components are polynomials in the combinatorial parameter $\zeta$ with integer coefficients, equal to the polynomials conjectured by Bazhanov–Mangazeev, Mangazeev–Bazhanov and Razumov–Stroganov — that is, identified with specializations of Rosengren's $T_n^{(\mathbf{k})}$ family (equivalently the $s_n,\bar{s}_n,p_n,q_n$ polynomials, or tau functions of Painlevé VI) — with all normalization constants fixed; prove the square-norm sum rule $|\Psi_-|^2=(4/3)^n\,\zeta^{n(n+1)}\,s_n(\zeta^{-2})\,s_{-n-1}(\zeta^{-2})$ (in the normalization and $\zeta$-parametrization recalled in the Background) and the factorization identities such as $s_{2n+1}(y^2)=s_{2n+1}(0)\,p_n(y)p_n(-y)$.

(3) Positivity and combinatorial interpretation. Prove that the polynomials in (2) have non-negative integer coefficients (currently only weaker positivity statements are known for $p_n,q_n$ via three-colour-model identities) and construct a combinatorial model or bijection explaining them, extending the alternating-sign-matrix combinatorics of the trigonometric case to the elliptic regime.

The literature notes that integral representations for the eight-vertex eigenvector components, analogous to the quantum-Knizhnik–Zamolodchikov integral formulas for the six-vertex model, would presumably allow one to complete the proofs of many properties of the ground state conjectured by Bazhanov–Mangazeev, Razumov–Stroganov and Zinn-Justin, and might also be useful for proving Hagendorf–Rosengren's Assumption 3.4 — the existence of a Q-operator solution of Baxter's TQ relations analytic in the crossing parameter $\eta$ in a neighbourhood of $\pi/3$ — which would in turn make their nearest-neighbour correlation formula (Theorem 2.3) unconditional. That correlation formula is a downstream consequence of the program of deliverables (1)–(3) rather than a separate required deliverable; the required completion boundary is exactly deliverables (1)–(3), so proving their Assumption 3.4 without establishing them does not resolve the task.

References

1. H. Rosengren, Special polynomials related to the supersymmetric eight-vertex model: A summary, Comm. Math. Phys. 340 (2015) 1143–1170. DOI: 10.1007/s00220-015-2439-0, arXiv:1503.02833, https://arxiv.org/abs/1503.02833
2. V. V. Bazhanov, V. V. Mangazeev, Eight-vertex model and non-stationary Lamé equation, J. Phys. A: Math. Gen. 38 (2005) L145–L153. DOI: 10.1088/0305-4470/38/8/l01, https://doi.org/10.1088/0305-4470/38/8/l01
3. V. V. Bazhanov, V. V. Mangazeev, The eight-vertex model and Painlevé VI, J. Phys. A: Math. Gen. 39 (2006) 12235–12243. DOI: 10.1088/0305-4470/39/39/s15, https://doi.org/10.1088/0305-4470/39/39/s15
4. V. V. Mangazeev, V. V. Bazhanov, The eight-vertex model and Painlevé VI equation II: eigenvector results, J. Phys. A: Math. Theor. 43 (2010) 085206. DOI: 10.1088/1751-8113/43/8/085206, arXiv:0912.2163, https://arxiv.org/abs/0912.2163
5. A. V. Razumov, Yu. G. Stroganov, A possible combinatorial point for the XYZ spin chain, Theor. Math. Phys. 164 (2010) 977–991. DOI: 10.1007/s11232-010-0078-3, arXiv:0911.5030, https://arxiv.org/abs/0911.5030
6. Yu. Stroganov, The importance of being odd, J. Phys. A: Math. Gen. 34 (2001) L179–L185. DOI: 10.1088/0305-4470/34/13/104, https://doi.org/10.1088/0305-4470/34/13/104
7. P. Zinn-Justin, Sum rule for the eight-vertex model on its combinatorial line, in: Symmetries, Integrable Systems and Representations, Springer Proc. Math. Stat. 40, Springer, London (2013), pp. 599–637. DOI: 10.1007/978-1-4471-4863-0_26, arXiv:1202.4420, https://arxiv.org/abs/1202.4420
8. H. Rosengren, Special polynomials related to the supersymmetric eight-vertex model. I. Behaviour at cusps, arXiv:1305.0666 (2013), https://arxiv.org/abs/1305.0666
9. H. Rosengren, Special polynomials related to the supersymmetric eight-vertex model. II. Schrödinger equation, arXiv:1312.5879 (2013), https://arxiv.org/abs/1312.5879
10. H. Rosengren, Special polynomials related to the supersymmetric eight-vertex model. III. Painlevé VI equation, arXiv:1405.5318 (2014), https://arxiv.org/abs/1405.5318
11. C. Hagendorf, P. Fendley, The eight-vertex model and lattice supersymmetry, J. Stat. Phys. 146 (2012) 1122–1155. DOI: 10.1007/s10955-012-0430-0, arXiv:1109.4090, https://arxiv.org/abs/1109.4090
12. C. Hagendorf, J. Liénardy, On the transfer matrix of the supersymmetric eight-vertex model. I. Periodic boundary conditions, J. Stat. Mech. (2018) 033106. DOI: 10.1088/1742-5468/aab01d, arXiv:1711.04397, https://arxiv.org/abs/1711.04397
13. C. Hagendorf, J. Liénardy, On the transfer matrix of the supersymmetric eight-vertex model. II. Open boundary conditions, J. Stat. Mech. (2020) 033104. DOI: 10.1088/1742-5468/ab7748, arXiv:1911.03348, https://arxiv.org/abs/1911.03348
14. L. Hietala, A combinatorial description of certain polynomials related to the XYZ spin chain, SIGMA 16 (2020) 101. DOI: 10.3842/SIGMA.2020.101, arXiv:2004.09924, https://arxiv.org/abs/2004.09924
15. L. Hietala, A combinatorial description of certain polynomials related to the XYZ spin chain. II. The polynomials $p_n$, SIGMA 18 (2022) 036. DOI: 10.3842/SIGMA.2022.036, arXiv:2104.04651, https://arxiv.org/abs/2104.04651
16. S. Brasseur, C. Hagendorf, Sum rules for the supersymmetric eight-vertex model, J. Stat. Mech. (2021) 023102. DOI: 10.1088/1742-5468/abda28, arXiv:2009.14077, https://arxiv.org/abs/2009.14077
17. C. Hagendorf, H. Rosengren, Nearest-neighbour correlation functions for the supersymmetric XYZ spin chain and Painlevé VI, Comm. Math. Phys. 405 (2024), article 97. DOI: 10.1007/s00220-024-04977-w, arXiv:2208.13533, https://arxiv.org/abs/2208.13533
18. L. Cantini, A. Sportiello, Proof of the Razumov–Stroganov conjecture, J. Combin. Theory Ser. A 118 (2011) 1549–1574. DOI: 10.1016/j.jcta.2011.01.007, https://doi.org/10.1016/j.jcta.2011.01.007
19. H. Rosengren, Elliptic pfaffians and solvable lattice models, J. Stat. Mech. (2016) 083106. DOI: 10.1088/1742-5468/2016/08/083106, arXiv:1605.02915, https://arxiv.org/abs/1605.02915
20. X. Zhang, Hypotheses regarding Baxter's $T-Q$ relation for the periodic XYZ chain, arXiv:2404.12615 (2024), https://arxiv.org/abs/2404.12615
21. X. Zhang, A. Klümper, V. Popkov, Pedestrian's way to Baxter's Bethe ansatz for the periodic XYZ chain, Phys. Rev. B 109 (2024) 115411. DOI: 10.1103/PhysRevB.109.115411, arXiv:2312.00161, https://arxiv.org/abs/2312.00161
22. P. Fendley, S. Gehrmann, E. Vernier, F. Verstraete, XYZ integrability the easy way, J. Stat. Phys. 193 (2026) 79. DOI: 10.1007/s10955-026-03645-8, arXiv:2511.04674, https://arxiv.org/abs/2511.04674
23. G. Niccoli, V. Terras, The open XYZ spin 1/2 chain: separation of variables and scalar products for boundary fields related by a constraint, SciPost Phys. 19 (2025) 090. DOI: 10.21468/SciPostPhys.19.4.090, arXiv:2402.04112, https://arxiv.org/abs/2402.04112
24. Z. Xin, J. Cao, W.-L. Yang, Y. Wang, Exact physical quantities of the XYZ spin chain in the thermodynamic limit, J. High Energy Phys. 02 (2025) 086. DOI: 10.1007/JHEP02(2025)086, arXiv:2411.12200, https://arxiv.org/abs/2411.12200