---
title: Lin-Shu-Kalnajs Modes in Stellar Discs
url: https://www.emergentmind.com/topics/lin-shu-kalnajs-modes
type: topic
---

# Lin-Shu-Kalnajs Modes in Stellar Discs

Searching arXiv for recent and foundational treatments of Lin–Shu–Kalnajs modes and related stellar-disc dynamics.
Relevant arXiv sources already identified in the provided corpus span foundational lecture notes, shearing-sheet analyses, unified linear-response theory, and modern matrix methods: [1202.3403], [1906.11696], [2312.02618], [2507.16950], and [2606.31523].
Lin–Shu–Kalnajs (LSK) modes are normal-mode solutions of the linearized collisionless Boltzmann–Poisson system for self-gravitating stellar discs, usually formulated for razor-thin, nearly axisymmetric equilibria and most commonly analyzed in the tightly wound, or WKB, limit. In that limit, spiral perturbations are represented by Fourier components of the form $e^{i(kR+m\phi-\omega t)}$, and the disc response is governed by the Lin–Shu–Kalnajs dispersion relation with its characteristic reduction factor $\mathcal{F}(s,\chi)$, which encodes the effect of finite stellar random motions. The LSK framework occupies a central place in spiral-density-wave theory because it provides both a local dispersion relation and a route to global normal-mode calculations; however, it is also now understood to be only one asymptotic sector of a broader linear-response problem that includes winding, swing amplification, corotation physics, and resonant absorption [1202.3403].

## 1. Normal-mode formulation in collisionless stellar discs

The formal starting point is the linearized collisionless Boltzmann equation coupled to Poisson’s equation for an equilibrium axisymmetric disc with distribution function $f_0(\mathbf{J})$ in action–angle variables $(\mathbf{J},\boldsymbol\theta)$:
\[
\frac{\partial f_1}{\partial t}+\{f_1,H_0\}=-\{f_0,\Phi_1\},
\qquad
\nabla^2\Phi_1=4\pi G\int f_1\,d^3\mathbf{v}.
\]
Seeking normal modes,
\[
f_1(\mathbf{J},\boldsymbol\theta,t)
=\sum_{\mathbf n}\hat f_{\mathbf n}(\mathbf J)\,
e^{\,i(\mathbf n\cdot\boldsymbol\theta-\omega t)},
\qquad
\Phi_1(\mathbf x,t)=\hat\Phi(\mathbf x)e^{-i\omega t},
\]
leads to an integral eigenvalue problem for the complex frequency $\omega$ and the mode shape $\hat\Phi$ [1202.3403].

A closely related matrix formulation is obtained by expanding the perturbation in angle harmonics. In the four-dimensional phase space $(J_1,J_2;w_1,w_2)$, with $J_2=L$, the linearized equation for a perturbation $F_1\propto e^{-\omega t+m\theta}$ becomes
\[
\frac{\partial F_1}{\partial t}
+\Omega_1\frac{\partial F_1}{\partial w_1}
+\Omega_2\frac{\partial F_1}{\partial w_2}
=
\frac{\partial\Phi_1}{\partial w_1}\frac{\partial F_0}{\partial J_1}
+\frac{\partial\Phi_1}{\partial w_2}\frac{\partial F_0}{\partial J_2},
\]
and, after harmonic decomposition, yields
\[
\bigl[\omega-l\Omega_1(J_1,J_2)-m\Omega_2(J_1,J_2)\bigr]\,
{\cal F}_l(J_1,J_2)
=
-F'_{0,l}(J_1,J_2)\,\Phi_l(J_1,J_2),
\]
with
\[
F'_{0,l}(J_1,J_2)\equiv
l\,\frac{\partial F_0}{\partial J_1}
+
m\,\frac{\partial F_0}{\partial J_2}.
\]
Together with the Poisson integral, this produces a standard eigenvalue problem
\[
\omega\,\mathbf{x}=\mathbf{A}\,\mathbf{x},
\qquad
\det|\mathbf{A}-\omega\mathbf{I}|=0,
\]
which is the modern matrix form of the global normal-mode calculation [2606.31523].

Within this formal structure, LSK modes are not defined merely by the existence of an eigenvalue problem, but by the additional asymptotic assumptions that reduce the global response to a local spiral-wave dispersion relation.

## 2. WKB reduction and the Lin–Shu–Kalnajs dispersion relation

The standard derivation imposes four assumptions: the spiral wave is tightly wound, so that $|kR|\gg1$; the unperturbed disc is razor-thin with surface density $\Sigma_0(R)$; the epicycle approximation holds; and attention is restricted to a single Fourier component proportional to $e^{i(kR+m\phi-\omega t)}$ [1202.3403].

In this regime, stellar orbits are decomposed into a guiding centre at $R_g$ rotating at $\Omega(R_g)$ plus small radial oscillations of frequency $\kappa(R_g)$. After Fourier–Laplace transformation in $(\phi,t)$ and averaging over epicyclic phases, the perturbed surface density can be written in terms of the potential amplitude through a reduction factor:
\[
s\equiv\frac{\omega-m\Omega}{\kappa},
\qquad
\chi\equiv\frac{k^2\sigma_R^2}{\kappa^2},
\]
\[
\mathcal{F}(s,\chi)\simeq
e^{-\chi}\sum_{\ell=-\infty}^{\infty}
\frac{I_\ell(\chi)}{1-s^2/\ell^2},
\]
or, equivalently,
\[
\mathcal{F}(s,\chi)
=
\frac{1-s^2}{\sin(\pi s)}
\int_0^\pi d\tau\;
e^{-\chi(1+\cos\tau)}\,\sin(s\tau)\sin\tau.
\]
For a razor-thin disc, Poisson’s equation gives
\[
\hat\Phi=-\frac{2\pi G}{|k|}\,\delta\Sigma,
\]
and combination of response and self-gravity yields the LSK dispersion relation
\[
1-\frac{2\pi G\,\Sigma_0\,|k|}{\kappa^2-(\omega-m\Omega)^2}\,
\mathcal{F}(s,\chi)=0,
\]
or
\[
(\omega-m\Omega)^2
=
\kappa^2-2\pi G\,\Sigma_0\,|k|\,\mathcal{F}(s,\chi)
\]
[1202.3403].

A closely related representation, derived from the Laplace-transformed Volterra kernel in a razor-thin stellar disc, is
\[
1
=
\frac{2\pi G\,|k|}{\kappa^2}\,\Sigma_0\;
e^{-\chi}
\sum_{n=-\infty}^{\infty}
\frac{n^2}{n^2-s^2}\,\frac{I_n(\chi)}{\chi},
\]
with
\[
k\equiv\sqrt{k_R^2+\frac{m^2}{R^2}}.
\]
This formulation makes explicit the equivalence between the local WKB approach and the more general linear-response theory [2507.16950].

The physical content of $\mathcal{F}(s,\chi)$ is the suppression of self-gravitating response by finite radial velocity dispersion. In the cold-disc limit, long waves approach $\mathcal{F}\to1$; at large $\chi$, short waves are stabilized as $\mathcal{F}\to0$ [1202.3403].

## 3. Propagation, resonances, and discrete global modes

For real $\omega$, the LSK dispersion relation defines neutrally stable wave propagation and permits the calculation of the radial group velocity
\[
v_g=\frac{\partial\omega}{\partial k}
=
\frac{\bigl[\pi G\Sigma_0\,\mathrm{sgn}(k)\,\mathcal{F}
+\pi G\Sigma_0\,|k|\,\partial_k\mathcal{F}\bigr]}
{\omega-m\Omega}
\]
in the reduction-factor formulation [1202.3403]. In the axisymmetric-sheet WKB limit, the corresponding expression is
\[
v_g
=
\frac{\sigma_R^2\,k-\pi G\Sigma\,\mathrm{sgn}(k)}{\omega},
\]
which is Toomre’s group-velocity result [1906.11696].

The local dispersion relation implies a resonance structure. The Lindblad resonances occur when
\[
\omega-m\Omega=\pm\kappa,
\]
and corotation occurs when
\[
\omega=m\Omega.
\]
Only between the Lindblad resonances can real waves propagate; outside that interval the local solutions are evanescent [1202.3403]. More generally, the reduction factor controls the growth or damping rate near resonances, and allowing a small imaginary part in the eigenfrequency leads to a Landau prescription for damping or growth in the lower half-$\omega$ plane [2507.16950].

To convert the local WKB picture into discrete global modes, one imposes boundary conditions at turning radii $R_{\rm in}$ and $R_{\rm out}$ where $k_R^2(R)$ changes sign. The quantization condition is
\[
\int_{R_{\rm in}(\omega)}^{R_{\rm out}(\omega)} k_R(R,\omega)\,dR
=
\pi\Bigl(n+\tfrac12\Bigr),
\qquad n\in\mathbb Z_{\ge0},
\]
with Maslov phase shifts of $\pi/2$ introduced when matching across corotation and Lindblad turning points [2507.16950]. The resulting radial dependence is
\[
\delta\Phi_n(R)\propto
\frac{1}{\sqrt{k_R(R)}}\,
\exp\Bigl[i\!\int^R k_R(R')\,dR'\Bigr].
\]

This WKB quantization gives a discrete sequence of pattern speeds $\Omega_p=\omega/m$ and mode numbers $n$, but its validity depends on the same tight-winding and slowly varying assumptions that underlie the local dispersion relation.

## 4. Winding, swing amplification, and the limits of the LSK picture

The main limitation of the LSK framework is that it assumes a single fixed $(k,\omega)$ for a normal mode, whereas differential rotation shears non-axisymmetric disturbances. In a differentially rotating disc,
\[
k(t)=k_0+m\,\frac{d\Omega}{dR}\,t,
\]
so an initially open spiral tightens with time and eventually leaves the WKB regime or is damped near resonances [1202.3403]. In the shearing sheet, the local radial wavenumber evolves as
\[
k_x(t)=2Ak_y t,
\]
with $A$ Oort’s first constant [1906.11696].

This winding dynamics is the basis of swing amplification. In the shearing sheet, a leading disturbance can be transiently amplified by factors of ten or more before it winds up, and this explains why cool stellar discs are responsive systems that amplify ambient noise [1202.3403]. The shearing-wave framework extends this by treating tight-winding Lin–Shu density waves and swing-amplified material patterns as limits of a common non-steady problem, with growth driven by the term $2Akm/R$ and interpreted either as swing amplification or as the Lynden-Bell–Kalnajs donkey effect at corotation [2312.02618].

The shearing-sheet analysis also alters the interpretation of wave propagation near corotation. Rather than a single packet bounded by a forbidden zone, two wavepackets emerge inside corotation: one each side of the inner Lindblad resonance [1906.11696]. Disturbances do not avoid an annulus around corotation, and as wavepackets approach Lindblad resonances they are subject to Landau absorption, resonant trapping, or partial reflection from the velocity-space distortion left by previous waves [1906.11696].

For this reason, the LSK dispersion relation and even global normal-mode calculations provide a very incomplete understanding of the dynamics of stellar discs [1202.3403]. A plausible implication is that LSK modes are best regarded as asymptotic organizing structures rather than exhaustive descriptions of long-lived spiral morphology.

## 5. Relation to unified linear-response theory

Modern work places LSK modes inside a broader response framework for razor-thin stellar discs. In this formulation, the perturbed potential and surface density are expanded in a biorthogonal basis, and the modal amplitudes satisfy a Volterra integral equation,
\[
B^p(t)=B^p_{\rm kin}(t)
+\int_0^t dt'\sum_{p'}\mathcal M^{pp'}(t-t')
\bigl[B^{p'}(t')+B^{p'}_{\rm ext}(t')\bigr],
\]
with kernel
\[
\mathcal M^{pp'}(\tau)
=
-\frac{(2\pi)^2}{E}\sum_{\mathbf n}\int dJ\;
\frac{\mathbf n\cdot\partial f_0/\partial J}{\mathbf n\cdot\Omega(J)}\,
e^{-i\mathbf n\cdot\Omega(J)\tau}\,
[\phi^{(p)}_{\mathbf n}(J)]^*\phi^{(p')}_{\mathbf n}(J)
\]
[2507.16950].

Within this theory, the Volterra, Landau, and van Kampen representations are equivalent, and several classic results emerge as limiting cases: Lindblad–Kalnajs density waves, swing amplification, Lin–Shu–Kalnajs modes, and groove instabilities [2507.16950]. The short-wavelength asymptotic expansion recovers the WKB kernel and hence the LSK dispersion relation precisely in the limit $|k_R|\gg m/R$ [2507.16950].

A complementary unification is provided by the shearing-wave cubic dispersion relation
\[
\bigl(\omega_e-m\Omega\bigr)^3
=
\bigl(\omega_e-m\Omega\bigr)\,\omega_0^2
+i\,\omega_c^3,
\]
which, in the steady-mode limit with $dk/dt=0$ and the imaginary term neglected, reduces to
\[
(\omega-m\Omega)^2=\omega_0^2
\]
and thereby recovers the Lin–Shu–Kalnajs dispersion relation in 3D [2312.02618]. The significance of these unified treatments is not that they replace LSK modes, but that they delimit the regime in which the classical normal-mode language is asymptotically valid.

## 6. Modern computation of global LSK-type modes

The contemporary numerical treatment of collisionless-disc modes generalizes the classical action–angle matrix method. For distribution functions with a sharp edge at $L=0$,
\[
F_0(J,L)=\mathfrak{F}(J,L)\,\Theta(L),
\qquad
\frac{\partial F_0}{\partial L}
=
\frac{\partial\mathfrak{F}}{\partial L}\,\Theta(L)
+\mathfrak{F}(J,0)\,\delta(L),
\]
the generalized linear-matrix equation acquires a boundary-integral term but no increase in matrix size [2606.31523]. This is important for unidirectional discs with sharp angular-momentum boundaries.

The same formulation naturally incorporates gravitational softening by replacing the Poisson kernel denominator $|{\bf x}-{\bf x}'|$ with
\[
\sqrt{|{\bf x}-{\bf x}'|^2+\beta^2},
\]
so the softening length $\beta$ enters directly through the matrix elements [2606.31523]. This enables direct comparison with softened $N$-body calculations.

For the Kuzmin–Toomre disk with Kalnajs DF index $m_K=6$, the five most unstable $m=2$ modes were computed with a sharp-edge calculation using grid $N_R=61$, $N_L=21$, and $\ell_{\max}=7$, giving
\[
(\Omega_p,\gamma)=
(0.801,0.650),\,
(0.384,0.159),\,
(0.291,0.087),\,
(0.236,0.052),\,
(0.198,0.037).
\]
A convergence study in $\ell_{\max}$ yielded
\[
\ell_{\max}=7:\;(\Omega_p,\gamma)=(0.801,0.650),\quad
\ell_{\max}=15:\;(0.814,0.655),\quad
\ell_{\max}\to\infty:\;(0.817,0.657),
\]
while comparison to an independent log-spiral method gave the “exact” values
\[
(\Omega_p,\gamma)=(0.813,0.653),
\]
so the extrapolated matrix method agreed to $\sim0.003$ in each component [2606.31523].

These calculations illustrate that the classical LSK formalism remains operational as a computational framework, even though its physical interpretation must now be embedded in a more general theory of transient, resonant, and shearing spiral structure.

Source: https://www.emergentmind.com/topics/lin-shu-kalnajs-modes