3.2: Phase Flows in Classical Mechanics
( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand\Dalpha
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[1], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dbeta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[2], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dgamma
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[3], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Ddelta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[4], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Depsilon
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[5], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dvarepsilon
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[6], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dzeta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[7], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Deta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[8], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dtheta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[9], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dvartheta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[10], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Diota
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[11], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dkappa
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[12], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dlambda
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[13], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Dvarpi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[14], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\DGamma
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[15], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\DDelta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[16], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\DTheta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[17], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vmu
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[18], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vnu
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[19], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vxi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[20], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vom
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[21], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vpi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[22], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vvarpi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[23], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vrho
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[24], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vvarrho
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[25], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vsigma
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[26], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vvarsigma
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[27], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vtau
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[28], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vupsilon
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[29], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vphi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[30], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vvarphi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[31], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vchi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[32], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vpsi
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[33], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\Vomega
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[34], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\VGamma
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[35], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\VDelta
Callstack: at (Template:MathJaxArovas), /content/body/div/p[1]/span[36], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\newcommand\BI{\mib I}}
\)
\newcommand { M}
\newcommand { m}
}
\( \newcommand\tcb{\textcolor{blue}\)
\( \newcommand\tcr{\textcolor{red}\)
1$#1_$
\newcommand\SZ{\textsf Z}} \( \newcommand\kFd{k\ns_{\RF\dar}\)
\newcommand\mutB{\tilde\mu}\ns_\ssr{B}
\( \newcommand\xhihOZ
Callstack: at (Template:MathJaxArovas), /content/body/div/span[1], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
\( \newcommand\labar
Callstack: at (Template:MathJaxArovas), /content/body/div/span[2], line 1, column 1 at template() at (Bookshelves/Thermodynamics_and_Statistical_Mechanics/Book:_Thermodynamics_and_Statistical_Mechanics_(Arovas)/03:_Ergodicity_and_the_Approach_to_Equilibrium/3.02:_Phase_Flows_in_Classical_Mechanics), /content/body/p/span, line 1, column 23
Hamiltonian evolution
The master equation provides us with a semi-phenomenological description of a dynamical system’s relaxation to equilibrium. It explicitly breaks time reversal symmetry. Yet the microscopic laws of Nature are (approximately) time-reversal symmetric. How can a system which obeys Hamilton’s equations of motion come to equilibrium?
Let’s start our investigation by reviewing the basics of Hamiltonian dynamics. Recall the Lagrangian L=L(q,˙q,t)=T−V. The Euler-Lagrange equations of motion for the action S[q(t)]=∫dtL are
˙p∗σ=ddt(∂L∂˙q∗σ)=∂L∂q∗σ ,
where p∗σ is the canonical momentum conjugate to the generalized coordinate q∗σ:
p†σ=∂L∂˙qσ .
The Hamiltonian, H(q,p) is obtained by a Legendre transformation,
H(q,p)=r∑σ=1p†σ˙q†σ−L .
Note that
dH=r∑σ=1(p†σd˙q†σ+˙q†σdp†σ−∂L∂qσdq†σ−∂L∂˙qσd˙q†σ)−∂L∂tdt =r∑σ=1(˙q†σdp†σ−∂L∂qσdq†σ)−∂L∂tdt .
Thus, we obtain Hamilton’s equations of motion,
∂H∂pσ=˙q†σ,∂H∂qσ=−∂L∂qσ=−˙p†σ
and
dHdt=∂H∂t=−∂L∂t .
Define the rank 2r vector φ by its components,
φ†i={q†iif\ 1≤i≤rp†i−rif\ r≤i≤2r .
Then we may write Hamilton’s equations compactly as
˙φi=Jij∂H∂φj ,
where
J=(0†r×r1†r×r−1†r×r0†r×r)
is a rank 2r matrix. Note that Jt=−J, J is antisymmetric, and that J2=−1†2r×2r.
Dynamical systems and the evolution of phase space volumes
Consider a general dynamical system,
dφdt=V(φ) ,
where φ(t) is a point in an n-dimensional phase space. Consider now a compact2 region R∗0 in phase space, and consider its evolution under the dynamics. That is, R∗0 consists of a set of points {φ|φ∈R∗0}, and if we regard each φ∈R∗0 as an initial condition, we can define the time-dependent set R(t) as the set of points φ(t) that were in R∗0 at time t=0:
R(t)={φ(t)|φ(0)∈R∗0} .
Now consider the volume Ω(t) of the set R(t). We have
Ω(t)=∫R(t)dμ
where
dμ=dφ∗1dφ∗2⋯dφ∗n ,
for an n-dimensional phase space. We then have
Ω(t+dt)=∫R(t+dt)dμ′=∫R(t)dμ|∂φ∗i(t+dt)∂φ∗j(t)| ,
where
|∂φ∗i(t+dt)∂φ∗j(t)|≡∂(φ′1,…,φ′n)∂(φ∗1,…,φ∗n)
is a determinant, which is the Jacobean of the transformation from the set of coordinates {φ∗i=φ∗i(t)} to the coordinates {φ′i=φ∗i(t+dt)}. But according to the dynamics, we have
φ∗i(t+dt)=φ∗i(t)+V∗i(φ(t))dt+O(dt2)
and therefore
∂φ∗i(t+dt)∂φ∗j(t)=δ∗ij+∂V∗i∂φ∗jdt+O(dt2) .
We now make use of the equality
lndetM=TrlnM ,
for any matrix M, which gives us3, for small ε,
det(1+εA)=expTrln(1+εA)=1+εTrA+12ε2((TrA)2−Tr(A2))+…
Thus,
Ω(t+dt)=Ω(t)+∫R(t)dμ∇⋅Vdt+O(dt2) ,
which says
dΩdt=∫R(t)dμ∇⋅V=∫∂R(t)dSˆn⋅V
Here, the divergence is the phase space divergence,
∇⋅V=n∑i=1∂V∗i∂φ∗i ,
and we have used the divergence theorem to convert the volume integral of the divergence to a surface integral of ˆn⋅V, where ˆn is the surface normal and dS is the differential element of surface area, and ∂R denotes the boundary of the region R. We see that if ∇⋅V=0 everywhere in phase space, then Ω(t) is a constant, and phase space volumes are preserved by the evolution of the system.
For an alternative derivation, consider a function ϱ(φ,t) which is defined to be the density of some collection of points in phase space at phase space position φ and time t. This must satisfy the continuity equation,
∂ϱ∂t+∇⋅(ϱV)=0 .
This is called the continuity equation. It says that ‘nobody gets lost’. If we integrate it over a region of phase space R, we have
ddt∫Rdμϱ=−∫Rdμ∇⋅(ϱV)=−∫∂RdSˆn⋅(ϱV) .
It is perhaps helpful to think of ϱ as a charge density, in which case J=ϱV is the current density. The above equation then says
dQ∗Rdt=−∫∂RdSˆn⋅J ,
where Q∗R is the total charge contained inside the region R. In other words, the rate of increase or decrease of the charge within the region R is equal to the total integrated current flowing in or out of R at its boundary.

The Leibniz rule lets us write the continuity equation as
∂ϱ∂t+V⋅∇ϱ+ϱ∇⋅V=0 .
But now suppose that the phase flow is divergenceless, ∇⋅V=0. Then we have
DϱDt≡(∂∂t+V⋅∇)ϱ=0 .
The combination inside the brackets above is known as the convective derivative. It tells us the total rate of change of ϱ for an observer moving with the phase flow. That is
ddtϱ(φ(t),t)=∂ϱ∂φ∗idφ∗idt+∂ϱ∂t=n∑i=1V∗i∂ρ∂φ∗i+∂ϱ∂t=DϱDt .
If Dϱ/Dt=0, the local density remains the same during the evolution of the system. If we consider the ‘characteristic function’
ϱ(φ,t=0)={1if φ∈R∗00otherwise
then the vanishing of the convective derivative means that the image of the set R∗0 under time evolution will always have the same volume.
Hamiltonian evolution in classical mechanics is volume preserving. The equations of motion are
˙qi=+∂H∂pi,˙pi=−∂H∂qi
A point in phase space is specified by r positions q∗i and r momenta p∗i, hence the dimension of phase space is n=2r:
φ=(qp),V=(˙q˙p)=( ∂H/∂p−∂H/∂q) .
Hamilton’s equations of motion guarantee that the phase space flow is divergenceless:
∇⋅V=r∑i=1{∂˙qi∂qi+∂˙pi∂pi}=r∑i=1{∂∂qi(∂H∂pi)+∂∂pi(−∂H∂qi)}=0 .
Thus, we have that the convective derivative vanishes, viz.
DϱDt≡∂ϱ∂t+V⋅∇ϱ=0 ,
for any distribution ϱ(φ,t) on phase space. Thus, the value of the density ϱ(φ(t),t) is constant, which tells us that the phase flow is incompressible. In particular, phase space volumes are preserved.
Liouville’s equation and the microcanonical distribution
Let ϱ(φ)=ϱ(q,p) be a distribution on phase space. Assuming the evolution is Hamiltonian, we can write
∂ϱ∂t=−˙φ⋅∇ϱ=−r∑k=1(˙q∗k∂∂q∗k+˙p∗k∂∂p∗k)ϱ=−iˆLϱ ,
where ˆL is a differential operator known as the Liouvillian:
ˆL=−ir∑k=1{∂H∂p∗k∂∂q∗k−∂H∂q∗k∂∂p∗k} .
Equation ???, known as Liouville’s equation, bears an obvious resemblance to the Schrödinger equation from quantum mechanics.
Suppose that Λ∗a(φ) is conserved by the dynamics of the system. Typical conserved quantities include the components of the total linear momentum (if there is translational invariance), the components of the total angular momentum (if there is rotational invariance), and the Hamiltonian itself (if the Lagrangian is not explicitly time-dependent). Now consider a distribution ϱ(φ,t)=ϱ(Λ∗1,Λ∗2,…,Λ∗k) which is a function only of these various conserved quantities. Then from the chain rule, we have
˙φ⋅∇ϱ=∑a∂ϱ∂Λ∗a˙φ⋅∇Λ∗a=0 ,
since for each a we have
dΛ∗adt=r∑σ=1(∂Λ∗a∂q∗σ˙q∗σ+∂Λ∗a∂p∗σ˙p∗σ)=˙φ⋅∇Λ∗a=0 .
We conclude that any distribution ϱ(φ,t)=ϱ(Λ∗1,Λ∗2,…,Λ∗k) which is a function solely of conserved dynamical quantities is a stationary solution to Liouville’s equation.
Clearly the microcanonical distribution,
ϱ∗E(φ)=δ(E−H(φ))D(E)=δ(E−H(φ))∫dμδ(E−H(φ)) ,
is a fixed point solution of Liouville’s equation.