Rigorous analytic status of the theta-function Zamolodchikov-algebra constructions for the Belavin model: the $N$-dependence of the case-2 alternative algebra

Shigemoto's construction of the Belavin $\mathbb{Z}_N\times\mathbb{Z}_N$ vertex model builds candidate Zamolodchikov-algebra state vectors from theta functions with characteristics; quasi-periodicity closure leaves exactly two cases, and for the alternative case-2 construction (components $\theta_{a/N,0}(u;N\tau)$, weight normalization $\gamma=2/N$) the algebra is verified only numerically by the REDUCE computer-algebra system: it holds at $N=2$ and fails for $N\ge 3$, with no analytic explanation, while even the standard case-1 construction is supported only by quasi-periodicity matching plus an undelivered completeness claim. The problem asks for a rigorous analytic determination: prove the case-1 algebra and the completeness of its vector set for all $N$, and determine exactly the set of $N$ for which the case-2 algebra holds.

Background

Exactly solvable ("integrable") two-dimensional statistical models are lattice models whose local Boltzmann weights satisfy the Yang–Baxter equation, a functional equation that guarantees commuting transfer matrices and hence exact solvability. The Belavin model (A. A. Belavin, 1981) is the $\mathbb{Z}_N\times\mathbb{Z}_N$-symmetric vertex model — an edge-interaction model whose edge states take $N$ values and whose weights carry a pair of clock-type $\mathbb{Z}_N$ symmetries — and it provides the general elliptic solution of the Yang–Baxter equation, with weights built from theta functions with characteristics: $\theta_{a,b}(u;\tau)=\sum_{n\in\mathbb{Z}}\exp\{\pi i(n+a)^2\tau+2\pi i(n+a)(u+b)\}$. These functions are quasi-periodic rather than periodic: shifting $u$ by $1$ or by $\tau$ multiplies them by explicit elementary factors. A classical route to such solutions for vertex models, which the source itself credits to Cherednik, is the Zamolodchikov algebra: one seeks a vector $A(u)$ and a weight matrix $S(u)$ obeying $A(u_1)\otimes A(u_2+\eta)=S(u_1-u_2)\,A(u_1+\eta)\otimes A(u_2)$, where $\eta$ is the crossing parameter; if vectors satisfying this algebra form a complete (linearly independent, spanning) set, the Yang–Baxter equation for $S$ follows. Shigemoto's review (arXiv:1603.01079) pursues this route for the Belavin model with state vectors whose components are $\theta_{a/N,Bk}(u;C\tau)$, $a,k\in\mathbb{Z}_N$, where $B$ scales the second characteristic and $C$ the modulus, and with weights $S(u)=\sum_{\alpha_1,\alpha_2}w^{\alpha_1,\alpha_2}(u)\,Z^{\alpha_1}X^{\alpha_2}\otimes X^{-\alpha_2}Z^{-\alpha_1}$, where $Z$ and $X$ are the clock and shift matrices ($ZX=\omega XZ$, $\omega=e^{2\pi i/N}$) generating the "cyclic representation of SU(2)" that the paper advocates as the common origin of integrability, and $w^{\alpha_1,\alpha_2}(u)=\frac{1}{N}\,\frac{\sigma^{\alpha_1,\alpha_2}(u+\eta/N)\,\sigma^{0,0}(\gamma\eta)}{\sigma^{\alpha_1,\alpha_2}(\eta/N)\,\sigma^{0,0}(u+\gamma\eta)}$ with $\sigma^{\alpha_1,\alpha_2}(u)=\theta_{\alpha_2/N+1/2,\,\alpha_1/N+1/2}(u;\tau)$, where $\gamma$ is a normalization constant left free at this stage and fixed afterwards by quasi-periodicity matching. Requiring the family of vectors to be closed under quasi-periodic shifts $u\to u+\xi_1+\xi_2\tau$, with shift prefactors independent of $k$, leaves exactly two possibilities: case 1, $B=C=N/(N-1)$, whose state vector has components $\theta_{a/N,Nk/(N-1)}(u;N\tau/(N-1))$ and forces $\gamma=1$; and case 2, $B=C=N$, whose state vector has components $\theta_{a/N,0}(u;N\tau)$ and forces $\gamma=2/N$. The paper establishes only that both sides of the algebra have matching quasi-periodicity (which is what fixes $\gamma$ in each case) and asserts completeness without proof. For case 2 it reports a computer-algebra check (using the REDUCE system) that the "another Zamolodchikov algebra" $\tilde A(u_1)\otimes\tilde A(u_2+\eta)=\tilde S(u_1-u_2)\,\tilde A(u_1+\eta)\otimes\tilde A(u_2)$ with the $\gamma=2/N$ weight is satisfied for $N=2$ but not for $N\ge 3$, and offers no analytic explanation. The paper also connects the case-2 vectors, when $N$ is a perfect square, to the classical (Mumford) representation of the finite Heisenberg group — the group generated by the shift operator $S_b f(u)=f(u+b)$ and the quasi-periodic multiplication operator $T_a f(u)=e^{\pi i a^2\tau+2\pi i a u}f(u+a\tau)$ — on theta functions, but that identification inherits the same unproven status of the underlying algebra.

Problem Statement

Replace the quasi-periodicity heuristics and REDUCE numerics of Shigemoto's Belavin-model construction by a rigorous analytic theory, with the explicit definitions of arXiv:1603.01079 (case-1 weight $S(u)$ of Eqs. (4.2)–(4.4) and case-1 state vector of Eqs. (4.29)–(4.39) in section 4.3, case-2 state vector $\tilde A(u)$ of section 4.5 and case-2 weight $\tilde S(u)$ of Eq. (4.50), Zamolodchikov algebras Eqs. (4.22)/(4.40) and (4.53)) or an explicitly demonstrated equivalent reformulation:

1. Prove or disprove that the case-1 state vector, with components $\theta_{a/N,Nk/(N-1)}(u;N\tau/(N-1))$ and the $\gamma=1$ weight, satisfies the Zamolodchikov algebra for every integer $N\ge 2$, and prove or disprove the completeness (linear independence and spanning) of the constructed vector set that is needed to deduce the Yang–Baxter equation from the algebra — specifying explicitly the family of vectors $A_k(u)$ (the range of the label $k$) and the target space that the completeness claim concerns.

2. Determine, with proof, the set of integers $N\ge 2$ for which the case-2 state vector, with components $\theta_{a/N,0}(u;N\tau)$ and the $\gamma=2/N$ weight, satisfies the Zamolodchikov algebra: give an analytic proof of the numerically observed validity at $N=2$, and settle whether the algebra fails for every $N\ge 3$ or holds for further values of $N$.

These two parts form one objective: an exact, proof-based characterization of when the theta-function state-vector construction yields a Zamolodchikov algebra for the Belavin model, so that the integrability argument no longer rests on numerical spot checks.

References

1. Kazuyasu Shigemoto, The Elliptic Function in Statistical Integrable Models, arXiv:1603.01079 (2016; v2 2019; originally Tezukayama Academic Review Vol. 17, pp. 15–32, 2011), https://arxiv.org/abs/1603.01079
2. Kazuyasu Shigemoto, Integrability Condition in the Statistical Model and the Addition Formula of g = 2 Hyperelliptic Function, arXiv:1603.02131 (2016), https://arxiv.org/abs/1603.02131
3. Masahito Hayashi, Kazuyasu Shigemoto, Takuya Tsukioka, Common Hirota Form Bäcklund Transformation for the Unified Soliton System, J. Phys. Commun. 4 (2020) 015014, DOI 10.1088/2399-6528/ab6941, arXiv:1910.06572, https://arxiv.org/abs/1910.06572
4. Masahito Hayashi, Kazuyasu Shigemoto, Takuya Tsukioka, Elliptic Solutions for Higher Order KdV Equations, J. Phys. Commun. 4 (2020) 045013, DOI 10.1088/2399-6528/ab88df, arXiv:2003.00005, https://arxiv.org/abs/2003.00005
5. Kazuyasu Shigemoto, Full Addition Formulae of Genus Two Hyperelliptic Functions by the Duplication Method, arXiv:2010.15505 (2020), https://arxiv.org/abs/2010.15505
6. Masahito Hayashi, Kazuyasu Shigemoto, Takuya Tsukioka, The Half-period Addition Formulae for Genus Two Hyperelliptic $\wp$ Functions and the Sp(4,R) Lie Group Structure, J. Phys. Commun. 6 (2022) 085004, DOI 10.1088/2399-6528/ac8521, arXiv:2206.05097, https://arxiv.org/abs/2206.05097
7. A. A. Belavin, Dynamical symmetry of integrable quantum systems, Nuclear Physics B 180 (1981) 189–200, DOI 10.1016/0550-3213(81)90414-4, https://doi.org/10.1016/0550-3213(81)90414-4