Formation of Suprathermal Electron Populations in the Expanding, Turbulent Solar Wind

This study utilizes the first fully kinetic particle-in-cell simulation of an expanding, turbulent solar wind to demonstrate that the combined effects of expansion-driven cooling and Alfvénic turbulence generate suprathermal electron populations with parallel power-law tails, suggesting their origin lies in parallel electric fields or resonant wave-particle interactions rather than simple velocity-space redistribution.

The Gist

Imagine the solar wind not as a gentle breeze, but as a chaotic, expanding river of invisible particles rushing away from the Sun. In this river, the electrons (tiny, fast particles) usually behave like a calm crowd, but often they suddenly develop "suprathermal" tails—groups of electrons that get kicked up to incredibly high speeds, forming a power-law distribution. Scientists have long wondered: How do these high-speed electrons get their energy in a space that is too empty for particles to bump into each other like billiard balls?

This paper acts as a high-speed, 3D movie simulation to answer that question. Here is what the researchers found, explained simply:

The Setup: A Stretching, Turbulent Box

The scientists built a virtual "box" representing a chunk of the solar wind. They set it up with two main ingredients:

1. Expansion: Like a balloon being blown up, the box stretches sideways (perpendicular to the magnetic field) but stays the same length forward and backward.
2. Turbulence: They stirred the pot with magnetic waves (Alfvénic turbulence), creating a chaotic, swirling environment similar to what exists in space.

They used a supercomputer to watch how the electrons and ions (heavier particles) reacted to this stretching and swirling.

The Stretching Effect: Cooling the Sideways Motion

As the box stretched sideways, something interesting happened to the electrons. Imagine a figure skater spinning; if they stretch their arms out, they slow down. Similarly, as the magnetic field stretched, the electrons' motion perpendicular to the field (sideways) cooled down and slowed. However, their motion parallel to the field (forward and backward) stayed roughly the same.

This created a lopsided situation: the electrons were "cold" sideways but "hot" forward. In physics terms, this pushed the plasma toward a tipping point called the firehose instability. Think of it like a garden hose that is pressurized too much; if the water pressure gets too high compared to the hose's strength, the hose starts to whip around uncontrollably. Here, the "whipping" is a magnetic instability that tries to fix the lopsidedness.

The Surprise: High-Speed Tails Form Forward

The researchers expected the instability to just shuffle particles around, making the distribution more even. Instead, they saw something more dramatic:

• The Perpendicular Side: The electrons got a bit "heated up" sideways due to the turbulence, forming a small group of fast movers.
• The Parallel Side (The Big Discovery): Even though the turbulence was mostly pushing things sideways, a massive group of electrons suddenly accelerated forward (parallel to the magnetic field). They formed a distinct "tail" of super-fast particles, following a mathematical pattern known as a power law.

Crucially, these high-speed tails formed before the firehose instability fully kicked in to regulate the system. This suggests that the instability isn't the cause of the high speeds, but rather a reaction to them.

The Mechanism: Direct Acceleration, Not Just Shuffling

The paper argues that these electrons didn't just get pushed from the side to the front (like shuffling cards). Instead, they were likely directly accelerated in the forward direction.

The Analogy:
Imagine a crowded dance floor (the plasma).

• Old Theory: The instability acts like a bouncer who grabs people dancing wildly in one spot and shoves them to a different spot to make the room even.
• This Paper's Finding: It's more like a DJ playing a specific beat that makes a specific group of people suddenly sprint forward in a straight line, creating a "tail" of runners, while the rest of the crowd stays put. The "DJ" here is likely the interaction between the particles and specific electric fields or waves moving along the magnetic field lines.

The Conclusion

The study provides the first direct evidence that in the expanding, turbulent solar wind:

1. Expansion cools the sideways motion, creating a lopsided state.
2. Turbulence and specific wave interactions directly accelerate electrons forward, creating high-energy "tails."
3. The firehose instability eventually steps in to stop the system from getting too lopsided, but it preserves the high-speed tails that were already formed.

In short, the solar wind doesn't just shuffle its electrons around; it actively cooks up high-speed populations in the direction of the magnetic field, a process driven by the unique combination of cosmic expansion and magnetic turbulence.

Technical Summary

Technical Summary: Formation of Suprathermal Electron Populations in the Expanding, Turbulent Solar Wind

Problem Statement
Nonthermal features, specifically suprathermal electron populations (strahl and halo) characterized by power-law tails, are ubiquitously observed in solar wind velocity distribution functions (VDFs). However, the mechanisms generating and regulating these features in the collisionless, turbulent, and expanding solar-wind plasma remain unclear. While previous studies have examined ion firehose instabilities in expanding-box simulations or electron-scale instabilities in homogeneous conditions, the coupled interaction between turbulence, expansion, and firehose instabilities for both ions and electrons has remained essentially unexplored. A specific open question concerns how suprathermal tails develop in the parallel direction (where the strahl is observed) given that highly Alfvénic turbulence typically drives perpendicular heating.

Methodology
The authors present the first fully kinetic, three-dimensional particle-in-cell (PIC) simulation of an expanding, turbulent plasma under heliospheric conditions, utilizing the relativistic PIC code Zeltron. The simulation employs a "kinetic expanding box" framework where the computational domain expands perpendicular to a guide magnetic field ($B_0$) while the parallel length remains fixed. This setup mimics the solar wind's radial expansion, causing the magnetic field strength and density to decrease as $a_{\perp}^{-2}(t)$.

Key numerical parameters and setup details include:

• Turbulence Driver: A solenoidal, charge-independent forcing is applied to both ions and electrons to maintain Alfvénic turbulence throughout the run. The forcing mimics an anisotropic turbulent cascade.
• Physical Regime: The simulation runs in a weakly compressive, highly Alfvénic regime. The initial ion plasma beta is $\beta_{i0} = 1$, and the initial ion Alfvén speed is $v_{Ai0} = 0.02c$.
• Numerical Constraints: To ensure computational feasibility, the ion-to-electron mass ratio is reduced to $m_i/m_e = 25$, and the light-to-Alfvén speed ratio is $c/v_{Ai} \gtrsim 50$. The authors note that while these reduce scale separation, the dimensionless parameters governing expansion-driven instabilities (parallel plasma-$\beta$ and temperature anisotropy) are initialized with solar-wind-like values.
• Evolution: The system evolves to a quasi-steady turbulent state before expansion begins at $t_{exp} = 0.5\tau_{exp}$. The expansion continues for a full timescale $\tau_{exp} = 10\tau_{A0}$.

Key Results

1. Global Evolution and Instability Onset:

   • Expansion drives the plasma toward the firehose-instability threshold by adiabatically cooling the perpendicular temperature while leaving the parallel temperature largely unchanged.
   • The electron firehose instability (EFI) is triggered at $t \approx 0.75\tau_{exp}$, marked by a rapid increase in the $T_{\perp e}/T_{\parallel e}$ ratio. The ion firehose instability develops more gradually, reaching its threshold around $t \approx 1.2\tau_{exp}$.
   • Magnetic field fluctuations increase in strength during expansion, consistent with Alfvénic turbulence scaling, though the background field decrease is slightly slowed by these fluctuations.

2. Formation of Suprathermal Electron Populations:

   • Perpendicular Direction: Turbulent heating generates suprathermal tails in the perpendicular direction prior to the onset of expansion. These tails remain relatively stable throughout the expansion phase.
   • Parallel Direction: A pronounced suprathermal tail develops in the parallel direction, well-fitted by a Kappa distribution ($\kappa = 5$) by the end of the simulation. Crucially, this parallel energization begins abruptly around $t \approx 0.67\tau_{exp}$, before the global onset of the EFI ($t_{EFI} = 0.75\tau_{exp}$).
   • Anisotropy Regulation: As the EFI triggers, it acts as a diffusive mechanism, redistributing particles from the parallel to the perpendicular direction, thereby regulating the temperature anisotropy and returning the plasma to a state of marginal stability. However, the suprathermal tails persist despite this regulation.

Significance and Claims
The paper claims to provide the first direct evidence of nonthermal electron features emerging within a unified kinetic framework that links expansion, turbulence, and instabilities. The primary findings and their implications are:

• Mechanism of Parallel Energization: The preferential development of suprathermal tails in the parallel direction, occurring prior to the global EFI onset, suggests that the energization mechanism is not merely a consequence of global instability or simple velocity-space redistribution (pitch-angle scattering). Instead, the authors propose the involvement of parallel electric fields or resonant wave–particle interactions that directly accelerate particles in velocity space.
• Decoupling of Formation and Regulation: The results indicate that the formation of suprathermal tails is largely decoupled from the state of marginal stability. The EFI serves to regulate the system's anisotropy after the tails have formed, rather than being the primary driver of their creation.
• Limitations and Context: The authors acknowledge that the perpendicular suprathermal populations present at the start of expansion (due to the pre-expansion turbulent heating) may influence the parallel tail formation. They posit that while the specific initial conditions are a limitation, the mechanism of direct parallel acceleration likely operates in more realistic, simultaneous evolution scenarios found in the solar wind.

In conclusion, the study demonstrates that in an expanding, turbulent solar wind, the interplay of these factors can generate and sustain parallel suprathermal electron populations, offering a potential explanation for the origin of the strahl and halo components observed in situ.