Relaxation of weakly collisional plasma: Continuous spectra, discrete eigenmodes, and the decay of echoes

The relaxation of a weakly collisional plasma, which is of fundamental interest to laboratory astrophysical plasmas, can be described by the self-consistent Boltzmann-Poisson equations with the Lenard-Bernstein collision operator. We perform a perturbative (linear and second-order) analysis of the Boltzmann-Poisson equations and obtain exact analytic solutions which resolve some longstanding controversies regarding the impact of weak collisions on the continuous spectra, the discrete Landau eigenmodes, and the decay of plasma echoes. We retain both damping and diffusion terms in the collision operator throughout our treatment. We find that the linear response is a temporal convolution of two types of contribution: a continuum that depends on the continuous velocities of particles (crucial for the plasma echo), and another, consisting of discrete modes that are coherent modes of oscillation of the entire system. The discrete modes are exponentially damped over time due to collective effects or wave-particle interactions (Landau damping), as well as collisional dissipation. The continuum is also damped by collisions but somewhat differently than the discrete modes. Up to a collision time, which is the inverse of the collision frequency ν_{c}, the continuum decay is driven by the diffusion of particle velocities and is cubic exponential, occurring over a timescale ∼ν_{c}^{-1/3}. After a collision time, however, the continuum decay is driven by the collisional damping of particle velocities and diffusion of their positions and occurs exponentially over a timescale ∼ν_{c}. This slow exponential decay causes perturbations to damp the most on scales comparable to the mean free path but very slowly on larger scales. This establishes the local thermal equilibrium, which is the essence of the fluid limit. The long-term decay of the linear response is driven by the discrete modes on scales smaller than the mean free path but, on larger scales, is governed by a combination of the slowly decaying continuum and the least damped discrete mode. This slow exponential decay implies that the echo, which results from the interference of the continuum response to two subsequent pulses, is detectable even on scales comparable to the mean free path, as long as the second pulse is introduced within a few phase-mixing timescales after the first.