World Library  
Flag as Inappropriate
Email this Article

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 traveling through phase-space is constant with time.

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


  • Liouville equations 1
  • Other formulations 2
    • Poisson bracket 2.1
    • Ergodic theory 2.2
    • Symplectic geometry 2.3
    • Quantum Liouville equation 2.4
  • Remarks 3
  • See also 4
  • References 5

Liouville equations

Evolution of an ensemble of classical systems in phase space (top). Each system consists of one massive particle in a one-dimensional potential well (red curve, lower figure). Whereas the motion of an individual member of the ensemble is given by Hamilton's equations, Liouville's equations describe the flow of the whole distribution. The motion is analogous to a dye in an incompressible fluid.

The Liouville equation describes the time evolution of the phase space distribution function. Although the equation is usually referred to as the "Liouville equation", Josiah Willard Gibbs was the first to recognize the importance of this equation as the fundamental equation of statistical mechanics.[1][2] It is referred to as the Liouville equation because its derivation for non-canonical systems utilises an identity first derived by Liouville in 1838.[3] Consider a Hamiltonian dynamical system with canonical coordinates q_i and conjugate momenta p_i, where i=1,\dots,n. Then the phase space distribution \rho(p,q) determines the probability \rho(p,q)\,d^nq\,d^n p that the system will be found in the infinitesimal phase space volume d^nq\,d^n p. The Liouville equation governs the evolution of \rho(p,q;t) in time t:

\frac{d\rho}{dt}= \frac{\partial\rho}{\partial t} +\sum_{i=1}^n\left(\frac{\partial\rho}{\partial q_i}\dot{q}_i +\frac{\partial\rho}{\partial p_i}\dot{p}_i\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 proof of Liouville's theorem uses the n-dimensional divergence theorem. This proof is based on the fact that the evolution of \rho obeys an n-dimensional version of the continuity equation:

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

That is, the tuplet (\rho, \rho\dot{q}_i,\rho\dot{p}_i) is a conserved current. Notice that the difference between this and Liouville's equation are the terms

\rho\sum_{i=1}^n\left( \frac{\partial\dot{q}_i}{\partial q_i} +\frac{\partial\dot{p}_i}{\partial p_i}\right) =\rho\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,

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 \rho/dt, is zero follows from the equation of continuity by noting that the 'velocity field' (\dot p , \dot q) 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 – p_i say – 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 invariant under time translations, and the generator (or Noether charge) of the symmetry is the Hamiltonian.

Other formulations

Poisson bracket

The theorem is often restated in terms of the Poisson bracket as

\frac{\partial\rho}{\partial t}=-\{\,\rho,H\,\}

or in terms of the Liouville operator or Liouvillian,

\mathrm{i}\hat{\mathbf{L}}=\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]=\{\cdot,H\}


\frac{\partial \rho }{\partial t}+{\mathrm{i}\hat{\mathbf{L}}}\rho =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 6n-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 phase space is represented as a symplectic manifold. The theorem then states that the natural volume form on a symplectic manifold is invariant under the Hamiltonian flows. The symplectic structure is represented as a 2-form, given as a sum of wedge products of dpi with dqi. The volume form is the top exterior power of the symplectic 2-form, and is just another representation of the measure on the phase space described above. One formulation of the theorem states that the Lie derivative of this volume form is zero along every Hamiltonian vector field.

In fact, the symplectic structure itself is preserved, not only its top exterior power.

Quantum Liouville equation

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, the Von Neumann equation. 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[4][5]

\frac{\partial \rho}{\partial t}=\frac{1}{i \hbar}[H,\rho]

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{d}{dt}\langle A\rangle = \frac{1}{i \hbar}\langle [A,H] \rangle

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.


See also


  • Modern Physics, by R. Murugeshan, S. Chand publications
  • Liouville's theorem in curved space-time : Gravitation § 22.6, by Misner,Thorne and Wheeler, Freeman
  1. ^ J. W. Gibbs, "On the Fundamental Formula of Statistical Mechanics, with Applications to Astronomy and Thermodynamics." Proceedings of the American Association for the Advancement of Science, 33, 57-58 (1884). Reproduced in The Scientific Papers of J. Willard Gibbs, Vol II (1906), pp. 16.
  2. ^  
  3. ^ [J. Liouville, Journ. de Math., 3, 349(1838)].
  4. ^ , by Breuer and Petruccione, p110The theory of open quantum systems.
  5. ^ , by Schwabl, p16Statistical mechanics.
  6. ^ [for a particularly clear derivation see "The Principles of Statistical Mechanics" by R.C. Tolman , Dover(1979), p48-51].
  7. ^ Nearly identical to proof in this WorldHeritage article. Assumes (without proof) the n-dimensional continuity equation. Retrieved 01/06/2014.
  8. ^ A rigorous proof based on how the Jacobian volume element transforms under Hamiltonian mechanics. Retrieved 01/06/2014.
  9. ^ Uses the n-dimensional divergence theorem (without proof) Retrieved 01/06/2014.
  10. ^ Proves Liouville's theorem using the language of modern differential geometry Retrieved 10/01/2015.
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.

Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.