Universal non-thermal power-law distribution functions from the self-consistent evolution of collisionless electrostatic plasmas

TL;DR

This paper demonstrates that the universal v^{-5} power-law tail in collisionless plasma distribution functions naturally arises from self-consistent electrostatic relaxation, providing a theoretical explanation for observed non-thermal tails in space plasmas.

Contribution

It introduces a self-consistent quasilinear framework showing how collisionless plasmas develop universal power-law tails, specifically v^{-5}, from electrostatic relaxation processes.

Findings

The v^{-5} tail emerges in 3D velocity distributions.

Energy distribution develops an E^{-2} tail.

Self-consistency is essential for tail formation.

Abstract

Distribution functions of collisionless systems are known to show non-thermal power law tails. Interestingly, collisionless plasmas in various physical scenarios, (e.g., the ion population of the solar wind) feature a $v^{-5}$ tail in the velocity ($v$) distribution, whose origin has been a long-standing mystery. We show this power law tail to be a natural outcome of the self-consistent collisionless relaxation of driven electrostatic plasmas. We perform a quasilinear analysis of the perturbed Vlasov-Poisson equations to show that the coarse-grained mean distribution function (DF), $f_0$, follows a quasilinear diffusion equation with a diffusion coefficient $D(v)$ that depends on $v$ through the plasma dielectric constant. If the plasma is isotropically forced on scales much larger than the Debye length with a white noise-like electric field, then $D(v)\sim v^4$ for $\sigma<v<\omega_{\mathrm{P}}/k$, with $\sigma$ the thermal velocity, $\omega_{\mathrm{P}}$ the plasma frequency and $k$ the maximum wavenumber of the perturbation; the corresponding $f_0$, in the quasi-steady state, develops a $v^{-\left(d+2\right)}$ tail in $d$ dimensions ($v^{-5}$ tail in 3D), while the energy ($E$) distribution develops an $E^{-2}$ tail irrespective of the dimensionality of space. Any redness of the noise only alters the scaling in the high $v$ end. Non-resonant particles moving slower than the phase-velocity of the plasma waves ($\omega_{\mathrm{P}}/k$) experience a Debye-screened electric field, and significantly less (power law suppressed) acceleration than the near-resonant particles. Thus, a Maxwellian DF develops a power law tail. The Maxwellian core ($v<\sigma$) eventually also heats up, but over a much longer timescale than that over which the tail forms. We definitively show that self-consistency (ignored in test-particle treatments) is crucial for the development of the universal $v^{-5}$ tail.