Aperiodic magnetic turbulence produced by relativistic ion beams

Magnetic-field generation by a relativistic ion beam propagating through an electron-ion plasma along a homogeneous magnetic field is investigated with 2.5D high-resolution particle-in-cell (PIC) simulations. The studies test predictions of a strong amplification of short-wavelength modes of magnetic turbulence upstream of nonrelativistic and relativistic parallel shocks associated with supernova remnants, jets of active galactic nuclei, and gamma-ray bursts. We find good agreement in the properties of the turbulence observed in our simulations compared with the dispersion relation calculated for linear waves with arbitrary orientation of ${\vec k}$. Depending on the parameters, the backreaction on the ion beam leads to filamentation of the ambient plasma and the beam, which in turn influences the properties of the magnetic turbulence. For mildly- and ultra-relativistic beams, the instability saturates at field amplitudes a few times larger than the homogeneous magnetic field strength. This result matches our recent studies of nonrelativistically drifting, hot cosmic-ray particles upstream of supernova-remnant shocks which indicated only a moderate magnetic-field amplification by nonresonant instabilities. We also demonstrate that the aperiodic turbulence generated by the beam can provide efficient particle scattering with a rate compatible with Bohm diffusion. Representing the ion beam as a constant external current, i.e. excluding a backreaction of the magnetic turbulence on the beam, we observe non-resonant parallel modes with wavelength and growth rate as predicted by analytic calculations. In this unrealistic setup the magnetic field is amplified to amplitudes far exceeding the homogeneous field, as observed in recent MHD and PIC simulations.

💡 Research Summary

This paper investigates magnetic‑field generation by a relativistic ion beam propagating through an electron‑ion plasma that is permeated by a homogeneous background magnetic field. Using high‑resolution 2.5‑D particle‑in‑cell (PIC) simulations, the authors test the long‑standing prediction that short‑wavelength, non‑resonant modes can be strongly amplified upstream of parallel shocks in supernova remnants (SNRs), active‑galactic‑nucleus (AGN) jets, and gamma‑ray bursts (GRBs).

The study explores two distinct modeling approaches. In the physically realistic case, the ion beam is represented by actual particles, allowing the back‑reaction of the generated turbulence on the beam. In the second, more idealized scenario, the beam is treated as a constant external current, thereby suppressing any feedback.

Key findings from the realistic simulations are:

1. Agreement with Linear Theory – The dispersion relation derived for linear waves of arbitrary wave‑vector orientation matches the growth rates and wavelengths observed in the simulations. Both mildly relativistic (γ≈2) and ultra‑relativistic (γ≈10) beams excite aperiodic modes whose characteristics are quantitatively consistent with analytic predictions.

2. Filamentation and Saturation – The back‑reaction induces strong filamentation of both the ambient plasma and the beam. These filaments generate localized current sheets that limit the magnetic‑field amplification. The field grows to only a few times the initial homogeneous field B₀ before saturating, in stark contrast to the order‑of‑magnitude amplification sometimes reported in magnetohydrodynamic (MHD) studies that neglect beam feedback.

3. Unrealistic External‑Current Limit – When the beam is forced to remain a fixed current, the non‑resonant parallel modes develop exactly as predicted by Bell‑type theory, achieving magnetic amplitudes tens of times larger than B₀. This demonstrates that the extreme amplification reported in some earlier works is an artifact of the imposed current‑source approximation.

4. Particle Scattering Efficiency – The aperiodic turbulence produced in the realistic runs provides efficient pitch‑angle scattering. Test‑particle analysis shows that the scattering rate approaches the Bohm limit (D≈(1/3) r_g c), implying that high‑energy cosmic‑ray particles can be confined and accelerated effectively in the upstream region of relativistic shocks.

The authors conclude that while the non‑resonant instability can, in principle, generate very strong magnetic fields, the inclusion of beam back‑reaction dramatically reduces the achievable amplification. Nevertheless, the resulting turbulence is sufficient to yield Bohm‑type diffusion, a crucial ingredient for models of cosmic‑ray acceleration in SNRs, AGN jets, and GRBs. The work therefore reconciles earlier analytic expectations with fully kinetic simulations and highlights the importance of self‑consistent beam‑plasma coupling in astrophysical shock environments.