Liouville-equation(2)

The equation of motion for the distribution function $ w _ {N} ( p , q ; t ) $ with respect to the momenta $ p = ( \mathbf p _ {1} \dots \mathbf p _ {N} ) $ and coordinates $ q = ( \mathbf r _ {1} \dots \mathbf r _ {N} ) $ of the $ N $- particle classical system

$$ \frac{\partial w _ {N} }{\partial t } = \ \{ H , w _ {N} \} _ { \mathop{\rm cl} . } \equiv \ \sum _ { i= } 1 ^ { N } \left ( \frac{\partial H }{\partial \mathbf r _ {i} } \frac{\partial w _ {N} }{\partial \mathbf p _ {i} } - \frac{\partial w _ {N} }{\partial \mathbf r _ {i} } \frac{\partial H }{\partial \mathbf p _ {i} } \right ) , $$

where $ H $ is the Hamiltonian of the system, and the curly brackets denote the classical Poisson brackets.

In the phase space $ ( p , q ) $ the distribution $ w _ {N} ( p , q ; t ) $ is associated with the density of the phase points (each of which corresponds to a definite mechanical state of the given system of $ N $ material points). By virtue of the fact that the trajectories of the motion of these points do not intersect, because of the uniqueness of the solutions of the mechanical equations of motion, and the fact that the phase volume is conserved, by the Liouville theorem (cf. Liouville theorems 4)), the ensemble of these points forms in the phase space a kind of incompressible liquid. The complete derivative of its density $ w _ {N} $ with respect to time is equal to zero:

$$ \frac{d w _ {N} }{d t } = \ \frac{\partial w _ {N} }{\partial t } + \sum _ { i= } 1 ^ { N } \left ( \frac{\partial w _ {N} }{\partial \mathbf r _ {i} } \frac{d \mathbf r _ {i} }{d t } + \frac{\partial w _ {N} }{\partial \mathbf p _ {i} } \frac{d \mathbf p _ {i} }{d t } \right ) = 0 . $$

This leads to Liouville's equation if, in accordance with the Hamilton equations, one expresses the derivatives of the coordinates and momenta in terms of the corresponding partial derivatives of the Hamiltonian.

Liouville's equation is used not only in the examination of general questions of statistical mechanics connected with the determination of the microscopic and macroscopic structures of the state of a system of several bodies, processes of tending towards equilibrium, problems of "mixing" in phase space, ergodicity, etc., but also in specific investigations, since Liouville's equation is the original equation for the construction of the Bogolyubov chain of equations, and consequently also for a different type of kinetic equation, by means of which problems of applied physics are solved.

In the case of quantum systems the role of Liouville's equation is played by the equation of motion for the statistical operator $ \rho ( t) $( the density matrix), which in the Schrödinger representation has the form

$$ \frac{\partial \rho }{\partial t } = \ \{ H , \rho \} _ { \mathop{\rm qu} . } \equiv \ \frac{1}{i \hbar } [ H \rho - \rho H ] , $$

where $ H $ is the Hamilton operator and the curly brackets denote the Poisson quantum bracket. This quantum Liouville equation is a consequence of the structure of the mixed state (described by a given statistical operator), in which each of the pure quantum-mechanical states comprising it evolves according to the Schrödinger equation.

References

[1]	H. Goldstein, "Classical mechanics" , Addison-Wesley (1957)
[2]	G.E. Uhlenbeck, G.V. Ford, "Lectures in statistical mechanics" , Amer. Math. Soc. (1963)
[3]	N.N. Bogolyubov, "Collected works" , 2 , Kiev (1970) (In Russian)

Comments

The Poisson quantum bracket is often written as $ \{ H , \rho \} _ { \mathop{\rm qu} . } = [ H , \rho ] / ( i \hbar ) $, with $ [ H , \rho ] = H \rho - \rho H $ the commutator bracket.

References

[a1]	C. Cercignani, "The Boltzmann equation and its applications" , Springer (1988)
[a2]	F.R. [F.R. Gantmakher] Gantmacher, "Lectures in analytical mechanics" , MIR (1975) (Translated from Russian)