ProofComplete

Proof of Liouville's Theorem

ProofVolume Preservation by Hamiltonian Flows

We prove that Hamiltonian flows preserve the Liouville volume form Ω=ωnn!\Omega = \frac{\omega^n}{n!}.

Setup: Let (M,ω)(M, \omega) be a 2n2n-dimensional symplectic manifold, H:MRH: M \to \mathbb{R} a Hamiltonian function, and ϕt\phi_t the flow of the Hamiltonian vector field XHX_H defined by ιXHω=dH\iota_{X_H} \omega = dH.

Step 1: Hamiltonian vector fields preserve ω\omega

By Cartan's formula:

LXHω=d(ιXHω)+ιXH(dω)\mathcal{L}_{X_H} \omega = d(\iota_{X_H} \omega) + \iota_{X_H}(d\omega)

Since ιXHω=dH\iota_{X_H} \omega = dH and dω=0d\omega = 0 (closedness), we have:

LXHω=d(dH)+0=0\mathcal{L}_{X_H} \omega = d(dH) + 0 = 0

Thus XHX_H preserves ω\omega: LXHω=0\mathcal{L}_{X_H} \omega = 0.

Step 2: Preservation extends to ωn\omega^n

We need to show LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n) = 0. Using the derivation property of the Lie derivative:

LXH(ωn)=LXH(ωωn1)\mathcal{L}_{X_H}(\omega^n) = \mathcal{L}_{X_H}(\omega \wedge \omega^{n-1})

Since LXH\mathcal{L}_{X_H} is a derivation on forms:

=(LXHω)ωn1+ωLXH(ωn1)= (\mathcal{L}_{X_H} \omega) \wedge \omega^{n-1} + \omega \wedge \mathcal{L}_{X_H}(\omega^{n-1})

By Step 1, LXHω=0\mathcal{L}_{X_H} \omega = 0, so:

=0ωn1+ωLXH(ωn1)=ωLXH(ωn1)= 0 \wedge \omega^{n-1} + \omega \wedge \mathcal{L}_{X_H}(\omega^{n-1}) = \omega \wedge \mathcal{L}_{X_H}(\omega^{n-1})

Proceeding inductively, if LXH(ωn1)=0\mathcal{L}_{X_H}(\omega^{n-1}) = 0, then LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n) = 0.

The base case is LXHω=0\mathcal{L}_{X_H} \omega = 0 from Step 1. Therefore LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n) = 0 for all nn.

Step 3: Integration yields flow preservation

The Lie derivative LXH\mathcal{L}_{X_H} measures the infinitesimal change along the flow:

LXH(ωn)=limt0ϕt(ωn)ωnt\mathcal{L}_{X_H}(\omega^n) = \lim_{t \to 0} \frac{\phi_t^*(\omega^n) - \omega^n}{t}

If LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n) = 0, then:

ddtϕt(ωn)=0\frac{d}{dt}\phi_t^*(\omega^n) = 0

Integrating from 00 to tt:

ϕt(ωn)=ϕ0(ωn)=ωn\phi_t^*(\omega^n) = \phi_0^*(\omega^n) = \omega^n

Since Ω=ωnn!\Omega = \frac{\omega^n}{n!}, we have ϕtΩ=Ω\phi_t^* \Omega = \Omega, completing the proof.

Remark

The key insight is that preservation of ω\omega by XHX_H automatically implies preservation of ωn\omega^n. This is a special feature of Hamiltonian vector fields: they preserve not just volume, but the entire symplectic structure.

DefinitionDivergence-Free Condition

In local coordinates (qi,pi)(q_i, p_i), the Hamiltonian vector field is:

XH=i=1n(HpiqiHqipi)X_H = \sum_{i=1}^n \left(\frac{\partial H}{\partial p_i} \frac{\partial}{\partial q_i} - \frac{\partial H}{\partial q_i} \frac{\partial}{\partial p_i}\right)

The divergence with respect to the volume form Ω=dq1dqndp1dpn\Omega = dq_1 \wedge \cdots \wedge dq_n \wedge dp_1 \wedge \cdots \wedge dp_n is:

div(XH)=i=1n(2Hqipi2Hpiqi)=0\text{div}(X_H) = \sum_{i=1}^n \left(\frac{\partial^2 H}{\partial q_i \partial p_i} - \frac{\partial^2 H}{\partial p_i \partial q_i}\right) = 0

The mixed partials cancel by equality of mixed partials, so XHX_H is divergence-free.

ExampleExplicit Computation for Harmonic Oscillator

For H=12(p2+q2)H = \frac{1}{2}(p^2 + q^2) on R2\mathbb{R}^2:

XH=pqqpX_H = p \frac{\partial}{\partial q} - q \frac{\partial}{\partial p}

The flow is ϕt(q,p)=(qcost+psint,qsint+pcost)\phi_t(q, p) = (q\cos t + p\sin t, -q\sin t + p\cos t), a rotation.

For any region AA:

ϕt(A)dqdp=Aϕt(dqdp)\int_{\phi_t(A)} dq \wedge dp = \int_A \phi_t^*(dq \wedge dp)

Computing the pullback:

ϕt(dqdp)=(costdq+sintdp)(sintdq+costdp)\phi_t^*(dq \wedge dp) = (\cos t \, dq + \sin t \, dp) \wedge (-\sin t \, dq + \cos t \, dp) =(cos2t+sin2t)dqdp=dqdp= (\cos^2 t + \sin^2 t) dq \wedge dp = dq \wedge dp

Thus area is preserved.

Remark

Alternative proof via Jacobian: If ϕt\phi_t is the Hamiltonian flow, its Jacobian matrix JϕtJ_{\phi_t} satisfies:

det(Jϕt)=1\det(J_{\phi_t}) = 1

for all tt. This follows from the fact that ϕtSp(2n,R)\phi_t \in \text{Sp}(2n, \mathbb{R}) at each point, and symplectic matrices have determinant 1.

The volume element transforms as:

ϕt(Ω)=det(Jϕt)Ω=Ω\phi_t^*(\Omega) = \det(J_{\phi_t}) \cdot \Omega = \Omega

DefinitionConservation of Phase Space Density

In statistical mechanics, the phase space density ρ(q,p,t)\rho(q, p, t) evolves according to Liouville's equation:

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

This is equivalent to ddt(ρϕt)=0\frac{d}{dt}(\rho \circ \phi_t) = 0, i.e., ρ\rho is constant along Hamiltonian flow. Combined with volume preservation, this gives the conservation of the measure ρΩ\rho \, \Omega.