L'équation de Liouville décrit l'évolution temporelle de la densité de probabilité ρ dans l'espace des phases. Cette densité de probabilité est définie comme la probabilité pour que l'état du système soit représenté par un point à l'intérieur du volume Γ considéré.
En mécanique classique
On utilise les coordonnées généralisées (q,p) où N est la dimension du système. La densité de probabilité est définie par la probabilité ρ(p,q)dNqdNp de rencontrer l'état du système dans le volume infinitésimal dNqdNp.
Lorsqu'on calcule l'évolution temporelle cette densité de probabilité ρ(p,q), on obtient :
dtdρ=∂t∂ρ+i=1∑N[∂qi∂ρq˙i+∂pi∂ρp˙i]=0
On peut utiliser les équations canoniques de Hamilton en les remplaçant dans l'équation précédente :