Hyperbolic Calogero-Sutherland System

• The hyperbolic Calogero–Sutherland system is an integrable many-body model defined by inverse hyperbolic-square interactions and organized by root systems such as Aₙ₋₁, Cₙ, and BCₙ.
• Its versatility is demonstrated through various formulations including spinless, spin, supersymmetric, and boundary-coupled models, making it central in both classical and quantum integrability.
• Lax pair constructions, Hamiltonian reductions, and connections to special-function theory emphasize its role in linking differential equations with advanced algebraic structures.

Searching arXiv for recent and foundational papers on the hyperbolic Calogero–Sutherland system and close variants. The hyperbolic Calogero–Sutherland system is a class of classical and quantum integrable many-body systems whose defining interactions are inverse hyperbolic-square kernels, typically of $1/\sinh^2$ type, and whose realizations are organized by root systems such as $A_{n-1}$, $C_n$, and $BC_n$. In current work the family encompasses spinless and spin models, supersymmetric extensions, boundary-coupled $BC_n$ systems, Hamiltonian and Hitchin reductions, Lax and $r$-matrix formulations, and exact quantum wave functions described by Baxter operators and Heckman–Opdam hypergeometric functions (Belousov et al., 26 Aug 2025, Pusztai, 2012, Diejen et al., 2023).

1. Standard hyperbolic Hamiltonians

For the $A_{n-1}$ quantum model, a standard nonrelativistic Hamiltonian is

$$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$

together with a commuting family obtained after gauge transformation by the weight $w(x)=|2\sinh \pi x|^g$ (Belousov et al., 26 Aug 2025). In the same setting, the basic kernel $\mathcal K(x)=(2\cosh \pi x)^{-g}$ appears as the hyperbolic limit of the Ruijsenaars double-sine kernels, so the hyperbolic interaction is tied from the outset to a bispectral integral-operator structure rather than only to a Schrödinger operator (Belousov et al., 26 Aug 2025).

In classical gauge-theoretic realizations, the same hyperbolic interaction emerges after diagonalizing a positive Hermitian matrix variable $A_{n-1}$0, eliminating off-diagonal gauge fields, and passing to logarithmic coordinates $A_{n-1}$1. The resulting bosonic action is

$A_{n-1}$2

which is the standard hyperbolic Calogero–Sutherland form of $A_{n-1}$3 type (Fedoruk et al., 2019). From the nonrelativistic limit of the trigonometric Ruijsenaars–Schneider model one also recovers the classical Lax matrix

$A_{n-1}$4

so the hyperbolic model sits naturally as a limit of a relativistic system with the same spectral data (Beketov et al., 2015).

2. Root-system realizations and boundary terms

Beyond the $A_{n-1}$5 case, the hyperbolic Calogero–Sutherland family is organized by non-$A_{n-1}$6 root systems through additional sum-type and boundary interactions. A standard $A_{n-1}$7 Hamiltonian is

$A_{n-1}$8

on the Weyl chamber

$A_{n-1}$9

with three independent couplings $C_n$0 in the repulsive regime (Pusztai, 2012). In this formulation the pair terms $C_n$1 encode the roots $C_n$2, while the one-body terms encode $C_n$3 and $C_n$4.

The $C_n$5 framework also supports more structured particle content. One Hamiltonian-reduction model interprets the system as two species of particles on the half-line, with same-species interactions of $C_n$6 type, cross-species interactions of $C_n$7 type, and species-dependent boundary terms involving $C_n$8, $C_n$9, and $BC_n$0 (Ayadi et al., 2011). In the three-coupling $BC_n$1 model, setting $BC_n$2 gives the hyperbolic $BC_n$3 Sutherland model, while $BC_n$4 gives the $BC_n$5 model (Pusztai, 2011). This root-system description is the mechanism by which hyperbolic Calogero–Sutherland dynamics acquires walls, mirror interactions, and boundary fields.

3. Reduction, Lax representations, and $BC_n$6-matrices

A central structural feature of the hyperbolic Calogero–Sutherland system is its derivation by reduction from higher-dimensional free dynamics. For the $BC_n$7 case, symplectic reduction of free geodesic motion on $BC_n$8 with a shifted momentum map produces the physical phase space, the hyperbolic Hamiltonian, a Lax matrix $BC_n$9, and a dynamical $BC_n$0-matrix. The reduced bracket takes the form

$BC_n$1

with $BC_n$2 depending only on the configuration variables $BC_n$3, not on the momenta $BC_n$4, and not explicitly on the coupling parameters $BC_n$5 used in the reduction (Pusztai, 2012). The corresponding Lax equation is equivalent to Hamilton’s equations, so integrability is encoded directly at the matrix level (Pusztai, 2012).

A different but equally explicit construction appears in the generalized two-spin hyperbolic model. There the Hamiltonian is

$BC_n$6

where $BC_n$7 and $BC_n$8 are two skew-symmetric spin matrices on commuting copies of $BC_n$9 (Kharchev et al., 2017). The model is realized as a Hitchin system on a singular curve formed by gluing two copies of $r$0; the reduced phase space is

$r$1

and complete integrability is established by an explicit Lax pair and a spectral-parameter-dependent classical $r$2-matrix on that singular curve (Kharchev et al., 2017). When one spin variable vanishes, the construction reduces to the standard hyperbolic spin Calogero–Sutherland model; for $r$3, it reproduces the $r$4 potential (Kharchev et al., 2017).

4. Duality, Ruijsenaars relatives, and action-angle structure

Hyperbolic Calogero–Sutherland systems are part of a larger reduction-theoretic hierarchy in which Sutherland and Ruijsenaars-type models appear as dual coordinate charts on the same reduced space. For the $r$5 root system, one common symplectic reduction yields both the hyperbolic $r$6 Sutherland model

$r$7

and a rational $r$8 Ruijsenaars–Schneider–van Diejen model, with global Darboux coordinates on both sides and a symplectomorphism under which the action variables of one system become the position variables of the other (Pusztai, 2011). The same pattern holds for the three-coupling $r$9 case: the hyperbolic $A_{n-1}$0 Sutherland model and the rational $A_{n-1}$1 Ruijsenaars–Schneider–van Diejen model arise from one reduction picture, and the positive eigenvalues of the Sutherland Lax matrix provide the action variables on the Sutherland side (Pusztai, 2011).

A relativistic companion to this picture is obtained by Poisson–Lie reduction on the Heisenberg double of $A_{n-1}$2. The reduced Hamiltonian is a $A_{n-1}$3-symmetric Ruijsenaars-type model with multiplicative hyperbolic interactions and one-body boundary terms, and in the cotangent bundle or nonrelativistic limit it reduces to the standard three-parameter hyperbolic $A_{n-1}$4 Sutherland Hamiltonian (Marshall, 2013). This places the hyperbolic Calogero–Sutherland system inside a hierarchy

$A_{n-1}$5

so the hyperbolic system functions both as an autonomous integrable model and as a nonrelativistic shadow of a relativistic one (Marshall, 2013).

5. Quantum wave functions and special-function theory

At the quantum level, one major development is the proof that two distinct integral representations of the $A_{n-1}$6 hyperbolic Calogero–Sutherland wave function are actually equivalent. One representation is Euler-type, built recursively from raising operators $A_{n-1}$7; the other is Mellin–Barnes type. The equality

$A_{n-1}$8

is established through the analysis of two families of Baxter operators, one acting in spatial variables and one in spectral variables, together with dominated-convergence arguments for the nonrelativistic limit from the hyperbolic Ruijsenaars system (Belousov et al., 26 Aug 2025). The same work identifies the resulting wave function with the renormalized Heckman–Opdam $A_{n-1}$9 hypergeometric function, fixes its normalization at the origin, and derives asymptotic formulas in the Weyl chamber (Belousov et al., 26 Aug 2025). This closes the gap between Baxter-operator constructions, bispectrality, and the standard Heckman–Opdam framework.

A separate line of work treats a Morse-deformed hyperbolic Calogero–Sutherland operator,

$$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$0

as a partial confluent limit of the $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$1 Heckman–Opdam hypergeometric function (Diejen et al., 2023). In that construction, the $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$2 Harish–Chandra series is translated and its couplings rescaled so that the $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$3 recurrence degenerates to the Morse-deformed hyperbolic Calogero–Sutherland recurrence (Diejen et al., 2023). The same paper shows that the resulting wave functions satisfy dual difference equations in the spectral variable and extend analytically in that variable, so their dependence on $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$4 is holomorphic rather than merely meromorphic (Diejen et al., 2023).

6. Spin and supersymmetric extensions

Supersymmetric hyperbolic Calogero–Sutherland systems are obtained by gauging matrix superfield models. For $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$5, gauging a $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$6 matrix system yields a model whose bosonic core is exactly the standard $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$7 hyperbolic Calogero–Sutherland system, while the fermionic sector consists of matrix fermions rather than a minimal particle-by-particle superpartner set (Fedoruk et al., 2019). The subsequent Hamiltonian analysis produces classical and quantum supercharges and an explicit Lax pair. A notable structural difference then appears: classically, the interaction terms in the supercharges are proportional to off-diagonal fermions, whereas quantum mechanically the supersymmetry generators admit a consistent invariant subsector without off-diagonal fermion operators (Fedoruk, 2019).

For $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$8, the bosonic sector becomes spinful. One gauged matrix construction uses $$\mathcal H = -\sum_{j=1}^n \partial_{x_j}^2 + \sum_{1\le j\ne k\le n}\frac{\pi^2 g(g-1)}{\sinh^2\!\pi(x_j-x_k)},$$9-spinor semi-dynamical variables and odd matrix fields so that the bosonic core is the $w(x)=|2\sinh \pi x|^g$0-spin hyperbolic Calogero–Sutherland system, with spin matrices

$w(x)=|2\sinh \pi x|^g$1

obeying the $w(x)=|2\sinh \pi x|^g$2 algebra (Fedoruk, 2020). The model admits explicit $w(x)=|2\sinh \pi x|^g$3 supercharges, a Lax representation, and an invariant reduction to the spinless hyperbolic Calogero–Sutherland system by setting one $w(x)=|2\sinh \pi x|^g$4 component of the spinor variables to zero (Fedoruk, 2020). More generally, Hamiltonian constructions starting from the classical $w(x)=|2\sinh \pi x|^g$5, $w(x)=|2\sinh \pi x|^g$6, $w(x)=|2\sinh \pi x|^g$7, and $w(x)=|2\sinh \pi x|^g$8 series produce $w(x)=|2\sinh \pi x|^g$9 and $\mathcal K(x)=(2\cosh \pi x)^{-g}$0 hyperbolic or trigonometric Calogero–Sutherland cousins whose bosonic cores are the standard systems; within that framework, the hyperbolic and trigonometric models appear to saturate at $\mathcal K(x)=(2\cosh \pi x)^{-g}$1 (2002.03929).

7. Related constructions and terminological boundaries

The family also appears under closely related nomenclature. In one formulation the same hyperbolic many-body system is treated as the hyperbolic Calogero–Moser system; the dynamics are derived from pole reduction of the semi-discrete KP equation, and the main result is a Lagrangian 1-form structure in both discrete and continuous time (Jairuk et al., 2016). The discrete Lagrangians involve $\mathcal K(x)=(2\cosh \pi x)^{-g}$2, the temporal Lax matrices satisfy compatibility relations yielding a discrete closure relation, and two successive continuum limits produce a continuous hierarchy together with the continuous closure relation $\mathcal K(x)=(2\cosh \pi x)^{-g}$3 (Jairuk et al., 2016). In this usage, the distinction between “Calogero–Moser” and “Calogero–Sutherland” is terminological rather than structural (Jairuk et al., 2016).

The adjective “hyperbolic” is also used in distinct senses across adjacent literature. One direction replaces line dynamics by motion on the hyperboloid $\mathcal K(x)=(2\cosh \pi x)^{-g}$4, with Hamiltonian

$\mathcal K(x)=(2\cosh \pi x)^{-g}$5

where $\mathcal K(x)=(2\cosh \pi x)^{-g}$6 is the angular Hamiltonian of a generalized rational Calogero model; the resulting hyperbolic Calogero-oscillator and Calogero-Coulomb systems on $\mathcal K(x)=(2\cosh \pi x)^{-g}$7 are maximally superintegrable (Hakobyan et al., 2014). A different construction, the “hyperbolic Kac–Moody Calogero model,” is based not on $\mathcal K(x)=(2\cosh \pi x)^{-g}$8 pair potentials on a line but on the real roots of the hyperbolic Lie algebra $\mathcal K(x)=(2\cosh \pi x)^{-g}$9, with a modular potential on the upper half-plane and Dunkl operators whose commutativity is obstructed by rank-2 hyperbolic subsystems (Lechtenfeld et al., 2022). That paper is careful not to claim nonintegrability outright; it states only that the usual Dunkl-based sufficient criterion fails (Lechtenfeld et al., 2022). The standard hyperbolic Calogero–Sutherland system is therefore best understood as the integrable $A_{n-1}$00-interaction family, while neighboring “hyperbolic” constructions may refer either to ambient hyperbolic geometry or to hyperbolic Weyl groups.