---
title: Calogero-Sutherland Model Overview
url: https://www.emergentmind.com/topics/calogero-sutherland-model
type: topic
---

# Calogero-Sutherland Model Overview

The Calogero–Sutherland model is a class of integrable one-dimensional quantum many-body systems with inverse-square interactions, realized most prominently in a rational form on the line and a trigonometric form on the circle. In a standard circle formulation with \(N\) particles on a circumference \(L\), the Hamiltonian is
\[
H_{\rm CS}
\;=\;
\sum_{j=1}^N\!-\frac{\partial^2}{\partial x_j^2}
\;+\;
g\sum_{i<j}\frac{(\pi/L)^2}{\sin^2\!\bigl[\pi(x_i-x_j)/L\bigr]},
\qquad
g=\lambda(\lambda-1),
\]
while a line version uses the pair potential \(1/(x_i-x_j)^2\) directly. Across these realizations, the model is distinguished by exact Jastrow-type ground states, solvability in terms of Jack polynomials and related algebraic structures, effective descriptions by extended conformal field theory and quantum hydrodynamics, and a large family of deformations including finite-range, spin, supersymmetric, PT-symmetric, and two-species variants [2501.03762] [2507.17397] [2404.10011].

## 1. Canonical formulations and coupling structures

The trigonometric Calogero–Sutherland Hamiltonian on the circle describes particles interacting through the periodic inverse-square potential \(\sin^{-2}\!\bigl[\pi(x_i-x_j)/L\bigr]\). In one formulation the model is written for \(N\) spinless fermions on a circle of circumference \(L\), with coupling \(g=\lambda(\lambda-1)\) and thermodynamic density \(n_0=N/L\) [2501.03762]. A line formulation replaces the chord distance by the direct separation,
\[
H \;=\; -\frac{\hbar^2}{2m}\sum_{i=1}^N\frac{\partial^2}{\partial x_i^2}
\;+\;\frac{\hbar^2\,\lambda(\lambda-1)}{m}\sum_{1\le i<j\le N}\frac{1}{(x_i - x_j)^2},
\]
and on a ring one may write the interaction in terms of
\[
\mathrm{crd}(x_{ij})
=\frac{L}{\pi}\,\sin\!\Bigl(\frac{\pi(x_i-x_j)}{L}\Bigr).
\]
In this language the dimensionless parameter \(\lambda\) governs the physics and is related to the Dyson index by \(\beta=2\lambda\) in random-matrix theory [2507.17397].

A broader formulation uses root systems. For a positive-root set \(R_+\), the bosonic hyperbolic/trigonometric Calogero–Sutherland Hamiltonian is
\[
H_{\rm bos}
=
\frac12\sum_{i=1}^n p_i^2
+
\sum_{\alpha\in R_+} g_\alpha^2\,f^2(\alpha\!\cdot\! x),
\]
with
\[
f(x)=
\begin{cases}
\dfrac1{x}, & \text{rational},\\[4pt]
\dfrac1{\sinh x}, & \text{hyperbolic},\\[4pt]
\dfrac1{\sin x}, & \text{trigonometric}.
\end{cases}
\]
For \(A_{n-1}\) this becomes
\[
H_{A_{n-1}}
=
\frac12\sum_{i=1}^n p_i^2
+
g^2\sum_{1\le i<j\le n} f^2(x_i-x_j),
\]
while \(B_n\), \(C_n\), and \(D_n\) add \(x_i+x_j\) terms and, where applicable, single-particle \(f^2(x_i)\) terms with a second coupling \(g'\) [2002.03929].

This multiplicity of forms is not merely a change of notation. It organizes the model into rational, trigonometric, and hyperbolic branches, and into root-system families whose algebraic structures differ while preserving the inverse-square character of the interactions.

## 2. Exact ground states, spectra, and Jack-polynomial structure

One of the defining properties of the model is the existence of closed-form ground states. On the infinite line,
\[
\Psi_0(x_1,\dots,x_N)
=
C_N(\lambda)\;\prod_{1\le i<j\le N}\bigl|x_i-x_j\bigr|^\lambda,
\]
and on a ring of length \(L\),
\[
\Psi_0
=
\widetilde C_N(\lambda,L)\;\prod_{i<j}\Bigl|\sin\frac{\pi(x_i-x_j)}{L}\Bigr|^\lambda.
\]
For the ring model the ground-state energy per particle is
\[
\frac{E_0}{N}
=
\frac{\hbar^2\pi^2\,\lambda^2}{6m}\,n^2,
\qquad
n=\frac{N}{L}.
\]
These formulas also exhibit the model’s direct connection to \(\beta\)-ensemble joint densities through \(|\Psi_0|^2\propto \prod_{i<j}|x_i-x_j|^{2\lambda}\) [2507.17397].

In the trigonometric \(A_{n-1}\) system, eigenfunctions may be written as
\[
\Psi_{\{m_i\}}(x)
=
\Delta(x)^g\,
P_{\{m_i\}}^{(1/g)}\!\bigl(e^{2ix_1},\dots,e^{2ix_n}\bigr)\,
e^{i\sum_k m_k x_k},
\]
where
\[
\Delta(x)=\prod_{i<j}\sin(x_i-x_j),
\]
and \(P_{\{m\}}^{(\alpha)}\) are symmetric Jack polynomials with \(\alpha=1/g\). The corresponding spectrum is
\[
E_{\{m_i\}}
=
\frac12\sum_{i=1}^n m_i^2,
\qquad
m_1<m_2<\cdots<m_n,
\]
subject to an overall shift depending on \(g\) [2002.03929].

The Jack-polynomial formulation becomes especially explicit after factoring out the ground state. For \(L=\pi\) and \(\beta=k^2\), one defines
\[
H_J
=
\sum_{i=1}^N\bigl(z_i\partial_{z_i}\bigr)^2
+
\beta\sum_{i<j}\frac{z_i+z_j}{z_i-z_j}\,(z_i\partial_{z_i}-z_j\partial_{z_j}),
\qquad
z_i=e^{2ix_i},
\]
and the Jack polynomial \(J_\lambda(z_1,\dots,z_N)\) satisfies
\[
H_J\,J_\lambda = E_\lambda\,J_\lambda,
\qquad
E_\lambda
=
\sum_{i=1}^N\lambda_i^2
+
\beta\sum_{i=1}^N(N-2i+1)\lambda_i.
\]
A hidden Virasoro structure permits a recursive construction of Jack states from singular vectors, with rectangular partitions playing a fundamental role in the skew hierarchy [1107.4234].

A recurrent theme in later work is that the Jack description is not only a solution method but also a bridge to conformal blocks, Yangian symmetry, and fractional-quantum-Hall wavefunctions.

## 3. Effective field theory, \(W_{1+\infty}\), and hydrodynamics

In the thermodynamic limit \(N,L\to\infty\) at fixed density \(n_0=N/L\), the low-energy theory reorganizes into chiral and anti-chiral sectors. One formulation identifies the effective theory as an extended conformal field theory with
\[
W_{1+\infty}\otimes\overline{W}_{1+\infty},
\qquad
(c,\bar c)=(1,1),
\]
and plasma velocity
\[
v=\pi\lambda n_0.
\]
In this second-quantized description the Hamiltonian becomes a bilinear combination of \(W_{1+\infty}\) generators, and the chiral density operator is proportional to the spin-one current [2501.03762].

A closely related effective Hamiltonian is written as
\[
H_{\rm CS}
=
\sum_{\ell\ge0} \epsilon_\ell
\bigl[
W^{(\ell+1)}_0(R)+W^{(\ell+1)}_0(L)
\bigr],
\]
with \(W^{(0)}_m=a_m\), \(W^{(1)}_m=\frac12\sum_r :a_r a_{m-r}:\), and \(W^{(2)}_m\) cubic in boson modes. Highest-weight states are labelled by the particle-number shift \(\Delta N\) and the pumped chiral charge \(\Delta D\), with
\[
Q=\sqrt{\lambda}\,\Delta N + \frac{\Delta D}{\sqrt{\lambda}},
\qquad
\bar Q=\sqrt{\lambda}\,\Delta N - \frac{\Delta D}{\sqrt{\lambda}},
\]
and exact finite-size energy and momentum
\[
E=(2\pi n_0)^2\Bigl\{Q^2+\sum_j k_j\Bigr\}+O(1/N),
\qquad
P=(2\pi n_0)\,Q
\]
[2305.06744].

The chiral density dynamics is governed by the quantum Benjamin–Ono equation. In one form, for the right-moving density field \(u(x)\),
\[
\partial_t u + v_F\,\partial_x u + (2\pi\lambda)^{-1}\partial_x(u^2)
+ \frac{1}{2\pi}H[\partial_x^2 u] = 0,
\]
where \(H[f](x)\equiv \mathrm{PV}\int (dy/\pi)\,f(y)/(x-y)\) is the Hilbert transform [2305.06744]. An operator-level derivation shows that this result does not rely on a semiclassical assumption. In that approach, the chiral density satisfies
\[
\partial_t\,\rho(x,t)
=
\frac{\pi\,\xi^2}{2m}\;\partial_x\!\bigl[\rho^2(x,t)\bigr]
+
\frac{(\xi-1)\,\pi\,\xi}{2m}\;\partial_x^2
\bigl(\mathcal H\{\rho\}(x,t)\bigr),
\]
which the source identifies as precisely the quantum Benjamin–Ono equation [2404.10011].

At the classical hydrodynamic level, the continuum limit yields a continuity equation, an Euler equation with a nonlocal enthalpy, and a bidirectional Benjamin–Ono equation. The latter appears as a real reduction of the modified KP hierarchy, and a chiral reduction gives a Chiral Non-linear Equation whose large-density degeneration is the conventional Benjamin–Ono equation [0810.5327]. A common simplification is to identify the thermodynamic Calogero–Sutherland theory with an ordinary free compactified boson; the sources instead emphasize that the Hilbert spaces may be isomorphic while the time evolution is deformed by the interaction-dependent effective Hamiltonian [2501.03762].

## 4. Dynamic structure factor and correlation physics

The dynamic structure factor provides one of the sharpest diagnostics of the model’s integrable dynamics. In the chiral effective theory and in the repulsive regime \(\lambda>1\), the zero-temperature structure factor takes the form
\[
S_{\rm CS}(q,\omega)
=
\mathcal{A}(\lambda)\;q\;
\delta\!\Bigl(\omega-qv\bigl[1-\tfrac{\lambda-1}{2\lambda\,n_0L}q^2\bigr]\Bigr),
\]
with \(v=\pi\lambda n_0\). In the strict thermodynamic limit \(L\to\infty\), the interaction-induced shift vanishes and the resonance sits at \(\omega=qv\). The corresponding free compactified boson has
\[
S_0(q,\omega)=\pi\,q\,\delta(\omega-vq),
\]
so the Hilbert spaces coincide while the resonance position differs away from the strict thermodynamic limit [2501.03762].

In the first-quantized treatment, the same response is described by a narrow rectangle,
\[
\omega\in[\omega_-,\omega_+],
\qquad
\omega_\pm=vq\pm\frac{q^2}{2m},
\]
with approximate form
\[
S_{\rm 1st}(q,\omega)
\approx
\frac{\Theta(\omega-\omega_-)\,\Theta(\omega_+-\omega)}{\omega_+-\omega_-}\,q.
\]
In the overlap regime \(N\gg1\), \(\lambda>1\), this rectangle collapses to the same delta peak at \(\omega=vq\), which the source presents as a precise matching between first- and second-quantized descriptions [2501.03762].

A broader account of dynamic correlations goes beyond the conventional Luttinger-liquid regime by combining sum rules with Monte Carlo sampling of the exact phase-space representation. At small \(q\), the excitation is always phononic,
\[
\omega(q)\simeq c\,q,
\qquad
c=\lambda v_F,
\qquad
v_F=\frac{\pi\hbar n}{m}.
\]
For weak coupling \(\lambda\ll1\), only the upper branch \(\omega_+(q)\) is populated and a Bogoliubov-type spectrum emerges. At \(\lambda=1\), the Tonks–Girardeau or ideal Fermi limit yields
\[
S(q,\omega)=
\begin{cases}
\mathrm{const.}, & \omega_-(q)<\omega<\omega_+(q),\\
0, & \text{otherwise,}
\end{cases}
\qquad
\omega_\pm(q)=c\,q\pm\frac{q^2}{2m}.
\]
For strong coupling \(\lambda>1\), the support again lies in \(\omega_-(q)<\omega<\omega_+(q)\), but the dominant behavior shifts to a singular lower edge and, in the \(\lambda\to\infty\) limit, to a quasi-crystalline Brillouin-zone structure [2507.17397].

The same source makes explicit that the exact ground state reproduces the universal Luttinger-liquid form with
\[
K=\frac1\lambda,
\qquad
c=\lambda v_F.
\]
This suggests an exact solvable realization of the full Luttinger-liquid universality class rather than a merely asymptotic low-energy approximation [2507.17397].

## 5. Finite-range, truncated, and two-species generalizations

A finite-range generalization on the line introduces particles in a harmonic trap with inverse-square two-body and three-body interactions truncated to \(r\) neighbors. The Hamiltonian is
\[
H
=
\sum_{i=1}^N\Bigl[-\tfrac12\,\partial_{x_i}^2 + \tfrac12\,\omega^2 x_i^2\Bigr]
+
\lambda(\lambda-1)\sum_{\substack{i<j\\|i-j|\le r}}\frac{1}{(x_i-x_j)^2}
+
\lambda^2
\sum_{\substack{i<j<k\\|i-j|\le r,\;|j-k|\le r}}
\frac{(x_j-x_i)(x_j-x_k)}{(x_j-x_i)^2(x_j-x_k)^2}.
\]
Its exact ground state is
\[
\Psi_0(x_1,\dots,x_N)
=
C_{N,\lambda,r}^{-1/2}
\exp\!\Bigl[-\tfrac{\omega}{2}\sum_{i=1}^N x_i^2\Bigr]
\prod_{\substack{i<j\\|i-j|\le r}} |x_i-x_j|^\lambda,
\]
with ground-state energy
\[
E_0(N,\lambda,r)
=
\frac{\omega}{2}\Bigl[N+\lambda\,r\,(2N-r-1)\Bigr].
\]
The limits \(r=N-1\), \(r=1\), and \(r=-1,\lambda=1\) recover, respectively, the full Calogero–Sutherland model, the Jain–Khare nearest- plus next-nearest-neighbor model, and the Tonks–Girardeau gas [1606.00895].

A circle version truncates the inverse-square interaction to an \(r\)-neighbor shell on each side and introduces an attractive three-body term. Its exact ground state is
\[
\psi_0(x_1,\dots,x_N)
=
{\mathcal N}_0\;
\prod_{\substack{i<j\\|i-j|\le r}}
\sin^{\beta}\!\bigl(\theta_{ij}\bigr),
\qquad
\theta_{ij}=\frac{\pi}{L}(x_i-x_j),
\]
with an exact ground-state energy that reduces to the full-range Sutherland result for \(r\ge c\) and to the Jain–Khare case at \(r=1\). A part of the excitation spectrum can be constructed in terms of elementary symmetric functions \(e_m(z)\) [1609.07928].

The finite-range models alter correlation properties in a controlled way. Numerically, increasing either \(\lambda>0\) or the interaction range \(r\) enhances spatial antibunching, broadens the density cloud, suppresses off-diagonal long-range order in the one-body reduced density matrix, and broadens the momentum distribution tails [1606.00895]. The source explicitly interprets truncation as a screening of the inverse-square force; although full long-range integrability is broken, the model remains quasi-exactly solvable because the ground state and a large class of collective excitations are still known in closed form [1606.00895].

A distinct extension is the two-species or deformed Calogero–Sutherland Hamiltonian,
\[
H_{N,M}(x,y;g)
=
-\sum_{j=1}^N\partial_{x_j}^2
-g\sum_{k=1}^M\partial_{y_k}^2
+\cdots,
\]
whose exact eigenfunctions are built from super-Jack polynomials \(P_\lambda^{(N,M)}(z,w;g)\). In the interpretation given by the source, the \(x\)-particles and \(y\)-particles may be viewed as electrons and quasi-holes [1603.06157].

Finally, a rational extension of the truncated model adds a radial term \(V_{\rm new}(\rho)\) while leaving the energy eigenvalues unchanged:
\[
E_{n,s}^{(\rm ext)}
=
\omega\bigl(2n+s+1+\lambda(2N-r-1)\bigr).
\]
The eigenfunctions are modified from classical Laguerre to exceptional \(X_1\) or, more generally, \(X_m\) Laguerre polynomials. The cases \(m=0\), \(r=1\), and \(r=N-1\) recover, respectively, the conventional truncated model, the extended Jain–Khare model, and the extended full-range Calogero–Sutherland model [1807.05163].

## 6. Spin, supersymmetric, and PT-symmetric variants

Spin degrees of freedom admit several nontrivial integrable realizations. A classical hyperbolic model with two interacting spin variables \(S\) and \(T\) has Hamiltonian
\[
H
=
\frac12\sum_{i=1}^N v_i^2
-
\sum_{1\le i<j\le N}
\frac{
S_{ij}^2 + T_{ij}^2\cosh(u_i-u_j) -2S_{ij}T_{ij}
}{
\sinh^2(u_i-u_j)
},
\]
with \(S_{ij}=-S_{ji}\), \(T_{ij}=-T_{ji}\). The source constructs a Lax pair, a classical \(r\)-matrix with spectral parameter, and a Hitchin-system realization, and proves complete integrability [1706.08793].

At the quantum level, the trigonometric spin-Calogero–Sutherland model can be analyzed by Bethe ansatz using Yangian symmetry. In additive coordinates,
\[
\widetilde H
=
-\tfrac12\sum_{i=1}^N\partial_{x_i}^2
+
\beta\sum_{1\le i<j\le N}
\frac{\beta\mp P_{ij}}{4\,\sin^2\!\bigl(\tfrac{x_i-x_j}{2}\bigr)},
\]
with \(P_{ij}\) the spin permutation operator. From a transfer matrix with diagonal twist one obtains a Bethe algebra that refines the usual commuting Hamiltonians and yields a new eigenbasis generalizing the Yangian Gelfand–Tsetlin basis of Takemura–Uglov [2308.16865].

Supersymmetric extensions are likewise extensive. For root systems \(A_n\), \(B_n\), \(C_n\), and \(D_n\), \({\cal N}=2\) and \({\cal N}=4\) hyperbolic/trigonometric supersymmetric models can be constructed so that the bosonic core is the standard Calogero–Sutherland system. In the \({\cal N}=4\) construction, fermionic bilinears \(I_{ij}\) enter the supercharges and Hamiltonian, and the supercharges satisfy
\[
\{Q^a,Q^b\}=0,
\qquad
\{\bar Q_a,\bar Q_b\}=0,
\qquad
\{Q^a,\bar Q_b\}=-2i\,\delta^a{}_b\,H.
\]
The source emphasizes that the full supersymmetric Hamiltonians are determined explicitly, while the bosonic limit preserves the standard integrable model [2002.03929].

A gauged-matrix approach yields \({\cal N}=2\) and \({\cal N}=4\) hyperbolic Calogero–Sutherland systems whose bosonic sectors reduce to
\[
H_{\rm bose}^{(2)}
=
\frac12\sum_{a=1}^n p_a^2
+
g^2\sum_{1\le a<b\le n}\frac{1}{\sinh^2(q_a-q_b)},
\]
and to a \(U(2)\)-spin hyperbolic model in the \({\cal N}=4\) case [1902.08023]. The \({\cal N}=4\) \(U(2)\)-spin model admits explicit supercharges, a Lax representation, and a consistent reduction to the spinless hyperbolic Calogero–Sutherland model [2007.11424]. The \({\cal N}=2\) hyperbolic model has both classical and quantum supercharges; the source stresses that the quantum supercharges, unlike the classical ones, can be restricted to an invariant subsector without off-diagonal fermion operators [1910.07348].

A non-Hermitian but PT-symmetric deformation acts by a constant complex shift \(q_i\to q_i+i\epsilon_i\), equivalently \(x_i\to x_i e^{-2\epsilon_i}\), and regularizes all two-body singularities while preserving PT symmetry. The deformed Hamiltonian is isospectral to the undeformed one, but the regularization adds an infinite tower of previously non-normalizable states. For integer coupling, additional degeneracy appears and a nonlinear conserved supersymmetry charge enlarges the ring of Liouville charges, while the Dunkl-operator integrability structure is maintained [1903.06481].

## 7. Conformal blocks, fractional quantum Hall structures, and universality

The Calogero–Sutherland model has a deep conformal-field-theoretic realization. Conformal blocks with only second-order degenerate Virasoro or \(W\)-algebra fields satisfy second-order differential equations that, after \(u(1)\) dressing and multiplication by appropriate Jastrow factors, become the Schrödinger equation of the trigonometric Calogero–Sutherland Hamiltonian,
\[
H_{\rm CS}^g\,\Psi(z_1,\dots,z_N)=E\,\Psi(z_1,\dots,z_N).
\]
The same framework exhibits a duality \(g\leftrightarrow 1/g\), organizes excited states by two partitions in the Virasoro case and by \(k\) partitions in \(W A_{k-1}\), and maps the Calogero–Sutherland integrals of motion to the operators \(I_r\) acting in \(u(1)\otimes\)Virasoro or \(W\)-algebra modules [1110.1101].

At special rational couplings
\[
g=\frac{k+r}{k+1},
\]
the relevant conformal blocks reduce to Jack polynomials with \((k,r)\)-clustering. The cited cases include the Laughlin state for \(k=1\), \(r=m-1\), the Moore–Read state for \(k=2\), \(r=2\), and the Read–Rezayi series for general \(k\), \(r=2\) [1110.1101]. This is a direct instance in which Calogero–Sutherland eigenfunctions become fractional-quantum-Hall trial states rather than merely analogous structures.

The deformed two-species model makes this relation more explicit. Its collective-field description is formulated in terms of a chiral boson CFT with anyonic vertex operators, and the source identifies the two species naturally with electrons and quasi-hole excitations in Wen’s effective field theory of the fractional quantum Hall effect. Super-Jack polynomials then provide simple explicit formulas for an orthonormal CFT basis originally proposed in that context [1603.06157].

The model’s broader significance extends beyond conformal theory. One source describes the effective field theory of the Calogero–Sutherland model as a universality class of one-dimensional quantum hydrodynamic fluids and develops a ribbon geometry in which the current profile becomes
\[
J_x(y)=J_0\,[1-(2y/L_y)^2],
\]
a Poiseuille-type parabola that is nonetheless dissipationless because all one-dimensional layers drift with the same velocity [2305.06744]. Another source places the model simultaneously in cold-atom waveguides, disordered metallic grains, and random-matrix theory, emphasizing that its exact ground state, sound velocity, Luttinger parameter, excitation spectrum, and dynamic correlations can be obtained analytically or by stochastic sampling of an exact phase-space integral [2507.17397].

Taken together, these developments establish the Calogero–Sutherland model not simply as a solvable inverse-square many-body problem, but as a structural node linking Jack and super-Jack polynomials, \(W_{1+\infty}\) symmetry, hydrodynamic Benjamin–Ono dynamics, Yangian and Dunkl-operator integrability, and fractional-quantum-Hall conformal blocks.

Source: https://www.emergentmind.com/topics/calogero-sutherland-model