Kinetic Theory of Plasmas, Vlasov Equation & Landau Damping

1. Phase Space Representation

When particles possess a spread in velocities that cannot be captured by a single mean fluid velocity $\vec{u}(\vec{r}, t)$, fluid theory fails. We must describe the system microscopically in a six-dimensional phase space $(\vec{r}, \vec{v}) = (x, y, z, v_x, v_y, v_z)$.

The state of species $\alpha$ is defined by its distribution function $f_\alpha(\vec{r}, \vec{v}, t)$:

$$dN_\alpha = f_\alpha(\vec{r}, \vec{v}, t) d^3r d^3v = f_\alpha(\vec{r}, \vec{v}, t) dx dy dz dv_x dv_y dv_z$$

where $dN_\alpha$ is the expected number of particles located in the spatial volume element $d^3r$ around $\vec{r}$ having velocities within the velocity volume element $d^3v$ around $\vec{v}$ at time $t$.

The distribution function must be strictly non-negative: $f_\alpha(\vec{r}, \vec{v}, t) \ge 0$, and integrating over all velocity space recovers the local macroscopic number density:

$$\int_{-\infty}^\infty \int_{-\infty}^\infty \int_{-\infty}^\infty f_\alpha(\vec{r}, \vec{v}, t) dv_x dv_y dv_z = n_\alpha(\vec{r}, t)$$

1. Liouville's Theorem in Hamiltonian Phase Space

In Hamiltonian mechanics, particle trajectories are governed by Hamilton's canonical equations $\dot{\vec{r}} = \nabla_p H, \dot{\vec{p}} = -\nabla_r H$. The phase space velocity field $\vec{w} = (\dot{\vec{r}}, \dot{\vec{p}})$ is strictly divergence-free:

$$\nabla_r \cdot \dot{\vec{r}} + \nabla_p \cdot \dot{\vec{p}} = \nabla_r \cdot (\nabla_p H) - \nabla_p \cdot (\nabla_r H) = 0$$

According to Liouville's Theorem, the flow of representative points in phase space behaves as an incompressible fluid:

$$\frac{df}{dt} = 0$$

The phase space density $f$ remains constant along any dynamical trajectory!

2. The Boltzmann Transport Equation

Expanding the total convective time derivative along a trajectory governed by acceleration $\vec{a} = \frac{q}{m}(\vec{E} + \vec{v}\times\vec{B})$:

$$\frac{df}{dt} = \frac{\partial f}{\partial t} + \dot{\vec{r}}\cdot\nabla_r f + \dot{\vec{v}}\cdot\nabla_v f = \frac{\partial f}{\partial t} + \vec{v}\cdot\nabla f + \frac{q}{m}(\vec{E} + \vec{v}\times\vec{B})\cdot\nabla_v f$$

In the presence of discrete, short-range binary collisions between particles:

$$\frac{\partial f}{\partial t} + \vec{v}\cdot\nabla f + \frac{q}{m}(\vec{E} + \vec{v}\times\vec{B})\cdot\nabla_v f = \left( \frac{\partial f}{\partial t} \right)_{\text{coll}}$$

This is the Boltzmann transport equation.

1. The Vlasov Equation

In high-temperature, low-density plasmas where the plasma parameter $N_D = \frac{4}{3}\pi n_e \lambda_D^3 \gg 1$, long-range collective Coulomb fields dominate over discrete binary collisions. Setting the collision integral to zero yields the Vlasov equation (or collisionless Boltzmann equation):

$$\frac{\partial f_\alpha}{\partial t} + \vec{v}\cdot\nabla f_\alpha + \frac{q_\alpha}{m_\alpha}\left( \vec{E}(\vec{r}, t) + \vec{v}\times\vec{B}(\vec{r}, t) \right)\cdot\frac{\partial f_\alpha}{\partial \vec{v}} = 0$$

Here, $\vec{E}$ and $\vec{B}$ are not external fields, but the self-consistent, macroscopic average electromagnetic fields generated by the collective charge density $\rho_q$ and current density $\vec{J}$ of all particles in the plasma:

$$\rho_q(\vec{r}, t) = \sum_\alpha q_\alpha \int f_\alpha(\vec{r}, \vec{v}, t) d^3v, \quad \vec{J}(\vec{r}, t) = \sum_\alpha q_\alpha \int \vec{v} f_\alpha(\vec{r}, \vec{v}, t) d^3v$$

2. The Self-Consistent Vlasov-Maxwell System

Coupled with Maxwell's equations:

$$\nabla \cdot \vec{E} = \frac{\rho_q}{\varepsilon_0}, \quad \nabla \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}$$ $$\nabla \cdot \vec{B} = 0, \quad \nabla \times \vec{B} = \mu_0 \vec{J} + \frac{1}{c^2}\frac{\partial \vec{E}}{\partial t}$$

The Vlasov-Maxwell system provides a complete, exact, collisionless description of plasma physics, capturing all kinetic phenomena including wave-particle resonances and Landau damping.

1. The General Moment Theorem

Let $\chi(\vec{v})$ be an arbitrary function of velocity. Multiplying the Vlasov equation by $\chi(\vec{v})$ and integrating over all velocity space $d^3v$:

$$\int \chi(\vec{v}) \frac{\partial f}{\partial t} d^3v + \int \chi(\vec{v}) \vec{v}\cdot\nabla f d^3v + \frac{q}{m}\int \chi(\vec{v}) (\vec{E} + \vec{v}\times\vec{B})\cdot\nabla_v f d^3v = 0$$

Using integration by parts on the third term and noting that $f \to 0$ as $v \to \pm\infty$:

$$\frac{\partial}{\partial t}\left( n \langle \chi \rangle \right) + \nabla \cdot \left( n \langle \vec{v} \chi \rangle \right) - \frac{q}{m} n \left\langle (\vec{E} + \vec{v}\times\vec{B})\cdot\nabla_v \chi \right\rangle = 0$$

2. Systematic Derivation of Fluid Equations

1. Zeroth Moment ($\chi = 1$): Since $\nabla_v(1) = 0$: $$\frac{\partial n}{\partial t} + \nabla \cdot (n \vec{u}) = 0 \quad \text{(Continuity Equation)}$$
2. First Moment ($\chi = m \vec{v}$): Since $\nabla_v(m\vec{v}) = m \mathbf{I}$: $$m n \left[ \frac{\partial \vec{u}}{\partial t} + (\vec{u}\cdot\nabla)\vec{u} \right] = q n (\vec{E} + \vec{u}\times\vec{B}) - \nabla \cdot \mathbf{P} \quad \text{(Momentum Equation)}$$
3. Second Moment ($\chi = \frac{1}{2}m v^2$): Yields the thermal energy transport equation containing the heat flux vector $\vec{Q} = \frac{1}{2}m \int (\vec{v} - \vec{u})|\vec{v} - \vec{u}|^2 f d^3v$.

1. Linearization of 1D Vlasov Equation

Consider 1D electrostatic perturbations in an unmagnetized plasma with stationary ion background. Decompose the electron distribution function:

$$f(x, v, t) = f_0(v) + f_1(x, v, t)$$

where $f_0(v)$ is an unperturbed homogeneous equilibrium (e.g., Maxwellian) and $f_1$ is small ($|f_1| \ll f_0$). Neglecting the second-order term $E_1 \frac{\partial f_1}{\partial v}$:

$$\frac{\partial f_1}{\partial t} + v \frac{\partial f_1}{\partial x} - \frac{e}{m_e} E_1 \frac{\partial f_0}{\partial v} = 0$$

coupled to Poisson's equation:

$$\frac{\partial E_1}{\partial x} = -\frac{e}{\varepsilon_0}\int_{-\infty}^\infty f_1(x, v, t) dv$$

2. Failure of the Standard Fourier Transform

If one attempts to solve this initial-value problem using a standard Fourier transform in time ($e^{-i\omega t}$):

$$-i(\omega - k v) f_1 = \frac{e}{m_e} E_1 \frac{\partial f_0}{\partial v} \implies f_1(k, v, \omega) = \frac{i e E_1}{m_e} \frac{\partial f_0 / \partial v}{\omega - k v}$$

Substituting $f_1$ into Poisson's equation yields the dielectric permittivity integral:

$$\varepsilon(k, \omega) = 1 + \frac{e^2}{\varepsilon_0 m_e k} \int_{-\infty}^\infty \frac{\partial f_0 / \partial v}{\omega - k v} dv = 1 - \frac{\omega_{pe}^2}{k^2} \int_{-\infty}^\infty \frac{f_0'(v)}{v - \omega/k} dv$$

Here lies the profound crisis that perplexed early plasma physicists: if $\omega$ is real, the integrand has a singularity pole on the real axis at $v = \omega / k$ (the phase velocity). Standard Fourier integration fails because it cannot handle the initial conditions or resolve the singularity properly.

1. Landau's Laplace Transform Method

In 1946, Soviet physicist Lev Landau resolved the singularity by treating the perturbation as an initial-value problem with a one-sided Laplace transform in time:

$$\tilde{f}_1(k, v, p) = \int_0^\infty f_1(k, v, t) e^{-p t} dt, \quad \text{Re}(p) > 0$$

Inverting the Laplace transform requires integrating along the Bromwich contour in the complex $p$-plane:

$$E_1(k, t) = \frac{1}{2\pi i} \int_{\sigma - i\infty}^{\sigma + i\infty} \tilde{E}_1(k, p) e^{p t} dp$$

To perform the analytic continuation into the damping half-plane ($\text{Re}(p) < 0$), Landau proved that the integration contour in velocity space must be deformed so that the pole at $v = \omega / k = -i p / k$ always remains above the contour. This defines the famous Landau contour $C$:

$$\int_C \frac{f_0'(v)}{v - \omega/k} dv = \mathcal{P}\int_{-\infty}^\infty \frac{f_0'(v)}{v - \omega/k} dv + i \pi f_0'\left( \frac{\omega}{k} \right)$$

where $\mathcal{P}$ denotes the Cauchy principal value, and the term $+i\pi f_0'(\omega/k)$ is the residue from wrapping around the pole from below.

2. The Landau Damping Rate

Writing the complex frequency as $\omega = \omega_r + i\gamma$, where $\gamma$ is the damping rate ($|\gamma| \ll \omega_r$). Expanding $\varepsilon(k, \omega_r + i\gamma) = 0$ in Taylor series:

$$\varepsilon_r(k, \omega_r) + i\gamma \frac{\partial \varepsilon_r}{\partial \omega} + i\varepsilon_i(k, \omega_r) = 0$$

The real part yields the Bohm-Gross frequency: $\omega_r^2 \approx \omega_{pe}^2 + 3 k^2 v_{\text{th}}^2$. The imaginary part yields the Landau damping rate $\gamma_L$:

$$\gamma_L = -\frac{\varepsilon_i}{\partial \varepsilon_r / \partial \omega} = \frac{\pi}{2} \frac{\omega_r \omega_{pe}^2}{k^2} \left[ \frac{\partial f_0}{\partial v} \right]_{v = \omega_r/k}$$

For a Maxwellian distribution $f_0(v) = \frac{1}{\sqrt{2\pi}v_{\text{th}}}e^{-v^2 / (2v_{\text{th}}^2)}$, the derivative at the phase velocity is strictly negative: $\left[ \frac{\partial f_0}{\partial v} \right]_{v_{ph}} < 0$. Consequently:

$$\gamma_L < 0$$

The wave is exponentially damped in time:

$$E_1(t) \propto e^{-|\gamma_L| t} \cos(\omega_r t)$$

Evaluating explicitly for a Maxwellian distribution:

$$\gamma_L = -\sqrt{\frac{\pi}{8}} \frac{\omega_{pe}}{(k\lambda_D)^3} \exp\left( -\frac{1}{2(k\lambda_D)^2} - \frac{3}{2} \right)$$

This is the discovery of Landau damping: irreversible wave dissipation occurring in a completely collisionless, conservative Hamiltonian plasma!

1. Physical Mechanism: Surfing the Wave

Landau damping is not caused by dissipative collisions, but by a resonant energy exchange between the electrostatic wave and particles moving near the wave's phase velocity:

$$v \approx v_{ph} = \frac{\omega}{k}$$

Consider a surfer catching an ocean wave:

• Particles slightly slower than the wave ($v < v_{ph}$): The wave pushes them forward, accelerating them and transferring energy from the wave to the particles.
• Particles slightly faster than the wave ($v > v_{ph}$): They push forward against the wave potential, decelerating and transferring energy from the particles to the wave.

In a thermal Maxwellian distribution, the distribution function decreases monotonically with speed ($\partial f_0 / \partial v < 0$). Therefore, there are always more slightly slower particles than slightly faster particles:

$$N(v < v_{ph}) > N(v > v_{ph})$$

On net, more particles gain kinetic energy from the wave than give energy to it. As a result, the wave loses electrostatic energy, damping exponentially away!

2. The Two-Stream Instability (Inverse Landau Damping)

If the plasma contains an electron beam moving at drift speed $v_d > 0$, the distribution function develops a bump on the tail. In the velocity range where the distribution function slope is positive:

$$\left[ \frac{\partial f_0}{\partial v} \right]_{v = v_{ph}} > 0$$

Now there are more faster particles than slower particles! The resonant particles transfer net kinetic energy to the wave, causing the wave amplitude to grow exponentially ($\gamma > 0$). This is inverse Landau damping, which drives the classic two-stream instability (or bump-on-tail instability).