JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

Phase Space, Flow, and Liouville's Theorem

Zoom out from single trajectories to the whole flow. A blob of states distorts wildly yet preserves its volume — the incompressibility that founds statistical mechanics and raises deep puzzles about time.

The Hamiltonian flow

Hamilton's equations do more than move one point. At every location in phase space they assign a velocity vector (\dot q,\dot p), so they define a velocity field covering the whole space. The trajectories are that field's flow lines — the Hamiltonian flow. Instead of asking where one particle goes, we now ask how the entire fluid of possible states moves.

Liouville's theorem: volume is sacred

Now the central result. Take a small region of phase space — a blob of nearby initial conditions, i.e. an ensemble — and let it flow. Its shape can stretch into thin filaments and fold like taffy, but its total phase-space volume never changes. This is Liouville's theorem.

The proof is a two-line calculation. The flow is volume-preserving exactly when its velocity field is divergence-free. Compute that divergence for Hamilton's equations and the mixed partial derivatives of H cancel identically:

\nabla\cdot\dot{\mathbf z} = \sum_i\!\left(\frac{\partial \dot q_i}{\partial q_i} + \frac{\partial \dot p_i}{\partial p_i}\right) = \sum_i\!\left(\frac{\partial^2 H}{\partial q_i\,\partial p_i} - \frac{\partial^2 H}{\partial p_i\,\partial q_i}\right) = 0

The Hamiltonian flow is divergence-free — hence volume-preserving.

So the flow is incompressible, exactly like an ideal fluid. Equivalently, the density of representative points measured moving along with the flow stays constant. That statement — the phase-space density is conserved along every trajectory — is the dynamical form of the theorem.

\frac{d\rho}{dt} = \frac{\partial \rho}{\partial t} + \{\rho, H\} = 0

The convective density is constant along the flow (the Poisson bracket appears — next guide).

The symplectic structure beneath it

Volume preservation is only the surface. The deeper truth is that phase space carries a symplectic form — a way of measuring oriented areas of (q_i,p_i) pairs — and the Hamiltonian flow preserves it exactly. Writing the state as a column vector \mathbf z=(q_1,\dots,p_n), Hamilton's equations take a single elegant matrix form.

\dot{\mathbf z} = \Omega\,\nabla_{\!\mathbf z} H, \qquad \Omega = \begin{pmatrix} 0 & I_n \\ -I_n & 0 \end{pmatrix}

The symplectic (canonical) matrix form of Hamilton's equations.

Preserving \Omega is a stronger, finer condition than preserving volume: in more than one degree of freedom it constrains each (q_i,p_i) area separately (the Poincaré invariants), not just the grand product. This symplectic rigidity is the true fingerprint of Hamiltonian dynamics and the reason its numerical integrators are designed to respect it.

Why this founds statistical mechanics — and puzzles about time

Liouville's theorem is the load-bearing wall of statistical mechanics. Because any density that depends on phase-space coordinates only through H automatically satisfies \{\rho,H\}=0, such a distribution is stationary — it does not evolve. That single fact justifies the equilibrium ensembles: the microcanonical ensemble as a uniform density on an energy surface, and by extension the canonical ensemble.