Source-linked AI summary
Revisiting kinetic electrostatic electron non-linear (KEEN) waves in the presence of dynamical ions
Rostislav-Paul Wilhelm, Philipp Krah, Kai Schneider, Fabio Bacchini, Virginie Grandgirard
TL;DR
Conventional grid solvers struggle with the fine phase-space structures and long timescales of KEEN waves, especially when ion dynamics are included. The paper extends CMM-NuFI with separate electron and ion maps coupled through Poisson, finding that ion motion leaves early electron dynamics essentially unchanged but produces later divergence while remaining computationally tractable.
Problem
Fine KEEN phase-space structures and long evolution times make conventional grid-based simulations prohibitively expensive, limiting extensions to dynamical ions.
Method
The paper generalizes CMM-NuFI to multiple species by maintaining independent characteristic submaps for electrons and ions coupled through a shared Poisson solve.
Results
Electron dynamics are essentially unaffected by ion motion during and shortly after the drive, but later ion feedback produces additional electron vortices, a shifted primary vortex, and drifting ρ-harmonics.
Takeaways & Limitations
The multi-map scheme makes long-time, fully kinetic KEEN simulations with mobile ions at realistic mass ratios computationally tractable.
Takeaways & Limitations
The ponderomotive-force-as-potential simplification is exact for hydrogen ions and the ion-static limit but requires revision for heavier or more highly charged ions.
Abstract
from arXiv · showhide
We revisit the kinetic electrostatic electron nonlinear (KEEN) waves studied by Afeyan et. al in 2014 using a hybrid flow-mapping strategy that combines the characteristic mapping method (CMM) with numerical flow iteration (NuFi). The study extends the classical setup to dynamical ions, comparing their impact on the long-time KEEN dynamics with the static-ion case. To this end, we extend the CMM-NuFI framework with a multi-map strategy, assigning one map to the ion and one to the electron characteristic flow. The resulting problem exhibits a wide separation of spatial and temporal scales, driven by the fine structures generated by the ponderomotive force and by the large ion-to-electron mass ratio, which renders it computationally prohibitive for conventional grid-based methods. Our multi-map CMM-NuFI method efficiently resolves these disparate scales, enabling long-time simulations of KEEN dynamics with fully dynamical ions.
1 Introduction
KEEN waves are driven, self-organized electron phase-space structures whose fine-scale, long-lived dynamics make conventional grid simulations costly. The paper addresses this challenge by extending flow-map methods to dynamical ions.
- KEEN wave background: KEEN waves are self-organized, non-stationary electron phase-space structures in the driven, ion-static Vlasov–Poisson system.They are excited by the ponderomotive force of two crossing laser beams and are not linear-dispersion or BGK equilibria.
- KEEN wave background: Strong, sufficiently long driving produces phase-locked, multi-harmonic KEEN structures that persist after the external drive is switched off.Unlike electron-acoustic waves or BGK modes, KEEN waves can form anywhere in the (ω, k) plane under the stated drive conditions.
- Numerical challenge: Small phase-space vorticlets shed, merge, and retrap around the driven phase velocity, generating filaments whose width scales with the square root of local field amplitude.These fine structures create stringent resolution demands for long-time simulations.
- Numerical challenge: Grid-based simulations require up to 8192 × 16384 points and 10^4–10^5 CPU-hours per run, while low-order ρ modes remain unconverged at those resolutions.The non-uniform grid also introduces visible subshocks in transition regions between cell sizes.
- Flow-map approach: The hybrid CMM-NuFi method combines NuFi’s conservation and resolution properties with CMM’s reduced cost growth by evolving characteristic-map compositions instead of storing f on a grid.NuFi iterates backward characteristic compositions, whereas CMM stores low-order polynomial submaps spanning many elementary time steps.
- Paper contribution: The paper extends the canonical KEEN setup to dynamical ions and generalizes CMM-NuFi with independent characteristic submaps for each species coupled through shared density and field calculations.The study compares dynamical and static ions using electron density, ρ-harmonics, and energy evolution.
2 KEEN waves with dynamical ions
The canonical model is extended to electrons and dynamical ions coupled through the total charge density and electric field. For hydrogen ions, the ponderomotive drive acts much more weakly on ions because of the small electron-to-ion mass ratio.
- Two-species model: The two-species Vlasov–Poisson model contains electron and ion distribution functions with charges, masses, charge number Zi, and mass ratio µ = me/mi.The equations are normalized using the electron Debye length, electron plasma frequency, and electron thermal speed.
- Two-species model: Electrons and ions are coupled through the total charge density and the shared self-consistent electric field.The species contribute their densities to the charge density entering the Poisson equation.
- Drive and scaling: For ions, the ponderomotive acceleration is scaled by qi/mi = −µZi, making the direct drive weaker whenever µZi ≪1.The opposite sign reflects opposite charges, while the small magnitude reflects the heavier ion mass.
- Ion-static limit: The ion-static system is recovered as µ →0, equivalently qi/mi →0, with ni ≡1 and Zi = 1.In this limit, ions remain a fixed, uniform neutralizing background.
- Modeling caveat: Writing the ponderomotive force as a potential contribution to the electric field is exact for hydrogen ions and the ion-static limit, but not for heavier or more highly charged ions.Those species would require a distinct external-force treatment and a revised coupling in CMM-NuFi.
- Canonical setup: For hydrogen ions with Zi = 1 and µ = 1/1836, the ion velocity domain is approximately 43 times narrower than the electron domain.This reflects the physical scale separation that the method must resolve for both species.
3 Multi-map CMM-NuFi for multiple species
The multi-species CMM-NuFi method maintains separate characteristic-map representations for each species while sharing density accumulation and the Poisson solve. This preserves the flow-map framework while accommodating species-dependent scales and charge-to-mass ratios.
- Species evolution: Each species applies the same NuFi and map-composition procedures, with species dependence entering through its charge-to-mass ratio.The electron and ion chains are advanced independently rather than through a shared flow map.
- Species coupling: The species are coupled only when their independently traced densities are accumulated and inserted into the shared Poisson solve for the electric field.The field history is then used by both species in the next time step.
- Time stepping: The algorithm samples each species’ backward maps to accumulate its density before solving Poisson for the updated electric field.The density accumulation and field solve are the shared portions of the time step.
- Multi-map representation: The multi-species method keeps one independent list of characteristic submaps, NuFi counter, and coarse map grid for each species.Species may use different remapping frequencies and map-grid resolutions in principle.
- Time stepping: At each time step, one NuFi step advances every species, and after Nremap iterations its chain is collapsed into a new coarse submap.The submap remapping follows the single-species hybrid algorithm.
- Scale adaptation: Slowly evolving ions can in principle use coarser, less frequently updated maps because their submaps remain well approximated over electron-scale remapping intervals.The reported section-4 simulations nevertheless use identical parameters for both species.
4 Results
The dynamical-ion and ion-static runs agree during the drive and shortly afterward, but diverge on longer ion timescales. Ion motion produces late-time changes in electron phase-space structure and ρ-harmonics while contributing little to the energy budget.
- Simulation setup: The multi-species CMM-NuFi simulation evolves the canonical KEEN wave to t = 2000 with independently maintained electron and ion characteristic maps.The run uses τ = 0.05, remapping every Δt = 20, and separate maps for the two species.
- Electron dynamics: At t = 1000ω^-1, ion dynamics shift the main electron vortex and generate additional small-scale phase-space vortices at larger v.These structures suggest that including ion dynamics produces more particle traps.
- Ion and energy evolution: The ion distribution remains close to a Maxwellian with a small-amplitude phase-locked perturbation, while electron kinetic energy stabilizes about 7% above its initial level.Ion kinetic energy and self-induced electric-field energy remain near negligible in the reported energy evolution.
- ρ-harmonics: The five largest electron-density ρ-harmonics agree during the drive but show a slow relative drift at late times between ion-static and ion-dynamic runs.The drift quantifies the cumulative effect of the ion response on the sustained KEEN structure.
5 Conclusion
The multi-map CMM-NuFi scheme extends KEEN-wave simulations to fully kinetic dynamical ions at realistic mass ratios. Ion motion leaves early electron dynamics largely unchanged but produces measurable cumulative effects at later times.
- The multi-map scheme maintains independent characteristic submaps for each species, coupled through a shared Poisson solve, to simulate two-species Vlasov–Poisson dynamics.
- Ion motion scarcely affects the electron dynamics during and shortly after the drive, but later produces additional small-scale vortices, a shifted primary vortex, and drifting ρ-harmonics.
- The ion distribution remains near Maxwellian, while the energy budget shows that the ponderomotive force drives mainly the electrons and ion kinetic and field energies contribute marginally.
- Although ion mobility does not qualitatively alter the canonical KEEN wave on these time scales, it leaves a measurable cumulative influence on electron phase space.