05.02.07 · symplectic / hamiltonian

Liouville's volume theorem

shipped3 tiersLean: none

Anchor (Master): Liouville 1838 (originator); Arnold *Mathematical Methods of Classical Mechanics* §16 + Appendix 1; Abraham-Marsden *Foundations of Mechanics* §3.5

Intuition [Beginner]

Imagine a swarm of particles drifting in phase space — every particle is one possible state of a mechanical system, with its own position and its own momentum. As time passes the particles move under Hamilton's equations. The shape of the cloud changes, sometimes stretching, sometimes folding, but Liouville's theorem says one number stays the same: the total volume of phase space the cloud fills.

The cloud is incompressible. It can flatten into a long thin filament, twist into a spiral, wrap itself around an attractor — yet wherever it goes, the multidimensional volume it occupies at every moment matches the volume it started with.

This is what makes Hamiltonian mechanics special. Friction-driven systems shrink phase volume, settling onto attractors. Hamiltonian systems do not. The total swarm cannot be compressed.

Visual [Beginner]

A schematic of a phase-space region drifting under a Hamiltonian flow: the same region is shown at three later times, deformed in shape but with identical area highlighted on each frame.

The picture marks the deformation of the boundary along the flow and the conserved area enclosed.

Worked example [Beginner]

Take the harmonic oscillator on a single particle: position $q$ and momentum $p$ in the plane, with energy $H (q, p) = \frac{1}{2} (q^{2} + p^{2})$ . The Hamiltonian flow rotates phase space at unit angular speed: each point $(q, p)$ travels around the origin on a circle of radius $q^{2} + p^{2}$ .

Pick a small disc of radius $r$ centred at $(1, 0)$ . Its area is $π r^{2}$ . After time $t$ , the centre of the disc has rotated to $(cos t, sin t)$ . Every point in the disc has rotated by the same angle, so the disc itself rotates rigidly to a new disc of radius $r$ centred at $(cos t, sin t)$ . The shape is unchanged, the area is still $π r^{2}$ .

Now pick a region between two energy levels, say $H = 1$ and $H = 4$ — an annulus with inner radius $2$ and outer radius $8$ . The flow rotates each circle at its own rate, but the annulus as a set rotates rigidly back to itself, with area $4 π - 2 π = 2 π$ at every moment. What this tells us: phase volume between energy levels is conserved by the flow, even when the flow stretches and twists individual orbits.

Check your understanding [Beginner]

Formal definition [Intermediate+]

Let $(M^{2 n}, ω)$ be a symplectic manifold and $H : M \to R$ a smooth function. The Hamiltonian vector field $X_{H}$ is the unique vector field satisfying $ι_{X_{H}} ω = d H$ 05.02.01. Sign convention: $ι_{X_{H}} ω = d H$ , matching the convention of 05.02.01 and the rest of chapter 05.02.

The symplectic volume form is the top-degree form

Ω = \frac{ω ^{n}}{n !},

where $ω^{n} = ω \land ω \land \dots \land ω$ ( $n$ factors). Non-degeneracy of $ω$ guarantees $Ω$ is nowhere-vanishing, hence a volume form. In Darboux coordinates $(q^{1}, p_{1}, \dots, q^{n}, p_{n})$ where $ω = \sum d q^{i} \land d p_{i}$ , one computes $Ω = d q^{1} \land d p_{1} \land \dots \land d q^{n} \land d p_{n}$ , the standard Lebesgue volume of $R^{2 n}$ .

A flow $φ_{t} : M \to M$ preserves the volume form $Ω$ if $φ_{t}^{*} Ω = Ω$ for every $t$ in the domain of the flow.

Key theorem with proof [Intermediate+]

Theorem (Liouville's volume theorem). Let $(M^{2 n}, ω)$ be a symplectic manifold and $H : M \to R$ a smooth function with Hamiltonian vector field $X_{H}$ and flow $φ_{t}^{H}$ . Then $φ_{t}^{H}$ preserves the symplectic volume form: $(\varphi_t^H)^ \Omega = \Omega $f or a l l$ t$ in the domain of the flow.*

Proof. Let $Ω = ω^{n} / n!$ . By the standard identity for the Lie derivative of a flow, $(φ_{t}^{H})^{*} Ω = Ω$ for all $t$ if and only if $L_{X_{H}} Ω = 0$ . So it suffices to compute $L_{X_{H}} Ω$ .

The Lie derivative is a graded derivation of the wedge product:

L_{X_{H}} (ω^{n}) = n ω^{n - 1} \land L_{X_{H}} ω .

Hence $L_{X_{H}} Ω = (1/ n!) \cdot n ω^{n - 1} \land L_{X_{H}} ω = (ω^{n - 1} / (n - 1)!) \land L_{X_{H}} ω$ . So $L_{X_{H}} Ω = 0$ once $L_{X_{H}} ω = 0$ .

Compute $L_{X_{H}} ω$ using Cartan's formula $L_{X} = d ι_{X} + ι_{X} d$ :

L_{X_{H}} ω = d ι_{X_{H}} ω + ι_{X_{H}} d ω .

Closedness of $ω$ (a symplectic form is closed by definition 05.01.02) kills the second term: $d ω = 0$ . The first term is $d ι_{X_{H}} ω = d (d H) = 0$ by the defining equation of $X_{H}$ and $d^{2} = 0$ . Hence $L_{X_{H}} ω = 0$ , so $L_{X_{H}} Ω = 0$ , so $(φ_{t}^{H})^{*} Ω = Ω$ . $□$

Bridge. The volume preservation proved here builds toward 05.02.08 (Poincaré recurrence theorem), where the same volume form $Ω = ω^{n} / n!$ appears again as a finite invariant measure on a compact energy surface, and recurrence follows from the pigeonhole principle applied to the volume of iterates. The construction also appears again in 05.09.01 (KAM theorem), where the persistence of invariant tori under perturbation is measured against the Liouville volume on phase space — a positive-measure subset of tori survives because the perturbed flow remains volume-preserving by the same Cartan-formula argument. Putting these together, the foundational reason equilibrium statistical mechanics works at all is exactly that Hamiltonian flows preserve $Ω$ : any density $ρ$ that is a function of conserved quantities — the canonical $e^{- β H}$ , the microcanonical $δ (H - E)$ — is automatically stationary under the flow, and the central insight is that volume preservation plus energy conservation produce the equilibrium ensembles without further input. The bridge between the analytic identity $L_{X_{H}} ω = 0$ and the geometric conclusion $(φ_{t}^{H})^{*} Ω = Ω$ is the foundational reason classical statistical mechanics is consistent with Hamiltonian dynamics.

Counterexamples to common slips

Volume preservation is not length preservation. A Hamiltonian flow can stretch directions arbitrarily; only the product of stretches across all $2 n$ directions is conserved. The harmonic oscillator preserves length and volume because rotations preserve both, but a generic Hamiltonian flow stretches some directions and shrinks others by exactly compensating amounts.
Gradient flows are not volume-preserving. A flow generated by $- \nabla V$ on Euclidean space contracts phase volume at rate $Δ V$ wherever $V$ is strictly convex. The Hamiltonian flow generated by $V$ on $T^{*} M$ is volume-preserving; the gradient flow of the same $V$ on $M$ is not.
Closedness is essential. If one drops $d ω = 0$ and uses a non-degenerate but non-closed two-form, the analogous "Hamiltonian" vector field has $L_{X} ω = ι_{X} d ω \neq = 0$ in general, and the volume $ω^{n} / n!$ is not preserved.
Time-dependent Hamiltonians still work. When $H = H (q, p, t)$ depends on time, the same calculation applies pointwise in $t$ : $L_{X_{H_{t}}} ω = 0$ for each $t$ , and the time-dependent flow remains volume-preserving. Energy is no longer conserved, but volume still is.

Exercises [Intermediate+]

Exercise 2 (easy, symbolic).

Compute $ω^{n} / n!$ in Darboux coordinates on $R^{2 n}$ where $ω = \sum_{i} d q^{i} \land d p_{i}$ , and confirm that the result is the standard Lebesgue volume form.

Hint

Expand $ω^{n}$ by multilinearity; only terms containing every $d q^{i}$ and every $d p_{i}$ exactly once survive antisymmetry.

Answer

$ω^{n} = (\sum_{i} d q^{i} \land d p_{i})^{n}$ . Each surviving term in the expansion picks one factor $d q^{i_{k}} \land d p_{i_{k}}$ from each copy of $ω$ with all $i_{k}$ distinct; $n!$ orderings of the same set of indices give the same wedge product up to sign, and the sign works out so that the sum is $n! \cdot d q^{1} \land d p_{1} \land \dots \land d q^{n} \land d p_{n}$ . Hence $ω^{n} / n! = d q^{1} \land d p_{1} \land \dots \land d q^{n} \land d p_{n}$ , the standard Lebesgue volume form on $R^{2 n}$ . Rubric: full credit for the combinatorial $n!$ and the explicit identification with Lebesgue volume.

Exercise 3 (medium, short-answer).

Derive the Liouville equation for a phase-space density $ρ (q, p, t)$ transported by a Hamiltonian flow:

\frac{\partial ρ}{\partial t} + {ρ, H} = 0.

Hint

Phase-space density along a flow line is constant: $ρ (φ_{t}^{H} (x), t) = ρ (x, 0)$ for all $x$ if $ρ$ is "advected" by the flow. Differentiate in $t$ .

Answer

By volume preservation, the density $ρ$ at the moving point is constant along each orbit when $ρ$ is advected (no source, no sink). Differentiating $ρ (φ_{t}^{H} (x), t) = const$ in $t$ at $t = 0$ gives $\partial_{t} ρ + X_{H} (ρ) = 0$ . The Hamiltonian vector field acts on functions as $X_{H} (ρ) = {ρ, H}$ in the convention ${f, g} = ω (X_{f}, X_{g}) = X_{g} (f)$ used elsewhere in chapter 05.02 05.02.02. Hence $\partial_{t} ρ + {ρ, H} = 0$ , the Liouville equation. Rubric: full credit for the advection argument and the identification of $X_{H} (ρ)$ with ${ρ, H}$ .

Exercise 4 (medium, short-answer).

Show that any density of the form $ρ = f (H)$ for a smooth $f : R \to R_{\geq 0}$ is a stationary solution of the Liouville equation. Why does this make $f (H) \cdot Ω$ an equilibrium measure?

Hint

Compute ${f (H), H}$ using the chain rule for Poisson brackets.

Answer

${f (H), H} = f^{'} (H) {H, H} = 0$ because the Poisson bracket of any function with itself vanishes by antisymmetry. Hence $\partial_{t} ρ + {ρ, H} = 0$ reduces to $\partial_{t} ρ = 0$ : any density that is a function of $H$ alone is time-independent under the Hamiltonian flow. Multiplying by the invariant volume $Ω$ , the measure $μ = f (H) Ω$ satisfies $(φ_{t}^{H})_{*} μ = μ$ , hence is a flow-invariant measure on phase space. The canonical ensemble $f (H) = e^{- β H} / Z$ and the microcanonical ensemble $f (H) = δ (H - E) / Z (E)$ are the two standard examples. Rubric: full credit for the bracket calculation, the stationarity conclusion, and naming at least one equilibrium ensemble.

Exercise 5 (medium, short-answer).

Let $V : R^{n} \to R$ be smooth and consider the gradient flow $\overset{x}{˙} = - \nabla V (x)$ on $R^{n}$ . Compute the divergence of the gradient vector field and conclude that gradient flows preserve volume only in a degenerate case.

Hint

The divergence of $- \nabla V$ is $- Δ V$ .

Answer

$\nabla \cdot (- \nabla V) = - Δ V = - \sum_{i} \partial^{2} V / \partial x_{i}^{2}$ . Volume is preserved at a point if and only if $Δ V = 0$ there. So the gradient flow is volume-preserving on all of $R^{n}$ if and only if $V$ is harmonic, in which case $V$ has no isolated minima and the gradient flow has no attracting equilibria — a degenerate case from the dynamics standpoint. For any $V$ with a strict local minimum, $Δ V > 0$ on a neighbourhood of the minimum and phase volume contracts toward the minimum. This contrasts sharply with Hamiltonian flows, which are volume-preserving regardless of the specific $H$ . Rubric: full credit for $- Δ V$ and the harmonic-only conclusion.

Exercise 6 (hard, short-answer).

Show that on a closed (compact, boundaryless) symplectic manifold, the symplectic volume of a smooth Hamiltonian flow is finite and the iterates of any positive-volume measurable set must eventually overlap. State the corollary precisely.

Hint

Compactness gives finite total volume. Apply the pigeonhole principle to the orbits of a measurable set under a discrete-time iterate of the flow.

Answer

On a closed manifold $M$ the volume $vol (M) = \int_{M} Ω < \infty$ because $Ω$ is a smooth top-degree form on a compact set. Fix the time-one map $φ = φ_{1}^{H}$ , which preserves $Ω$ . Let $A \subset M$ have positive volume. The iterates $φ^{k} (A)$ all have the same volume $vol (A) > 0$ . If they were pairwise disjoint, the total volume of $⨆_{k} φ^{k} (A)$ would be $\sum_{k} vol (A) = \infty$ , exceeding $vol (M)$ . Hence $φ^{j} (A) \cap φ^{k} (A) \neq = \emptyset$ for some $j < k$ , and applying $φ^{- j}$ gives $A \cap φ^{k - j} (A) \neq = \emptyset$ . Corollary: every positive-volume set returns to itself in positive time. This is the volume-preservation step of Poincaré recurrence 05.02.08. Rubric: full credit for the pigeonhole argument and the explicit return statement.

Exercise 7 (hard, short-answer).

A volume-preserving vector field $X$ on a Riemannian manifold $(M, g)$ need not be Hamiltonian. Construct an example on $R^{4}$ of a vector field with vanishing Euclidean divergence that is not the Hamiltonian vector field of any smooth function.

Hint

Hamiltonian vector fields satisfy a global cohomological condition: $ι_{X} ω$ is exact. A divergence-free vector field whose contraction with $ω$ is closed but not exact is volume-preserving without being Hamiltonian.

Answer

The condition $ι_{X} ω = d H$ is the Hamiltonian condition; the relaxed condition $d (ι_{X} ω) = 0$ characterises symplectic (locally Hamiltonian) vector fields. Volume preservation is yet weaker: $L_{X} (ω^{n}) = 0$ . On $R^{4}$ with standard $ω = d q^{1} \land d p_{1} + d q^{2} \land d p_{2}$ , every symplectic vector field is global Hamiltonian because $H^{1} (R^{4}) = 0$ , so this example must use a different ambient form. Take instead $T^{4} = R^{4} / Z^{4}$ with the same $ω$ . The constant vector field $X = \partial_{q^{1}}$ has $ι_{X} ω = d p_{1}$ , which is closed but not exact on $T^{4}$ (the cohomology class $[d p_{1}]$ is non-zero in $H^{1} (T^{4})$ ). Hence $L_{X} ω = d (d p_{1}) = 0$ , so $X$ preserves $ω$ and therefore preserves $ω^{2} /2! = d q^{1} \land d p_{1} \land d q^{2} \land d p_{2}$ , the volume. But there is no global $H : T^{4} \to R$ with $d H = d p_{1}$ , so $X$ is symplectic but not globally Hamiltonian. The Euclidean-divergence calculation also gives zero in coordinates, confirming volume preservation directly. Rubric: full credit for the topological obstruction and the explicit example on $T^{4}$ .

Lean formalization [Intermediate+]

[object Promise]

The proof would proceed by establishing $L_{X_{H}} ω = 0$ via Cartan's formula plus closedness of $ω$ , then using the graded-Leibniz rule for Lie derivatives over wedge products to extend to $ω^{n} / n!$ . Each step is a candidate Mathlib contribution, listed in the lean_mathlib_gap field above.

Advanced results [Master]

The Hamiltonian-flow case of Liouville's theorem fits inside a wider family of measure-preservation results, each arising from a vanishing-divergence condition. Restating the theorem in the language most useful for downstream applications, and listing the parallel results, makes the place of Liouville's theorem in classical mechanics, ergodic theory, and statistical physics precise.

Statement and equivalent forms. On $(M^{2 n}, ω)$ with smooth Hamiltonian $H$ , the following are equivalent:

$(φ_{t}^{H})^{*} (ω^{n} / n!) = ω^{n} / n!$ for all $t$ in the domain of the flow.
$L_{X_{H}} (ω^{n} / n!) = 0$ pointwise.
$L_{X_{H}} ω = 0$ pointwise — the flow is symplectic, not merely volume-preserving.
In Darboux coordinates, $\sum_{i} (\partial X_{H}^{q^{i}} / \partial q^{i} + \partial X_{H}^{p_{i}} / \partial p_{i}) = 0$ — the Euclidean divergence of $X_{H}$ vanishes.

The implication chain is (1) ⇔ (2) by Lie-derivative-of-flow, (2) ⇐ (3) by the graded-Leibniz expansion $L_{X_{H}} (ω^{n}) = n ω^{n - 1} \land L_{X_{H}} ω$ , (3) by Cartan's formula plus closedness of $ω$ , and (3) ⇒ (4) by direct coordinate calculation.

The implication (3) ⇒ volume preservation is strict: a symplectic flow preserves $ω$ itself, not just $ω^{n}$ . Symplectic flows form a much smaller class than volume-preserving flows when $n \geq 2$ — a generic divergence-free vector field on $R^{4}$ does not preserve $ω = d q^{1} \land d p_{1} + d q^{2} \land d p_{2}$ . Liouville's theorem says Hamiltonian flows are at the strictest end of this hierarchy.

Liouville equation for phase-space densities. For a transported density $ρ (q, p, t)$ , the equation $\partial_{t} ρ + {ρ, H} = 0$ governs the evolution. Equilibrium solutions are functions of conserved quantities; on a symplectic manifold whose Hamiltonian has only $H$ as a smooth conserved quantity, these are functions of $H$ alone. The canonical ensemble $ρ_{β} = e^{- β H} / Z$ and the microcanonical ensemble $ρ_{E} = δ (H - E) / Z (E)$ are the two standard equilibrium measures. The compatibility of equilibrium statistical mechanics with Hamiltonian dynamics rests on this single computation: ${f (H), H} = 0$ .

Volume-rigid versus length-flexible. Symplectic geometry is volume-rigid in the sense of Liouville's theorem: every symplectomorphism preserves $ω^{n} / n!$ . But it is length-flexible in the sense of Gromov non-squeezing 05.07.01: a symplectomorphism can stretch one coordinate direction arbitrarily as long as it shrinks the conjugate direction by the same factor. Riemannian geometry is the opposite: length-rigid (isometries preserve every distance) and volume-rigid as a consequence. Symplectic geometry is the only $C^{\infty}$ -rigid geometry on phase space that allows the volume-but-not-length combination — a fact captured precisely by Gromov's symplectic-capacity theory.

Generalisations to non-Hamiltonian flows. Any divergence-free vector field on a Riemannian manifold preserves the Riemannian volume. The Liouville equation $\partial_{t} ρ + \nabla \cdot (ρX) = 0$ for densities transported by such a flow reduces to $\partial_{t} ρ + X (ρ) = 0$ when $\nabla \cdot X = 0$ . Examples: incompressible fluid flow (Euler equations on a Riemannian manifold), magnetic fields (Maxwell on closed manifolds), the geodesic flow on the unit cotangent bundle preserves the contact volume $α \land (d α)^{n - 1}$ which is the restriction of $ω^{n}$ to the level set. Hamiltonian flows are the special case where the divergence-free condition arises from a closed two-form rather than from an additional input.

Ergodic-theoretic consequences. On a closed symplectic manifold, $Ω = ω^{n} / n!$ is a finite invariant measure, so the Hamiltonian flow is a measure-preserving system in the sense of Birkhoff. The Birkhoff ergodic theorem then asserts that time averages along orbits converge to spatial averages with respect to a flow-invariant probability measure (which is $Ω/ vol (M)$ when $H$ is the only conserved quantity, or a measure restricted to a smaller invariant set otherwise). The Poincaré recurrence theorem 05.02.08 is the volume-preservation corollary established by the pigeonhole argument of Exercise 6. KAM theory 05.09.01 uses Liouville volume to measure the surviving-tori set: a positive proportion of phase volume remains foliated by invariant tori under small perturbation, and the precise quantitative statement is a Lebesgue-measure-of-tori bound with respect to $Ω$ .

Originator's setting. Liouville's 1838 paper Sur la théorie de la variation des constantes arbitraires worked in the setting of Hamilton's canonical equations in $R^{2 n}$ — the modern coordinate-free symplectic-manifold framing came a century later with Cartan, Souriau, and Abraham-Marsden. Liouville's calculation was the divergence form $\nabla \cdot X_{H} = 0$ , derived directly from the equality of mixed partial derivatives of $H$ . The cohomology-and-Cartan-formula framing presented above is logically equivalent and globally clean, but Liouville's coordinate proof is essentially the same computation written without the language of differential forms.

Synthesis. Hamiltonian flows on a symplectic manifold preserve a canonical volume form $Ω = ω^{n} / n!$ , derived from the symplectic form by purely combinatorial wedge-product algebra. The proof reduces in three lines to the closedness of $ω$ and the defining equation $ι_{X_{H}} ω = d H$ — no input from the specific $H$ enters. Read in the opposite direction, every volume-preserving Hamiltonian flow has its volume form computable from the symplectic structure alone, with no choice of measure required: the symplectic structure determines the equilibrium measure up to a function of conserved quantities. Putting these together, classical statistical mechanics is consistent with Hamiltonian dynamics not by accident but by the same identity $L_{X_{H}} ω = 0$ : any equilibrium measure of the form $f (H) Ω$ is automatically flow-invariant, recovering the canonical and microcanonical ensembles as the simplest possible solutions of the equilibrium condition. The bridge between the analytic identity $L_{X_{H}} ω = 0$ and the geometric conclusion that Hamiltonian flows preserve $Ω$ is the foundational reason every Poincaré-recurrence, KAM, and ergodic-theory result on phase space rests on a single closedness computation. The central insight is that volume preservation is a corollary of symplectic preservation, and symplectic preservation is a corollary of $d ω = 0$ — the entire structural apparatus follows from one identity.

Full proof set [Master]

Theorem (Liouville's volume theorem). Restated: on a symplectic manifold $(M^{2 n}, ω)$ with smooth $H$ and Hamiltonian flow $φ_{t}^{H}$ , $(φ_{t}^{H})^{*} (ω^{n} / n!) = ω^{n} / n!$ for all $t$ in the flow domain.

Proof. Given in the Intermediate Key theorem section above. The chain of identities is $L_{X_{H}} ω = d ι_{X_{H}} ω + ι_{X_{H}} d ω = d (d H) + 0 = 0$ , then $L_{X_{H}} (ω^{n}) = n ω^{n - 1} \land L_{X_{H}} ω = 0$ , then equivalence with flow invariance via the standard $\frac{d}{d t} (φ_{t}^{*} α) = φ_{t}^{*} (L_{X} α)$ identity. $□$

Lemma (Lie-derivative-of-flow identity). For any smooth vector field $X$ with flow $φ_{t}$ and any smooth differential form $α$ , $\frac{d}{dt}(\varphi_t^ \alpha) = \varphi_t^(\mathcal{L}_X \alpha)$.

Proof. By the semigroup property $φ_{t + s} = φ_{t} \circ φ_{s}$ , $\frac{d}{d t} (φ_{t}^{*} α) ∣_{t = t_{0}} = \frac{d}{d s} ∣_{s = 0} (φ_{t_{0} + s}^{*} α) = φ_{t_{0}}^{*} \frac{d}{d s} ∣_{s = 0} (φ_{s}^{*} α) = φ_{t_{0}}^{*} (L_{X} α)$ , where the last equality is the definition of the Lie derivative. $□$

Lemma (graded Leibniz for Lie derivative). For any vector field $X$ and any differential forms $α, β$ , $L_{X} (α \land β) = L_{X} α \land β + α \land L_{X} β$ .

Proof. The Lie derivative is a derivation of the algebra of differential forms, so it satisfies the (un-graded) Leibniz rule on wedge products. The standard derivation comes from differentiating $φ_{t}^{*} (α \land β) = φ_{t}^{*} α \land φ_{t}^{*} β$ at $t = 0$ . Iterating gives $L_{X} (ω^{n}) = n ω^{n - 1} \land L_{X} ω$ when applied to $ω \land ω \land \dots \land ω$ . $□$

Theorem (Liouville equation for phase-space densities). Let $ρ_{t} : M \to R$ be a smooth time-dependent density on $(M^{2 n}, ω)$ , transported by $φ_{t}^{H}$ in the sense $ρ_{t} (φ_{t}^{H} (x)) = ρ_{0} (x)$ . Then $\partial_{t} ρ + {ρ, H} = 0$ .

Proof. Differentiate $ρ_{t} (φ_{t}^{H} (x)) = ρ_{0} (x)$ in $t$ at $t = 0$ for fixed $x$ : $\partial_{t} ρ_{t} (x) ∣_{t = 0} + d ρ_{0} (X_{H} (x)) = 0$ . The directional derivative $d ρ (X_{H}) = X_{H} (ρ) = {ρ, H}$ in the convention ${f, g} = ω (X_{f}, X_{g}) = X_{g} (f)$ of 05.02.02. Hence $\partial_{t} ρ + {ρ, H} = 0$ at every point and every time. $□$

Theorem (equilibrium measures are functions of $H$ ). For any smooth $f : R \to R$ , the density $ρ = f (H)$ satisfies $\partial_{t} ρ + {ρ, H} = 0$ , hence is a stationary solution of the Liouville equation.

Proof. ${f (H), H} = f^{'} (H) {H, H} = 0$ because ${H, H} = ω (X_{H}, X_{H}) = 0$ by antisymmetry of $ω$ . Combined with $\partial_{t} f (H) = 0$ (since $f (H)$ is time-independent as a function on $M$ ), the Liouville equation holds. $□$

Connections [Master]

Hamiltonian vector field 05.02.01. Liouville's theorem is the canonical example of a structural theorem about Hamiltonian flows that follows from the defining equation $ι_{X_{H}} ω = d H$ together with closedness of $ω$ . The construction sets the template: any structural property of Hamiltonian flows is one Cartan-formula calculation away from the symplectic data.
Symplectic manifold 05.01.02. The closedness condition $d ω = 0$ — half of the definition of a symplectic structure — is exactly what makes the Liouville calculation work. Non-closed almost-symplectic forms do not preserve their associated volume forms under the analogous "Hamiltonian" flows, so Liouville's theorem is structurally tied to closedness.
Poisson bracket 05.02.02. The Liouville equation $\partial_{t} ρ + {ρ, H} = 0$ is the bracket form of the volume preservation theorem. Equilibrium measures satisfy ${ρ, H} = 0$ , identifying flow-invariant measures with functions of conserved quantities.
Poincaré recurrence theorem 05.02.08. Direct corollary on a finite-volume (typically compact-energy-surface) phase space: every positive-volume measurable set returns close to itself in positive time. The pigeonhole argument of Exercise 6 is the volume-preservation step of the recurrence proof.
KAM theorem 05.09.01. Liouville volume measures the surviving-tori set under small Hamiltonian perturbation. Quantitative KAM bounds are stated in terms of the Lebesgue measure of the Diophantine-frequency set with respect to $ω^{n} / n!$ ; volume preservation is what makes "most" tori a measure-theoretic statement.
Action-angle coordinates 05.02.04. On an integrable Hamiltonian system, the Liouville volume in action-angle coordinates is $\prod_{i} d I^{i} \land d θ_{i}$ , and the flow $\dot{θ} = ω (I)$ , $\dot{I} = 0$ preserves this product directly; volume preservation is built into the action-angle structure at the coordinate level.
Geodesic flow as a Hamiltonian flow 05.02.06. The geodesic flow on $T^{*} M$ preserves the Liouville volume of $T^{*} M$ , and the restricted flow on the unit cotangent bundle preserves the contact volume $α \land (d α)^{n - 1}$ . The ergodic theory of geodesic flow on negatively-curved manifolds rests on this.
Symplectic capacity / non-squeezing [05.07.01, 05.07.02]. Volume preservation is necessary but not sufficient for symplectic flows: the "much stronger" rigidity captured by symplectic capacities — and the Gromov non-squeezing theorem — distinguishes symplectic from merely volume-preserving maps, marking the difference between volume rigidity and length flexibility.

Historical & philosophical context [Master]

Joseph Liouville published Sur la théorie de la variation des constantes arbitraires in the Journal de Mathématiques Pures et Appliquées in 1838 ^{[Liouville 1838]}. The paper concerned the integration of Hamilton's canonical equations through variation of parameters, and as a corollary established the volume-preservation property of the canonical flow in $R^{2 n}$ . Liouville's argument used the divergence form: he computed $\sum_{i} (\partial \overset{q}{˙}^{i} / \partial q^{i} + \partial \overset{p}{˙}_{i} / \partial p_{i}) = 0$ directly from Hamilton's equations and the equality of mixed partial derivatives, and concluded that the Jacobian of the canonical flow is identically one. The differential-forms reformulation $L_{X_{H}} (ω^{n} / n!) = 0$ followed Cartan's Leçons sur les invariants intégraux (1922), which placed integral invariants at the centre of Hamiltonian mechanics in the modern coordinate-free language.

The connection to statistical mechanics was made explicit by Boltzmann and Gibbs in the 1870s and 1900s. Boltzmann's 1872 H-theorem and Gibbs's 1902 Elementary Principles in Statistical Mechanics both rested on volume preservation as the structural reason that equilibrium ensembles — canonical, microcanonical, grand canonical — are functions of conserved quantities. Gibbs's framing of statistical mechanics through phase-space density evolution under the Liouville equation became the standard formulation, and the Liouville equation continues to bear his name in some texts.

Poincaré 1890 Sur le problème des trois corps drew the recurrence corollary from Liouville's theorem in the context of celestial mechanics: in a closed phase region of finite volume, almost every initial condition returns arbitrarily close to itself infinitely often. The recurrence theorem became the flashpoint of the Boltzmann-Loschmidt-Zermelo controversy on the foundation of irreversibility — a debate resolved on the kinetic-theory side by recognising recurrence times as astronomically large for macroscopic systems while remaining finite, and on the formal side by the development of mixing and the SRB-measure framework in the 20th century.

Arnold Mathematical Methods of Classical Mechanics §16 ^[Arnold] gives the modern symplectic-manifold treatment at the level of advanced graduate physics, deriving Liouville's theorem from the Cartan-formula argument and using it as input for Poincaré recurrence and the Poincaré-Cartan integral invariant. Abraham-Marsden Foundations of Mechanics §3.5 ^{[Abraham-Marsden]} gives the canonical mathematical-physics treatment, with the volume form $ω^{n} / n!$ identified as the canonical Lebesgue measure on $T^{*} M$ in Darboux coordinates. Goldstein Classical Mechanics §9.9 ^[Goldstein] is the standard graduate-physics reference, covering the divergence form alongside the differential-forms statement.

Bibliography [Master]

[object Promise]

Prerequisites

05.02.01
05.01.02

Used in

05.02.08
05.09.01

Tier anchors

beginner: Strogatz *Nonlinear Dynamics and Chaos* §6 informal phase-space density; Goldstein *Classical Mechanics* §9.9 informal
intermediate: Goldstein *Classical Mechanics* §9.9; Arnold *Mathematical Methods of Classical Mechanics* §16
master: Liouville 1838 (originator); Arnold *Mathematical Methods of Classical Mechanics* §16 + Appendix 1; Abraham-Marsden *Foundations of Mechanics* §3.5

References

TODO_REF
Liouville 1838 — Sur la théorie de la variation des constantes arbitraires · J. Math. Pures Appl. (1) 3, 342–349; originator paper
TODO_REF
Arnold — Mathematical Methods of Classical Mechanics · §16 phase-volume preservation; Appendix 1 Riemannian curvature contrast
TODO_REF
Goldstein — Classical Mechanics · §9.9 Liouville's theorem and the equation of motion for phase-space density
TODO_REF
Abraham-Marsden — Foundations of Mechanics · §3.5 symplectic volume and Liouville's theorem on cotangent bundles

Reviewer

TBD

Estimated time

beginner: 12m
intermediate: 32m
master: 60m