Solar wind density turbulence and solar flare electron transport from the Sun to the Earth

Solar flare accelerated electron beams propagating away from the Sun can interact with the turbulent interplanetary media, producing plasma waves and type III radio emission. These electron beams are detected near the Earth with a double power-law energy spectrum. We simulate electron beam propagation from the Sun to the Earth in the weak turbulent regime taking into account the self-consistent generation of plasma waves and subsequent wave interaction with density fluctuations from low frequency MHD turbulence. The rate at which plasma waves are induced by an unstable electron beam is reduced by background density fluctuations, most acutely when fluctuations have large amplitudes or small wavelengths. This suppression of plasma waves alters the wave distribution which changes the electron beam transport. Assuming a 5/3 Kolmogorov-type power density spectrum of fluctuations often observed near the Earth, we investigate the corresponding energy spectrum of the electron beam after it has propagated 1 AU. We find a direct correlation between the spectrum of the double power-law below the break energy and the turbulent intensity of the background plasma. For an initial spectral index of 3.5, we find a range of spectra below the break energy between 1.6-2.1, with higher levels of turbulence corresponding to higher spectral indices.

💡 Research Summary

The paper investigates how density turbulence in the solar wind modifies the transport of electron beams accelerated during solar flares as they travel from the Sun to the Earth. Observations show that near‑Earth electron spectra often exhibit a double power‑law shape, with a break at a few tens of keV, a feature that standard beam‑propagation models have struggled to reproduce. The authors address this discrepancy by coupling a self‑consistent treatment of Langmuir‑wave generation with the effects of background density fluctuations that arise from low‑frequency magnetohydrodynamic (MHD) turbulence.

A one‑dimensional kinetic model is employed. The electron distribution function f(v,x) evolves according to the quasilinear diffusion equation, while the spectral energy density of Langmuir waves W(k,x) follows a wave‑propagation equation that includes growth driven by the beam and damping due to Landau resonance. The background electron density n₀ is modulated by a stochastic component δn(x) whose power spectrum follows a Kolmogorov law, P(k) ∝ k⁻⁵ᐟ³, consistent with in‑situ measurements near 1 AU. The turbulence intensity is parameterized by the root‑mean‑square (RMS) value of δn/n₀, and the characteristic wavelength range is varied to explore its impact on wave growth.

Simulations start with an electron beam having an initial power‑law index of 3.5 and a limited pitch‑angle spread. In a homogeneous plasma (δn = 0) the beam rapidly excites strong Langmuir turbulence, leading to efficient energy loss and a steepening of the spectrum well before the beam reaches 1 AU. Introducing density fluctuations changes the picture dramatically. When the fluctuation amplitude is large (δn/n₀ ≈ 10⁻³) or the dominant wavelengths are short, the resonance condition ω ≈ k v is intermittently broken, suppressing the growth rate of Langmuir waves. Consequently, the wave energy density remains low, the beam suffers less scattering and deceleration, and a larger fraction of the original high‑energy electrons survive to 1 AU.

The key result is a quantitative link between the turbulence level and the low‑energy spectral index after propagation. For the chosen initial index of 3.5, the simulated spectra below the break energy fall in the range 1.6 – 2.1, with higher turbulence producing higher indices (i.e., flatter spectra). This reproduces the variety of double‑power‑law spectra reported by spacecraft, suggesting that the observed spread is largely a consequence of differing solar‑wind turbulence conditions along individual field lines.

In addition, the Langmuir waves generated in the simulations have frequencies below ~3 kHz, matching the characteristics of type III radio bursts observed in the interplanetary medium. This provides a self‑consistent explanation for both the particle and radio signatures of flare‑associated electron beams.

Overall, the study demonstrates that weak, Kolmogorov‑type density turbulence can substantially modulate beam‑driven plasma wave activity, thereby shaping the electron energy distribution that arrives at Earth. By incorporating realistic turbulence spectra, the model bridges the gap between theory and observation, offering a robust framework for interpreting solar energetic electron events and improving space‑weather forecasting.