Liouville's Theorem

Liouville's Theorem: Conservation of Phase Space Density

Liouville's Theorem is a fundamental result in classical mechanics that connects the deterministic laws of motion (Hamilton's equations) with statistical descriptions of large systems. It states that phase space behaves like an incompressible fluid: as a system evolves, the "cloud" of possible states may distort, but its volume remains constant.

1. The Statement

Let \( \rho(q, p, t) \) be the probability density of finding a system at coordinates \( (q, p) \) in phase space at time \( t \). Liouville's Theorem asserts that the density of states in the immediate neighborhood of a moving point in phase space remains constant:

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

Where \( \{ \rho, H \} \) is the Poisson bracket. Geometrically, this implies the conservation of phase space volume.

2. Detailed Proof

Step 1: The Continuity Equation

We treat the ensemble of systems as a fluid flowing through \( 2N \)-dimensional phase space. Since the number of systems is conserved, the density \( \rho \) must satisfy the continuity equation:

\[ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{V}) = 0 \]

Here, \( \mathbf{V} = (\dot{q}_1, \dots, \dot{q}_N, \dot{p}_1, \dots, \dot{p}_N) \) is the generalized velocity vector.

Step 2: Divergence of the Flow

We expand the divergence term \( \nabla \cdot (\rho \mathbf{V}) \):

\[ \nabla \cdot (\rho \mathbf{V}) = \sum_{i=1}^N \left[ \frac{\partial}{\partial q_i}(\rho \dot{q}_i) + \frac{\partial}{\partial p_i}(\rho \dot{p}_i) \right] \] \[ = \sum_{i=1}^N \left[ \left( \dot{q}_i \frac{\partial \rho}{\partial q_i} + \dot{p}_i \frac{\partial \rho}{\partial p_i} \right) + \rho \left( \frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}_i}{\partial p_i} \right) \right] \]

Step 3: Applying Hamilton's Equations

We substitute Hamilton's equations of motion, \( \dot{q}_i = \frac{\partial H}{\partial p_i} \) and \( \dot{p}_i = -\frac{\partial H}{\partial q_i} \), into the divergence term multiplying \( \rho \):

\[ \frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}_i}{\partial p_i} = \frac{\partial}{\partial q_i}\left(\frac{\partial H}{\partial p_i}\right) + \frac{\partial}{\partial p_i}\left(-\frac{\partial H}{\partial q_i}\right) \] \[ = \frac{\partial^2 H}{\partial q_i \partial p_i} - \frac{\partial^2 H}{\partial p_i \partial q_i} = 0 \]

The velocity field in phase space has zero divergence (it is solenoidal).

Step 4: The Total Derivative

The continuity equation thus simplifies to:

\[ \frac{\partial \rho}{\partial t} + \sum_{i=1}^N \left( \dot{q}_i \frac{\partial \rho}{\partial q_i} + \dot{p}_i \frac{\partial \rho}{\partial p_i} \right) = 0 \]

The left side is exactly the total time derivative (convective derivative) \( \frac{d\rho}{dt} \). Therefore, \( \frac{d\rho}{dt} = 0 \). Q.E.D.

3. Applications in Physics

• Statistical Mechanics: It provides the justification for the Microcanonical Ensemble. Since \( \rho \) is constant along trajectories, a distribution that is uniform over the energy surface (\( \rho = \text{const} \)) is a valid equilibrium solution.
• Beam Physics (Particle Accelerators): It leads to the conservation of "Emittance." It implies that you cannot simply squeeze a particle beam into a smaller spatial volume without increasing its velocity spread (heating it up).
• Stellar Dynamics: Used to model the evolution of large clusters of stars (like globular clusters) where individual stellar collisions are rare, treating the cluster as a collisionless fluid in phase space.