Cyclic reformation of subcritical perpendicular fast magnetosonic shocks due to oblique Whistler waves

Two-dimensional PIC simulations reveal that subcritical perpendicular fast magnetosonic shocks undergo cyclic reformation driven by an oblique lower-hybrid gradient drift instability, where the coupling between backward-propagating oblique Whistler waves and forward-propagating ion acoustic waves modulates the magnetic field and forces the shock to collapse and reform without external perturbations.

The Gist

The Big Picture: A Cosmic Traffic Jam

Imagine a massive, invisible highway in space where a super-fast river of gas (the solar wind) is rushing toward a wall. In our everyday world, if a car hits a wall, it crumples, and the energy turns into heat and noise. In space, there is no air to crash into, so the "crash" happens differently.

This crash is called a shock wave. It's a boundary where the fast-moving space gas suddenly slows down and gets hot. Scientists have known for a long time that these shock waves aren't perfectly smooth; they wiggle, ripple, and sometimes even collapse and rebuild themselves.

This paper is about figuring out why a specific type of shock wave (one that is moving fast but not too fast) keeps rebuilding itself, and what invisible force is pushing the buttons.

The Characters in the Story

1. The Shock Wave: Think of this as a moving wall of compressed air and magnetic fields. It's the "front line" of the crash.
2. The Electrons: These are tiny, light, and fast particles. They are like a swarm of hyperactive bees.
3. The Ions: These are heavier, slower particles (like nitrogen atoms in this study). They are like the bees' lazy, heavy cousins.
4. The Magnetic Field: Imagine this as invisible rubber bands stretching through space. They hold the plasma together and guide the particles.

The Mystery: Why Does the Wall Rebuild Itself?

In previous studies, scientists saw these shock walls collapse and reform. They thought it was because the wall was "bumpy" or because an outside force pushed it.

In this new study, the researchers set up a computer simulation (a virtual universe) to watch what happens when the magnetic field is tilted at a specific angle (45 degrees). They discovered something surprising: The shock wave isn't being pushed from the outside. It is being pushed from the inside by a specific type of invisible wave.

The Culprit: The "Whistler" Wave

The paper identifies a specific type of wave called an oblique Whistler wave.

• The Analogy: Imagine a long, stretched-out rubber band (the magnetic field). If you flick one end, a wave travels down it. Now, imagine that wave is also spinning and twisting as it moves. That's a Whistler wave.
• The "Oblique" part: Usually, these waves travel straight across the magnetic field. In this study, they are traveling at a diagonal angle (oblique), like a skier cutting across a slope rather than going straight down.

The Mechanism: The "Reactive" Dance

Here is the step-by-step process of how the shock wave gets rebuilt, explained simply:

1. The Setup: The shock wave creates a region where the magnetic field is squeezed tight (the "overshoot"). This squeezes the electrons, making them drift sideways very fast, like water rushing through a narrow pipe.
2. The Instability: This fast-drifting electron current becomes unstable. It's like a crowd of people running in a hallway; if they run too fast and in the wrong direction, they start bumping into each other and creating chaos.
3. The Wave Growth: This chaos creates the oblique Whistler wave. Think of it as a ripple forming on the surface of a pond, but this ripple is made of magnetic energy.
4. The Collision: This magnetic ripple moves backward (relative to the electrons) and crashes into the forward-moving ions.
5. The Reformation: The magnetic ripple is so strong that it acts like a piston (a mechanical pusher). It smashes the existing shock wave, collapsing it. But immediately after, the pressure builds up again, and a new shock wave forms right next to it.

The Key Difference: In the past, scientists thought the shock wave collapsed because it hit a bump in the road (an external perturbation). This paper shows that the shock wave collapses because of an internal engine (the Whistler wave) that it creates itself. It's like a car that doesn't just hit a pothole; the car's own engine suddenly revs up, flips the car over, and then the car rights itself and keeps driving.

Why Does This Matter?

• Space Weather: Earth's magnetic field acts as a shield against the solar wind. The "bow shock" is the front line of that shield. Understanding how these shocks wiggle and rebuild helps us predict how space weather might affect our satellites and power grids.
• The "Reactive" Nature: The paper proves that this process is "reactive," meaning the whole crowd of particles is involved, not just a few lucky ones. This is different from other types of instabilities where only a few fast particles cause the trouble.
• The Growth Rate: The waves in the simulation grew much faster than simple math predicted. The authors suspect this is because the electrons in the shock are "hot" in one direction and "cold" in another (like a stretched balloon), which makes the instability explode even faster.

The Takeaway

This research is like finding the hidden gear inside a clock that makes the hands jump forward. The scientists discovered that tilted magnetic fields allow a specific type of wave (the oblique Whistler) to form. This wave acts as a self-correcting mechanism, constantly smashing the shock wave down and rebuilding it, keeping the cosmic traffic jam in a state of constant, rhythmic motion.

In short: The shock wave isn't just a static wall; it's a living, breathing entity that uses its own internal magnetic waves to constantly tear itself down and build itself back up.

Technical Summary

1. Problem Statement

The paper investigates the stability and evolution mechanisms of subcritical perpendicular fast magnetosonic (FMS) shocks in collisionless plasmas.

• Context: While supercritical shocks are known to undergo cyclic reformation driven by reflected ion beams, subcritical shocks (Mach number $M < 2.7$) reflect fewer ions and were previously thought to be more stable or to reform only due to external perturbations or magnetic tension in specific geometries.
• Gap in Knowledge: Previous 2D Particle-in-Cell (PIC) simulations showed that the orientation of the background magnetic field ($B_0$) relative to the simulation plane significantly alters shock dynamics.
  • If $B_0$ is resolved in the plane, magnetic tension drives oscillations.
  • If $B_0$ is perpendicular to the plane, lower-hybrid drift instabilities arise.
• Objective: This study explores a configuration where $B_0$ is oriented at $45^\circ$ relative to the simulation plane normal (but remains perpendicular to the shock normal). The authors aim to identify the specific instability driving wave growth and shock reformation in this intermediate geometry, where traditional lower-hybrid drift modes are damped.

2. Methodology

The authors utilized two-dimensional Particle-in-Cell (PIC) simulations using the EPOCH code.

• Plasma Composition: Fully ionized Nitrogen ($Z=7$) and electrons.
• Initial Conditions:
  • A dense plasma slab expands into a uniform ambient plasma, driving a shock.
  • Magnetic Field: Uniform $B_0$ oriented at $45^\circ$ to the simulation plane normal ($z$-axis) and perpendicular to the shock propagation direction ($x$-axis).
  • Parameters: $M \approx 1.7$ (subcritical), $\beta \approx 0.56$, and a specific electron drift velocity ($v_D$) driven by the diamagnetic current.
• Simulation Suite:
  • Main Simulation (2D): Tracks the shock evolution in a box of $180\lambda_e \times 36\lambda_e$ to observe non-linear saturation and reformation.
  • Control Simulations (1D & 2D): Used to isolate wave modes, verify dispersion relations, and test the effect of box size (wavenumber resolution) on instability growth.
  • Perturbed Shock Simulation: Introduced a non-uniform ion density upstream to test if external corrugation suppresses the intrinsic instability.

3. Key Contributions

The paper makes several novel theoretical and observational contributions:

1. Identification of the Instability: The authors identify the driver of shock reformation in this geometry as the oblique lower-hybrid gradient drift instability. This instability arises from a reactive coupling between:
   • A backward-propagating oblique Whistler wave (in the electron rest frame).
   • A forward-propagating ion acoustic wave.
2. Mechanism of Reformation: Unlike previous models where reformation was driven by reflected ion beams (supercritical) or external perturbations, this study demonstrates that the oblique Whistler wave itself forces the shock to collapse into a "magnetic piston" and reform.
3. Growth Rate Discrepancy: The simulations reveal that the exponential growth rate of these waves is approximately 8 times faster than predicted by linear analytic theory (specifically the model by [34] for Harris-type sheaths). The authors attribute this to thermal anisotropy of electrons caused by magnetic compression at the shock.
4. Stability Criteria: The study establishes that if the upstream plasma is non-uniform (corrugated), the oblique Whistler instability is suppressed, and the shock reverts to oscillations driven by magnetic tension rather than intrinsic wave growth.

4. Key Results

• Wave Growth and Saturation:
  • Oblique Whistler waves grew in the magnetic overshoot region ahead of the density ramp.
  • The waves saturated and induced a periodic collapse and reformation of the shock on timescales of $\omega_{lh}^{-1}$ (lower-hybrid frequency), which is much faster than the ion cyclotron timescale ($\omega_{ci}^{-1}$) typical of supercritical shocks.
• Shock Dynamics:
  • The electrostatic modulation of ion density by the Whistler wave forces the shock front to collapse.
  • The shock transforms into a magnetic piston (mediated by Hall electric fields) before reforming into a new FMS shock.
  • This cycle creates a "standing wave" pattern in the downstream density if the shock is perturbed, but a traveling wave pattern if the shock is initially planar.
• Dispersion and Phase Velocity:
  • The phase velocity of the oblique Whistler wave in the simulation box frame was found to be supersonic relative to the ion acoustic speed ($v_{ph} \gg c_s$).
  • This high speed prevents resonant wave-particle interactions, confirming the instability is reactive (fluid-like) rather than resonant.
  • The instability only grows for wavenumbers where the phase speed is slower than the electron drift speed; higher harmonics were found to be stable.
• Effect of Upstream Perturbation:
  • When the upstream plasma density was modulated (simulating a corrugated shock front), the growth of the oblique Whistler wave was suppressed. The shock oscillations were then driven by magnetic tension, and the reformation frequency dropped significantly.

5. Significance

• Space and Astrophysical Relevance: The findings are highly relevant to the Earth's bow shock and other astrophysical collisionless shocks. The MMS mission has observed surface oscillations on quasi-perpendicular shocks; this paper provides a mechanism (oblique Whistler waves) that could explain these oscillations without requiring supercritical conditions or external triggers.
• Theoretical Advancement: It bridges the gap between magnetic reconnection physics (Harris sheet instabilities) and shock physics. It demonstrates that the oblique lower-hybrid drift instability is a universal mechanism for subcritical shock reformation when the magnetic field has a component resolved in the simulation plane.
• Implications for Modeling: The results suggest that the thermal anisotropy of electrons (a common feature in shocks) plays a critical role in accelerating instability growth beyond linear predictions. This implies that standard linear theories may underestimate the turbulence and dissipation rates in collisionless shocks.
• Future Directions: The authors propose that this mechanism might also suppress reformation in supercritical shocks if the Whistler waves grow fast enough to saturate before the shock collapses, potentially explaining "shock rippling" observed by MMS.

In summary, the paper establishes that oblique Whistler waves, driven by the diamagnetic current and a reactive instability mechanism, are the primary agents of cyclic reformation in subcritical perpendicular shocks when the magnetic field is oblique to the simulation plane, operating on timescales significantly faster than previously understood ion-driven mechanisms.