Plasma turbulence driven by wave-hole interaction

Technical Summary: Plasma Turbulence Driven by Wave-Hole Interaction

Problem Statement
In turbulent plasmas, nonlinear wave-like oscillations typically drive an energy cascade from large to small scales, leading to dissipation and particle heating. In weakly collisional systems, this cascade extends into the kinetic regime where electromagnetic energy converts to particle energization, often mediated by electrostatic fluctuations. While plasmas are inherently inhomogeneous, and waves interact with coherent structures such as Debye-scale density holes (formed via particle trapping described by Bernstein–Greene–Kruskal theory), the specific dynamics of wave-obstacle interactions and their role in transitioning to phase-space turbulence remain poorly understood. Specifically, it is unclear how these inhomogeneities modify the turbulent cascade and structure formation in both physical and velocity space.

Methodology
The authors investigate the kinetic dynamics of high-frequency electrostatic waves propagating through an inhomogeneous electron plasma using a novel 2D–2V time-splitting Vlasov–Poisson (VP) solver. The study employs normalized units based on the electron Debye length, plasma frequency, and thermal velocity.

• Numerical Setup: The system consists of a homogeneous proton background and an equilibrium electron distribution function (PDF) featuring a radially symmetric density hole ($f_h$). This configuration represents a stationary solution of the 2D–2V VP system.
• Perturbation: The equilibrium is perturbed by unidirectionally propagating electron-acoustic waves (EAWs), a class of fluctuations known for weak Landau damping.
• Simulation Runs: Four distinct numerical runs are performed:
  • Run I: Control simulation without density holes, featuring a large-scale density oscillation.
  • Run II: Includes periodic density holes (acting as obstacles) but no initial large-scale density perturbation ($A_0 = 0$).
  • Run III: Combines density holes with an initial density perturbation ($A_0 = 0.2$).
  • Run IV: Repeats Run III with a doubled domain size in the $y$-direction ($L_y$).
• Analysis Techniques: The study utilizes Fourier-Hermite spectral analysis to characterize turbulence in both fluid and kinetic frameworks. The Hermite transform is used to decompose PDF fluctuations at different velocity scales, isolating non-thermal components and calculating the enstrophy (the second Casimir invariant) to track the energy cascade in velocity space.

Key Results

• Phase-Space Turbulence Transition: The presence of density holes fundamentally alters the turbulent dynamics. While the nonlinear cascade initially develops in the 1D–1V phase subspace ($x-v_x$) driven by the EAWs, the interaction with the holes rapidly spreads the cascade to the full 2D–2V phase space. The holes act as physical obstacles, scattering the strongly anisotropic one-dimensional turbulence and inducing a full 2D–2V enstrophy cascade.
• Wave-Hole Interaction: In physical space, the electron deficit near the holes modifies wavefront velocities, causing waves to slow down and fold around the obstacles. This generates closed structures with small-scale features in both spatial directions. The electron bulk velocity field resembles a laminar wake, with enhanced vorticity localized near the density depletions.
• Velocity Space Structures: The interaction generates complex non-thermal features in velocity space. Specifically, collimated electron beams with well-defined $v_y/v_x$ ratios emerge. These beams are a result of particles being scattered by the lattice of density holes. The opening angle of these beams depends on the lattice geometry; doubling the domain size ($L_y$) increases the vertical spacing between holes, resulting in beams with larger minimum arrival angles.
• Spectral Scaling:
  • Velocity Space: The Hermite spectrum of enstrophy ($\Delta C_2$) shows a dominant cascade along the $v_x$ direction for low modes, consistent with EAW-driven 1D turbulence. However, at intermediate scales ($m \gtrsim 30$), the flow becomes nearly isotropic, following a power law scaling of $m^{-2}$ (shell-integrated $m^{-1}$). This behavior aligns with Batchelor-type theory for enstrophy cascades driven by large-scale electric field fluctuations.
  • Physical Space: The Fourier spectrum of density fluctuations exhibits an approximately isotropic flux due to the hole lattice distributing the cascade along both $x$ and $y$ directions. At intermediate scales ($k \gtrsim 1$), the density spectrum follows a $k^{-3}$ scaling, and the electric field spectrum follows a $k^{-5}$ scaling, consistent with Batchelor theory for 2D–2V VP turbulence.
• Heating Mechanism: The enstrophy density ($\Omega$) exhibits stronger gradients along the wave propagation direction ($x$) and correlates spatially with the electron temperature ($T_e$). Regions farther from equilibrium (higher deviation from Maxwellian) show stronger heating, suggesting that the enstrophy cascade acts as a key heating mechanism in weakly collisional plasmas.

Significance and Claims
The paper claims to be the first to reproduce and analyze a reticulated pattern of stable electron density gaps in 2D–2V phase space with the explicit purpose of identifying turbulent cascading mechanisms in kinetic plasma regimes. The authors propose a new paradigm for energy transfer in weakly collisional plasmas, where medium inhomogeneity induces strong distortion of wavefronts, enhancing wave coupling and promoting turbulence development even orthogonally to the primary wave propagation.

The study suggests that the continuous wave-obstacle interaction sustains a kinetic cascade that efficiently transfers fluctuations toward dissipative scales, providing a pathway for irreversible plasma heating and entropy production. These findings have direct implications for understanding space plasmas, such as the solar wind, where small-scale inhomogeneities may modulate the propagation of Langmuir and other high-frequency waves, potentially leaving statistical signatures of wave modulation and particle heating detectable by spacecraft. The authors note that future work could explore the role of specific inhomogeneity parameters (hole size and separation) and extend the approach to magnetic kinetic turbulence.