--- title: Periodic Elliptic Ruijsenaars Chain url: https://www.emergentmind.com/topics/periodic-elliptic-ruijsenaars-chain type: topic --- # Periodic Elliptic Ruijsenaars Chain Searching arXiv for recent 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 [2012.05664] [2508.12656]. ## 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 \[ L^a_{ij}(z)=\phi(z,q_i^{a-1}-q_j^a+\eta)\, \frac{\prod_{l=1}^{N}(q_j^a-q_l^{a-1}-\eta)} {(-\eta)\prod_{l\neq j}(q_j^a-q_l^a)}\,e^{p_j^a/c}, \] where \(\phi\) is the elliptic Kronecker function and \(\eta\) is the deformation parameter. The ordered product \[ T(z)=L^1(z)L^2(z)\cdots L^n(z) \] is the monodromy matrix. The Hamiltonian is extracted from the highest-order pole expansion at \(z=0\), \[ \exp(H/c)=\underset{z=0}{\mathrm{res}}\, z^{n-1}T(z), \] which yields \[ H=c\sum_{a=1}^n \log h_{a-1,a}, \qquad h_{a-1,a}= \sum_{j=1}^N \frac{\prod_{l=1}^{N}(q_j^a-q_l^{a-1}-\eta)} {(-\eta)\prod_{l\neq j}(q_j^a-q_l^a)}\,e^{p_j^a/c}. \] This ring topology is not auxiliary: periodicity is built into the definition of the chain, the neighboring-site couplings, and the monodromy construction [2602.08143]. 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 \(N=2\), after center-of-mass reduction, each site is described by a single pair \(({\bf q}_a,{\bf p}_a)\), and the resulting \(2\times 2\) Lax matrices furnish the periodic elliptic Ruijsenaars-Toda chain. The further limit \(\eta\to 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 [2602.08143]. ## 2. Classical \(r\)-matrix structure and integrability A central structural result is the explicit quadratic dynamical \(r\)-matrix Poisson algebra for the site Lax matrices. The bracket has the schematic form \[ c\{L_1^a(z),L_2^b(w)\} = \delta^{ab}\big(L_1^bL_2^b r_{12}^b-r_{12}^{b-1}L_1^bL_2^b +L_1^b s_{12}^{+,b}L_2^b-L_2^b s_{12}^{-,b}L_1^b\big) +\delta^{a,b-1}L_1^{b-1}L_2^b s_{12}^{-,b} -\delta^{a,b+1}L_1^{b+1}L_2^b s_{12}^{+,b+1}, \] with dynamical matrices \(r_{12}^a(z,w)\) and \(s_{12}^{\pm,a}\) built from the elliptic Kronecker function and the first Eisenstein function. This local algebra induces a Poisson bracket for the monodromy matrix \(T(z)\). The resulting trace identities imply \[ \{\operatorname{tr} T^k(z),\operatorname{tr} T^l(w)\}=0, \] which provides Liouville integrability of the periodic elliptic Ruijsenaars chain [2508.12656]. 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\)-matrix structure of the field analogue of the elliptic Calogero-Moser model. This connects the periodic elliptic Ruijsenaars chain to \(1+1\)-dimensional integrable field theory through a precise continuum procedure rather than a purely heuristic analogy [2508.12656]. A related construction starts from a trace of a chain product of Lax matrices, \[ t(z)=\operatorname{tr}(L_1(z)L_2(z)\cdots L_n(z)), \] 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 \(2\)D Toda solutions. This identifies the periodic chain as both a finite-dimensional integrable lattice model and a spatial discretization of a field-theoretic extension [2107.01697]. ## 3. Spectral curves, Toda correspondence, and deformed dynamics Elliptic solutions of the \(2\)D Toda lattice hierarchy provide one of the most direct routes to the elliptic Ruijsenaars-Schneider hierarchy. For tau-functions of the form \[ \tau(x,t,\bar t)=\exp\!\left(-\sum_k t_k\bar t_k\right) \prod_{i=1}^N \sigma\bigl(x-x_i(t,\bar t)\bigr), \] the zeros \(x_i\) move as particles of the elliptic Ruijsenaars-Schneider model. The higher \(t_m\) and \(\bar t_m\) 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 \[ X(z)=\sum_{m\geq 1} H_m z^{-m}, \] while near the second one \[ X(z)=N\eta+\sum_{m\geq 1}\widetilde H_m z^{m}. \] This identifies the full \(2\)D Toda hierarchy with the full hierarchy of commuting Ruijsenaars-Schneider Hamiltonians rather than only the first physical flow [2103.00214]. A distinct but related deformation arises from the Toda lattice with constraint of type \(B\). In that setting, elliptic pole solutions obey the deformed equations of motion \[ \ddot{x}_i+\sum_{j\ne i}\big[\zeta(x_{ij}+\eta)+\zeta(x_{ij}-\eta)-2\zeta(x_{ij})\big] +g(U_i^- - U_i^+)=0, \] with \[ U_i^+=\prod_{j\ne i} \frac{\sigma(x_{ij}+2\eta)\sigma(x_{ij})} {\sigma(x_{ij}+\eta)\sigma(x_{ij}-\eta)}, \qquad U_i^- = U_i^+\ \text{with}\ \eta\mapsto -\eta. \] At \(g=0\), this reduces to the standard periodic elliptic Ruijsenaars-Schneider model. The corresponding commutation representation is a Manakov triple, \[ \dot L+[L,M]=R\,(L-kI), \] and the spectral curve \[ \det(kI-L(k,\lambda))=0 \] admits a holomorphic involution \[ (z,\lambda)\mapsto (z^{-1},\,2N\eta-\lambda). \] Its genus is \(h=2N\), 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 [2302.12085]. 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 \(B\) deformation it survives a non-isospectral evolution because the tracelessness of \(R\) preserves the spectrum encoded by the characteristic equation [2103.00214] [2302.12085]. ## 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 \[ D^{(r)}(p)= \sum_{\substack{I\subset\{1,\dots,n\}\\ |I|=r}} \prod_{i\in I,\,j\notin I} \frac{\theta(x_i/x_j;p)}{\theta(t x_i/x_j;p)} \prod_{i\in I} T_{q,x_i}, \] and the spectral problem is to solve \[ D^{(r)}(p)\Psi(x;p)=\lambda^{(r)}(p)\Psi(x;p). \] A perturbative construction produces two kinds of eigenfunctions: elliptic deformations of the Macdonald polynomials \(P_\lambda(x;p)\), and asymptotically free eigenfunctions \(f(x;s;p)\). 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 [2012.05664]. A complementary proposal is furnished by non-stationary Ruijsenaars functions. These functions are explicit series over \(N\)-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 \(\kappa\to 1\) limit. The same work introduces \(\mathcal T\)-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 [2006.07171]. Dual spectral structures also appear. Symmetric elliptic polynomials \(E_\lambda(x)\) are eigenfunctions of the dual elliptic Ruijsenaars-Schneider Hamiltonians acting on the mother-function variable \(y_i\), while their orthogonal complements \(P_R(x)\) 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 [2103.02508]. The operator algebra extends further through \(Q\)-operators. A one-parameter commuting family \(Q_c\) 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 \(Q_c\) at special singularity loci [2503.18057]. ## 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 \(R\)-matrices. The relevant Hilbert space is \[ H=(\mathbb C^M)^{\otimes N}, \] and the commuting difference operators take the form \[ D_k=\sum_{1\le i_1<\dots