Algebraic completion of the Q-operator programme for the XXZ chain at roots of unity: FM-string creation operators, the hidden Onsager-type symmetry, Bethe completeness, and FM-string thermodynamics

The spin-1/2 Heisenberg XXZ chain at root-of-unity anisotropy $\eta=\mathrm{i}\pi\ell_1/\ell_2$ is exactly solvable, yet its spectrum there is only partially understood: six-vertex transfer-matrix eigenvalues come in exponentially large degenerate blocks organized by Fabricius–McCoy (FM) strings and 'descendant towers'. Miao, Lamers and Pasquier (SciPost Phys. 11, 067 (2021)) constructed the truncated Baxter Q operator at roots of unity, proved the associated truncated fusion and Wronskian relations, and left the structural core as conjectures: derivatives of semicyclic transfer matrices create and annihilate FM strings (Conjectures 1–4, supported by numerics only), the imaginary parts of FM-string centres take fixed values at commensurate twist (Conjecture 5), and an unidentified algebraic symmetry 'compatible with the Q operator' should make primitive ('highest-weight') and descendant states intrinsic, with rigorous Bethe/Q completeness and the thermodynamic-Bethe-ansatz role of FM strings left open. Rational Q-systems, non-commuting transfer matrices and affine Temperley–Lieb arguments have addressed parts of this picture, but the creation-operator conjectures, the Q-operator-compatible symmetry algebra, a rigorous completeness proof, and the derivation of FM-string thermodynamics all remain open. The problem asks for this missing algebraic and thermodynamic completion.

Background

The spin-1/2 Heisenberg XXZ chain is the paradigmatic quantum integrable model: a one-dimensional lattice of spin-1/2 degrees of freedom with Hamiltonian

$$H=\sum_{j=1}^{N}\left(\sigma^x_j\sigma^x_{j+1}+\sigma^y_j\sigma^y_{j+1}+\Delta\,\sigma^z_j\sigma^z_{j+1}\right),\qquad\Delta=\cosh\eta,$$

with twisted boundary conditions, where $\sigma^{x,y,z}_j$ are Pauli matrices. It is solved by the Bethe ansatz: eigenstates are parametrized by Bethe roots $u_m$ satisfying algebraic Bethe equations, and the conserved charges are generated by the commuting family of six-vertex transfer matrices $T_s(u)$ (row-to-row transfer matrices of the six-vertex statistical model, with spectral parameter $u$ and auxiliary spin $s$); the Hamiltonian is obtained from the fundamental transfer matrix by a logarithmic derivative, up to rescaling and an additive constant. Baxter's Q operator is an operator commuting with all $T_s(u)$ whose eigenvalues are polynomials whose zeros are exactly the Bethe roots; it underlies the functional $TQ$ relation (the identity linking the eigenvalues of $T$ and $Q$ that encodes the Bethe equations) and is the deepest layer of the solution.

When $\eta$ is a root of unity, $\eta=\mathrm{i}\pi\,\ell_1/\ell_2$ with $\ell_1,\ell_2$ coprime (equivalently $q=e^{\eta}$ with $q^{2\ell_2}=1$), the spectrum develops a singular structure that is still not fully understood. Transfer-matrix eigenvalues then appear in exponentially large degenerate blocks. The degenerate partners form 'descendant towers': a lowest member, called the primitive state, together with all states obtained from it by adjoining Fabricius–McCoy (FM) strings — bound states of exactly $\ell_2$ Bethe roots $u_0+\mathrm{i}\pi(n-1)/\ell_2$, $n=1,\ldots,\ell_2$, with zero string deviation, which scatter trivially off all other roots and carry no local charge generated by the fundamental transfer matrix $T_{1/2}(u)$ (apart from a possible momentum $\pi$), although they do carry the quasilocal Z charges. Two symmetry structures are known to be involved, and both fall short. The $\mathfrak{sl}_2$ loop algebra (the Lie algebra of polynomial maps into $\mathfrak{sl}_2$) commutes with the transfer matrices but acts only on part of the Hilbert space and does not commute with the Q operator, so its highest-weight and descendant notions do not reproduce the observed towers. The quasilocal Z charges — conserved charges that are sums of local densities with exponentially decaying overlap (quasilocality), existing only at roots of unity, which explain the non-vanishing high-temperature spin Drude weight (the zero-frequency ballistic part of the spin-transport conductivity) and are generated by transfer matrices built from semicyclic (non-highest-weight, cyclic up to one matrix entry) auxiliary representations — control the thermodynamics through string-charge duality, the bijection between Bethe-string densities and conserved-charge eigenvalues that underlies generalized hydrodynamics (the exact macroscopic theory of integrable systems built on their infinitely many conserved charges); but their operator relation to the degeneracies was unclear.

Miao, Lamers and Pasquier (2021) constructed the Q operator at arbitrary anisotropy from a two-parameter transfer matrix with a complex-spin auxiliary space, and showed that at roots of unity a truncation makes the Q operator an explicit finite-dimensional trace over the auxiliary space. This yielded proofs of the truncated fusion and Wronskian relations (functional relations between the Q and P eigenvalue polynomials), a proof of a previously conjectured interpolation formula, the quantization condition identifying FM-string centres with zeros of the product of the regular parts of the Q and P eigenvalues (the Drinfeld polynomial of the earlier loop-algebra literature, i.e. the polynomial characterising the highest-weight representations of that loop algebra), and the picture of FM strings behaving like free fermions within each descendant tower. The operator content nevertheless remained conjectural. Conjectures 1–4 assert that the linearization of the semicyclic transfer matrix in its cyclic parameter at auxiliary spin $s=(\ell_2-1)/2$,

$$\mathbf{B}^{\mathrm{FM}}(u)=\left.\partial_\beta\,\mathbf{T}^{\mathrm{sc}}_s(u,\phi,\beta)\right|_{\beta=0},$$

creates an FM string on any eigenstate of the Q operator, with an analogous operator $\mathbf{C}^{\mathrm{FM}}(u)$ annihilating one; these were verified only in numerical examples, and the commutation relations of these operators away from eigenstates ('off shell') are unknown. (An earlier creation operator proposed by Fabricius and McCoy yields, in the examples checked, eigenvectors of $H$ and of the fundamental transfer matrix, but not of the Q operator.) Conjecture 5 asserts that at commensurate twist — when the boundary phase is tuned so that FM strings solve the Bethe equations — FM-string centres have imaginary part $\pi/(2\ell_2)$ for odd $\ell_1$ and $0$ for even $\ell_1$. The authors singled out the central open question: a deeper algebraic structure compatible with the Q operator, which would make the terms 'highest-weight' and 'descendant' intrinsic — conjecturally an Onsager algebra, the infinite-dimensional Lie algebra with generators $A_m$, $G_l$ obeying $[A_l,A_m]=4G_{l-m}$, $[G_l,A_m]=2A_{l+m}-2A_{m-l}$, $[G_l,G_m]=0$, introduced in Onsager's solution of the two-dimensional Ising model — together with a rigorous proof of completeness of the Bethe/Q description at roots of unity and a derivation of the FM-string contributions within the thermodynamic Bethe ansatz (TBA), the finite-temperature formalism that encodes thermodynamics in the densities of Bethe strings in Takahashi's classification, whose 'last two' string types are special at roots of unity.

Since 2021 the finite-size picture has advanced considerably — an algorithmic determination of all physical solutions and exact FM-string centres (Hou–Jiang–Miao), explicit Onsager generators from a quadratic algebra of non-commuting transfer matrices (Vernier–Miao–Yamazaki), and rigorous degeneracy bounds from the representation theory of the affine Temperley–Lieb algebra (the diagrammatic algebra underlying the XXZ chain; Hu–Gerken–Posske) — yet none of these proves the conjectured creation/annihilation operators, identifies the algebra compatible with the Q operator, proves completeness rigorously, or derives the FM-string thermodynamics. That missing algebraic and thermodynamic completion is the open core of the problem.

Problem Statement

Objective: replace the numerically supported, partly conjectural picture of the root-of-unity spectrum of the periodic spin-1/2 XXZ chain by a proven algebraic theory, and establish its thermodynamic consequences. The setting is the twisted periodic spin-1/2 XXZ chain at anisotropy $\eta=\mathrm{i}\pi\,\ell_1/\ell_2$ (with $\ell_1,\ell_2$ coprime; commensurate twist means the boundary phase $\phi$ is tuned so that FM strings solve the Bethe equations), where the six-vertex transfer matrix and the truncated Baxter Q operator of Miao, Lamers and Pasquier exhibit exponentially degenerate spectra organized by Fabricius–McCoy (FM) strings and descendant towers. Four linked deliverables — all understood for all roots of unity, in both sign sectors $\varepsilon=q^{\ell_2}=\pm1$, and for all commensurate twists — jointly constitute the resolution.

(1) Operator algebra. Prove Conjectures 1–4 of the source paper: the semicyclic-transfer-matrix derivatives

$$\mathbf{B}^{\mathrm{FM}}(u)=\left.\partial_\beta\,\mathbf{T}^{\mathrm{sc}}_s(u,\phi,\beta)\right|_{\beta=0,\;2s=\ell_2-1}$$

and their transposed-Lax twins $\mathbf{C}^{\mathrm{FM}}(u)=\left.\partial_\gamma\,\widehat{\mathbf{T}}^{\mathrm{sc}}_s(u,\phi,\gamma)\right|_{\gamma=0,\;2s=\ell_2-1}$ (the semicyclic transfer matrix built from the transposed Lax operator, whose cyclic parameter $\gamma$ deforms the auxiliary matrix entry as $\mathbf{S}^-\mapsto\mathbf{S}^-+\gamma\,|\ell_2-1\rangle\langle 0|$) create and annihilate FM strings $\{u+\delta_\varepsilon,\,u+\delta_\varepsilon+\mathrm{i}\pi/\ell_2,\,\ldots,\,u+\delta_\varepsilon+\mathrm{i}\pi(\ell_2-1)/\ell_2\}$ on eigenstates of the truncated Q operator — where $u$ is the operator's spectral parameter and, following the convention of Conjectures 3–4, $\delta_\varepsilon=0$ for $\varepsilon=+1$ while $\delta_\varepsilon=\mathrm{i}\pi/(2\ell_2)$ for $\varepsilon=-1$ — and determine the off-shell (non-eigenstate) commutation relations of these operators.

(2) Intrinsic symmetry. Identify the deeper algebraic structure compatible with the Q operator — the structure whose representation theory organizes the degenerate spectrum so that the labels 'primitive' (highest-weight) and 'descendant' become intrinsic — and establish its relation to the Onsager algebra and to the semicyclic transfer matrices.

(3) Completeness. Building on (1)–(2) or independently, prove completeness of the Bethe/Q-operator description at roots of unity: a rigorous classification of all eigenstates of the twisted periodic chain into primitive states and descendant towers — including degenerate primitive states, Bethe roots sent to $u=\pm\infty$ (which drop out of the Bethe equations), and all parameter-dependent exceptional cases — whose counting matches the Hilbert-space dimension in every magnetization sector.

(4) Thermodynamics. Derive the FM-string contributions within the thermodynamic Bethe ansatz, their precise relation to the quasilocal Z charges and string-charge duality for arbitrary roots of unity and commensurate twists, and prove that the finite-size FM-string picture yields the correct thermodynamic root densities (including the last two Takahashi string types).

These deliverables form one scientific objective: (1) supplies the operator algebra, (2) its representation-theoretic meaning, (3) the resulting finite-size spectral classification, and (4) its thermodynamic limit. The route is not prescribed: any rigorous construction with the stated proved content — for instance via cyclic-representation Q operators or the non-commuting transfer-matrix algebra — is admissible, provided it covers the stated generality rather than special points or fixed system sizes.

References

1. Y. Miao, J. Lamers, V. Pasquier, 'On the Q operator and the spectrum of the XXZ model at root of unity', SciPost Phys. 11, 067 (2021). DOI: 10.21468/SciPostPhys.11.3.067; arXiv:2012.10224; https://scipost.org/10.21468/SciPostPhys.11.3.067
2. Y. Miao, 'Conjectures on hidden Onsager algebra symmetries in interacting quantum lattice models', SciPost Phys. 11, 066 (2021). DOI: 10.21468/SciPostPhys.11.3.066; arXiv:2103.14569; https://scipost.org/10.21468/SciPostPhys.11.3.066
3. J. Hou, Y. Jiang, Y. Miao, 'Rational Q-systems at root of unity I. Closed chains', SciPost Phys. 16, 129 (2024). DOI: 10.21468/SciPostPhys.16.5.129; arXiv:2310.14966; https://scipost.org/10.21468/SciPostPhys.16.5.129
4. E. Vernier, Y. Miao, M. Yamazaki, 'Lattice non-invertible symmetry from non-commuting transfer matrices', arXiv:2606.25660 (2026); https://arxiv.org/abs/2606.25660
5. R. Weston, 'Cyclic representations of $U_q(\widehat{\mathfrak{sl}}_2)$ and its Borel subalgebras at roots of unity and Q-operators', SciPost Phys. Core 9, 016 (2026). DOI: 10.21468/SciPostPhysCore.9.1.016; arXiv:2412.14811; https://scipost.org/10.21468/SciPostPhysCore.9.1.016
6. Y. Hu, F. Gerken, T. Posske, 'Hidden Twisted Sectors and Exponential Degeneracy in Root-of-Unity XXZ Heisenberg Chains', arXiv:2602.15098 (2026); https://arxiv.org/abs/2602.15098
7. Y. Miao, 'Generalised Onsager Algebra in Quantum Lattice Models', SciPost Phys. 13, 070 (2022). DOI: 10.21468/SciPostPhys.13.3.070; arXiv:2203.16594; https://scipost.org/10.21468/SciPostPhys.13.3.070
8. Y. Jiang, Y. Miao, 'Spectrum-preserving deformations of integrable spin chains with open boundaries', Phys. Rev. B 110, L121406 (2024). DOI: 10.1103/PhysRevB.110.L121406; arXiv:2406.17171; https://doi.org/10.1103/PhysRevB.110.L121406
9. Y. Jiang, Y.-C. Liu, Y. Miao, Z.-X. Tan, 'Rational Q-systems for integrable spin chains without U(1) symmetry', SciPost Phys. 20, 175 (2026). DOI: 10.21468/SciPostPhys.20.6.175; arXiv:2512.01551; https://scipost.org/10.21468/SciPostPhys.20.6.175
10. J. Hou, Y. Jiang, R.-D. Zhu, 'Spin-s rational Q-system', SciPost Phys. 16, 113 (2024). DOI: 10.21468/SciPostPhys.16.4.113; arXiv:2303.07640; https://scipost.org/10.21468/SciPostPhys.16.4.113
11. K. Fabricius, B. M. McCoy, 'Evaluation Parameters and Bethe Roots for the Six-Vertex Model at Roots of Unity', Birkhäuser Boston, pp. 119–144 (2002). DOI: 10.1007/978-1-4612-0087-1_6; arXiv:cond-mat/0108057; https://doi.org/10.1007/978-1-4612-0087-1_6
12. C. Korff, 'Auxiliary matrices for the six-vertex model at $q^N=1$: II. Bethe roots, complete strings and the Drinfeld polynomial', J. Phys. A 37, 385 (2004). DOI: 10.1088/0305-4470/37/2/009; arXiv:math-ph/0305035; https://doi.org/10.1088/0305-4470/37/2/009
13. T. Deguchi, K. Fabricius, B. M. McCoy, 'The $\mathfrak{sl}_2$ Loop Algebra Symmetry of the Six-Vertex Model at Roots of Unity', J. Stat. Phys. 102, 701 (2001). DOI: 10.1023/A:1004894701900; arXiv:cond-mat/9912141; https://doi.org/10.1023/A:1004894701900
14. T. Deguchi, 'Regular XXZ Bethe states at roots of unity as highest weight vectors of the $\mathfrak{sl}_2$ loop algebra', J. Phys. A 40, 7473 (2007). DOI: 10.1088/1751-8113/40/27/005; arXiv:cond-mat/0503564; https://doi.org/10.1088/1751-8113/40/27/005
15. T. Prosen, E. Ilievski, 'Families of Quasilocal Conservation Laws and Quantum Spin Transport', Phys. Rev. Lett. 111, 057203 (2013). DOI: 10.1103/PhysRevLett.111.057203; arXiv:1306.4498; https://doi.org/10.1103/PhysRevLett.111.057203
16. L. Zadnik, M. Medenjak, T. Prosen, 'Quasilocal conservation laws from semicyclic irreducible representations of $U_q(\mathfrak{sl}_2)$ in XXZ spin-1/2 chains', Nucl. Phys. B 902, 339 (2016). DOI: 10.1016/j.nuclphysb.2015.11.023; arXiv:1510.08302; https://doi.org/10.1016/j.nuclphysb.2015.11.023
17. E. Ilievski, E. Quinn, J. De Nardis, M. Brockmann, 'String-charge duality in integrable lattice models', J. Stat. Mech. (2016) 063101. DOI: 10.1088/1742-5468/2016/06/063101; arXiv:1512.04454; https://doi.org/10.1088/1742-5468/2016/06/063101
18. E. Vernier, E. O'Brien, P. Fendley, 'Onsager symmetries in U(1)-invariant clock models', J. Stat. Mech. (2019) 043107. DOI: 10.1088/1742-5468/ab11c0; arXiv:1812.09091; https://doi.org/10.1088/1742-5468/ab11c0
19. M. Takahashi, 'Thermodynamics of One-Dimensional Solvable Models', Cambridge University Press (1999). DOI: 10.1017/CBO9780511524332; https://doi.org/10.1017/CBO9780511524332
20. C. Rylands, E. Vernier, P. Calabrese, 'Dynamical symmetry restoration in the Heisenberg spin chain', J. Stat. Mech. (2024) 123102. DOI: 10.1088/1742-5468/ad97b3; arXiv:2409.08735; https://doi.org/10.1088/1742-5468/ad97b3
21. V. Popkov, X. Zhang, A. Klümper, 'Phantom Bethe excitations and spin helix eigenstates in integrable periodic and open spin chains', Phys. Rev. B 104, L081410 (2021). DOI: 10.1103/PhysRevB.104.L081410; arXiv:2102.03295; https://doi.org/10.1103/PhysRevB.104.L081410