Give every possible state one point: position across, momentum up. A living system is that point, sliding. A cloud of systems is a flowing region — and Liouville proved the astonishing thing: as the cloud flows, it may stretch and shear beyond recognition, but its area never changes. The flow is incompressible; the point that was you is still one point; the curves never cross. Down the center, states go in, the Hamiltonian flow carries them, the conserved area comes out. The blue team builds and defends; the red team tries to break it.
source Liouville, J. (1838), Note sur la théorie de la variation des constantes arbitraires, J. de Math. Pures Appl. sér. 1, 3, 342–349; Gibbs (1902), Elementary Principles in Statistical Mechanics — archive.org/details/elementaryprinci00gibbrich. AMBER: Liouville's 1838 note is pure ODE theory — he never wrote "phase space"; Gibbs and Boltzmann gave it the mechanical reading. Rendered, not quoted.
Take the harmonic oscillator in natural units (m=1, k=1): energy H = ½p² + ½q². Hamilton's equations give the flow field:
q̇ = ∂H/∂p = p and ṗ = −∂H/∂q = −q.
Its divergence is identically zero: ∂q̇/∂q + ∂ṗ/∂p = ∂²H/∂q∂p − ∂²H/∂p∂q = 0 — mixed partials cancel. Zero divergence is incompressibility, and incompressibility is Liouville's theorem: the flow moves points around but never compresses or rarefies the cloud. Live divergence at sampled points:
| (q, p) | div(flow) |
|---|
A single trajectory obeying Hamilton's equations lives in the-hamiltonian. Phase space is the arena that holds all of them at once — Liouville 1838, Gibbs 1902.
Once you know the volume of a state-cloud is conserved, you can count states without following any one of them. That single fact — volume is an invariant of the flow — is the load-bearing beam under statistical mechanics: the microcanonical measure, entropy, the whole ensemble picture. Each sphere is the next one's premise.
The blue team's live check: evolve a patch of initial conditions by the symplectic map, re-measure its area by the shoelace formula, and confirm preservation to 1e-6 — plus unit Jacobian and zero divergence. If red adds damping, this badge is where it shows.
Seed the flow with a region, not a point: a ring of initial conditions (q, p) centred in phase space, each one a whole system waiting to evolve. Its starting area is measured once by the shoelace formula.
This region — its centre, radius, and vertex count — is exactly what the engine below takes and pushes along Hamilton's flow. Watch it stretch and shear while its area holds.
The patch is advanced by a symplectic integrator (semi-implicit Euler), whose one-step Jacobian determinant is exactly 1 — so area is preserved by construction, not by luck. Plain Euler would spiral out; this does not. The concentric rings are constant-energy contours; every point rides one and stays on it.
What the machine proves, live: the patch's area after — symplectic steps equals its starting area to within 1e-6; the one-step Jacobian is 1; the flow's divergence is 0; each trajectory stays on its energy contour and no two ever cross. That conserved area is the output — Liouville's theorem, measured.
The blue team's witness (left) confirms preservation live; the red team (right) adds damping to make the area shrink.
And it is volume, not shape: the region shears without bound. After enough time it filaments so finely that any coarse-grained view sees it "fill" the accessible space — the origin of the arrow of time in Gibbs' picture, which fine-grained Liouville flatly denies. Both are true at different resolutions; that tension is real.
"Liouville's theorem says energy is conserved." Cut. It says phase-space volume is conserved. Energy conservation is a separate fact (time-invariance of H); a driven Hamiltonian can conserve volume while energy changes.
"Any integrator that keeps area is fine." Cut. Plain (explicit) Euler has Jacobian 1+½h²ω²·… ≠ 1 and spirals outward; only a symplectic map preserves the 2-form exactly — the engine uses semi-implicit Euler, det = 1.
"Trajectories can merge if two states end up equal." Kept, corrected. The flow is a bijection (det ≠ 0), so distinct states stay distinct forever — determinism is non-crossing.
The red team's move: slip a damping term into the flow — ṗ becomes −q − γp — while still claiming Liouville holds. The map's determinant drops below 1; the patch spirals inward and its area shrinks. The blue team's witness (window 7) is watching.
Add friction and the cloud collapses — the witness recomputes the area, sees it shrink, and turns red. Nothing is faked; the dissipation is real and it is caught.