Periodic Elliptic Ruijsenaars Chain

• Periodic elliptic Ruijsenaars chain is an integrable system with doubly periodic elliptic dependence, formulated through both quantum difference operators and classical lattice Lax matrices.
• It incorporates a quadratic dynamical r-matrix structure that guarantees Liouville integrability and connects to Toda hierarchies, spectral curves, and gauge theory.
• Advanced quantization techniques yield commuting Hamiltonians and dual spectral structures, paving the way for spin generalizations and compact phase-space completions.

Searching arXiv for papers on the periodic elliptic Ruijsenaars chain and closely related structures. The periodic elliptic Ruijsenaars chain denotes an integrable system with doubly periodic elliptic dependence that appears in two closely related formulations in the literature represented here: as a quantum family of commuting elliptic difference operators, and as a classical periodic lattice chain built from site Lax matrices and a monodromy matrix. In the quantum setting, the commuting Hamiltonians are elliptic Ruijsenaars difference operators acting on functions of several variables; in the classical setting, one studies a ring of $n$ sites with monodromy $T(z)=L^1(z)L^2(z)\cdots L^n(z)$. Across these formulations, the subject is organized by Lax and Manakov representations, spectral curves, explicit eigenfunction constructions, $R$-matrix identities, and correspondences with Toda hierarchies, long-range spin chains, compact phase spaces, gauge theory, and field analogues (Langmann et al., 2020, Murinov et al., 18 Aug 2025).

1. Classical periodic chain

In the periodic lattice formulation, the phase space consists of $n$ sites, each carrying $N$ canonical pairs $(q_i^a,p_i^a)$, with Poisson brackets

$$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$

and site indices taken modulo $n$, so that $q_i^{n+1}=q_i^1$ and similarly for the other variables. The local Lax matrix at site $a$ is

$T(z)=L^1(z)L^2(z)\cdots L^n(z)$0

where $T(z)=L^1(z)L^2(z)\cdots L^n(z)$1 is the elliptic Kronecker function and $T(z)=L^1(z)L^2(z)\cdots L^n(z)$2 is the deformation parameter. The ordered product

$T(z)=L^1(z)L^2(z)\cdots L^n(z)$3

is the monodromy matrix. The Hamiltonian is extracted from the highest-order pole expansion at $T(z)=L^1(z)L^2(z)\cdots L^n(z)$4,

$T(z)=L^1(z)L^2(z)\cdots L^n(z)$5

which yields

$T(z)=L^1(z)L^2(z)\cdots L^n(z)$6

This ring topology is not auxiliary: periodicity is built into the definition of the chain, the neighboring-site couplings, and the monodromy construction (Murinov et al., 8 Feb 2026).

The same framework specializes to lower-rank or reduced models. In particular, the elliptic Ruijsenaars-Toda chain and the elliptic Toda chain are obtained as particular cases of the elliptic Ruijsenaars chain. For $T(z)=L^1(z)L^2(z)\cdots L^n(z)$7, after center-of-mass reduction, each site is described by a single pair $T(z)=L^1(z)L^2(z)\cdots L^n(z)$8, and the resulting $T(z)=L^1(z)L^2(z)\cdots L^n(z)$9 Lax matrices furnish the periodic elliptic Ruijsenaars-Toda chain. The further limit $R$0 produces the elliptic Toda chain introduced by Krichever. This places the periodic elliptic Ruijsenaars chain in a reduction scheme where relativistic and non-relativistic nearest-neighbor chains arise from the same elliptic parent system (Murinov et al., 8 Feb 2026).

2. Classical $R$1-matrix structure and integrability

A central structural result is the explicit quadratic dynamical $R$2-matrix Poisson algebra for the site Lax matrices. The bracket has the schematic form

$R$3

with dynamical matrices $R$4 and $R$5 built from the elliptic Kronecker function and the first Eisenstein function. This local algebra induces a Poisson bracket for the monodromy matrix $R$6. The resulting trace identities imply

$R$7

which provides Liouville integrability of the periodic elliptic Ruijsenaars chain (Murinov et al., 18 Aug 2025).

The same analysis yields a continuous non-relativistic limit. In that limit, the discrete site index becomes a spatial coordinate, the lattice Poisson algebra passes to a non-ultralocal Maillet-type bracket, and one reproduces the $R$8-matrix structure of the field analogue of the elliptic Calogero-Moser model. This connects the periodic elliptic Ruijsenaars chain to $R$9-dimensional integrable field theory through a precise continuum procedure rather than a purely heuristic analogy (Murinov et al., 18 Aug 2025).

A related construction starts from a trace of a chain product of Lax matrices,

$n$0

which generates commuting Hamiltonians and leads to a lattice field analogue of the Ruijsenaars-Schneider model with continuous time. In this formulation the model is gauge equivalent to a classical elliptic spin chain, and the lattice version coincides with the system derived from elliptic families of $n$1D Toda solutions. This identifies the periodic chain as both a finite-dimensional integrable lattice model and a spatial discretization of a field-theoretic extension (Zabrodin et al., 2021).

3. Spectral curves, Toda correspondence, and deformed dynamics

Elliptic solutions of the $n$2D Toda lattice hierarchy provide one of the most direct routes to the elliptic Ruijsenaars-Schneider hierarchy. For tau-functions of the form

$n$3

the zeros $n$4 move as particles of the elliptic Ruijsenaars-Schneider model. The higher $n$5 and $n$6 flows are governed by Hamiltonians extracted from expansions of the spectral curve of the Ruijsenaars Lax matrix at two marked points. Near one marked point one obtains

$n$7

while near the second one

$n$8

This identifies the full $n$9D Toda hierarchy with the full hierarchy of commuting Ruijsenaars-Schneider Hamiltonians rather than only the first physical flow (Prokofev et al., 2021).

A distinct but related deformation arises from the Toda lattice with constraint of type $N$0. In that setting, elliptic pole solutions obey the deformed equations of motion

$N$1

with

$N$2

At $N$3, this reduces to the standard periodic elliptic Ruijsenaars-Schneider model. The corresponding commutation representation is a Manakov triple,

$N$4

and the spectral curve

$N$5

admits a holomorphic involution

$N$6

Its genus is $N$7, and the pole-ansatz wave function is identified with the Baker-Akhiezer function on the spectral curve. The same paper also suggests a field-theory extension through more general elliptic families (Prokofev et al., 2023).

These Toda correspondences clarify what the spectral curve controls. It is not only a generating device for integrals of motion: in the hierarchy picture it organizes the positive and negative flows, while in the type $N$8 deformation it survives a non-isospectral evolution because the tracelessness of $N$9 preserves the spectrum encoded by the characteristic equation (Prokofev et al., 2021, Prokofev et al., 2023).

4. Quantum difference operators and spectral theory

On the quantum side, the periodic elliptic Ruijsenaars system is generated by a commuting family of elliptic difference operators. One standard form is

$(q_i^a,p_i^a)$0

and the spectral problem is to solve

$(q_i^a,p_i^a)$1

A perturbative construction produces two kinds of eigenfunctions: elliptic deformations of the Macdonald polynomials $(q_i^a,p_i^a)$2, and asymptotically free eigenfunctions $(q_i^a,p_i^a)$3. These are given by convergent infinite series in suitable domains, and in the domain where the elliptic Ruijsenaars operators define a relativistic quantum mechanical system, the elliptic deformations of the Macdonald polynomials provide a family of orthogonal functions with respect to the pertinent scalar product (Langmann et al., 2020).

A complementary proposal is furnished by non-stationary Ruijsenaars functions. These functions are explicit series over $(q_i^a,p_i^a)$4-tuples of partitions with coefficients expressed through Nekrasov factors, and they are conjectured to provide explicit eigenfunctions of the elliptic Ruijsenaars model after a balancing and $(q_i^a,p_i^a)$5 limit. The same work introduces $(q_i^a,p_i^a)$6-difference operators, proves diagonal action in the trigonometric case, and conjectures the analogous diagonal action in the general elliptic case. The series convergence is established in a nontrivial parameter domain, so these objects are genuine special functions rather than only formal expansions (Langmann et al., 2020).

Dual spectral structures also appear. Symmetric elliptic polynomials $(q_i^a,p_i^a)$7 are eigenfunctions of the dual elliptic Ruijsenaars-Schneider Hamiltonians acting on the mother-function variable $(q_i^a,p_i^a)$8, while their orthogonal complements $(q_i^a,p_i^a)$9 are eigenfunctions of the elliptic reduction of the Koroteev-Shakirov Hamiltonians. In that formulation, the duality between the two operator families is mediated by an orthogonality transformation that is natural in the full space of time variables and becomes less transparent after the Miwa transform (Mironov et al., 2021).

The operator algebra extends further through $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$0-operators. A one-parameter commuting family $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$1 is realized as integral operators with kernels built from the elliptic gamma function. Their commutativity is equivalent to an elliptic hypergeometric integral transformation conjectured by Gadde et al., and they act diagonally on elliptic Macdonald polynomials. The Noumi-Sano operators appear as discrete degenerations, obtained from residues of $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$2 at special singularity loci (Rains et al., 23 Mar 2025).

5. Spin generalizations and freezing to long-range chains

Matrix-valued or spin versions of the periodic elliptic Ruijsenaars system replace scalar coefficients by Baxter-Belavin $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$3-matrices. The relevant Hilbert space is

$$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$4

and the commuting difference operators take the form

$$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$5

For $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$6, these operators reduce to the scalar elliptic Macdonald-Ruijsenaars operators. Their commutativity is equivalent to a family of matrix-valued identities, and the proof relies on the quantum Yang-Baxter equation, the associative Yang-Baxter equation, and unitarity of the elliptic $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$7-matrix (Matushko et al., 2022).

The corresponding $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$8-matrix identities were later isolated and proved in a form that directly generalizes Ruijsenaars’ scalar kernel identities. In the scalar case $$\{p_i^a,q_j^b\}=\delta^{ab}\delta_{ij}, \qquad \{p_i^a,p_j^b\}=\{q_i^a,q_j^b\}=0,$$9, the identities reduce to the elliptic function identities used by Ruijsenaars to construct integral solutions of the quantum spinless model. In the matrix case, they are described as the first step toward constructing solutions of the quantum eigenvalue problem for the anisotropic spin Ruijsenaars model (Matushko et al., 2022).

The freezing mechanism converts these many-body operators into long-range spin chains. For the elliptic $n$0-deformed anisotropic construction, the Polychronakos freezing trick is applied to elliptic spin Ruijsenaars-Macdonald operators, and the required identities are equivalent to the equilibrium conditions of the underlying classical spinless Ruijsenaars-Schneider system. In the $n$1 case the construction provides a $n$2-deformation of the anisotropic XXZ Haldane-Shastry model; trigonometric degenerations recover Uglov’s $n$3-deformed XXZ Haldane-Shastry chain and the standard Haldane-Shastry model (Matushko et al., 2022).

A later development focuses on the elliptic case at fixed elliptic parameter and shows that the spinless elliptic Ruijsenaars-Schneider system possesses a modular family of classical equilibrium configurations. These equilibria typically have constant but nonzero momenta. Using deformation quantisation, one can freeze elliptic spin-Ruijsenaars systems at any classical equilibrium while preserving quantum integrability. The resulting modular families include the Heisenberg, Inozemtsev and Haldane-Shastry chains together with their $n$4-like $n$5-deformations, as well as the antiperiodic Haldane-Shastry chain of Fukui-Kawakami, its elliptic generalisation of Sechin-Zotov, and their completely anisotropic $n$6-deformations due to Matushko-Zotov (Klabbers et al., 17 Jul 2025).

6. Compactifications, reductions, and related elliptic chains

A geometric compactification of the center-of-mass reduced trigonometric and elliptic Ruijsenaars-Schneider systems realizes the completed phase space as $n$7 equipped with the Fubini-Study symplectic structure. The resulting compact real forms are labelled by an integer $n$8 relatively prime to $n$9 and by a coupling parameter $q_i^{n+1}=q_i^1$0 varying in a punctured interval around $q_i^{n+1}=q_i^1$1. The local phase space $q_i^{n+1}=q_i^1$2 embeds symplectically into $q_i^{n+1}=q_i^1$3, and the elliptic Hamiltonian and Lax matrix extend smoothly to the compact phase space. This construction generalizes earlier compactifications that imposed the restriction $q_i^{n+1}=q_i^1$4 (Feher et al., 2016).

Within the periodic chain formalism, the elliptic Ruijsenaars-Toda chain and the elliptic Toda chain appear as particular cases, and both admit classical $q_i^{n+1}=q_i^1$5-matrix structures derived from the parent elliptic Ruijsenaars chain. For the Ruijsenaars-Toda case, the by-product is an explicit gauge equivalence to the discrete Landau-Lifshitz model of XYZ type. For the elliptic Toda chain, the gauge-equivalent spin description is an XYZ chain with special values of the Casimir functions at each site. The same gauge transformation identifies the reduced monodromy with a product of Sklyanin-type Lax operators, so the relation to the XYZ chain is structural rather than merely spectral (Murinov et al., 8 Feb 2026).

Taken together, these results show that the periodic elliptic Ruijsenaars chain is not tied to a single phase-space realization. It admits both noncompact many-body descriptions and compact projective-space completions, and it functions as a master elliptic model from which Ruijsenaars-Toda, elliptic Toda, and XYZ-type chains can be obtained by reduction, limiting procedures, and gauge equivalence (Feher et al., 2016, Murinov et al., 8 Feb 2026).

7. Exact quantization, finite-volume comparisons, and conjectural gauge-theoretic realizations

For the $q_i^{n+1}=q_i^1$6-particle quantum elliptic Ruijsenaars-Schneider model, exact quantization conditions have been proposed from the correspondence with five-dimensional $q_i^{n+1}=q_i^1$7 $q_i^{n+1}=q_i^1$8 gauge theory in the Nekrasov-Shatashvili limit. Two natural sets of quantization conditions arise, related by electro-magnetic duality: a $q_i^{n+1}=q_i^1$9-type quantization that fixes the periods $a$0, and an $a$1-type quantization in terms of the dual periods obtained by derivatives of the full nonperturbatively completed twisted effective superpotential. The paper emphasizes that naive all-order WKB quantization is insufficient because it develops poles at special values of $a$2; non-perturbative corrections in $a$3 are required to obtain pole-free quantization conditions compatible with modular duality. The eigenfunction problem is recast through Separation of Variables and the Baxter equation, whose entire solutions occur only at quantized values (Hatsuda et al., 2018).

A separate line of work compares the spectrum of the elliptic Ruijsenaars-Schneider model with the finite-size spectrum of sine-Gordon theory in the two-particle sector. The analytic and numerical conclusion is that the Bethe-Yang quantization conditions are exact in the rational and trigonometric cases, but acquire finite-size corrections in the hyperbolic and elliptic cases, both for the relativistic model and its non-relativistic limit. The same work stresses that these corrections are different in quantum field theories and in many-body systems. Thus the infinite-volume correspondence based on scattering phases does not extend straightforwardly to the finite-volume periodic elliptic setting (Bajnok et al., 20 Nov 2025).

Recent gauge-theoretic work establishes that cohomological and $a$4-theoretic Coulomb branches of $a$5d $a$6 necklace quiver gauge theories reproduce the rational and hyperbolic spin Ruijsenaars-Schneider models. For the elliptic case, the paper states a conjecture: the Poisson algebras of elliptic Coulomb branches should similarly reproduce the elliptic spin Ruijsenaars-Schneider model, presumably through an elliptic $a$7-operator algebra with elliptic analogues of the $a$8-matrices governing the rational and hyperbolic cases. This conjectural extension places the periodic elliptic Ruijsenaars chain within a broader program relating monopole operators, Coulomb branches, and classical integrable systems (Arutyunov et al., 3 Mar 2026).

The current picture is therefore structurally coherent but not closed. Exact quantization is available in a gauge-theoretic formulation, finite-volume behavior is now sharply distinguished from field-theoretic analogues, and several elliptic spin and Coulomb-branch correspondences remain explicitly conjectural rather than fully constructed (Hatsuda et al., 2018, Bajnok et al., 20 Nov 2025, Arutyunov et al., 3 Mar 2026).