Theta-function solution space of the level-one elliptic qKZB equation: wheel-condition classification and the first-excited-state conjecture in the three-colour model

Finch, Weston and Zinn-Justin (J. Phys. A 49 (2016) 064001) constructed the first explicit solutions of the level-one quantum Knizhnik-Zamolodchikov-Bernard (qKZB) equation for the elliptic solid-on-solid (SOS) face model: multiple contour integrals whose components are theta functions in the dynamical and spectral parameters, with prescribed pseudo-period lattice and degree matrix, satisfying a wheel condition. Two linked structural questions about these solutions remain open. First, the linear space of theta functions with the same pseudo-periods and degree that satisfy the wheel condition has never been classified in the elliptic setting, although its trigonometric analogue - Kasatani's classification via the double affine Hecke algebra and nonsymmetric Macdonald polynomials - underpins the exact ground-state and sum-rule theory of the XXZ/six-vertex combinatorial point. Second, on the combinatorial line where the SOS model becomes Baxter's three-colour model, the constructed states are conjectured, from evidence at system sizes $L=4$ and $L=6$, to be first excited states of the homogeneous three-colour Hamiltonian; no proof or refutation for general even $L$ exists. The problem asks for the classification (dimension and explicit basis of the wheel-condition theta space, including whether the constructed solutions span it) and for a determination of the spectral role of the states for general even chain length - the two being linked through the recurrence relations and wheel condition on which the eigenvector property rests.

Background

The Knizhnik-Zamolodchikov (KZ) equations are differential equations satisfied by conformal blocks in two-dimensional conformal field theory; their difference analogues, the quantum KZ (qKZ) equations, are satisfied by matrix elements of quantum integrable models. Felder's qKZB equations generalize both by using an elliptic solution of the dynamical Yang-Baxter equation - an $R$-matrix that also depends on and shifts a 'dynamical' parameter. The elliptic solid-on-solid (SOS) or face model of Andrews, Baxter and Forrester supplies the canonical such $R$-matrix: the model lives on a periodic chain of $L$ sites carrying integer heights $a_i$ with $|a_{i-1}-a_i|=1$, and its face (Boltzmann) weights are ratios of the Jacobi theta function $\vartheta_1(u|\tau)$, where $\tau$ sets the elliptic nome $p=e^{i\pi\tau}$ and $\eta$ is the crossing parameter. At level $\ell$ (a representation-theoretic parameter), the qKZB system is the pair of difference equations consisting of $R$-matrix exchange relations in the spectral parameters $u_1,\dots,u_L$ together with a cyclic condition whose shift is $s=(\ell+2)\eta$.

Finch, Weston and Zinn-Justin (FWZ) constructed, for even $L=2n$ at level $\ell=1$, two solutions $\Psi^{(2)},\Psi^{(3)}$ of this system as explicit multiple contour integrals, and proved for each component $\Psi_a$ that: it is holomorphic in all variables $(u_0=\zeta,u_1,\dots,u_L)$ and quasi-periodic under an index-2 sublattice $\Lambda^e$ of $(\pi\mathbb{Z}+\pi\tau\mathbb{Z})^{L+1}$ - hence a theta function, the elliptic analogue of a polynomial - with prescribed degree matrix $\tfrac{1}{2}\alpha_{ij}(a)$, where $\alpha_{ij}$ is computed from the height pattern $a$ ($\alpha_{ij}=L+3$ on $(i,j)=(0,0)$; $a_j-a_{j-1}$ for $i=0$, $j\geq 1$; $a_i-a_{i-1}$ for $j=0$, $i\geq 1$; $L-1$ on the diagonal $i=j\geq 1$; $-1$ off-diagonal for $i,j\geq 1$, heights extended periodically); it obeys recursion relations relating systems of size $L$ and $L-2$; and it satisfies a wheel condition - $\Psi$ vanishes whenever three spectral parameters form an arithmetic progression with spacing $\eta$, i.e. $u_j=u_i-\eta$ and $u_k=u_j-\eta$ for some $i<j<k$. The ambient theta space of fixed lattice and degree has dimension $d_n=(n+1)n^{2(n-1)}=2,12,324,20480$ for $L=2,4,6,8$, so the wheel conditions cut out an a priori far smaller subspace whose structure is unknown.

On the 'combinatorial line' $\eta=2\pi/3$ the level-1 shift is $s=3\eta=2\pi$, on which the $R$-matrix is invariant, so the qKZB solutions become eigenvectors of the SOS transfer matrix (the row-to-row evolution operator of the lattice model) with periodic boundary conditions. For $\zeta\notin(\pi/3)\mathbb{Z}$ the heights may be reduced modulo 3 and the model maps onto Baxter's three-colour model (proper 3-colourings of the square lattice, the zero-field limit of the associated vertex model); for $\zeta\in(\pi/3)\mathbb{Z}$ it degenerates to the 2-height restricted solid-on-solid (RSOS) model. In the homogeneous limit $u_1=\cdots=u_L$, the transfer matrix generates a quantum Hamiltonian $H_{3C}=(d/du)\ln t_{3C}(u)$ at $u=0$, and for $L=4$ and $L=6$ the images $\phi^{-1}\Psi^{(2)}$, $\phi^{-1}\Psi^{(3)}$ of the constructed states under the mod-3 height map $\phi$ were found computationally to be first excited states - the lowest energy level strictly above the ground state - rather than ground states as in the trigonometric precedent. FWZ conjecture this for all larger even $L$. (The two constructed solutions are distinct theta functions in the elliptic SOS model and in the three-colour problem - odd elements of the full lattice $\Lambda$ exchange $\Psi^{(2)}$ and $\Psi^{(3)}$ - and they coincide, $\Psi^{(2)}=\Psi^{(3)}$, only in the degenerate 2-height RSOS case $\zeta\in(\pi/3)\mathbb{Z}$.)

The trigonometric case explains why these questions matter. For the XXZ spin chain at anisotropy $-1/2$ (the 'combinatorial point' of the six-vertex model), the inhomogeneous ground state satisfies the qKZ system, and Kasatani classified the polynomial qKZ solutions obeying the wheel condition as subrepresentations of the polynomial representation of the double affine Hecke algebra (DAHA, the algebra whose nonsymmetric Macdonald polynomials diagonalize the commuting Cherednik operators); Kasatani and Takeyama connected the classification to nonsymmetric Macdonald polynomials, and this structure underlies the exact ground-state formulas and sum rules of the six-vertex and loop models. FWZ explicitly posed the elliptic analogue - classify the theta-function space of fixed pseudo-periods and degree satisfying the wheel condition - noting that no elliptic counterpart of the DAHA machinery is available; they also left open the vertex-operator side, since the trace construction of qKZB solutions via vertex operators (Foda-Jimbo-Miwa-Miki-Nakayashiki; Felder-Varchenko-Tarasov) is ill-defined precisely at the critical levels $\ell=1$ and $\ell\to-2$, and no vertex-IRF (interaction-round-a-face) map identifies their states with the eight-vertex eigenvectors studied on the combinatorial line by Zinn-Justin. The classification is the expected route to the spectral conjecture, because the wheel condition and the recursion relations are exactly the tools that establish the eigenvector property; as of 2024 the vertex-language counterpart remains open (Hagendorf-Rosengren note that eight-vertex/XYZ eigenvectors are not known explicitly and cite FWZ as the SOS-side state of the art).

Problem Statement

Fix an even chain length $L=2n$, an admissible periodic height pattern $a=(a_1,\dots,a_L)$ of integers with $|a_i-a_{i-1}|=1$ and $a_{i+L}=a_i$, a crossing parameter $\eta$, and an elliptic nome $p=e^{i\pi\tau}$. Work in the variables $(u_0=\zeta,u_1,\dots,u_L)$ on $\mathbb{C}^{L+1}$, let $\Lambda=(\pi\mathbb{Z}+\pi\tau\mathbb{Z})^{L+1}$, and let $\Lambda^e$ be its index-2 sublattice consisting of vectors $(\pi(m_0+n_0\tau),\dots,\pi(m_L+n_L\tau))$ whose $\tau$-multiplicities satisfy $\sum_i n_i$ even. Let $\alpha_{ij}(a)$ be the matrix defined in the source ($\alpha_{ij}=L+3$ for $i=j=0$; $a_j-a_{j-1}$ for $i=0$, $j=1,\dots,L$; $a_i-a_{i-1}$ for $j=0$, $i=1,\dots,L$; $L-1$ for $i=j$ in $1,\dots,L$; $-1$ for $i\neq j$ in $1,\dots,L$), and say that $f$ is a theta function of pseudo-period lattice $\Lambda^e$ and degree matrix $\tfrac{1}{2}\alpha_{ij}(a)$ if, for every $\lambda\in\Lambda^e$, $f(u+\lambda)=c_\lambda\exp(-i\sum_{i,j=0}^{L}n_i\alpha_{ij}(a)u_j)f(u)$ for constants $c_\lambda$ depending on $\lambda$ (the multiplier conventions of the source's Corollary 3.4). Say that $f$ satisfies the wheel condition if $f$ vanishes whenever $u_j=u_i-\eta$ and $u_k=u_j-\eta$ for some $1\leq i<j<k\leq L$.

(I) Classification. For each even $L$ and admissible height pattern $a$, determine the linear space $V_{L,a}$ of theta functions of pseudo-period lattice $\Lambda^e$ and degree matrix $\tfrac{1}{2}\alpha_{ij}(a)$ that satisfy the wheel condition: give its dimension as an explicit function of $L$ and $a$, and an explicit basis. In particular, determine how $V_{L,a}$ compares with the span of the Finch-Weston-Zinn-Justin solution components $\Psi_a^{(2)}$ and $\Psi_a^{(3)}$, which are exchanged by odd elements of $\Lambda$ and hence distinct in the generic elliptic case (in the trigonometric case the corresponding classification problem was solved by Kasatani, leading to a connection with nonsymmetric Macdonald polynomials; no elliptic counterpart of that structure is known, and the classification may proceed by any means).

(II) Spectral role. On the combinatorial line $\eta=2\pi/3$ with $\zeta\notin(\pi/3)\mathbb{Z}$, let $t_{3C}(u)$ be the transfer matrix of the three-colour model restricted to the zero sector (the subsector $H^{3C}_{L;0}$ spanned by three-colour configurations whose integer height lifts satisfy $a_{i+L}=a_i$, on which the mod-3 height map $\phi$ induces an isomorphism onto its image in the SOS space $H_L$, namely the height-shift-invariant subspace $\{v\in H_L: Tv=v\}$, where $T$ shifts all heights by three) and $H_{3C}=(d/du)\ln t_{3C}(u)$ at $u=0$ the associated homogeneous Hamiltonian. Determine the position of the constructed states in the spectrum for general even $L$: specifically, prove or refute the source's conjecture that, in the homogeneous limit $u_1=\cdots=u_L$, the images $\phi^{-1}\Psi^{(2)}$ and $\phi^{-1}\Psi^{(3)}$ of the level-1 qKZB solutions under the mod-3 height map are first excited states of $H_{3C}$ - eigenstates whose energy is the smallest eigenvalue strictly greater than the ground-state energy - for all even $L$ (the cases $L=4$ and $L=6$ being already verified by the authors, so any counterexample would occur at an even $L\geq 8$).

These two deliverables form one objective - completing the structural theory of the level-one theta-function solutions of the elliptic qKZB equation - and are scientifically linked: the wheel condition and the size-recursion relations are precisely the tools that establish the transfer-matrix eigenvector property, and the classification of the ambient wheel-condition space (its dimension and degeneracy structure) is the natural route to the level-placement statement, as the trigonometric classification was for the six-vertex sum-rule program.

References

1. P. E. Finch, R. Weston, P. Zinn-Justin, Theta function solutions of the quantum Knizhnik-Zamolodchikov-Bernard equation for a face model, J. Phys. A: Math. Theor. 49 (2016) 064001, doi:10.1088/1751-8113/49/6/064001, https://doi.org/10.1088/1751-8113/49/6/064001 (preprint arXiv:1509.04594, https://arxiv.org/abs/1509.04594)
2. G. Felder, Conformal field theory and integrable systems associated to elliptic curves, in: Proceedings of the International Congress of Mathematicians (Zurich, 1994), Birkhauser, 1995, pp. 1247-1255, preprint arXiv:hep-th/9407154, https://arxiv.org/abs/hep-th/9407154
3. G. Felder, V. Tarasov, A. Varchenko, Solutions of the elliptic qKZB equations and Bethe ansatz I, in: Topics in Singularity Theory, Amer. Math. Soc. Transl. Ser. 2 180 (1997) 45-75, doi:10.1090/trans2/180/04, https://doi.org/10.1090/trans2/180/04 (preprint arXiv:q-alg/9606005, https://arxiv.org/abs/q-alg/9606005)
4. O. Foda, M. Jimbo, T. Miwa, K. Miki, A. Nakayashiki, Vertex operators in solvable lattice models, J. Math. Phys. 35 (1994) 13-46, doi:10.1063/1.530783, https://doi.org/10.1063/1.530783 (preprint arXiv:hep-th/9305100, https://arxiv.org/abs/hep-th/9305100)
5. G. Andrews, R. Baxter, P. Forrester, Eight-vertex SOS model and generalized Rogers-Ramanujan-type identities, J. Statist. Phys. 35 (1984) 193-266, doi:10.1007/BF01014383, https://doi.org/10.1007/BF01014383
6. R. J. Baxter, Three-colorings of the square lattice: A hard squares model, J. Math. Phys. 11 (1970) 3116-3124, doi:10.1063/1.1665102, https://doi.org/10.1063/1.1665102
7. M. Kasatani, Subrepresentations in the polynomial representation of the double affine Hecke algebra of type GL_n at t^{k+1}q^{r-1}=1, Int. Math. Res. Not. 2005 (2005), no. 28, 1717-1742, doi:10.1155/IMRN.2005.1717, https://doi.org/10.1155/IMRN.2005.1717 (preprint arXiv:math/0501272, https://arxiv.org/abs/math/0501272)
8. M. Kasatani, Y. Takeyama, The quantum Knizhnik-Zamolodchikov equation and non-symmetric Macdonald polynomials, Funkcial. Ekvac. 50 (2007) 491-509, doi:10.1619/fesi.50.491, https://doi.org/10.1619/fesi.50.491 (preprint arXiv:math/0608773, https://arxiv.org/abs/math/0608773)
9. 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 (2013) 599-637, doi:10.1007/978-1-4471-4863-0_26, https://doi.org/10.1007/978-1-4471-4863-0_26 (preprint arXiv:1202.4420, https://arxiv.org/abs/1202.4420)
10. C. Hagendorf, H. Rosengren, Nearest-neighbour correlation functions for the supersymmetric XYZ spin chain and Painleve VI, Commun. Math. Phys. 405 (2024), article 97, doi:10.1007/s00220-024-04977-w, https://doi.org/10.1007/s00220-024-04977-w (preprint arXiv:2208.13533, https://arxiv.org/abs/2208.13533)
11. R. Klabbers, J. Lamers, Modular families of elliptic long-range spin chains from freezing, Lett. Math. Phys. 116 (2026), no. 3, article 58, doi:10.1007/s11005-026-02075-7, https://doi.org/10.1007/s11005-026-02075-7 (preprint arXiv:2507.13104, https://arxiv.org/abs/2507.13104)
12. H. Rosengren, Elliptic pfaffians and solvable lattice models, J. Stat. Mech. (2016) 083106, doi:10.1088/1742-5468/2016/08/083106, https://doi.org/10.1088/1742-5468/2016/08/083106 (preprint arXiv:1605.02915, https://arxiv.org/abs/1605.02915)
13. H. Konno, Elliptic weight functions and elliptic q-KZ equation, preprint arXiv:1706.07630, https://arxiv.org/abs/1706.07630
14. H. De Clercq, N. Reshetikhin, J. Stokman, Graphical calculus for quantum vertex operators, II: q-KZB and coordinate Macdonald-Ruijsenaars equations, preprint arXiv:2411.17306, https://arxiv.org/abs/2411.17306