# Liouville's theorem (Hamiltonian)

﻿
Liouville's theorem (Hamiltonian)

In physics, Liouville's theorem, named after the French mathematician Joseph Liouville, is a key theorem in classical statistical and Hamiltonian mechanics. It asserts that the phase-space distribution function is constant along the trajectories of the system - that is that the density of system points in the vicinity of a given system point travelling through phase-space is constant with time.

There are also related mathematical results in symplectic topology and ergodic theory.

Liouville equation

The Liouville equation describes the time evolution of phase space distribution function (while measure is the correct term from mathematics, physicists generally call it a distribution). Consider a dynamical system with canonical coordinates $q^ i$ and conjugate momenta $p_i$, where $i=1,dots,d$. Then the phase space distribution $ho\left(p,q\right)$ determines the probability $ho\left(p,q\right),d^dq,d^dp$ that a particle will be found in the infinitesimal phase space volume $d^dq,d^dp$. The Liouville equation governs the evolution of $ho\left(p,q;t\right)$ in time $t$:

:$frac\left\{d ho\right\}\left\{dt\right\}=frac\left\{partial ho\right\}\left\{partial t\right\}+sum_\left\{i=1\right\}^dleft\left(frac\left\{partial ho\right\}\left\{partial q^i\right\}dot\left\{q\right\}^i+frac\left\{partial ho\right\}\left\{partial p_i\right\}dot\left\{p\right\}_i ight\right)=0.$

Time derivatives are denoted by dots, and are evaluated according to Hamilton's equations for the system. This equation demonstrates the conservation of density in phase space (which was Gibbs's name for the theorem). Liouville's theorem states that

:"The distribution function is constant along any trajectory in phase space."

A simple proof of the theorem is to observe that the evolution of $ho$ is "defined" by the continuity equation:

:$frac\left\{partial ho\right\}\left\{partial t\right\}+sum_\left\{i=1\right\}^dleft\left(frac\left\{partial\left( hodot\left\{q\right\}^i\right)\right\}\left\{partial q^i\right\}+frac\left\{partial\left( hodot\left\{p\right\}_i\right)\right\}\left\{partial p_i\right\} ight\right)=0.$

That is, the tuplet $\left( ho, hodot\left\{q\right\}^i, hodot\left\{p\right\}_i\right)$ is a conserved current. Notice that the difference between this and Liouville's equation are the terms

:$hosum_\left\{i=1\right\}^dleft\left(frac\left\{partialdot\left\{q\right\}^i\right\}\left\{partial q^i\right\}+frac\left\{partialdot\left\{p\right\}_i\right\}\left\{partial p_i\right\} ight\right)= hosum_\left\{i=1\right\}^dleft\left(frac\left\{partial^2 H\right\}\left\{partial q^i,partial p_i\right\}-frac\left\{partial^2 H\right\}\left\{partial p_i partial q^i\right\} ight\right)=0,$

where $H$ is the Hamiltonian, and Hamilton's equations have been used. That is, viewing the motion through phase space as a 'fluid flow' of system points, the theorem that the convective derivative of the density $d ho/dt$ is zero follows from the equation of continuity by noting that the 'velocity field' $\left(dot p , dot q\right)$ in phase space has zero divergence (which follows from Hamilton's relations).

Another illustration is to consider the trajectory of a cloud of points through phase space. It is straightforward to show that as the cloud stretches in one coordinate &ndash; $p_i$ say &ndash; it shrinks in the corresponding $q^i$ direction so that the product $Delta p_i Delta q^i$ remains constant.

Equivalently, the existence of a conserved current implies, via Noether's theorem, the existence of a symmetry. The symmetry is invariance under time translations, and the generator (or Noether charge) of the symmetry is the Hamiltonian.

Physical interpretation

The expected total number of particles is the integral over phase space of the distribution::$N=int d^dq,d^dp, ho\left(p,q\right).$A normalizing factor is conventionally included in the phase space measure but has here been omitted. In the simple case of a nonrelativistic particle moving in Euclidean space under a force field $mathbf\left\{F\right\}$ with coordinates $mathbf\left\{x\right\}$ and momenta $mathbf\left\{p\right\}$, Liouville's theorem can be written:$frac\left\{partial ho\right\}\left\{partial t\right\}+frac\left\{mathbf\left\{p\left\{m\right\}cdot abla_mathbf\left\{x\right\} ho+mathbf\left\{F\right\}cdot abla_mathbf\left\{p\right\} ho=0.$This is different from the Vlasov equation, or sometimes the Collisionless Boltzmann Equation, in astrophysics. The latter, which has a 6-D phase space, is used to describe the evolution of a large number of collisionless particles moving under the influence of gravity and/or electromagnetic field.

In classical statistical mechanics, the number of particles $N$ is very large, (typically of order Avogadro's number, for a laboratory-scale system). Setting $partial ho/partial t=0$ gives an equation for the stationary states of the system and can be used to find the density of microstates accessible in a given statistical ensemble. The stationary states equation is satisfied by $ho$ equal to any function of the Hamiltonian $H$: in particular, it is satisfied by the Maxwell-Boltzmann distribution $hopropto e^\left\{-H/kT\right\}$, where $T$ is the temperature and $k$ the Boltzmann constant.

Other formulations

Poisson bracket

The theorem is often restated in terms of the Poisson bracket as:$frac\left\{partial ho\right\}\left\{partial t\right\}=-\left\{, ho,H,\right\}$or in terms of the Liouville operator or Liouvillian,:$hat\left\{mathbf\left\{L=sum_\left\{i=1\right\}^\left\{d\right\}left \left[frac\left\{partial H\right\}\left\{partial p_\left\{ifrac\left\{partial\right\}\left\{partial q^\left\{i-frac\left\{partial H\right\}\left\{partial q^\left\{ifrac\left\{partial \right\}\left\{partial p_\left\{i ight\right] ,$as:$frac\left\{partial ho \right\}\left\{partial t\right\}+\left\{hat\left\{L ho =0.$

Ergodic theory

In ergodic theory and dynamical systems, motivated by the physical considerations given so far, there is a corresponding result also referred to as Liouville's theorem. In Hamiltonian mechanics, the phase space is a smooth manifold that comes naturally equipped with a smooth measure (locally, this measure is the 6"n"-dimensional Lebesgue measure). The theorem says this smooth measure is invariant under the Hamiltonian flow. More generally, one can describe the necessary and sufficient condition under which a smooth measure is invariant under a flow. The Hamiltonian case then becomes a corollary.

Symplectic geometry

In terms of symplectic geometry, the theorem states that the d-power of the symplectic structure (2-form, formed by summation of the wedge products of $Delta p_i$ and $Delta q_i$) has a zero Lie derivative for its Hamiltonian evolution. (The d-power of the symplectic structure is just the measure on the phase space described above.)

In fact, the symplectic structure itself is preserved (not only the d-power). For this reason, in this context, symplectic structure is also called Poincaré invariant. Hence the theorem about Poincaré invariant is a generalization of the Liouville's theorem.

Further generalization is also possible. In the frame of [http://daarb.narod.ru/wircq-eng.html#ihf invariant Hamiltonian formalism] , the theorem about existence of symplectic structure on invariant phase space turns out to be a deep generalization of the theorem about Poincaré invariant.

Quantum mechanics

The analog of Liouville equation in quantum mechanics describes the time evolution of a mixed state. Canonical quantization yields a quantum-mechanical version of this theorem. This procedure, often used to devise quantum analogues of classical systems, involves describing a classical system using Hamiltonian mechanics. Classical variables are then re-interpreted as quantum operators, while Poisson brackets are replaced by commutators. In this case, the resulting equation is:$frac\left\{partial\right\}\left\{partial t\right\} ho=-frac\left\{i\right\}\left\{hbar\right\} \left[H, ho\right]$where ρ is the density matrix.

When applied to the expectation value of an observable, the corresponding equation is given by Ehrenfest's theorem, and takes the form

:$frac\left\{d\right\}\left\{dt\right\}langle A angle = frac\left\{i\right\}\left\{hbar\right\}langle \left[H,A\right] angle$

where $A$ is an observable. Note the sign difference, which follows from the assumption that the operator is stationary and the state is time-dependent.

Remarks

*The Liouville equation is valid for both equilibrium and nonequilibrium systems. It is a fundamental equation of nonequilibrium statistical mechanics. Its approximation to collisional systems is called the Boltzmann equation.
*The Liouville equation is integral to the proof of the fluctuation theorem from which the second law of thermodynamics can be derived. It is also the key component of the derivation of Green-Kubo relations for linear transport coefficients such as shear viscosity, thermal conductivity or electrical conductivity.
* Virtually any textbook on Hamiltonian mechanics, advanced statistical mechanics, or symplectic geometry will derive the Liouville theorem.

ee also

* reversible reference system propagation algorithm (r-RESPA)

Wikimedia Foundation. 2010.

### Look at other dictionaries:

• Liouville's theorem — has various meanings, all mathematical results named after Joseph Liouville:*In complex analysis, see Liouville s theorem (complex analysis). *In conformal mappings, see Liouville s theorem (conformal mappings). *In Hamiltonian mechanics, see… …   Wikipedia

• Liouville's equation — For Liouville s equation in dynamical systems, see Liouville s theorem (Hamiltonian). In differential geometry, Liouville s equation, named after Joseph Liouville, is the equation satisfied by the conformal factor f of a metric f^2 (dx^2 + dy^2)… …   Wikipedia

• Hamiltonian mechanics — is a re formulation of classical mechanics that was introduced in 1833 by Irish mathematician William Rowan Hamilton. It arose from Lagrangian mechanics, a previous re formulation of classical mechanics introduced by Joseph Louis Lagrange in 1788 …   Wikipedia

• Liouville, Joseph — ▪ French mathematician born March 24, 1809, Saint Omer, France died September 8, 1882, Paris       French mathematician known for his work in analysis, differential geometry, and number theory and for his discovery of transcendental numbers i.e …   Universalium

• Hamiltonian system — In classical mechanics, a Hamiltonian system is a physical system in which forces are velocity invariant. Hamiltonian systems are studied in Hamiltonian mechanics. In mathematics, a Hamiltonian system is a system of differential equations which… …   Wikipedia

• Joseph Liouville — Infobox Scientist name =Joseph Liouville box width = image width =150px caption =Joseph Liouville birth date = March 24 1809 birth place = death date = September 8 1882 death place = residence = citizenship = nationality = French ethnicity =… …   Wikipedia

• Ehrenfest theorem — The Ehrenfest theorem, named after Paul Ehrenfest, relates the time derivative of the expectation value for a quantum mechanical operator to the commutator of that operator with the Hamiltonian of the system. It is:frac{d}{dt}langle A angle =… …   Wikipedia

• Poincaré recurrence theorem — In mathematics, the Poincaré recurrence theorem states that certain systems will, after a sufficiently long time, return to a state very close to the initial state. The Poincaré recurrence time is the length of time elapsed until the recurrence.… …   Wikipedia

• Microcanonical ensemble — Statistical mechanics Thermodynamics · …   Wikipedia

• Canonical transformation — In Hamiltonian mechanics, a canonical transformation is a change of canonical coordinates (mathbf{q}, mathbf{p}, t) ightarrow (mathbf{Q}, mathbf{P}, t) that preserves the form of Hamilton s equations, although it might not preserve the… …   Wikipedia