Induced Scattering of Strong Waves in Pair Plasmas

Abstract: We study induced (stimulated) scattering of linearly polarized, strong electromagnetic waves in pair plasmas, which is crucial for understanding the propagation of fast radio bursts (FRBs). Magnetars are the most promising progenitors of FRBs, and FRBs propagate through the magnetar wind and successfully escape before being significantly scattered. We revisit the steady-state solution of linearly polarized electromagnetic waves in pair plasmas with arbitrary amplitude, and demonstrate that the nonlinearity is characterized by the nonlinearity parameter $a_0ω<em>{pe}/ω_0$ rather than the dimensionless amplitude $a_0$, where $ω</em>{pe}$ is the electron plasma frequency and $ω<em>0$ is the wave frequency. We follow the time evolution of the steady-state solution for the linear regime $a_0ω</em>{pe}/ω<em>0 \ll 1$ by performing one-dimensional particle-in-cell simulations, and show that the conventional linear analysis of induced scattering assuming $a_0 \ll 1$ is applicable even for $a_0 &gt; 1$ when the Lorentz boost due to the plasma motion in the incident wave is considered. The saturation level is controlled by $a_0ω_0/ω</em>{pe}$, which corresponds to the ratio of the wave energy to the plasma energy, and the incident wave is hardly scattered for $a_0ω<em>0/ω</em>{pe} \gg 1$. We discuss the application of our results to FRBs.

• The paper demonstrates that induced scattering in pair plasmas is governed by the nonlinearity parameter a₀ωₚₑ/ω₀ rather than by the wave amplitude alone.
• A fully relativistic two-fluid model combined with PIC simulations reveals a transition from undistorted cosine profiles to sawtooth-like forms under strong coupling.
• The findings imply that fast radio bursts traverse magnetar winds with minimal scattering losses due to high incident wave-plasma energy ratios.

Induced Scattering of Strong Waves in Pair Plasmas: An Authoritative Essay

Motivation and Context

The investigation of induced (stimulated) scattering of strong electromagnetic waves in pair plasmas is critical for understanding wave-plasma interactions relevant to the propagation of fast radio bursts (FRBs) generated by magnetars. High-amplitude FRB pulses traverse magnetar winds with strength parameter $a_0 > 1$, and their escape depends on whether induced scattering (SBS/Compton) significantly impedes their propagation. Conventional treatments employ linear theory valid for $a_0 \ll 1$, but the pertinence of nonlinear effects when $a_0 > 1$ in tenuous pair plasmas remains contentious, inviting a rigorous kinetic and analytical re-examination.

Analytical Formulation and Steady-State Structure

The paper constructs a fully relativistic, cold two-fluid model and associated Maxwell equations for linearly polarized monochromatic plane waves in a pair plasma. All field and particle quantities are recast as functions of the wave phase, $\phi = \omega_0 t - k_0 x$, with the condition $\omega_0/k_0 > c$ ensuring superluminal propagation.

The analytical approach yields a set of self-consistent equations characterizing the plasma response, wave field profile, dispersion, and current. The plasma parameters—Lorentz factor, four-velocity, density—are derived as explicit functions of the normalized electric field $y = E_y/E_0$. Crucially, the nonlinear plasma feedback is shown to be governed not by $a_0$, but by the nonlinearity parameter $a_0\omega_{pe}/\omega_0$, encapsulating the ratio of the driven quiver energy to the plasma energy scale.

The steady-state solution reveals two asymptotic regimes:

• For $a_0\omega_{pe}/\omega_0 \ll 1$, plasma feedback is negligible, the solution converges to $y = \cos \phi$, and the waveform remains undistorted (test-particle limit).
• For $a_0 \ll 1$0, nonlinear effects dominate, yielding sawtooth-like profiles for $a_0 \ll 1$1 due to strong amplitude-dependent coupling with the plasma.

  Figure 1: The normalized wave electric field $a_0 \ll 1$2 across wave phase $a_0 \ll 1$3, for different $a_0 \ll 1$4.

The dispersion relation parameter $a_0 \ll 1$5, encoding frequency shifts and group velocity, is also shown to depend exclusively on $a_0 \ll 1$6, asymptoting to unity in the linear regime and decreasing for strong nonlinear coupling.

Figure 2: The plasma dispersion parameter $a_0 \ll 1$7 as a function of $a_0 \ll 1$8.

Induced Scattering: Parametric Instability Theory

Induced scattering is quantified as a parametric instability—stimulated Brillouin (SBS) or Compton scattering—whose linear growth rate and wavenumber are well-characterized analytically for $a_0 \ll 1$9. The analysis is extended by incorporating Lorentz boosts from the bulk plasma motion driven by intense incident waves, showing that the laboratory-frame growth rates and scattered wavenumbers depend both on $a_0 > 1$0 and the incident amplitude $a_0 > 1$1.

The transformation to the center-of-momentum frame allows the linear SBS theory to be extrapolated to $a_0 > 1$2 provided $a_0 > 1$3. In this regime, the maximal growth rates and scattered wavenumbers agree with Lorentz-boosted predictions, except in strong-coupling cases where the thermal speed $a_0 > 1$4 is below threshold.

Kinetic Simulations: Validation and Nonlinear Dynamics

One-dimensional particle-in-cell (PIC) simulations are employed, systematically exploring combinations $a_0 > 1$5 to keep $a_0 > 1$6 fixed (linear regime), varying both $a_0 > 1$7 and $a_0 > 1$8. The simulations confirm the theoretical predictions for growth rates and wavenumbers, including the $a_0 > 1$9-dependent Lorentz boost effect.

Figure 3: Initial spatial profiles of normalized laboratory density $\phi = \omega_0 t - k_0 x$0 for selected $\phi = \omega_0 t - k_0 x$1.

Figure 4: Temporal evolution of the Poynting flux spectrum, discriminating incident and backscattered components, for distinct parameter sets.

Figure 5: Growth dynamics of the Poynting flux at the fastest-growing mode, for variable $\phi = \omega_0 t - k_0 x$2 and $\phi = \omega_0 t - k_0 x$3, evidencing theoretical agreement.

Figure 6: Maximal SBS growth rate (top) and corresponding wavenumber (bottom) as functions of incident amplitude $\phi = \omega_0 t - k_0 x$4, for fixed nonlinearity parameter.

Nonlinear saturation of SBS is shown to depend sensitively on the dimensionless energy ratio $\phi = \omega_0 t - k_0 x$5. When this ratio is large ($\phi = \omega_0 t - k_0 x$6), the incident wave remains largely unscattered; only a minor fraction of the energy is transferred to the plasma over the simulation timescale. For small ratios ($\phi = \omega_0 t - k_0 x$7), SBS induces substantial dissipation and plateau formation in the velocity distribution, a precursor to the quenching of further resonant interaction.

Figure 7: Time-evolution of incident and scattered Poynting flux, showing saturation dynamics and the impact of $\phi = \omega_0 t - k_0 x$8.

Figure 8: Longitudinal four-velocity $\phi = \omega_0 t - k_0 x$9 distributions tracked over time, illustrating SBS-induced heating and distribution flattening.

Astrophysical Implications for FRB Propagation

Applying these results to FRB propagation in magnetar winds, the relevant parameters are $\omega_0/k_0 > c$0, $\omega_0/k_0 > c$1, and $\omega_0/k_0 > c$2 at $\omega_0/k_0 > c$3 cm for fiducial magnetar wind models. This places FRBs in the regime where linear SBS theory holds, Lorentz boosts are essential, instability growth rates are rapid, but the saturation level is minimal due to large wave-plasma energy ratios.

Practically, FRBs are predicted to traverse the magnetar wind with negligible energy loss via SBS, despite prompt particle heating and velocity distribution broadening. The complex interplay of scattered wave escape, filamentation instability, broadband pulse structure, and wind acceleration further diminishes nonlinear SBS effects in realistic astrophysical environments.

Theoretical and Future Directions

The determination that nonlinear plasma feedback can be neglected when $\omega_0/k_0 > c$4 advances the theoretical foundation for strong wave-plasma interaction models in both astrophysical and laboratory contexts. The identification of $\omega_0/k_0 > c$5 as the governing parameter, irrespective of amplitude $\omega_0/k_0 > c$6, is a robust result supported by both kinetic simulations and analytical theory. Saturation dynamics controlled by the energy ratio $\omega_0/k_0 > c$7 are likely universal for SBS-driven heating.

Extensions to magnetized plasmas, where the relevant parameter may shift to $\omega_0/k_0 > c$8, promise further refinement, and multi-dimensional, broadband, and open-boundary simulations are signposted as next steps in bridging theory with observational constraints and laboratory experiment design.

Conclusion

The paper rigorously shows that induced scattering of strong electromagnetic waves in unmagnetized pair plasmas is dictated by the nonlinearity parameter $\omega_0/k_0 > c$9, not amplitude $y = E_y/E_0$0, with linear theory valid even for $y = E_y/E_0$1 as long as $y = E_y/E_0$2. Kinetic simulations verify Lorentz boost effects and support the extrapolation of linear SBS arguments to strong-wave regimes. The incident wave energy ratio $y = E_y/E_0$3 governs nonlinear saturation, and large values inhibit dissipation, ensuring FRB transparency in magnetar winds for physically motivated parameter sets. Broader theoretical implications include extensions to magnetized environments and guidance for experimental plasma wave studies.