Exact finite-size correlation functions of the supersymmetric XYZ spin chain: Q-operator analyticity at the supersymmetric point, ground-state eigenvectors, further correlators, and the emergence of Painlevé VI

The spin-1/2 XYZ chain on the supersymmetric coupling line $J_xJ_y+J_xJ_z+J_yJ_z=0$ with periodic boundary conditions and odd length $L=2n+1$ is one of the rare lattice models with exact finite-size results: Stroganov's finite-chain ground-state energy, and, by Hagendorf and Rosengren (2024), nearest-neighbour ground-state correlation functions expressed exactly through the Bazhanov-Mangazeev polynomials $s_n,\bar s_n$, which Rosengren proved to be tau functions of the Painlevé VI equation. That result is conditional on the unproven Assumption 3.4 — analyticity in the crossing parameter $\eta$ near the supersymmetric point $\eta=\pi/3$ of a scalar Q-function satisfying Baxter's TQ relation — and the authors explicitly leave open five linked issues: (i) a proof of this assumption, or an alternative rigorous route to the TQ relation at $\eta=\pi/3$; (ii) explicit finite-$L$ ground-state eigenvector formulas (the eight-vertex analogue of the XXZ integral component formulas); (iii) exact finite-size expressions for further correlators, such as one-point functions and the emptiness-formation probability; (iv) the double-scaling limit $n,\tau/i\to\infty$ with $e^{\pi i\tau}\sim n^{-2/3}$, expected to degenerate Painlevé VI to Painlevé III towards the supersymmetric sine-Gordon model; and (v) a conceptual explanation of why Painlevé VI governs these finite-size observables. None of these five components has been resolved in the literature.

Background

The XYZ spin chain is the spin-1/2 anisotropic Heisenberg chain with Hamiltonian $H=-\frac{1}{2}\sum_{j=1}^{L}(J_x\sigma^x_j\sigma^x_{j+1}+J_y\sigma^y_j\sigma^y_{j+1}+J_z\sigma^z_j\sigma^z_{j+1})$ and periodic boundary conditions, where $\sigma^\alpha_j$ are Pauli matrices acting on site $j$. Via its commuting family of transfer matrices it is equivalent to Baxter's two-dimensional eight-vertex model, whose Boltzmann weights are elliptic functions of the spectral parameter $u$ with nome $\tau$ and crossing parameter $\eta$. Baxter's solution of the eight-vertex model rests on the Q-operator, an auxiliary operator whose eigenvalues $Q(u)$ satisfy the scalar TQ relation $\lambda(u)Q(u)=\phi(u-\eta)Q(u+2\eta)+\phi(u+\eta)Q(u-2\eta)$, where $\lambda(u)$ is the transfer-matrix eigenvalue and $\phi$ an explicit function.

Exact results for such models usually hold only in the infinite-lattice limit; exact finite-size results are rare. Stroganov (2001) observed that on the coupling line

$$J_xJ_y+J_xJ_z+J_yJ_z=0,\qquad J_x+J_y+J_z>0,$$

now called the supersymmetric line, the ground-state energy per site equals $-(J_x+J_y+J_z)/2$ exactly already at odd length $L=2n+1$ with periodic boundary conditions. Hagendorf and Fendley (2012) explained this through lattice supersymmetry — the Hamiltonian is the anticommutator of two nilpotent operators — and Hagendorf and Liénardy (2018) turned the observation into a rigorous proof. The XXZ specialization $J_x=J_y$, i.e. $\Delta=J_z=-1/2$, is the celebrated combinatorial point related to alternating sign matrices (square arrays with entries $0,\pm 1$ whose non-zero entries alternate in sign along every row and column) and the Razumov-Stroganov conjecture, proved by Cantini and Sportiello (2011); a conjectural combinatorial point for the full XYZ chain was proposed by Razumov and Stroganov (2010).

In the elliptic parametrization the supersymmetric line corresponds to $\eta=\pi/3$, a root-of-unity value at which the standard Q-operator machinery degenerates (the shift $u\mapsto u+2\eta$ in the TQ relation has finite order: three shifts equal the period $2\pi$ of $Q$). Bazhanov and Mangazeev (2005, 2006) studied the ground-state Q-operator eigenvalue in Stroganov's odd-length setting and found that at special values of the spectral parameter it is expressed through polynomials $s_n(z)$ and $\bar s_n(z)$ of $z=\zeta^{-2}$, where the couplings are parametrized as $J_x=1+\zeta$, $J_y=1-\zeta$, $J_z=(\zeta^2-1)/2$, with $s_0=s_1=1$ and $\bar s_0=1$, $\bar s_1=3$. Rosengren (2015) proved that these polynomials satisfy Toda-type recursions identifying them with tau functions of the Painlevé VI equation — the most general of the six Painlevé transcendents, a nonlinear second-order ODE whose tau function is a log-derivative of a solution. Conjecturally, the same polynomials describe the ground-state eigenvector itself: its square norm is essentially $s_ns_{-n-1}$, its component with all spins equal is given by $s_n$ and its component with all but one spin equal by $\bar s_n$; these relations are only partially proved (Zinn-Justin 2013; Brasseur and Hagendorf 2021).

Hagendorf and Rosengren (2024) extended this to correlation functions. For the odd-length supersymmetric chain, all nearest-neighbour correlators reduce to a single rational function $f_n$ of a normalized coupling discriminant, and their Theorem 2.3 expresses $f_n$ exactly through the ratio $\bar s_n(\zeta^{-2})\bar s_{-n-1}(\zeta^{-2})/(s_n(\zeta^{-2})s_{-n-1}(\zeta^{-2}))$, with a companion result (Theorem 2.4) interpreting $f_n$ as a Painlevé VI Hamiltonian evaluated on another Painlevé VI solution with shifted parameters. The proof is incomplete: differentiating the TQ relation with respect to $\eta$ at $\eta=\pi/3$ requires a Q-function analytic in $\eta$ near $\pi/3$, which is Assumption 3.4 and remains unproven. Baxter's first Q-operator construction degenerates at $\eta=\pi/3$ (an auxiliary operator is no longer invertible), his second works only for even $L$, and the constructions of Fabricius (2007) and Roan (2007), which do exist at $\eta=\pi/3$, are fixed there and carry no $\eta$-dependence. In the XXZ limit $\Delta=-1/2$ none of this is an obstruction: Stroganov (2001) computed the analogous correlators rigorously, integral formulas for all eigenvector components are known (Razumov, Stroganov and Zinn-Justin 2007), and Cantini (2012) computed the finite-size emptiness-formation probability, i.e. the probability that a block of consecutive spins is perfectly aligned. One-point functions are known in the infinite-lattice limit (Baxter and Kelland 1974; Jimbo, Miwa and Nakayashiki 1993) but not at finite size, and a double-scaling limit $n\to\infty$, $\tau/i\to\infty$ with $e^{\pi i\tau}\sim n^{-2/3}$ is believed to connect the supersymmetric XYZ chain to the supersymmetric sine-Gordon model and to polymers on a cylinder (Fendley and Saleur 1992; Bazhanov and Mangazeev 2005), with Painlevé VI expected to degenerate to Painlevé III; the only available theta-function integral formulas for eigenvectors of an eight-vertex-type model are those of Finch, Weston and Zinn-Justin (2016), and for the eight-vertex-solid-on-solid face model, not the vertex model itself.

Problem Statement

The problem, stated by Hagendorf and Rosengren (2024) in the open problems concluding their introduction and in their Assumption 3.4, is to complete the exact finite-size solution of the supersymmetric XYZ spin chain and to understand structurally its Painlevé VI description. It has five linked components: the authors identify the second as the likely key to the first and third, the fourth extends the third to the scaling regime, and the fifth explains what the others establish computationally.

(1) Rigorous foundations at the supersymmetric point. Prove, or disprove, Assumption 3.4 of Hagendorf-Rosengren: there exists a non-zero function $Q(u;\eta,\tau)$, analytic for $u\in\mathbb{C}$, $\mathrm{Im}\,\tau>0$, and $\eta$ in some neighbourhood of $\pi/3$ (that may depend on $u$ and $\tau$), satisfying Baxter's scalar TQ relation

$$\lambda(u)Q(u)=\phi(u-\eta)Q(u+2\eta)+\phi(u+\eta)Q(u-2\eta)$$

together with the quasi-periodicity $Q(u+2\pi)=Q(-u)=Q(u)$ and $Q(u+2\pi\tau)=e^{-2iL(u+\pi\tau)}Q(u)$ at odd length $L=2n+1$; or otherwise justify rigorously the differentiation of the TQ relation with respect to $\eta$ at $\eta=\pi/3$, thereby removing the conditional status of Theorem 2.3, the exact formula for the nearest-neighbour correlation functions.

(2) Explicit ground-state eigenvectors. Find explicit formulas — integral representations or comparable closed forms — for the finite-$L$ ground-state eigenvector of the supersymmetric XYZ chain at odd $L$, the eight-vertex analogue of the XXZ integral component formulas, and use them to prove the conjectural relations between the eigenvector's square norm and distinguished components and the polynomials $s_n,\bar s_n$.

(3) Further finite-size correlation functions. Derive exact expressions, at finite odd $L$ on the supersymmetric line, for correlation functions beyond nearest neighbours — at minimum the one-point functions and the emptiness-formation probability — and determine whether, like the nearest-neighbour correlators, they are expressible through Painlevé VI tau functions or related special polynomials.

(4) Double-scaling limit. Compute the limit of these finite-size correlation functions as $n\to\infty$ and $\tau/i\to\infty$ with $e^{\pi i\tau}\sim n^{-2/3}$, and demonstrate the expected degeneration of Painlevé VI to Painlevé III, believed to reflect the supersymmetric sine-Gordon model and polymers on a cylinder.

(5) Conceptual explanation. Provide a structural derivation — representation-theoretic, isomonodromic, Lax-pair-based or comparable — explaining why Painlevé VI, its tau functions and its Hamiltonian govern finite-size observables of the supersymmetric XYZ chain, including whether the observed link between quantum and classical Painlevé VI in this setting connects to the quantum Painlevé-Calogero correspondence.

A complete answer grounds the exact finite-size correlation functions rigorously (components 1-3), carries them to the scaling regime (component 4), and supplies the structural explanation (component 5). Substantial partial answers — for instance a proof of Assumption 3.4 alone — constitute meaningful progress on the objective but not its completion.

References

1. Christian Hagendorf and Hjalmar Rosengren, Nearest-Neighbour Correlation Functions for the Supersymmetric XYZ Spin Chain and Painlevé VI, Commun. Math. Phys. 405 (2024), no. 4, article 97, DOI 10.1007/s00220-024-04977-w; preprint arXiv:2208.13533, https://doi.org/10.1007/s00220-024-04977-w
2. V. V. Bazhanov and V. V. Mangazeev, Eight-vertex model and non-stationary Lamé equation, J. Phys. A 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 and V. V. Mangazeev, The eight-vertex model and Painlevé VI, J. Phys. A 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 and V. V. Bazhanov, The eight-vertex model and Painlevé VI equation II: eigenvector results, J. Phys. A 43 (2010) 085206, DOI 10.1088/1751-8113/43/8/085206, https://doi.org/10.1088/1751-8113/43/8/085206
5. Hjalmar Rosengren, Special Polynomials Related to the Supersymmetric Eight-Vertex Model: A Summary, Commun. Math. Phys. 340 (2015) 1143-1170, DOI 10.1007/s00220-015-2439-0, https://doi.org/10.1007/s00220-015-2439-0
6. Yu. G. Stroganov, The importance of being odd, J. Phys. A 34 (2001) L179-L185, DOI 10.1088/0305-4470/34/13/104, https://doi.org/10.1088/0305-4470/34/13/104
7. Yu. G. Stroganov, XXZ spin chain of the asymmetry parameter Δ = -1/2: evaluation of the simplest correlators, Theor. Math. Phys. 129 (2001) 1596-1608, DOI 10.1023/A:1012925110210, https://doi.org/10.1023/A:1012925110210
8. Christian Hagendorf and Paul Fendley, The Eight-Vertex Model and Lattice Supersymmetry, J. Stat. Phys. 146 (2012) 1122-1155, DOI 10.1007/s10955-012-0430-0, https://doi.org/10.1007/s10955-012-0430-0
9. Christian Hagendorf and Jean Liénardy, On the transfer matrix of the supersymmetric eight-vertex model. I. Periodic boundary conditions, J. Stat. Mech. Theory Exp. (2018) 033106, DOI 10.1088/1742-5468/aab01d, https://doi.org/10.1088/1742-5468/aab01d
10. Paul 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, 2013, pp. 599-637, DOI 10.1007/978-1-4471-4863-0_26, https://doi.org/10.1007/978-1-4471-4863-0_26
11. Sandrine Brasseur and Christian Hagendorf, Sum rules for the supersymmetric eight-vertex model, J. Stat. Mech. Theory Exp. (2021) 023102, DOI 10.1088/1742-5468/abda28, https://doi.org/10.1088/1742-5468/abda28
12. Luigi Cantini, Finite size emptiness formation probability of the XXZ spin chain at Δ = -1/2, J. Phys. A 45 (2012) 135207, DOI 10.1088/1751-8113/45/13/135207, https://doi.org/10.1088/1751-8113/45/13/135207
13. A. V. Razumov, Yu. G. Stroganov and P. Zinn-Justin, Polynomial solutions of qKZ equation and ground state of XXZ spin chain at Δ = -1/2, J. Phys. A 40 (2007) 11827-11847, DOI 10.1088/1751-8113/40/39/009, https://doi.org/10.1088/1751-8113/40/39/009
14. Peter E. Finch, Robert Weston and Paul Zinn-Justin, Theta function solutions of the quantum Knizhnik-Zamolodchikov-Bernard equation for a face model, J. Phys. A 49 (2016) 064001, DOI 10.1088/1751-8113/49/6/064001, https://doi.org/10.1088/1751-8113/49/6/064001
15. P. Fendley and H. Saleur, N = 2 supersymmetry, Painlevé III and exact scaling functions in 2D polymers, Nucl. Phys. B 388 (1992) 609-626, DOI 10.1016/0550-3213(92)90556-Q, https://doi.org/10.1016/0550-3213(92)90556-Q
16. M. Jimbo, T. Miwa and A. Nakayashiki, Difference equations for the correlation functions of the eight-vertex model, J. Phys. A 26 (1993) 2199-2209, DOI 10.1088/0305-4470/26/9/015, https://doi.org/10.1088/0305-4470/26/9/015
17. R. J. Baxter and S. B. Kelland, Spontaneous polarization of the eight-vertex model, J. Phys. C: Solid State Phys. 7 (1974) L403-L406, DOI 10.1088/0022-3719/7/22/003, https://doi.org/10.1088/0022-3719/7/22/003
18. Klaus Fabricius, A new Q-Matrix in the eight-vertex model, J. Phys. A 40 (2007) 4075-4086, DOI 10.1088/1751-8113/40/15/002, https://doi.org/10.1088/1751-8113/40/15/002
19. Shi-shyr Roan, The Q-operator and functional relations of the eight-vertex model at root-of-unity η = 2mK/N for odd N, J. Phys. A 40 (2007) 11019-11044, DOI 10.1088/1751-8113/40/36/004, https://doi.org/10.1088/1751-8113/40/36/004
20. Rodney J. Baxter, Partition function of the Eight-Vertex lattice model, Ann. Phys. 70 (1972) 193-228, DOI 10.1016/0003-4916(72)90335-1, https://doi.org/10.1016/0003-4916(72)90335-1
21. Rodney J. Baxter, Eight-vertex model in lattice statistics and one-dimensional anisotropic Heisenberg chain. I. Some fundamental eigenvectors, Ann. Phys. 76 (1973) 1-24, DOI 10.1016/0003-4916(73)90439-9, https://doi.org/10.1016/0003-4916(73)90439-9
22. A. V. Razumov and 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, https://doi.org/10.1007/s11232-010-0078-3
23. Luigi Cantini and Andrea 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
24. Zhirong Xin, Junpeng Cao, Wen-Li Yang and Yupeng Wang, Exact physical quantities of the XYZ spin chain in the thermodynamic limit, arXiv:2411.12200 (2024), https://arxiv.org/abs/2411.12200
25. G. Niccoli and V. Terras, On correlation functions of the open XYZ spin 1/2 chain with boundary fields related by a constraint, arXiv:2507.20299 (2025), https://arxiv.org/abs/2507.20299
26. Eric Rains and Hjalmar Rosengren, Q-operators for the Ruijsenaars model, arXiv:2503.18057 (2025), https://arxiv.org/abs/2503.18057
27. Paul Fendley, Sascha Gehrmann, Eric Vernier and Frank Verstraete, XYZ integrability the easy way, arXiv:2511.04674 (2025), https://arxiv.org/abs/2511.04674
28. A. M. Povolotsky and A. A. Trofimova, Exact loop densities in O(1) dense loop model on the cylinder of odd circumference and clusters in half-turn self-dual critical percolation, arXiv:2406.15133 (2024), https://arxiv.org/abs/2406.15133