Violent Relaxation: How a Collapsing Galaxy Reaches Equilibrium in One Free-Fall Time

What violent relaxation is: relaxation without collisions

In a galaxy or star cluster, stars almost never physically collide, and even close gravitational encounters are astronomically rare. The two-body relaxation time — the time for random star-on-star deflections to redistribute energy — scales as t_relax ≈ (0.1 N / ln N) × t_cross, where N is the number of stars and t_cross is the crossing time. For a galaxy with N ≈ 10^11 stars, this is around 10^18 years, far longer than the age of the Universe (1.38 × 10^10 yr). Galaxies should be dynamically 'frozen' — yet ellipticals are strikingly smooth and relaxed.

Lynden-Bell's 1967 insight resolved this. When a system is far from equilibrium and collapsing, its overall gravitational potential Φ is not static — it swings violently as mass rushes inward, overshoots, and re-expands. A star moving through a changing potential does not conserve its energy: dE/dt = ∂Φ/∂t. This collective mechanism redistributes energies across the whole system at once, in one dynamical time, with no collisions required.

The mechanism: energy exchange with a fluctuating potential

The dynamics are governed by the collisionless Boltzmann (Vlasov) equation, which states the fine-grained phase-space density f(x, v, t) is conserved along orbits: df/dt = ∂f/∂t + v·∇f − ∇Φ·∇_v f = 0, with Φ set self-consistently by Poisson's equation ∇²Φ = 4πGρ.

• Collapse: An out-of-equilibrium cloud free-falls; its potential well deepens on the timescale t_ff ≈ (Gρ)^(−1/2).
• Energy scrambling: Because ∂Φ/∂t ≠ 0, each star gains or loses energy depending on its phase during the collapse — a purely collective effect, independent of the star's mass.
• Damping: The potential oscillations die out within a few crossing times as the system virializes (2K + U = 0), leaving a quasi-stationary state.

Lynden-Bell derived the predicted equilibrium by maximizing a mixing entropy subject to conserved mass, momentum, energy, and phase-space volume. The result is a distribution resembling Fermi–Dirac statistics (an exclusion principle from Liouville's theorem forbidding overlapping phase elements), which reduces to the familiar isothermal Maxwell–Boltzmann form in the dilute, non-degenerate limit.

Key quantities: a worked example for a protogalaxy

The controlling number is the free-fall time, t_ff = √(3π / (32 G ρ)). Take a collapsing protogalaxy with mean density ρ ≈ 10^(−24) g/cm³ (roughly a few hydrogen atoms per cm³ smeared over the volume):

• t_ff ≈ 1 / √(Gρ) ≈ 2 × 10^15 s ≈ 60–100 million years. That is the entire duration of violent relaxation.
• Two-body t_relax for N = 10^11 and t_cross ≈ 10^8 yr: (0.1 × 10^11 / ln(10^11)) × 10^8 ≈ 4 × 10^17 yr — about 30 million times the age of the Universe.
• Ratio: violent relaxation is faster by a factor ~10^9, which is exactly why it dominates.

A crucial feature: because the energy change depends on ∂Φ/∂t and not on stellar mass m, violent relaxation produces no mass segregation and does not drive equipartition. Heavy and light stars end up with the same spatial distribution — a direct, testable contrast with slow collisional relaxation.

Where it appears and how we see its fingerprints

Violent relaxation is not observed as an event — it lasts only ~10^8 yr — but its end products are everywhere:

• Elliptical galaxies and bulges: Their smooth, centrally concentrated light follows the de Vaucouleurs R^(1/4) law (I(R) ∝ exp[−7.67((R/R_e)^(1/4) − 1)]). Numerical collapse simulations relaxing violently reproduce this profile, strongly linking the two.
• Galaxy mergers: When two disk galaxies merge, the transient deep potential wells drive fresh violent relaxation, scrambling the ordered disk orbits into a pressure-supported elliptical — the basis of the merger-origin picture for ellipticals.
• Dark-matter halos: N-body cosmological simulations show halos relaxing violently into the universal NFW profile (ρ ∝ 1/[(r/r_s)(1+r/r_s)²]).
• Young star clusters: Recent Gaia-based studies (e.g., 2024 work on the Lagoon Nebula cluster) find expanding subgroups consistent with violent relaxation after gas expulsion.

How it differs from its close cousins

Violent relaxation is easy to confuse with two neighboring processes:

• vs. Two-body relaxation: Two-body relaxation is collisional, driven by discrete encounters, scales with N, conserves each star's energy on average until slow diffusion, and pushes toward equipartition and mass segregation. Violent relaxation is collisionless, N-independent, changes individual energies wholesale, and is mass-blind.
• vs. Phase mixing: Phase mixing occurs even in a static potential: orbits of different periods shear apart, smoothing the coarse-grained distribution while each star's energy stays fixed. Violent relaxation additionally changes energies through ∂Φ/∂t. In practice the two act together during collapse.
• vs. Landau damping: The plasma-physics analog — collisionless damping of potential fluctuations via wave–particle energy exchange — is mathematically kin, since both obey Vlasov dynamics.

The unifying idea: violent relaxation preserves the fine-grained phase-space density exactly (Liouville), but tangles it into ever-finer filaments so the observable coarse-grained density relaxes to a smooth quasi-equilibrium.

Significance, open questions, and Lynden-Bell's legacy

Violent relaxation is a cornerstone of galactic dynamics: it explains why galaxies formed within a Hubble time can look dynamically old and relaxed, and why elliptical galaxies share regular, near-universal structure. Donald Lynden-Bell (1935–2018), who also co-predicted supermassive black holes powering quasars, laid the foundation in his landmark 1967 MNRAS paper 'Statistical mechanics of violent relaxation in stellar systems.'

Yet the theory remains genuinely incomplete — a well-known open problem:

• Incomplete relaxation: Real systems do not fully reach the predicted Lynden-Bell state; mixing halts before the maximum-entropy distribution is attained, depleting the high-energy tail. A true global maximum-entropy state for self-gravity does not even exist (the entropy is unbounded).
• Predictive limits: The final coarse-grained profile depends on the (messy) details of the initial conditions and the collapse, so the theory predicts form better than exact numbers.
• Ongoing work: Statistical-mechanics refinements (Chavanis and others) and high-resolution N-body/Vlasov simulations continue to probe why relaxation stalls — a debate still active in 2020s literature.

Violent relaxation vs. two-body (collisional) relaxation vs. phase mixing — the three routes a stellar system can take toward equilibrium.

Property	Violent relaxation	Two-body relaxation	Phase mixing
Driver	Global time-varying potential Φ(t)	Individual star–star gravitational encounters	Shearing of orbits with different periods
Timescale	~1–2 t_ff (crossing time)	t_relax ≈ (0.1 N / ln N) × t_cross	A few crossing times
Energy conservation	Individual E not conserved (dΦ/dt ≠ 0)	Individual E slowly changed by kicks	Individual E strictly conserved
Mass dependence	Independent of stellar mass	Favors equipartition (mass segregation)	None
For a galaxy (N~10^11)	~10^8 yr — effective	~10^18 yr — utterly negligible	~10^8 yr — complementary
End state	Quasi-equilibrium, near-Lynden-Bell DF	Full thermal equilibrium (never reached)	Smooth coarse-grained DF