Physics of Circular Polarized Ion-Scale Waves in Hybrid Simulations of Alfvénic Fluctuations

Abstract: Ion cyclotron waves (ICW) and fast magnetosonic/whistler waves (FMW) are fundamental electromagnetic modes at ion kinetic scales, yet their generation mechanisms and roles in plasma evolution remain poorly understood. We analyze a 2.5D hybrid simulation of broadband Alfvénic fluctuations, where the proton velocity distribution is modeled as a sum of two bi-Maxwellian components: a thermal core and a drifting beam. Using wavelet-based wave identification, bi-Maxwellian VDF fitting, and the PLUME linear dispersion solver, we find that ICW behave as linear modes. Growth is intermittent, occurring when core temperature anisotropy builds up, and is driven mainly by the core (the beam contributes negligibly). Poynting flux analysis shows that ICW are predominantly forward-propagating, with a net energy flux ratio of $+1$ across all frequencies, consistent with the initial condition. FMW present a stark contrast: PLUME solutions often yield very small (near-zero) linear growth/damping rates. The species decomposition breaks down when $|γ/ω_r| \gtrsim 0.368$, indicating that linear theory predicts these waves to be strongly damped and not describable by linear eigenmodes. Nevertheless, FMW are clearly observed in the wavelet helicity spectrogram, indicating that they are generated by nonlinear processes (e.g., parametric decay or phase steepening) and persist despite linear damping. The net energy flux ratio for FMW is close to $+1$ at low frequencies but decreases at higher frequencies, yet never reaches zero (net energy flow remains forward). These results demonstrate that ICW are linear, core-driven waves that transfer energy to the plasma, while FMW are heavily damped, nonlinearly generated waves.

• The paper combines wavelet polarization, proton distribution fits, and PLUME dispersion analysis to distinguish forward-propagating ion cyclotron waves (0.03–0.52 Ωp) from fast magnetosonic/whistler waves (0.09–2 Ωp).
• The paper finds that ion cyclotron waves are intermittent linear eigenmodes driven almost entirely by core temperature anisotropy, with growth reaching approximately 0.05 Ωp near kvA/Ωp≈0.65 before relaxing toward marginal stability.
• The paper shows that fast magnetosonic/whistler waves remain forward-propagating despite dominant linear damping, supporting nonlinear generation through processes such as parametric decay or phase steepening rather than local instability.

Overview

This paper analyzes the physics of circularly polarized ion-scale waves generated during the nonlinear evolution of broadband Alfvénic fluctuations, using a 2.5D hybrid particle-in-cell simulation performed with the CAMELIA code (2608.14151). The central question is whether ion cyclotron waves (ICW) and fast magnetosonic/whistler waves (FMW) observed at proton kinetic scales behave as linear eigenmodes of the local plasma state or arise from nonlinear processes. The authors combine three diagnostics applied to time series recorded at 16 fixed spatial probes: wavelet-based polarization analysis of the magnetic and electric fields, bi-Maxwellian (core-plus-beam) fits to the proton velocity distribution functions (VDFs), and linear dispersion solutions from the PLUME solver with species-resolved growth rates.

The simulation initializes a homogeneous, isotropic plasma ($\beta_{p,e}=0.5$, 8000 particles per cell) in a $128\,d_i$ domain and imposes an outward-propagating, left-handed circularly polarized Alfvénic pump with spectrum $E_B(k_\parallel)\propto k_\parallel^{-2}$ over $k_\parallel d_i\in[0.049,0.490]$. Fields are sampled at cadence $0.1\,\Omega_p^{-1}$ for $300\,\Omega_p^{-1}$.

Wave identification via wavelet polarization analysis

The authors transform fields into the plasma rest frame using the time-averaged bulk velocity, apply Morlet wavelet transforms, and compute reduced helicities $\sigma_B$ and $\sigma_E$ analogous to Stokes parameters. A wave event requires $|\sigma_B|$ and $|\sigma_E|>0.7$ (extending while $128\,d_i$0), a phase-speed fit with $128\,d_i$1, and duration of at least one wavelet e-folding length. Because the full electric field is available in the plasma frame, the analysis avoids the Doppler ambiguity that afflicts spacecraft measurements, where a left-handed fluctuation could be either an ICW or a Doppler-shifted FMW.

Two coherent populations emerge: ICW occupy $128\,d_i$2 and FMW occupy $128\,d_i$3, both spanning from the pump frequency to ion-kinetic scales. Dispersion points derived from measured phase speeds align well with the cold-plasma branches—left-handed points on the ICW branch, right-handed points on the FMW branch—validating the identification methodology.

Poynting flux analysis yields a net energy flux ratio $128\,d_i$4 that is consistently $128\,d_i$5 for ICW across all frequencies, i.e., purely forward-propagating, consistent with the initial outward pump. For FMW, $128\,d_i$6 at low frequencies ($128\,d_i$7) but decreases at higher frequencies, reaching a minimum of $128\,d_i$8 at $128\,d_i$9 without ever becoming zero or negative; net energy flow remains forward at all scales.

Core-beam decomposition of the proton VDFs

Each probe VDF is fit with two drifting bi-Maxwellians (core and beam) via constrained SLSQP minimization, with sequential initialization from the previous timestep. Fits fail or produce nonphysical parameters in fewer than about 5% of timesteps. Representative fits show beam densities of $E_B(k_\parallel)\propto k_\parallel^{-2}$0–$E_B(k_\parallel)\propto k_\parallel^{-2}$1, beam drifts $E_B(k_\parallel)\propto k_\parallel^{-2}$2–$E_B(k_\parallel)\propto k_\parallel^{-2}$3, and beam temperature anisotropies $E_B(k_\parallel)\propto k_\parallel^{-2}$4 ranging from 1.20 to 10.4 at the three highlighted times. The authors are explicit that the apparent constancy of core density is an artifact of normalization rather than physical invariance, and that spurious spikes in $E_B(k_\parallel)\propto k_\parallel^{-2}$5 occur when the beam is tenuous or cold parallel to the field, making those anisotropy values unreliable where the beam is poorly resolved.

Linear stability: ICW as core-driven eigenmodes

The fitted parameters feed PLUME, which solves the hot-plasma Vlasov–Maxwell dispersion relation with nearly parallel propagation ($E_B(k_\parallel)\propto k_\parallel^{-2}$6) and returns species-resolved growth/damping rates valid for $E_B(k_\parallel)\propto k_\parallel^{-2}$7. Growth is deemed dynamically significant when $E_B(k_\parallel)\propto k_\parallel^{-2}$8 over the simulation timescale.

For the ICW branch, positive total growth occurs only intermittently (e.g., near $E_B(k_\parallel)\propto k_\parallel^{-2}$9, 77, and $k_\parallel d_i\in[0.049,0.490]$0), with damping dominating elsewhere. Each growth episode follows a systematic buildup of core anisotropy $k_\parallel d_i\in[0.049,0.490]$1—attributed to phase steepening of the Alfvénic fluctuations—and during growth $k_\parallel d_i\in[0.049,0.490]$2 decreases, consistent with cyclotron-resonant extraction of free energy, followed by relaxation toward marginal stability. Species decomposition shows the core contributes essentially all of the growth; the beam contribution is negligible (typically under 5%) despite its substantial drift speed. This establishes ICW as genuine linear modes driven by the core population, not by the beam. In wavenumber space, most detected ICW events cluster near $k_\parallel d_i\in[0.049,0.490]$3, while the densest positive growth rates ($k_\parallel d_i\in[0.049,0.490]$4) sit near $k_\parallel d_i\in[0.049,0.490]$5; the offset reflects stronger damping at higher $k_\parallel d_i\in[0.049,0.490]$6, so observed waves preferentially populate the weakly damped portion of the spectrum. The first $k_\parallel d_i\in[0.049,0.490]$7 of PLUME output is excluded as an initial-transient artifact.

FMW: strong damping and nonlinear generation

The FMW branch behaves qualitatively differently. PLUME returns very small or near-zero growth rates at most times and frequencies, yet FMW are unambiguously present in the helicity spectrogram. Two independent lines of evidence indicate that linear theory fails for these waves:

• Breakdown of species decomposition: the sum $k_\parallel d_i\in[0.049,0.490]$8 ceases to equal $k_\parallel d_i\in[0.049,0.490]$9 whenever $0.1\,\Omega_p^{-1}$0, beyond which the anti-Hermitian susceptibility decomposition underlying species-resolved rates is invalid. FMW rates frequently exceed this threshold.
• No correspondence between growth and detection: across the entire $0.1\,\Omega_p^{-1}$1 range, damping dominates, and the wavenumber where FMW events peak shows no concentration of positive growth rates.

The authors conclude that the observed FMW are strongly damped structures generated by nonlinear processes—parametric decay or phase steepening producing forced compressible perturbations of the fast type—and persist despite linear damping. A consequence is that the driving species for FMW cannot be determined from this linear analysis. Notably, the observed FMW still track the cold-plasma dispersion curve $0.1\,\Omega_p^{-1}$2 even under strong damping, which the authors flag as curious and unresolved.

Limitations and open questions

Several caveats bear directly on the results. The bi-Maxwellian core-plus-beam model may not capture the full complexity of the evolving distributions, particularly for the FMW case; the authors suggest that a generalized dispersion treatment such as ALPS could resolve some discrepancies. Beam anisotropy values are unreliable where the beam is poorly resolved, and the first quarter of the simulation is excluded from quantitative stability analysis. The persistence of FMW along the cold-plasma dispersion relation despite strong damping remains unexplained. More broadly, the study leaves open how these nonlinearly generated, heavily damped FMW contribute quantitatively to proton heating relative to the core-driven ICW channel.

Conclusion

By coupling wavelet-based wave identification, bi-Maxwellian VDF fitting, and PLUME stability analysis within a single high-cross-helicity hybrid simulation, this work separates two coexisting ion-scale wave populations with distinct physics. ICW are confirmed as intermittent, forward-propagating linear modes driven almost entirely by core temperature anisotropy built up at phase-steepened fronts, relaxing toward marginal stability after each growth episode. FMW, by contrast, exhibit no sustained linear growth, violate the validity condition for species-resolved linear decomposition, and persist against dominant damping—consistent with nonlinear generation via parametric decay or phase steepening. The results sharpen the interpretation of in situ ion-scale wave observations by demonstrating that the presence of a coherent polarized signal does not, by itself, imply a locally unstable linear mode.