PIC simulations of nonrelativistic high-Mach-number oblique shocks propagating in a turbulent medium

This paper presents the first 2D3V particle-in-cell simulations demonstrating that pre-existing compressive turbulence in nonrelativistic oblique shocks enhances whistler-wave instabilities, resulting in a shorter, hotter electron foreshock and more efficient non-thermal electron acceleration.

The Gist

Imagine the universe is filled with invisible, high-speed "wind" made of charged particles (plasma). Sometimes, this wind hits a wall of magnetic fields and slams into a shockwave, much like a car crashing into a brick wall. In space, these crashes are called collisionless shocks. They are famous for being cosmic particle accelerators, blasting tiny electrons to near-light speeds.

For a long time, scientists thought these shocks happened in a perfectly smooth, empty vacuum. But in reality, the space ahead of these shocks is often turbulent—think of it as a calm river suddenly turning into a choppy, foamy rapid with swirling eddies and bumps.

This paper asks a simple question: What happens to the particle acceleration if the "wind" hitting the shock is already choppy and turbulent, rather than smooth?

Here is the story of what the researchers found, using some everyday analogies:

1. The Setup: The Smooth Road vs. The Bumpy Road

The scientists used a supercomputer to run a virtual experiment (a "Particle-in-Cell" simulation). They created two scenarios:

• Scenario A (The Smooth Road): A shockwave moves through a perfectly smooth, calm stream of particles.
• Scenario B (The Bumpy Road): A shockwave moves through a stream that is already 15% turbulent, full of density bumps and magnetic swirls (mimicking the real interstellar medium).

They focused on oblique shocks, which are like hitting a wall at an angle rather than head-on. This angle allows some particles to bounce back upstream, creating a "foreshock" region—a waiting area before the main crash.

2. The "Whistler" Waves: The Bouncy Ball Effect

In the smooth scenario, the shock creates a specific type of wave called a whistler wave. Imagine these waves as bouncy balls that knock the incoming electrons, giving them a little push to get them ready for the big acceleration.

• What happened in the turbulent scenario?
  The pre-existing turbulence acted like a giant mixer. It made these "bouncy balls" (whistler waves) much stronger and created larger, more chaotic structures.
  • The Result: The "bouncy balls" appeared earlier and grew larger (about 3.5 times bigger in size) in the turbulent simulation. It's like having a trampoline that is already being shaken by a storm; when you jump on it, the bounce is wilder and more unpredictable.

3. The "Foreshock" Shrinkage: A Shorter Waiting Room

Usually, the "foreshock" is a long region where reflected electrons bounce back and forth, getting heated up and scattered before they hit the main shock.

• The Finding: When the upstream medium was turbulent, this waiting room shrank. The electrons didn't travel as far upstream before being turned around.
• The Analogy: Imagine a hallway where people are bouncing off walls. If the walls are smooth, people bounce far down the hall. If the hallway is filled with obstacles (turbulence), people get bounced back much sooner. The result? The electrons in the turbulent scenario were hotter (more energetic) right from the start because they were being scattered more aggressively by the pre-existing chaos.

4. The Final Crash: More Energy, More Particles

The ultimate goal of these shocks is to accelerate particles to high energies.

• The Smooth Scenario: A small fraction of electrons got supercharged.
• The Turbulent Scenario: The results were significantly better.
  • More Particles: There were about 60% more high-energy electrons.
  • More Energy: These electrons carried nearly double the total energy compared to the smooth scenario.
  • Higher Speeds: The fastest electrons reached energies 40% higher than in the smooth case.

5. The "Cavities": Giant Bubbles of Heat

The turbulence helped create massive, bubble-like structures in the magnetic field (called nonlinear cavities).

• What are they? Think of them as giant, hollow bubbles made of magnetic force. Inside these bubbles, hot, fast electrons get trapped.
• The Effect: Because the turbulence made these bubbles bigger and stronger, they distorted the shockwave more violently when they finally merged with it. This created a more chaotic and powerful environment for acceleration.

The Bottom Line

The paper concludes that pre-existing turbulence is a game-changer. It doesn't just add a little noise; it fundamentally rewrites the rules of the crash. By making the "waiting room" (foreshock) shorter and hotter, and by creating larger, more powerful magnetic bubbles, turbulence makes the shockwave a much more efficient particle accelerator.

In simple terms: If you want to blast particles to high speeds in space, you don't want a smooth, calm approach. You want a bumpy, turbulent one. The chaos ahead of the crash actually helps the crash happen better.

Technical Summary

Technical Summary: PIC Simulations of Nonrelativistic High-Mach-Number Oblique Shocks in a Turbulent Medium

Problem Statement
Collisionless shocks are ubiquitous in astrophysical systems, such as supernova remnants (SNRs), and serve as primary sites for particle acceleration. While Diffusive Shock Acceleration (DSA) provides a foundational framework for understanding particle energization, it typically assumes a homogeneous upstream medium, neglecting the pre-existing turbulence common in astrophysical environments. A critical challenge in shock physics is the "electron injection problem": thermal electrons have dynamical scales much smaller than the shock width and cannot directly participate in DSA without a pre-acceleration mechanism. Previous studies have established that at oblique shocks, shock-reflected electrons drive instabilities (specifically whistler waves) in an "electron foreshock," facilitating stochastic shock drift acceleration (SSDA). However, the interplay between these kinetic-scale processes and pre-existing compressive turbulence remains poorly understood. This work addresses the gap in understanding how pre-existing turbulence influences the structure of oblique shocks, the development of foreshock instabilities, and the efficiency of electron acceleration.

Methodology
The authors employed fully relativistic 2D3V (two spatial, three velocity components) Particle-in-Cell (PIC) simulations using the code THATMPI. The study modeled nonrelativistic, high-Mach-number electron-ion shocks ($M_s \approx 36$, $M_A \approx 32$) with an obliquity angle of $\theta_0 = 60^\circ$.

Two distinct simulation runs were conducted to isolate the effects of turbulence:

1. Run H (Homogeneous): A standard simulation with a uniform upstream medium.
2. Run T (Turbulent): A simulation where the upstream medium contains pre-existing compressive density fluctuations with a relative amplitude of $\delta n/n_0 = 15\%$, consistent with interstellar medium estimates.

To incorporate turbulence, the authors utilized a method of continuous turbulent plasma injection. Prefabricated slabs of turbulent plasma were generated in separate periodic boxes, allowed to evolve until correlated electromagnetic and density fluctuations were established, and then injected into the far-upstream region of the shock simulation. The simulation box was extended in the transverse direction ($L_y \approx 12\lambda_{si}$) compared to previous perpendicular shock studies to fully capture the electron foreshock and large nonlinear structures driven by whistler instabilities. The simulations evolved for approximately 20 ion cyclotron times ($\Omega_i^{-1}$).

Key Contributions and Results
The study presents the first 2D3V PIC simulations of nonrelativistic oblique shocks interacting with pre-existing compressive turbulence. The key findings are:

• Shock Structure and Magnetic Amplification: Pre-existing turbulence modifies the global shock structure, resulting in a wider overshoot and a less pronounced undershoot compared to the homogeneous case. The Weibel instability is suppressed due to the out-of-plane magnetic field configuration; thus, magnetic field amplification is driven primarily by whistler waves ahead of the shock. In the turbulent run, magnetic field amplitudes in the foreshock region were approximately 25% higher than in the homogeneous run.
• Foreshock Dynamics and Whistler Instability: The presence of turbulence significantly alters the electron foreshock. The foreshock region in the turbulent simulation is shorter (reduced by $\sim$5 ion Larmor radii) and "hotter" (shock-reflected electrons possess higher temperatures).
  • Whistler Growth: The growth rate of the right-hand circularly polarized whistler waves was found to be lower in the turbulent case ($\gamma_T \approx 5.8 \times 10^{-3}\Omega_e$) compared to the homogeneous case ($\gamma_H \approx 8.2 \times 10^{-3}\Omega_e$). This reduction is attributed to the higher temperature of the reflected electrons and potentially reduced temperature anisotropy.
  • Nonlinear Structures: Nonlinear cavities driven by the whistler instability emerged earlier (by $\sim$2 ion cyclotron times) and reached larger maximum sizes ($7\lambda_{si}$ vs. $3\lambda_{si}$) in the turbulent simulation. These structures are characterized by high electron pressure and bipolar magnetic field features, inflating as they advect toward the shock.
• Electron Acceleration: The most significant outcome concerns particle acceleration. At the end of the simulations, the turbulent run produced a non-thermal electron population that was:
  • More numerous (0.23% of total electrons vs. 0.14% in the homogeneous run).
  • Carrying a larger fraction of the total energy (3.6% vs. 1.9%).
  • Reaching higher maximum energies (40% higher maximum Lorentz factor).

Significance and Claims
The paper claims that pre-existing compressive turbulence plays a beneficial role in electron acceleration at high-Mach-number oblique shocks. By enhancing magnetic field fluctuations and creating larger, earlier-emerging nonlinear structures, turbulence facilitates more efficient scattering and energization of electrons. The authors posit that these magnetic fluctuations, along with waves spawned in the shock ramp, are crucial for stochastic shock-drift acceleration.

The study acknowledges limitations inherent to 2D simulations, noting that 3D effects might alter the interplay between whistler and Weibel instabilities and the properties of the turbulence itself. However, the authors argue that their 2D approach captures the essential physics of oblique whistler waves and provides a reasonable description of the acceleration mechanisms. The results suggest that realistic models of cosmic-ray sources, such as SNRs, must account for the influence of upstream turbulence to accurately predict electron injection and acceleration efficiency. The findings also draw parallels to microscale coherent structures observed in the heliosphere, suggesting that such solitary structures may be a universal feature of collisionless shocks across various plasma environments.