Supernovae drive large-scale, incompressible turbulence through small-scale instabilities

This paper demonstrates that core-collapse supernovae sustain galactic turbulence by generating incompressible modes through baroclinic vorticity at unstable contact discontinuities, where small-scale surface instabilities create a $k^{-3/2}$ spectrum that is transferred to kiloparsec scales via an inverse cascade mechanism.

The Gist

The Big Picture: How the Galaxy Gets "Shaken"

Imagine our galaxy, the Milky Way, as a giant, swirling bowl of soup. This soup isn't just liquid; it's a mix of hot gas, cold gas, dust, and stars. For the galaxy to stay healthy and form new stars, this soup needs to be constantly stirred. If it sits still, it settles and stops working.

Scientists have long known that supernovae (exploding stars) are the spoons that stir this soup. When a star explodes, it blasts out huge amounts of energy. But there was a big mystery: How exactly does a single explosion turn into a galaxy-wide swirl?

Previous theories suggested that the turbulence (the swirling) happened because different explosions crashed into each other, or because the shockwaves hit bumpy patches of gas.

This paper says: "No, that's not the main story."

The author, James Beattie, argues that the turbulence is actually generated by a specific, tiny, unstable layer inside the explosion itself, before it even hits anything else.

---

The Core Discovery: The "Wrinkled Balloon" Effect

1. The Setup: A Hot Bubble in a Cold Room

Imagine you blow up a hot air balloon (the supernova remnant) inside a cold room (the interstellar medium).

• Inside the balloon: It's scorching hot.
• Outside: It's cool.
• The Skin: There is a thin, fragile membrane separating the hot inside from the cool outside.

2. The Instability: The "Wrinkling"

As the balloon expands, this thin skin becomes unstable. It doesn't stay smooth. It starts to wrinkle, fold, and ripple like a crumpled piece of aluminum foil.

• The Paper's Insight: The author found that these wrinkles aren't just random noise. They are a specific type of instability (like a drum skin vibrating) that creates a very specific kind of "spin" or vortex.
• The Mechanism: Because the pressure and density change so sharply at this wrinkled skin, it generates a force called baroclinicity. Think of this as a "twist generator." The more the skin wrinkles, the more it twists the gas around it.

3. The Analogy: The "Ripple in the Pond"

Imagine dropping a stone in a pond. Usually, we think the ripples come from the stone hitting the water.

• Old Theory: The ripples happen because the stone hits a second stone, or the water is already choppy.
• New Theory: The ripples happen because the surface of the water itself becomes unstable and starts folding over on itself immediately after the stone hits, creating its own internal spin that spreads out.

---

The "Magic" Math: Why It Matters

The paper uses some heavy math to prove two things:

1. The "Recipe" for Turbulence
The author found a perfect mathematical recipe connecting the wrinkles on the balloon skin to the swirls in the gas.

• The wrinkles follow a specific pattern (a power law).
• The resulting swirls follow a matching pattern.
• The Result: The paper proves that the "wrinkled skin" is the only thing needed to create the swirling motion. You don't need other explosions to crash into it. The explosion creates its own engine for turbulence.

2. The "Shedding" Process
Here is the most creative part of the analogy.

• Imagine the wrinkles on the balloon are like 2D ripples moving only along the surface of the skin.
• The paper shows that as the balloon expands, these 2D ripples get "stretched" by the wind blowing out from the center of the explosion.
• The Shedding: This stretching pulls the ripples off the surface and throws them into the surrounding space, turning them into 3D turbulence that fills the galaxy.
• The Timing: This works best when the balloon is "young" (just after the explosion). As the balloon gets old and slows down, it stops shedding these ripples efficiently.

---

Why This Changes Everything

The "Galactic Engine"
This paper suggests that the entire galaxy's turbulence is built from the bottom up.

1. A star explodes.
2. A thin, unstable layer forms around the explosion.
3. This layer wrinkles and folds (creating 2D turbulence).
4. The explosion's own wind stretches these folds and flings them outward (creating 3D turbulence).
5. These 3D swirls travel outward, merging with other swirls to create the massive, galaxy-wide currents we see.

The "Inverse Cascade"
Usually, in physics, big things break into small things (like a big wave breaking into foam). This paper suggests the opposite happens here: Small things build up into big things.
The tiny, unstable wrinkles on the surface of a single supernova remnant are the seeds that grow into the massive, galaxy-spanning currents.

Summary in One Sentence

Supernovae don't just stir the galaxy by crashing into things; they create a thin, wrinkled skin around themselves that acts like a factory, manufacturing tiny twists that get flung out to spin the entire galaxy.

This discovery helps us understand how galaxies mix their ingredients (like heavy metals and dust), how they generate magnetic fields, and how they transport energy across vast distances, all starting from the microscopic instability of a single exploding star's shell.

Technical Summary

1. Problem Statement

The origin of turbulence in the interstellar medium (ISM) of galaxies is a fundamental problem in astrophysics. While core-collapse supernovae (SNe) provide sufficient energy to sustain galactic turbulence, the specific mechanism by which they generate incompressible (solenoidal) turbulence remains debated.

• Prevailing Theory: It is widely believed that turbulence arises from the interaction of multiple supernova remnants (SNRs), curved shock waves, or collisions between hot superbubbles and inhomogeneities in the ISM.
• The Gap: Previous studies suggest that isolated SN blast waves alone are insufficient to excite even weak acoustic turbulence, let alone a Kolmogorov-style cascade.
• The Challenge: This paper challenges the prevailing view by proposing that isolated SNRs, specifically through instabilities in their cooling shells, are the primary source of incompressible turbulence, without requiring interactions with other SNRs or a pre-existing inhomogeneous medium.

2. Methodology

The study utilizes high-resolution numerical simulations and a novel statistical analysis of extracted SNRs.

• Numerical Simulations:

  • Code: RAMSES (Teyssier 2002).
  • Setup: A gravitohydrodynamical, stratified, multiphase model of a Milky Way-analogue galactic disk.
  • Physics: Ideal hydrodynamics (no explicit viscosity), time-dependent cooling network (modeling WNM, WIM, HIM phases), and a supernova explosion rate of $0.1$ per $10^4$ years.
  • Exclusions: Self-gravity, magnetic fields, cosmic rays, and large-scale galactic shear are excluded to isolate the impact of SN explosions on turbulence.
  • Resolution: $1024^3$ grid cells in a $1 \text{ kpc}^3$ domain, resolving scales down to $\sim 10$ pc.

• SNR Extraction and Analysis:

  • Cataloging: The author identified $\approx 100$ SNRs (specifically $n=87$) from the global simulation by isolating hot bubbles ($T \gg 10^6$ K) using a temperature mask, morphological filtering (binary opening), and a "friend-of-friends" algorithm.
  • Local Analysis: Instead of analyzing the global domain, the study focuses on $128^3$ grid-cell regions ($\ell_0 = 125$ pc) centered on each SNR to analyze local Fourier spectral properties.
  • Spectral Definitions:
    • Incompressible Velocity Spectrum ($P_{us}(k)$): Derived via Helmholtz decomposition.
    • Baroclinic-Vorticity Cospectrum ($P_{\omega B}(k)$): Measures the transfer of energy from baroclinicity ($\nabla \rho \times \nabla P / \rho^2$) to vorticity ($\omega$).
    • Spherical Harmonics: Used to analyze the corrugations of the contact discontinuity surface.

3. Key Contributions and Results

A. Baroclinicity as the Primary Driver

The study demonstrates that baroclinicity generated at the unstable contact discontinuity (the thin shell between hot and warm plasma) is the dominant source of incompressible turbulence.

• Volume vs. Contribution: Although the unstable shell occupies a tiny fraction of the volume, it accounts for $\approx 100\%$ of baroclinic production at the onset of turbulence and $\gtrsim 70\%$ in the steady state.
• Rejection of Alternatives: Shock-shock interactions and curved shocks contribute only $\approx 10\%$ to the total baroclinicity, refuting the idea that they are the primary drivers.
• Analytical Relation: The author derives a relation between the driving spectrum and the resulting turbulence:
  $$P_{\omega B}(k) \sim k^3 P_{us}^{3/2}(k)$$
  This relation holds perfectly in the data, confirming that the baroclinic term directly fuels the vorticity cascade.

B. The $k^{-3/2}$ Spectrum

The paper identifies a specific power-law spectrum for the incompressible velocity modes generated by the unstable layer:
$$P_{us}(k) \propto k^{-3/2}$$

• This spectrum is observed locally within individual SNRs and matches the spectrum found in global galactic disk simulations.
• The driving spectrum (baroclinicity) follows $P_{\omega B}(k) \propto k^{3/4}$, peaking at high wavenumbers ($k$), indicating that turbulence is driven by small-scale instabilities rather than a single large-scale driving scale.

C. Mechanism of Turbulence Shedding (Vortex Stretching)

A critical question is how 2D surface modes on the shell become 3D turbulence in the surrounding medium.

• Mechanism: Radial velocities from the SN explosion ($u_r$) stretch the surface vorticity ($\omega_{\theta, \phi}$) into the third dimension via vortex stretching.
• Timescale Criterion: The author derives a criterion for when modes are "shed" into the ISM versus confined to the shell:
  $$\frac{t_{nl}}{t_{3D}} \sim \frac{u_r}{u_t}$$
  Where $t_{nl}$ is the nonlinear timescale of surface modes and $t_{3D}$ is the vortex stretching timescale.
• Age Dependence:
  • Young SNRs ($R_c \lesssim 50$ pc, $t_{age} \approx t_{cool}$): $t_{nl}/t_{3D} > 1$. These are highly efficient at radiating turbulence into the surrounding medium.
  • Old SNRs: $t_{nl}/t_{3D} < 1$. Modes remain confined to the shell and cascade within the 2D layer.

D. Evidence of 2D Surface Turbulence

Analysis of the spherical harmonics of the shell corrugations ($C(\ell)$) reveals a power law:
$$C(\ell) \propto \ell^{-8/3}$$

• This matches the prediction for 2D Kraichnan turbulence ($\ell \sim k_\perp R_c$).
• This provides empirical evidence that the unstable shell supports a self-generated 2D turbulence cascade, which seeds the baroclinic instability.

4. Significance and Implications

1. New Phenomenology for SN-Driven Turbulence: The paper establishes that large-scale galactic turbulence does not require complex interactions between multiple SNRs or a turbulent ISM. Instead, it originates from small-scale instabilities (Vishniac instability, 2D surface turbulence) on the cooling shell of isolated SNRs.
2. Inverse Cascade Connection: The study links the local $k^{-3/2}$ spectrum of the SNR shell to the global galactic turbulence cascade. Through the inverse cascade mechanism (previously identified in Beattie et al. 2025c), the small-scale energy generated at the cooling radius is transported to kiloparsec scales, imprinting the $k^{-3/2}$ spectrum on the entire galaxy.
3. Resolution of the "Driving Scale" Problem: The findings suggest that SN-driven turbulence has no single driving scale; rather, it is driven by a spectrum peaking at high-$k$ (small scales) due to the fractal nature of the unstable shell.
4. Future Directions: The results highlight the importance of the Biermann battery effect (generation of magnetic fields from baroclinicity) in these shells, suggesting that future MHD simulations are needed to fully characterize the magnetic dynamo action in these environments.

Conclusion

James R. Beattie provides robust evidence that core-collapse supernovae drive galactic turbulence primarily through baroclinic vorticity generation at the unstable contact discontinuity of isolated SNRs. The process involves the formation of a fractal, 2D turbulent shell that generates a $k^{-3/2}$ incompressible spectrum, which is then ejected into the interstellar medium via vortex stretching in young SNRs and transported to large scales via an inverse cascade. This redefines the understanding of how small-scale hydrodynamic instabilities shape the large-scale dynamics of galaxies.