The elliptic Calogero-Inozemtsev model of BCn\mathrm{BC}_n type - nn particles on a torus interacting through Weierstrass ℘\wp-potentials at qi−qjq_i-q_j and mirror images qi+qjq_i+q_j, with one-body reflection terms at the half-periods and five complex couplings - admits an RR-matrix-valued generalization of Takasaki's 2n×2n2n\times 2n Lax pair built from the elliptic Shibukawa-Ueno RR-operator and Komori-Hikami KK-operators satisfying an associative reflection equation. Freezing this Lax pair at a classical equilibrium yields a proposed quantum long-range spin chain of BCn\mathrm{BC}_n type with fully anisotropic (XYZ, eight-vertex) interactions on spin-1/2 sites. The open problem posed in the source concerns this frozen chain: its integrability - a complete commuting family of conserved operators, commuting transfer matrices, or a Bethe-type solution - is conjectured but unproven, the quantum Lax relation being structural evidence rather than proof. The equilibrium positions ζi\zeta_i at which the freezing is performed solve nn transcendental equilibrium equations and are unknown for arbitrary couplings (only an implicit characterization via zeros of Heun-equation solutions exists), so the chain cannot be explicitly instantiated in general. The problem is to prove or refute integrability of the frozen BCn\mathrm{BC}_n elliptic XYZ chain.

Background

Integrable many-body systems of Calogero-Moser type describe nn particles on a line or circle whose pairwise interaction is an inverse-square potential. In the elliptic version the positions qiq_i live on the torus C/(Z+τZ)\mathbb{C}/(\mathbb{Z}+\tau\mathbb{Z}), Im τ>0\mathrm{Im}\,\tau>0, and the interaction is the Weierstrass ℘\wp-function, the doubly periodic meromorphic function with a double pole ℘(z)∼z−2\wp(z)\sim z^{-2} at the origin. The elliptic Calogero-Moser model corresponds to the root system An−1A_{n-1}. Its BCn\mathrm{BC}_n analogue, introduced by Inozemtsev, keeps the pair potential ℘(qi−qj)\wp(q_i-q_j) and adds mirror terms ℘(qi+qj)\wp(q_i+q_j) together with one-body terms ℘(qi+ωa)\wp(q_i+\omega_a) at the half-periods ω0=0\omega_0=0, ω1=1/2\omega_1=1/2, ω2=(1+τ)/2\omega_2=(1+\tau)/2, ω3=τ/2\omega_3=\tau/2:

H=12∑k=1npk2−g2∑i<jn(℘(qi−qj)+℘(qi+qj))−12∑a=03∑k=1nνa2 ℘(qk+ωa),H=\frac{1}{2}\sum_{k=1}^{n}p_{k}^{2}-g^{2}\sum_{i<j}^{n}\Big(\wp(q_{i}-q_{j})+\wp(q_{i}+q_{j})\Big)-\frac{1}{2}\sum_{a=0}^{3}\sum_{k=1}^{n}\nu_{a}^{2}\,\wp(q_{k}+\omega_{a}),

with five arbitrary complex couplings g,ν0,ν1,ν2,ν3g,\nu_0,\nu_1,\nu_2,\nu_3. Like all Calogero-Moser systems it is integrable and admits a Lax representation: the equations of motion are equivalent to a matrix equation L˙(z)=[L(z),M(z)]\dot L(z)=[L(z),M(z)] with spectral parameter zz, so that the invariants of L(z)L(z), in particular the coefficients of the spectral curve, are conserved. For the Inozemtsev model a 2n×2n2n\times 2n Lax pair is due to Takasaki.

A newer development replaces the scalar entries of LL by RR-matrices. An RR-matrix is a solution of the quantum Yang-Baxter equation acting on a tensor product; the relevant elliptic one is the Baxter-Belavin RR-matrix, or, in an infinite-dimensional form acting on functions of nn variables, the completely Z\mathbb{Z}-symmetric Shibukawa-Ueno RR-operator. The key algebraic input is the associative Yang-Baxter equation (AYBE), a cubic relation for an RR-matrix depending on two pairs of parameters (a spectral parameter and a Planck constant) which makes products of RR-operators associative. R-matrix-valued Lax pairs built on it, with entries given by RR-operators evaluated at particle differences, were introduced by Levin-Olshanetsky-Zotov and developed by Grekov-Zotov and Sechin-Zotov, who derived from them, by a freezing trick, anisotropic (XYZ-type) long-range spin chains of AA-type.

The freezing trick, due to Polychronakos, turns a quantum spin many-body Hamiltonian into a quantum spin chain by evaluating it at a classical equilibrium of the spinless classical model: at pi=0p_i=0, qi=ζiq_i=\zeta_i with ζ\zeta an equilibrium, the (quantum) Lax equation collapses to a relation among the restricted matrices, and the restricted Hamiltonian becomes the Hamiltonian of a spin chain with long-range, all-to-all interactions at the fixed sites ζi\zeta_i. The Haldane-Shastry chain, the paradigmatic integrable long-range spin chain with inverse-square interactions, is the rational/trigonometric ancestor of this family. Applied to the elliptic Calogero-Moser model, freezing gives Inozemtsev's elliptic spin chain; applied to RR-matrix-valued AA-type Lax pairs it gives the anisotropic Sechin-Zotov chain; the elliptic qq-deformed Haldane-Shastry-type chain of Matushko-Zotov arises likewise from elliptic spin Ruijsenaars-Macdonald operators, whose commutativity was proven from RR-matrix identities.

The source of this problem (Matushko-Mostovskii-Zotov) extends the whole construction to BCn\mathrm{BC}_n. It extends Kirillov's BB-type associative Yang-Baxter equations to relations depending on spectral parameters and Planck constants, builds an RR-matrix-valued version of Takasaki's 2n×2n2n\times 2n Lax pair from the Shibukawa-Ueno RR-operator together with elliptic KK-operators of Komori-Hikami satisfying a reflection equation (Sklyanin's boundary formalism augments an RR-matrix with a boundary matrix KK so that integrable models with boundaries can be built), and, via the Felder-Pasquier construction, rewrites the answer through Baxter's eight-vertex RR-matrix, giving an explicit 2×22\times 2 (spin-1/2) realization. Freezing at an equilibrium then proposes a BCn\mathrm{BC}_n elliptic XYZ long-range spin chain: a Hamiltonian on (C2)⊗n(\mathbb{C}^2)^{\otimes n} with two-body anisotropic couplings σ4−ak⊗σ4−al\sigma^k_{4-a}\otimes\sigma^l_{4-a} (a=1,2,3a=1,2,3; σ0\sigma_0 the identity, σ1,σ2,σ3\sigma_1,\sigma_2,\sigma_3 the Pauli matrices) built from the Kronecker elliptic function and the Eisenstein-Kronecker series E1,E2E_1,E_2 evaluated at ζk∓ζl\zeta_k\mp\zeta_l and half-period shifts, plus one-site anisotropic terms with coefficients νa\nu_a.

Two facts delimit the open problem. First, the equilibrium positions needed to instantiate the frozen chain are not known in general: for the AA-type elliptic model the equally spaced points ζj=j/n\zeta_j=j/n are equilibria by symmetry, but for BCn\mathrm{BC}_n with five couplings the equilibrium equations form a transcendental system whose general solution is unknown; the classical theory of Treibich-Verdier potentials (Gesztesy-Weikard) identifies the ζi\zeta_i with zeros of solutions of the Heun equation in elliptic form, an implicit characterization only. Second, integrability of the frozen chain is conjectural: the frozen Lax matrices satisfy a quantum Lax equation [H0′,L′(z)]=[L′(z),M′(z)][\mathcal{H}^{0\prime},\mathcal{L}^{\prime}(z)]=[\mathcal{L}^{\prime}(z),\mathcal{M}^{\prime}(z)], which is structural evidence, but, as the authors emphasize, the existence of a Lax equation is not sufficient for integrability. For comparison, the AA-type RR-matrix chain was recently proven integrable by Chalykh-Matushko using RR-matrix Dunkl operators (Dunkl operators are difference-differential operators that permute variables and encode root-system symmetries; suitable symmetric combinations of them commute), and the isotropic Inozemtsev chain was proven integrable by Chalykh via elliptic Dunkl operators; no analogous construction exists so far for the BCn\mathrm{BC}_n RR-matrix-valued case with reflection.

Problem Statement

The problem, as posed in the concluding section of Matushko-Mostovskii-Zotov (arXiv:2503.22659, published in Letters in Mathematical Physics, 2026), is to establish the proposed BCn\mathrm{BC}_n elliptic XYZ long-range spin chain as a quantum integrable model.

The chain is obtained by freezing the RR-matrix-valued Takasaki Lax pair of the BCn\mathrm{BC}_n Calogero-Inozemtsev system at a classical equilibrium position. For the classical Hamiltonian HH displayed in the Background, on a fixed elliptic curve and with five arbitrary complex couplings g,ν0,ν1,ν2,ν3g,\nu_0,\nu_1,\nu_2,\nu_3, the sites ζ1,…,ζn\zeta_1,\dots,\zeta_n are the positions of an equilibrium (with pi=0p_i=0) solving

g2∑k=1k≠in(℘′(ζi−ζk)+℘′(ζi+ζk))+12∑a=03νa2 ℘′(ζi+ωa)=0,i=1,…,n.g^{2}\sum_{\substack{k=1\\ k\neq i}}^{n}\Big(\wp'(\zeta_{i}-\zeta_{k})+\wp'(\zeta_{i}+\zeta_{k})\Big)+\frac{1}{2}\sum_{a=0}^{3}\nu_{a}^{2}\,\wp'(\zeta_{i}+\omega_{a})=0,\qquad i=1,\dots,n.

As the source notes, these equilibrium positions are unknown in the general case to the best of the authors' knowledge; it is known that the ζi\zeta_i coincide with positions of zeros of solutions to the Heun equation in elliptic form, an implicit characterization.

The problem is to prove or refute the conjecture that the quantum spin chain obtained by freezing the RR-matrix-valued Takasaki Lax pair at such an equilibrium position is integrable. Concretely, the chain acts on (C2)⊗n(\mathbb{C}^2)^{\otimes n} with Hamiltonian H0′\mathcal{H}^{0\prime} equal to eq. (7.5) of the source: two-body terms −g2∑k<l(E2(ζk−ζl)+E2(ζk+ζl))σ0k⊗σ0l-\frac{g}{2}\sum_{k<l}(E_2(\zeta_k-\zeta_l)+E_2(\zeta_k+\zeta_l))\sigma^k_0\otimes\sigma^l_0 and anisotropic terms g2∑k<l∑a=13[φa(ζk−ζl,ωa)(E1(ζk−ζl+ωa)−E1(ζk−ζl)−E1(ωa))+(ζk−ζl→ζk+ζl)]σ4−ak⊗σ4−al\frac{g}{2}\sum_{k<l}\sum_{a=1}^{3}\big[\varphi_a(\zeta_k-\zeta_l,\omega_a)(E_1(\zeta_k-\zeta_l+\omega_a)-E_1(\zeta_k-\zeta_l)-E_1(\omega_a))+( \zeta_k-\zeta_l\to\zeta_k+\zeta_l)\big]\sigma^k_{4-a}\otimes\sigma^l_{4-a}, plus one-site terms −12∑k=1n[ν0E2(ζk)σ0k+∑a=13νaφa(ζk−ωa,ωa)(E1(ζk+ωa)−E1(ζk)−E1(ωa))σ4−ak]-\frac{1}{2}\sum_{k=1}^{n}\big[\nu_0 E_2(\zeta_k)\sigma^k_0+\sum_{a=1}^{3}\nu_a\varphi_a(\zeta_k-\omega_a,\omega_a)(E_1(\zeta_k+\omega_a)-E_1(\zeta_k)-E_1(\omega_a))\sigma^k_{4-a}\big], where E1,E2E_1,E_2 are the Eisenstein-Kronecker series and φa\varphi_a the half-period-shifted Kronecker function, in the conventions of the source. Integrability means the existence of a complete set of mutually commuting, independent conserved operators containing this Hamiltonian; mechanisms through which such a family might be constructed include commuting transfer matrices for the chain with boundaries, a Bethe-ansatz solution of its spectrum, or a Dunkl/Cherednik-type operator construction whose freezing yields the commuting charges. The conjecture is for arbitrary nn and arbitrary five couplings, with the chain defined, as in the source, at some equilibrium position satisfying the equilibrium equations above; the quantum Lax equation [H0′,L′(z)]=[L′(z),M′(z)][\mathcal{H}^{0\prime},\mathcal{L}^{\prime}(z)]=[\mathcal{L}^{\prime}(z),\mathcal{M}^{\prime}(z)] satisfied by the frozen Lax matrices is evidence, not proof, and settling the conjecture requires the commutativity (or its failure) to be established.

Because the equilibrium positions are unknown for general couplings, the frozen Hamiltonian cannot be concretely instantiated in general, and without a commuting family the chain remains a proposal supported only by a quantum Lax relation; establishing, or refuting, its integrability is the open question posed by the source.

References

  1. M. Matushko, A. Mostovskii, A. Zotov, R-matrix valued Lax pair for elliptic Calogero-Inozemtsev system and associative Yang-Baxter equations of BC_n type, Letters in Mathematical Physics (2026), DOI: 10.1007/s11005-025-02043-7, arXiv:2503.22659, https://arxiv.org/abs/2503.22659
  2. O. Chalykh, M. Matushko, R-matrix Dunkl Operators and Spin Calogero-Moser System, Annales Henri Poincare (2026), DOI: 10.1007/s00023-026-01671-6, arXiv:2509.18989, https://arxiv.org/abs/2509.18989
  3. O. Chalykh, Integrability of the Inozemtsev spin chain, arXiv:2407.03276 (2024), https://arxiv.org/abs/2407.03276
  4. R. Klabbers, J. Lamers, Landscapes of integrable long-range spin chains, arXiv:2405.09718 (2024), https://arxiv.org/abs/2405.09718
  5. R. Klabbers, J. Lamers, Modular families of elliptic long-range spin chains from freezing, Letters in Mathematical Physics (2026), DOI: 10.1007/s11005-026-02075-7, arXiv:2507.13104, https://arxiv.org/abs/2507.13104
  6. A. Liashyk, N. Reshetikhin, I. Sechin, Quantum Integrable Systems on a Classical Integrable Background, Communications in Mathematical Physics (2025), DOI: 10.1007/s00220-025-05523-y, arXiv:2405.17865, https://arxiv.org/abs/2405.17865
  7. F. Gesztesy, R. Weikard, Treibich-Verdier potentials and the stationary (m)KdV hierarchy, Mathematische Zeitschrift 219 (1995) 451-476, DOI: 10.1007/BF02572375, https://doi.org/10.1007/BF02572375
  8. I. Sechin, A. Zotov, R-matrix-valued Lax pairs and long-range spin chains, Physics Letters B 781 (2018) 1-7, DOI: 10.1016/j.physletb.2018.03.062, arXiv:1801.08908, https://arxiv.org/abs/1801.08908
  9. M. Matushko, A. Zotov, Elliptic generalisation of integrable q-deformed anisotropic Haldane-Shastry long-range spin chain, Nonlinearity 36 (2023) 319-353, DOI: 10.1088/1361-6544/aca510, arXiv:2202.01177, https://arxiv.org/abs/2202.01177
  10. V. I. Inozemtsev, On the connection between the one-dimensional S=1/2 Heisenberg chain and Haldane-Shastry model, Journal of Statistical Physics 59 (1990) 1143-1155, DOI: 10.1007/BF01334745, https://doi.org/10.1007/BF01334745
  11. A. P. Polychronakos, Lattice Integrable Systems of Haldane-Shastry Type, Physical Review Letters 70 (1993) 2329-2331, DOI: 10.1103/PhysRevLett.70.2329, arXiv:hep-th/9210109, https://arxiv.org/abs/hep-th/9210109
  12. K. Takasaki, Elliptic Calogero-Moser Systems and Isomonodromic Deformations, Journal of Mathematical Physics 40 (1999) 5787-5821, DOI: 10.1063/1.533056, arXiv:math/9905101, https://arxiv.org/abs/math/9905101
  13. Y. Komori, K. Hikami, Elliptic K-matrix associated with Belavin's symmetric R-matrix, Nuclear Physics B 494 (1997) 687-701, DOI: 10.1016/S0550-3213(97)00177-6, https://doi.org/10.1016/S0550-3213(97)00177-6
  14. G. Felder, V. Pasquier, A simple construction of elliptic R-matrices, Letters in Mathematical Physics 32 (1994) 167-171, DOI: 10.1007/BF00739425, arXiv:hep-th/9402011, https://arxiv.org/abs/hep-th/9402011
  15. R. J. Baxter, Partition function of the Eight-Vertex lattice model, Annals of Physics 70 (1972) 193-228, DOI: 10.1016/0003-4916(72)90335-1, https://doi.org/10.1016/0003-4916(72)90335-1
  16. Y. Shibukawa, K. Ueno, Completely Z symmetric R matrix, Letters in Mathematical Physics 25 (1992) 239-248, DOI: 10.1007/BF00406551, https://doi.org/10.1007/BF00406551
  17. A. Enciso, F. Finkel, A. González-López, M. A. Rodríguez, Haldane-Shastry spin chains of BC_N type, Nuclear Physics B 707 (2005) 553-576, DOI: 10.1016/j.nuclphysb.2004.10.064, arXiv:hep-th/0406054, https://arxiv.org/abs/hep-th/0406054
  18. A. Mostovskii, A. Zotov, Classical elliptic BC_1 Ruijsenaars-van Diejen model: relation to Zhukovsky-Volterra gyrostat and 1-site classical XYZ model with boundaries, Theoretical and Mathematical Physics 226 (2) (2026) 189-216, DOI: 10.1134/S0040577926020029, arXiv:2601.06826, https://arxiv.org/abs/2601.06826