--- title: Hyperbolic Ruijsenaars Model url: https://www.emergentmind.com/topics/hyperbolic-ruijsenaars-model type: topic --- # Hyperbolic Ruijsenaars Model The hyperbolic Ruijsenaars model, usually called the hyperbolic Ruijsenaars–Schneider system, is a relativistic integrable deformation of Calogero–Sutherland-type many-body systems in which the quantum Hamiltonians are commuting finite-difference operators with hyperbolic coefficients. In the \(A\)-type setting it admits an \(N\)-particle modular-double formulation with periods \(\omega_1,\omega_2\), a distinguished two-particle reduction with explicitly known eigenfunctions, Baxter \(Q\)-operators, bispectrality, and exact integral transforms; it also appears in Liouville conformal field theory, quantum Teichmüller theory, Poisson–Lie reduction, and gauge-theoretic Coulomb-branch constructions [2508.12449] [2411.19659] [2603.24123]. ## 1. Definition and operator realizations A standard \(N\)-particle hyperbolic Ruijsenaars Hamiltonian acts on a wave function \(\psi(\{x_j\})\) by \[ \bigl(H^{(\omega_1)}_{\mathrm{RS}}\psi\bigr)(\{x_j\}) =\sum_{j=1}^{N}\left[\prod_{k\neq j}\frac{\sinh\!\Bigl(\tfrac{\pi}{\omega_2}(x_j-x_k+i g)\Bigr)}{\sinh\!\Bigl(\tfrac{\pi}{\omega_2}(x_j-x_k)\Bigr)}\right]\psi(\ldots,x_j+i\omega_1,\ldots), \] and it has a commuting partner obtained by interchanging \(\omega_1\) and \(\omega_2\): \[ \bigl(H^{(\omega_2)}_{\mathrm{RS}}\psi\bigr)(\{x_j\}) =\sum_{j=1}^{N}\left[\prod_{k\neq j}\frac{\sinh\!\Bigl(\tfrac{\pi}{\omega_1}(x_j-x_k+i g)\Bigr)}{\sinh\!\Bigl(\tfrac{\pi}{\omega_1}(x_j-x_k)\Bigr)}\right]\psi(\ldots,x_j+i\omega_2,\ldots). \] The two operators commute, and their coefficients are governed by the modular-double periods \(\omega_1,\omega_2\) and the coupling \(g\) [2508.12449]. In the Macdonald gauge, the same system is described by the commuting family \[ M_r(x; g|\omega) =\sum_{\substack{I\subset [N]\\ |I|=r}} \Biggl(\prod_{i\in I}\prod_{j\notin I} \frac{\sinh\!\big(\tfrac{\pi}{\omega_2}(x_i-x_j)\big)} {\sinh\!\big(\tfrac{\pi}{\omega_2}(x_i-x_j-g)\big)}\Biggr)\,T_{-\omega_1}^{I}, \] where \(T_a^{(i)}=e^{a\partial_{x_i}}\) and \(T_a^I=\prod_{i\in I}T_a^{(i)}\). These operators are gauge-equivalent to the Ruijsenaars operators \(H_r\) via a diagonal factor built from the double sine \(S_2\): \[ p^{-1}(x)\,M_r(x;g|\omega)\,p(x)=H_r(x;g|\omega), \qquad p(x)=\prod_{1\le i0,\qquad \mathrm{Re}\,\omega_2>0,\qquad 0<\mathrm{Re}\,g<\mathrm{Re}(\omega_1+\omega_2), \] with the stronger condition \(\mathrm{Re}\,g<\mathrm{Re}\,\omega_2\) often imposed for adjointness and contour-separation arguments [2303.06382]. The special-function backbone is the hyperbolic gamma or double-sine function. In \((\omega_1,\omega_2)\)-notation, \(\gamma^{(2)}(u;\omega_1,\omega_2)\) satisfies \[ \frac{\gamma^{(2)}(u+\omega_1)}{\gamma^{(2)}(u)}=2\sin\!\Bigl(\frac{\pi u}{\omega_2}\Bigr),\qquad \frac{\gamma^{(2)}(u+\omega_2)}{\gamma^{(2)}(u)}=2\sin\!\Bigl(\frac{\pi u}{\omega_1}\Bigr), \] together with the reflection relation \[ \gamma^{(2)}(u)\,\gamma^{(2)}(\omega_1+\omega_2-u)=1. \] In Liouville notation one writes \(b=\sqrt{\omega_1/\omega_2}\) and \(Q=b+b^{-1}\) [2508.12449]. ## 2. The two-particle hyperbolic model In the two-particle center-of-mass reduction, the hyperbolic Ruijsenaars operator takes a particularly explicit finite-difference form. In the \(b\)-parametrization used in the Liouville-theoretic treatment, \[ H=\frac{\sin(\pi b(x-g))}{\sin(\pi b x)}\,e^{-b\partial_x} + \frac{\sin(\pi b(x+g))}{\sin(\pi b x)}\,e^{b\partial_x}, \qquad (e^{\pm b\partial_x}\psi)(x)=\psi(x\pm b), \] so the shift step is \(b\). Its distinguished eigenfunction is \[ F^g_{\lambda}(x)=\frac{1}{S_b(g)}\int_{-i\infty}^{i\infty} S_b(\pm z+\lambda/2+g/2)\,S_b(\pm z-\lambda/2+g/2)\, e^{2i\pi z x}\,\frac{dz}{i}, \] and it satisfies the spectral equation \[ H\big[F^g_{\lambda}\big](x)=2\cos(\pi b\lambda)\,F^g_{\lambda}(x). \] Here \(x\in\mathbb{R}\), \(\lambda\in\mathbb{R}\), the integration variable lies on \(i\mathbb{R}\), and the spectrum is continuous with eigenvalues \(2\cos(\pi b\lambda)\) [2603.24123]. The same two-particle system can be written in modular-double variables as \[ H_h =\frac{\sh\!\bigl(\tfrac{\pi}{\omega_2}(x-ig)\bigr)}{\sh\!\bigl(\tfrac{\pi}{\omega_2}x\bigr)}\,e^{-i\omega_1\partial_x} +\frac{\sh\!\bigl(\tfrac{\pi}{\omega_2}(x+ig)\bigr)}{\sh\!\bigl(\tfrac{\pi}{\omega_2}x\bigr)}\,e^{+i\omega_1\partial_x}, \] together with the commuting partner \[ H_h' =\frac{\sh\!\bigl(\tfrac{\pi}{\omega_1}(x-ig)\bigr)}{\sh\!\bigl(\tfrac{\pi}{\omega_1}x\bigr)}\,e^{-i\omega_2\partial_x} +\frac{\sh\!\bigl(\tfrac{\pi}{\omega_1}(x+ig)\bigr)}{\sh\!\bigl(\tfrac{\pi}{\omega_1}x\bigr)}\,e^{+i\omega_2\partial_x}. \] Their joint eigenfunctions are analytic in \(x\) and can be expressed in Mellin–Barnes form through \(\gamma^{(2)}\) [2508.12449]. A characteristic structural property is self-duality: \[ F^g_{\lambda}(x)=F^{Q-g}_{x}(\lambda). \] This exchanges the coordinate and spectral variables while reflecting the coupling \(g\mapsto Q-g\) [2603.24123]. ## 3. Baxter operators, product formulas, and spectral transform For the two-particle model, the product formula \[ 4\,F^g_{\lambda}(x_1)\,F^g_{\lambda}(x_2) =S_b(g)\int_{-i\infty}^{i\infty} \frac{S_b(g\pm z)}{S_b(\pm z)}\, S_b\!\left(\frac{\bar g\pm z\pm x_1\pm x_2}{2}\right)\, F^g_{\lambda}(z)\,\frac{dz}{i}, \qquad \bar g=Q-g, \] leads directly to an integral Baxter operator \(Q^\rho\) with kernel built from the same \(S_b\)-factors. Its action on the eigenfunction is diagonal: \[ \big[Q^{\rho}F^g_{\lambda}\big](x) =\frac{4\,S_b(\pm \rho\pm\lambda/2+g/2)}{(S_b(g))^2}\,F^g_{\lambda}(x), \] and the operators satisfy \[ [Q^\rho,H]=0,\qquad [Q^\rho,Q^\lambda]=0. \] In that two-particle Liouville-derived framework, a standard scalar \(T\)–\(Q\) functional relation is not written explicitly; the characterization is by commutation and diagonal action on \(F^g_\lambda\) [2603.24123]. For the \(A_{N-1}\) hyperbolic system, Baxter operators form a commuting family of integral operators \[ (Q_N(\lambda)f)(x)=d_N(g|\omega)\int_{\mathbb{R}^N}dy\; e^{2\pi i\lambda(\sum_i x_i-\sum_i y_i)}\,K(x,y)\,\prod_{j=1}^N p(y_j)\,f(y), \] and they commute with the Macdonald operators. The Hallnäs–Ruijsenaars joint eigenfunction \(\Phi_N(x;\lambda)\), constructed recursively by the raising operators \(A_N(\lambda_N)\cdots A_2(\lambda_2)\), diagonalizes them: \[ Q_N(\mu)\,\Phi_N(x;\lambda) =\Big(\prod_{j=1}^N K(\mu-\lambda_j)\Big)\,\Phi_N(x;\lambda), \] while simultaneously satisfying the Macdonald spectral problem \[ M_r(x;g|\omega)\,\Phi_N(x;\lambda) =e_r\!\big(e^{\tfrac{2\pi}{\omega_1}\lambda_1},\dots,e^{\tfrac{2\pi}{\omega_1}\lambda_N}\big)\,\Phi_N(x;\lambda). \] The same wave function is therefore jointly diagonal for the commuting hyperbolic Macdonald and Baxter families [2303.06382]. Orthogonality and completeness were then established by exploiting Baxter diagonalization and coordinate–spectral duality. The resulting transform is an \(n\)-particle analogue of Fourier inversion: for symmetric functions, \[ [Tf](\lambda)=\int_{\mathbb{R}^n}dx\,\mu(x)\,\Psi_\lambda(-x)\,f(x), \qquad [T^\dagger\chi](x)=\int_{\mathbb{R}^n}d\lambda\,\hat\mu(\lambda)\,\Psi_\lambda(x)\,\chi(\lambda), \] with \[ [T^\dagger T f](x)=f(x). \] In distributional form, the wave functions satisfy orthogonality and completeness relations \[ (\Psi_\lambda(y),\Psi_\lambda(x))_\mu=\mu^{-1}(x)\delta(x,y), \qquad (\Psi_x(\lambda),\Psi_x(\xi))_{\hat\mu}=\hat\mu^{-1}(\lambda)\delta(\lambda,\xi), \] and the transform is unitary in four regimes: I: \(w_1,w_2\in\mathbb{R}\), \(g\in\mathbb{R}\); II: \(\bar w_1=w_2\), \(g\in\mathbb{R}\); III: \(w_1,w_2\in\mathbb{R}\), \(g=g^*\); IV: \(\bar w_1=w_2\), \(g=g^*\) [2307.16817] [2411.19659]. ## 4. Liouville CFT, quantum Teichmüller theory, and coupling reflection A central recent result is that the two-particle hyperbolic Ruijsenaars Hamiltonian, its Baxter operator, and the product formula for its eigenfunctions can be derived from the genus-one Moore–Seiberg identity in Liouville conformal field theory. The Liouville parameters are identified by \[ b\leftrightarrow b,\qquad Q=b+b^{-1},\qquad \alpha=\frac{Q}{2}+ip,\qquad \lambda=2p,\qquad g=\beta_5, \] and the Ruijsenaars wavefunction coincides, up to known normalization factors, with the Liouville modular one-point \(S\)-kernel. A degenerate insertion \(\alpha_1=-b/2\) reduces the Moore–Seiberg identity to a two-term finite-difference equation reproducing the Ruijsenaars Hamiltonian, while the condition \[ \alpha_2-\alpha_1=\beta_5=\beta_3 \] leads to the exact product formula. The authors further expect extension to \(N\)-particle hyperbolic Ruijsenaars systems via conformal Toda field theory and to supersymmetric variants via super Liouville theory [2603.24123]. A different geometric realization appears in \(GL(2)\) quantum Teichmüller theory on the punctured torus. There, the order-four modular element \[ S=\begin{pmatrix}0&-1\\ 1&0\end{pmatrix}\in SL(2,\mathbb{Z}) \] acts on the Hilbert space of the quantized cluster Poisson variety, and the matrix coefficient \[ \Phi_\mu(\lambda)=\langle \Psi_\lambda,\, I_S\,\Psi_\mu\rangle \] is a joint eigenfunction of the two-particle hyperbolic Macdonald operators. In this framework the Hallnäs–Ruijsenaars eigenfunction is proportional to the modular \(S\)-matrix coefficient, and \(GL(2)\) Macdonald polynomials arise as special values of its analytic continuation. The construction uses an \(SL(2,\mathbb{Z})\)-equivariant embedding of the \(GL(2)\) spherical DAHA into the quantized coordinate ring of the cluster Poisson variety [2402.14214]. Coupling reflection is encoded by a second Baxter family. Writing \[ g^*=\omega_1+\omega_2-g, \] one introduces \(Q^*(\lambda)\), conjugate to the original \(Q_n(\lambda)\) at reflected coupling. The two Baxter families commute with each other and with the Macdonald operators, and the wave function obeys the exact symmetry \[ \Psi_n(x_n; g|w) = \nu(x_n)^{-1}\,\tilde\nu(\Lambda)^{-1}\,\Psi_n(x_n; g^*|w). \] This refines the coordinate–spectral duality by adding a direct \(g\leftrightarrow g^*\) reflection symmetry at the level of wave functions [2308.07619]. ## 5. Classical structure, Lax formulations, and spin generalizations The classical hyperbolic Ruijsenaars model arises by Poisson reduction from the Heisenberg double of \(GL(N,\mathbb{C})\). In reduced variables \(z_i=e^{q_i}\), \(P_i=e^{p_i}\), the Lax matrix can be written as \[ L = \sum_{i,j=1}^{N} \frac{(1-\omega)\, z_i}{z_i - \omega\, z_j}\; \prod_{a\neq j}^{N}\frac{\omega\, z_j - z_a}{z_j - z_a}\;P_j^{-1}\,E_{ij}, \qquad \omega=e^{-\gamma}. \] Its spectral invariants \[ I_k=\mathrm{Tr}\,L^k \] commute, and the Poisson bracket of \(L\) has quadratic dynamical \(r\)-matrix form. Upon introduction of a spectral parameter, the corresponding baxterized \(r\)-matrices satisfy shifted Yang–Baxter-type equations; on the quantum side, trace-type formulae for commuting integrals are conjectured to generate the same commutative ring as the Macdonald operators [1902.06755]. A structurally related family is the hyperbolic spin Ruijsenaars–Schneider model. In one Poisson-reduction derivation, the spin degrees of freedom are rectangular matrices \[ a\in \mathrm{Mat}_{N,\ell}(\mathbb{C}),\qquad b\in \mathrm{Mat}_{\ell,N}(\mathbb{C}), \] the phase space is built from the Heisenberg double together with a deformed oscillator manifold, and the reduced Lax matrix has collective spin entries \((\mathbf a\mathbf c)_{ij}\). The commuting Hamiltonians remain \[ H_k=\mathrm{Tr}L^k, \] while the model carries Poisson–Lie symmetry of the spin group \(GL_\ell(\mathbb{C})\), providing superintegrability [1906.02619]. A gauge-theoretic realization identifies the hyperbolic spin model with the \(K\)-theoretic Coulomb branch of 3d \(\mathcal N=4\) affine \(A\) necklace quiver gauge theory. In that description the positions are \(x_i=\log Q_i^0\), the relativistic coupling is \(\gamma=-\log t\), the Hamiltonians are \[ H[n]=\operatorname{Tr}L^n, \] and the equations of motion reproduce the Krichever–Zabrodin hyperbolic spin RS system with potential \[ V(z)=\tfrac{1}{2}\coth \tfrac{z}{2}-\tfrac{1}{2}\coth \tfrac{z+\gamma}{2}. \] The same construction exhibits the classical limit of the quantum toroidal algebra of \(\mathfrak{gl}_\ell\), and the Hamiltonians lie in the center of this symmetry, making superintegrability manifest [2603.03048]. ## 6. Related families, lattice versions, and degenerations The hyperbolic Ruijsenaars model sits inside a larger family of Ruijsenaars–Schneider–van Diejen systems. In the two-parameter hyperbolic van Diejen model, \[ H_g=\sum_{a=1}^n\cosh(\theta_a)\,u_a(p;g), \] the interaction factor \(u_a\) contains both \(\lambda_a-\lambda_c\) and \(\lambda_a+\lambda_c\) terms as well as a one-body contribution depending on \(2\lambda_a\). Setting \(v=0\) removes the one-body term, and suppressing the \(BC\)-type sum interactions recovers the \(A\)-type hyperbolic RS Hamiltonian structure. The same work establishes global action–angle variables, self-duality, and a factorized scattering map for the van Diejen system, thereby extending Ruijsenaars-type phenomena beyond the translation-invariant \(A\)-type case [1701.08558]. A different \(BC_n\)-symmetric hyperbolic Ruijsenaars-type system is obtained by Hamiltonian reduction on the Heisenberg double of \(SU(n,n)\). Its reduced Hamiltonian is a three-parameter relativistic many-body system with canonical coordinates \((q_i,p_i)\), and its cotangent-bundle limit reproduces the standard three-parameter hyperbolic \(BC_n\) Sutherland Hamiltonian. The construction is explicitly distinct from van Diejen’s \(BC_n\) relativistic models [1311.4641]. There is also a lattice hyperbolic Ruijsenaars model with an exponential Morse term. On the lattice \(p+\Lambda\), its \(n\)-particle Hamiltonian is diagonalized by multivariate continuous dual \(q\)-Hahn polynomials, obtained as a parameter reduction of Macdonald–Koornwinder polynomials. This produces an exact spectral transform, a bispectral dual \(q\)-difference system, a full commuting family of self-adjoint operators, and an \(n\)-particle scattering operator whose phase factor decomposes into two-body RS terms and one-body Morse contributions [1508.03829]. Several limits connect the hyperbolic model to other integrable systems. In the two-particle setting, the complex rational degeneration \(\omega_1+\omega_2\to 0\), equivalently \(b\to i\), yields commuting finite-difference operators with rational coefficients in complex variables, Mellin–Barnes and Euler integral representations for eigenfunctions, and complex limits of the hyperbolic Baxter operators. A complementary degeneration occurs at \(\omega_1-\omega_2\to 0\), or \(b\to 1\). The nonrelativistic limit recovers hyperbolic Calogero–Sutherland operators, while additional complex limits produce commuting hypergeometric operators on the cylinder [2508.12449]. Source: https://www.emergentmind.com/topics/hyperbolic-ruijsenaars-model