Macroscopic fluid approximations fail whenever plasma behavior depends on non-Maxwellian particle distributions, wave-particle resonances, or collisionless damping. To capture the full microscopic reality of high-temperature plasmas, we must transition to the six-dimensional phase space (mathbf{r}, mathbf{v}) governed by Kinetic Theory. In this installment, we systematically derive the Vlasov-Maxwell System from the exact microscopic Klimontovich equation and demonstrate how the Liouville theorem guarantees the incompressibility of phase space flow.
1. The Klimontovich Density and the Exact Microscopic Equation
Consider an ensemble of (N) point particles of species (s) with mass (m_s) and charge (q_s). The exact microscopic particle density in 6D phase space ((mathbf{r}, mathbf{v})) is expressed by the Klimontovich distribution:
[
N_s(mathbf{r}, mathbf{v}, t) = sum_{i=1}^N delta(mathbf{r} – mathbf{r}_i(t)) , delta(mathbf{v} – mathbf{v}_i(t))
]
Since particles are neither created nor destroyed, the phase space density satisfies the microscopic continuity equation:
[
rac{partial N_s}{partial t} +
abla_mathbf{r} cdot (mathbf{v} N_s) +
abla_mathbf{v} cdot (mathbf{a}_i N_s) = 0
]
where the acceleration is determined by the exact, rapidly fluctuating microscopic electric and magnetic fields (mathbf{E}^m, mathbf{B}^m):
[
mathbf{a}_i = rac{q_s}{m_s} left( mathbf{E}^m(mathbf{r}_i, t) + mathbf{v}_i imes mathbf{B}^m(mathbf{r}_i, t)
ight)
]
Because (
abla_mathbf{v} cdot mathbf{v} = 0) and (
abla_mathbf{v} cdot (mathbf{v} imes mathbf{B}) = 0), the phase space velocity field (dot{mathbf{X}} = (mathbf{v}, mathbf{a})) is divergence-free. Thus, the exact Klimontovich equation becomes:
[
rac{partial N_s}{partial t} + mathbf{v} cdot
abla_mathbf{r} N_s + rac{q_s}{m_s}left( mathbf{E}^m + mathbf{v} imes mathbf{B}^m
ight) cdot
abla_mathbf{v} N_s = 0
]
2. Ensemble Averaging and Derivation of the Vlasov Equation
To eliminate wild micro-fluctuations on scales smaller than the Debye length, we perform an ensemble average (langle cdot
angle). We split each variable into a smooth mean field and a fluctuation:
[
N_s = f_s(mathbf{r}, mathbf{v}, t) + delta N_s, quad mathbf{E}^m = mathbf{E}(mathbf{r}, t) + delta mathbf{E}, quad mathbf{B}^m = mathbf{B}(mathbf{r}, t) + delta mathbf{B}
]
where (f_s(mathbf{r}, mathbf{v}, t) quiv langle N_s
angle) is the smooth single-particle distribution function normalized such that the number density is (n_s(mathbf{r}, t) = int f_s(mathbf{r}, mathbf{v}, t) , d^3v). Taking the ensemble average of the Klimontovich equation:
[
rac{partial f_s}{partial t} + mathbf{v} cdot
abla_mathbf{r} f_s + rac{q_s}{m_s}left( mathbf{E} + mathbf{v} imes mathbf{B}
ight) cdot
abla_mathbf{v} f_s = -rac{q_s}{m_s} leftlangle (delta mathbf{E} + mathbf{v} imes delta mathbf{B}) cdot
abla_mathbf{v} delta N_s
ight
angle
]
The right-hand side represents the Collision Integral (left(rac{partial f_s}{partial t}
ight)_{ ext{coll}}), which encapsulates discrete binary Coulomb interactions (the BBGKY hierarchy closure). In a collisionless plasma ((N_D gg 1)), discrete particle correlations are suppressed by a factor of (mathcal{O}(1/N_D)). Setting the collision integral to zero yields the fundamental Vlasov (Collisionless Boltzmann) Equation:
[
rac{d f_s}{dt} = rac{partial f_s}{partial t} + mathbf{v} cdot rac{partial f_s}{partial mathbf{r}} + rac{q_s}{m_s} left( mathbf{E} + mathbf{v} imes mathbf{B}
ight) cdot rac{partial f_s}{partial mathbf{v}} = 0
]
3. Self-Consistent Coupling: The Vlasov-Maxwell System
The macroscopic fields (mathbf{E}) and (mathbf{B}) are not external parameters; they are generated self-consistently by the charge and current moments of the distribution functions (f_s):
[
ho(mathbf{r}, t) = sum_s q_s int f_s(mathbf{r}, mathbf{v}, t) , d^3v, quad mathbf{J}(mathbf{r}, t) = sum_s q_s int mathbf{v} f_s(mathbf{r}, mathbf{v}, t) , d^3v
]
These source terms couple directly into Maxwell’s equations:
[
abla cdot mathbf{E} = rac{
ho}{arepsilon_0}, quad
abla imes mathbf{E} = -rac{partial mathbf{B}}{partial t}
]
[
abla cdot mathbf{B} = 0, quad
abla imes mathbf{B} = mu_0 mathbf{J} + mu_0 arepsilon_0 rac{partial mathbf{E}}{partial t}
]
Together, the Vlasov equations for all species coupled to Maxwell’s equations constitute the Vlasov-Maxwell System. Because (df_s/dt = 0), phase space volume is strictly conserved by Liouville’s theorem: the distribution function behaves as an incompressible fluid moving through 6-dimensional phase space.
