In 1946, the Soviet theoretical physicist Lev Landau made one of the most profound and counter-intuitive discoveries in modern theoretical physics: an electrostatic wave in a completely collisionless plasma can attenuate exponentially in time without any thermodynamic entropy production or Joule heating. Known as Landau Damping, this purely kinetic phenomenon is governed by wave-particle phase resonance and analytic continuation along the Bromwich-Landau contour in the complex frequency plane.
1. The Pole Singularity and Laplace Transform Treatment
Consider 1D longitudinal electrostatic oscillations described by the linearized Vlasov equation for electrons with an unperturbed distribution (f_0(v)):
[
rac{partial f_1}{partial t} + v rac{partial f_1}{partial x} – rac{e}{m_e} E_1 rac{partial f_0}{partial v} = 0
]
Fourier transforming in space ((e^{ikx})) and taking the Laplace transform in time (( ilde{f}_1(k, v, s) = int_0^infty f_1(k, v, t) e^{-st} dt)):
[
(s + ikv) ilde{f}_1 – rac{e}{m_e} ilde{E}_1 rac{df_0}{dv} = g(v)
]
where (g(v) = f_1(k, v, t=0)) is the initial perturbation. Solving for ( ilde{f}_1) and substituting into Gauss’s law (ik arepsilon_0 ilde{E}_1 = -e int_{-infty}^infty ilde{f}_1 dv):
[
ilde{E}_1(k, s) = rac{-rac{e}{ikarepsilon_0} int_{-infty}^infty rac{g(v)}{s + ikv} dv}{1 – rac{omega_{pe}^2}{k^2} int_{-infty}^infty rac{df_0/dv}{v – i s/k} dv} quiv rac{N(k, s)}{D(k, s)}
]
Setting (s = -iomega), the denominator defines the plasma dielectric function (D(k, omega)):
[
D(k, omega) = 1 – rac{omega_{pe}^2}{k^2} int_C rac{df_0/dv}{v – omega/k} dv = 0
]
When (omega) has an imaginary part ((omega = omega_r + igamma)), the integral possesses a pole at (v = omega/k). Landau demonstrated that to preserve causality, the integration path (C) along the real (v)-axis must be analytically continued by passing beneath the singularity in the complex velocity plane.
2. The Plemelj Formula and Derivation of the Damping Rate
Applying the Sokhotski-Plemelj theorem to resolve the resonant denominator as (gamma o 0^-):
[
lim_{psilon o 0^+} rac{1}{v – (omega_r/k + ipsilon)} = mathcal{P}left(rac{1}{v – omega_r/k}
ight) + ipi deltaleft(v – rac{omega_r}{k}
ight)
]
where (mathcal{P}) denotes the Cauchy principal value. Decomposing the dielectric function into real and imaginary parts (D = D_r + i D_i):
[
D_r(k, omega_r) = 1 – rac{omega_{pe}^2}{k^2} mathcal{P} int_{-infty}^infty rac{df_0/dv}{v – omega_r/k} dv approx 1 – rac{omega_{pe}^2}{omega_r^2} – rac{3 k^2 v_{th}^2 omega_{pe}^2}{omega_r^4}
]
[
D_i(k, omega_r) = -pi rac{omega_{pe}^2}{k^2} left. rac{df_0}{dv}
ight|_{v = omega_r/k}
]
Assuming weak damping ((|gamma| ll omega_r)), a Taylor expansion about (omega_r) yields (D(omega_r + igamma) approx D_r(omega_r) + igamma rac{partial D_r}{partial omega_r} + i D_i(omega_r) = 0). The Landau damping rate (gamma) is therefore:
[
gamma = -rac{D_i(k, omega_r)}{partial D_r / partial omega_r} = left. rac{pi omega_{pe}^2 omega_r}{2 k^2} rac{df_0}{dv}
ight|_{v = omega_r/k}
]
For a Maxwellian distribution (f_0(v) = rac{1}{sqrt{2pi} v_{th}} xp(-v^2 / 2 v_{th}^2)), the derivative at phase velocity (v_phi = omega_r/k) is strictly negative: (rac{df_0}{dv} = -rac{v}{v_{th}^2} f_0(v) < 0). Evaluating explicitly yields the classical Landau Damping Coefficient:
[
gamma = -sqrt{rac{pi}{8}} rac{omega_{pe}}{(k lambda_D)^3} xpleft( -rac{1}{2 k^2 lambda_D^2} – rac{3}{2}
ight)
]
The physical mechanism is intuitive: particles traveling slightly slower than the wave phase velocity ((v omega/k)) are decelerated, pumping energy into the wave. Because (rac{df_0}{dv} < 0) in a thermal distribution, there are always more slow particles absorbing energy than fast particles giving it up, resulting in net collisionless wave damping.
