Fast Magnetosonic Turbulence in Two-Dimensional Relativistic Plasmas

Abstract: We present fully kinetic simulations of driven 2D turbulence in a relativistic plasma, designed for the first time to induce a fast magnetosonic cascade. As the driving strength increases, turbulence transitions from a weak wave-dominated regime to strong shock-driven dynamics. Using spatiotemporal Fourier analysis, we identify fast modes, finding that the weak turbulence regime exhibits spectral properties that are in excellent agreement with theoretical expectations. Our results are relevant for the modeling of turbulence in high-energy astrophysical plasmas.

• The paper demonstrates that fast magnetosonic turbulence cascades persist in 2D relativistic pair plasmas using kinetic PIC simulations.
• It identifies a clear transition from wave-dominated to shock-dominated regimes, with spectral indices shifting from -1.35 to -2.2.
• Findings imply that fluid-like behavior persists at kinetic scales, challenging standard MHD models in high-energy astrophysical contexts.

Fast Magnetosonic Turbulence in 2D Relativistic Plasmas: Kinetic Regimes and Scaling Properties

Introduction

The paper "Fast Magnetosonic Turbulence in Two-Dimensional Relativistic Plasmas" (2604.04276) investigates the kinetic dynamics of fast magnetosonic (FM) turbulence within a collisionless, ultra-relativistic pair plasma in two dimensions. The study addresses a pivotal gap in plasma astrophysics, where compressible turbulence—crucial for energy transfer and dissipation in high-energy environments—remains poorly characterized beyond the magnetohydrodynamic (MHD) framework, especially under relativistic and collisionless conditions. Conventional MHD models do not account for kinetic and relativistic effects that fundamentally alter mode dynamics, damping, and cascade mechanisms, particularly for FM modes, which dominate compressible energy transfer in many astrophysical plasmas.

Simulation Methodology

The authors employ fully kinetic, particle-in-cell (PIC) simulations with the Zeltron code to drive turbulence in a 2D electron-positron plasma with temperatures $\theta \equiv T / m_e c^2 \gg 1$. The configuration ensures an initially homogeneous magnetic field perpendicular to the simulation plane, so only FM modes—with wavevector perpendicular to the field—can be excited. This geometry minimizes collisionless damping for $k_\perp \rho_e \lesssim 1$, enabling robust isolation of FM cascades. External compressive driving is implemented via a velocity fluctuation field constructed from large-scale, randomly evolving Fourier modes.

A systematic scan over a dimensionless driving parameter $F$ provides access to regimes ranging from weakly driven, wave-dominated turbulence (subsonic) to strongly driven, shock-dominated turbulence (supersonic). The simulations resolve inertial and kinetic scales with high spatial and particle resolution, and spectral analysis is performed using spatiotemporal Fourier transforms and Helmholtz decomposition of velocity and electric fields, augmented by explicit PIC noise subtraction.

Results: Wave-Shock Transition and Mode Identification

A central result is the demonstration of a clear transition between wave-dominated and shock-dominated turbulence as driving strength increases. The "supersonic fraction" $S$ serves as a quantitative metric, reflecting the fraction of local FM Mach numbers exceeding unity. At low driving ($F \ll 1$), $S \ll 1$ and the density morphology is coherent, whereas for strong driving ($F \gtrsim 1$), $S$ saturates and density exhibits irregular, shock-like structures.

Spatiotemporal spectral analysis reveals that, in the weak regime, spectral power is concentrated along the analytical FM dispersion relation, persisting from fluid to kinetic scales ($k_\perp \rho_e \ll 1$ to $k_\perp \rho_e \gtrsim 1$) with minimal nonlinear correction. Importantly, even at kinetic scales, FM power does not transition to Langmuir or electromagnetic branch as predicted by classic theory, indicating effective persistence of fluid-like behavior due to possible enhanced collisionality or suppressed damping mechanisms.

In the strong driving (supersonic) regime, spectral power distribution is broader with substantial low-frequency, non-wave contributions and a spectral index steepening toward Burgers-like values, characteristic of shock-dominated turbulence.

Scaling of Turbulence and Spectral Indices

A bold, contradictory result emerges regarding the spectral scaling: In the weakly driven (wave-dominated) regime, the spatial spectra of magnetic field, density, electric field, and velocity are nearly identical, with inertial-range power-law indices consistently near $k_\perp \rho_e \lesssim 1$0. For strong driving, these indices steepen to approximately $k_\perp \rho_e \lesssim 1$1. This result both affirms classical predictions for FM turbulence and challenges the expectation of kinetic-scale corrections and strong damping effects. The persistence of linear MHD FM dispersion and similar scaling across fields suggests nontrivial physics in energy transfer and dissipation, likely driven by the absence of common 2D dissipation channels (reconnection, parallel Landau damping, anisotropic stress) and the dominance of cyclotron resonance.

The authors highlight that theoretical predictions for 2D FM turbulence remain challenging, as standard weak turbulence theory breaks down due to secular term growth and singularities in the wave kinetic equation, unless weak dispersion is present. The observed scaling thus motivates further analytical studies.

Practical and Theoretical Implications

The findings establish that FM wave turbulence can persist in fully kinetic, relativistic plasmas without inevitable steepening into shocks—a result not previously demonstrated. This supports the robustness of the wave-cascade paradigm even under extreme plasma conditions and may necessitate updates to turbulence models in high-energy astrophysical contexts, including cosmic-ray scattering, pulsar wind nebulae, AGN jets, and potentially laboratory analogues such as electron-hole gases in 2D materials (e.g., graphene).

The suppression of shock formation and maintenance of a wave-dominated energy cascade implies that dissipation proceeds primarily via cyclotron damping, with stochastic heating and shocks only relevant at higher driving amplitudes. This distinction alters entropy production mechanisms and may impact particle acceleration and heating rates in astrophysical systems.

Conclusions

The paper provides the first kinetic demonstration of FM-dominated turbulence in a relativistic, collisionless pair plasma, elucidating the wave-shock transition, spectral scaling, and mode persistence across fluid and kinetic regimes. The major numerical result is the inertial-range spectral index of $k_\perp \rho_e \lesssim 1$2 for magnetic, velocity, electric, and density fluctuations in the weak regime, steepening to $k_\perp \rho_e \lesssim 1$3 in the shock-dominated case. Theoretical implications extend to the nature of wave turbulence and dissipation mechanisms, while practical consequences affect energy transport models for extreme astrophysical plasmas.

Future directions include detailed investigation of dissipation channels, entropy production, and extension to three dimensions and laboratory analogs. The persistence of FM wave turbulence and its scaling properties under relativistic kinetic conditions provides a foundation for advancing the theory of compressible plasma turbulence in astrophysical and laboratory settings.