🔭 astrophysics

Interplay between Escaping Cosmic Rays and Interstellar Medium: Driving of Galactic Winds and Shaping the Local Proton Spectrum

This paper utilizes spherically symmetric CR-hydrodynamical simulations to demonstrate that escaping cosmic rays can drive galactic winds and heat the interstellar medium, potentially explaining anomalous emission lines in old supernova remnants and offering an alternative interpretation for the local proton spectrum observed by CALET, AMS02, and Voyager 1.

Original authors: Jiro Shimoda, Katsuaki Asano, Shu-ichiro Inutsuka

Published 2026-04-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jiro Shimoda, Katsuaki Asano, Shu-ichiro Inutsuka

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine our galaxy, the Milky Way, not as a static swirl of stars, but as a bustling, chaotic city. In this city, Supernova Remnants (SNRs) are like massive, ancient construction sites where stars have exploded. These explosions don't just send out a shockwave of debris; they also launch a flood of invisible, high-speed particles called Cosmic Rays (CRs).

For a long time, scientists thought these cosmic rays just drifted through the galaxy like smoke in a gentle breeze, spreading out evenly. But this paper suggests a much more dramatic story: Cosmic rays are actually the wind that drives the galaxy's weather.

Here is the simple breakdown of what the authors, Jiro Shimoda and his team, discovered:

1. The "Traffic Jam" of Particles

Imagine you are at a crowded party. If everyone walks out the door at a normal pace (a high diffusion coefficient), they spread out quickly, and the room empties fast. But what if the exit is narrow, or the crowd is so dense they can't move freely? They get stuck, piling up right near the door.

The authors studied what happens when cosmic rays get "stuck" near their source (the supernova) because the space around them is turbulent and hard to navigate (a low diffusion coefficient).

The Result: Instead of drifting away, these particles pile up, creating a massive, invisible "pressure cooker." This pressure is so strong that it physically pushes the gas in space (the Interstellar Medium) away, creating a powerful wind.

2. The "Galactic Fountain" vs. The "Galactic Wind"

Usually, when a star explodes, it heats up the gas around it. But that gas is heavy, and gravity pulls it back down to the galaxy's disk, like water from a fountain falling back into the pool. This is called the "Galactic Fountain."

However, the authors found that if cosmic rays get trapped and build up enough pressure, they act like a jet engine. They don't just push the gas; they heat it and accelerate it so fast that it escapes the galaxy's gravity entirely.

The Analogy: Think of it like a pressure cooker releasing steam. If the steam (cosmic rays) builds up enough pressure, it blasts the lid (the gas) off the pot and shoots it into the sky.
The Impact: This "wind" carries gas and heavy metals out of the galaxy and into the vast space between galaxies. The amount of gas being blown away is roughly equal to the amount of new stars being born. This solves a mystery: How does the galaxy get rid of all that extra gas and metal? The cosmic rays are the exhaust pipe.

3. The "Glowing Shell" Mystery

The paper also explains a weird observation made by astronomers looking at old supernova remnants. They see two types of glowing gas:

H-alpha: A red glow from cooler gas.
[OIII]: A blue-green glow from very hot, highly energized gas.

Usually, you'd expect the hottest gas to be right next to the explosion. But in these old remnants, the hot gas is actually on the outside, and the cooler gas is on the inside. It's like a fire where the smoke is hotter than the flames!

The Paper's Explanation:
The trapped cosmic rays act like a microwave oven for the gas. As they push the gas outward, they generate waves (Alfvén waves) that heat the gas from the inside out.

The gas right at the edge of the expanding bubble gets hit by the most intense cosmic ray "wind," heating it up to extreme temperatures (creating the [OIII] glow).
The gas closer to the center doesn't get hit as hard, so it stays cooler (creating the H-alpha glow).
This "hot shell, cool core" structure perfectly matches what telescopes are seeing.

4. Why We See What We See on Earth

Finally, the authors looked at the cosmic rays hitting Earth right now. Standard theories say these particles come from all over the galaxy, mixed together like a smoothie.

But the authors propose a different idea: We are living in a "Local Bubble."

The Analogy: Imagine you are standing in a quiet cul-de-sac (the Local Bubble) surrounded by a few loud neighbors (nearby supernovae). Because the "traffic" (diffusion) in your neighborhood is bad, the noise from your loud neighbors drowns out the noise from the rest of the city.
The Result: The cosmic rays hitting Earth aren't a mix from the whole galaxy. They are mostly from just a few nearby explosions that happened a few million years ago. Because these particles got "stuck" and interacted with the local gas, their energy spectrum (the mix of fast vs. slow particles) looks different than we expected. This explains some recent, puzzling measurements from space telescopes like AMS-02 and CALET.

The Big Picture

This paper suggests that Cosmic Rays are the architects of our galaxy's environment. They aren't just passive passengers; they are the active force that:

Blows gas out of the galaxy (driving the Galactic Wind).
Heats up the gas in specific patterns (explaining the strange glowing shells).
Dictates what kind of radiation we detect on Earth (by trapping particles in our local neighborhood).

It's a shift from thinking of cosmic rays as just "space dust" to realizing they are the engine that shapes the life and death of the galaxy.

1. Problem Statement

The paper addresses two critical, interconnected issues in Galactic astrophysics:

The Galactic Wind Mechanism: While Galactic winds are known to transport mass and metals from the disk to the halo (impacting star formation and halo metallicity), the specific launching mechanism remains unclear. Standard Supernova Remnant (SNR) models often fail to drive sustained outflows because radiative cooling causes the gas to fall back to the disk (the "galactic fountain") before reaching significant heights.
The Local Cosmic Ray (CR) Spectrum: The origin and propagation of Galactic CRs are traditionally modeled using a diffusion coefficient ( $D \sim 10^{28} \text{ cm}^2 \text{ s}^{-1}$ ) averaged over the Galactic halo. However, recent observations (e.g., TeV halos, $^{60}\text{Fe}$ in ocean crusts) suggest that diffusion may be locally suppressed near sources. The authors investigate how this local suppression and the resulting CR-ISM interaction affect the CR spectrum observed at Earth.

2. Methodology

The authors employ spherically symmetric CR-hydrodynamical simulations to model the evolution of the Interstellar Medium (ISM) driven by escaping CRs from a single source.

Physical Framework:
- Coupled Equations: They solve the continuity, momentum, and energy conservation equations for the thermal gas, coupled with a transport equation for the CR momentum distribution function $N_{cr}(p)$ .
- CR Effects: The model includes CR pressure ( $P_{cr}$ ), adiabatic compression/expansion, collisional energy losses (ionization/Coulomb), and CR-induced heating via the dissipation of Alfvén waves ( $|V_A \partial_r \epsilon_{cr}|$ ).
- Thermal Physics: Radiative cooling (collisional ionization equilibrium) and thermal conduction are included to realistically model gas thermodynamics.
Initial Conditions:
- An idealized scenario where a uniform ISM ( $T_g = 7000 \text{ K}$ , $n_g \approx 0.2 \text{ cm}^{-3}$ ) is suddenly injected with a population of CRs ( $E_{cr} = 10^{50} \text{ erg}$ ) within a radius $r_0 = 5 \text{ pc}$ .
- The CR spectrum follows a power law ( $p^{-2.3}$ ).
Key Parameter: The study systematically varies the CR diffusion coefficient $D(p) = D_0 (p/m_pc)^\delta$ $D (p) = D_{0} (p / m_{p} c)^{δ}$ , focusing on the normalization $D_0$ $D_{0}$ .
- Standard Case: $D_0 = 10^{28} \text{ cm}^2 \text{ s}^{-1}$ .
- Suppressed Cases: $D_0 = 10^{27}$ and $10^{26} \text{ cm}^2 \text{ s}^{-1}$ (simulating local environments like the Local Bubble or TeV halos).
- Spectral Index: $\delta = 1/3$ (Kolmogorov turbulence).

3. Key Contributions & Results

A. Acceleration of the ISM and Galactic Winds

Diffusion Dependence: The efficiency of gas acceleration is inversely related to the diffusion coefficient.
- Low Diffusion ( $D_0 = 10^{26} \text{ cm}^2 \text{ s}^{-1}$ ): CRs remain confined longer, creating a steep pressure gradient. This drives a strong shock, accelerating the ISM to supersonic velocities ( $>10 \text{ km s}^{-1}$ ).
- High Diffusion ( $D_0 = 10^{28} \text{ cm}^2 \text{ s}^{-1}$ ): CRs escape rapidly, resulting in a weaker pressure gradient and less efficient gas acceleration.
Mass Loss Rate: For a single source, the model expels $\sim 10^3 M_\odot$ $\sim 1 0^{3} M_{⊙}$ of gas. Scaling this to the Galactic supernova rate ( $\sim 0.01 \text{ yr}^{-1}$ $\sim 0.01 yr^{- 1}$ ), the total mass loss rate is $\sim 10 M_\odot \text{ yr}^{-1}$ $\sim 10 M_{⊙} yr^{- 1}$ .
- Significance: This rate is comparable to the Galactic Star Formation Rate (SFR), suggesting CR-driven outflows are a viable mechanism for driving the Galactic wind and polluting the circumgalactic medium with metals.

B. CR Heating and Optical Line Emission

Heating Mechanism: The dissipation of Alfvén waves generated by CR streaming heats the gas, particularly at the diffusion front (the edge of the expanding cavity).
Temperature Structure:
- In the suppressed diffusion case, the gas inside the cavity reaches temperatures of $\sim 10^5 \text{ K}$ .
- A distinct "bump-tail" structure emerges: a high-temperature region (bump) at the outer edge and a cooler region (tail) inside.
Observational Signature: This structure explains recent observations of old SNRs (e.g., by Fesen et al. 2024) where:
- [OIII] $\lambda$ 5007 (requiring $T \sim 10^4\text{--}10^5 \text{ K}$ ) is bright at the outer filaments.
- H $\alpha$ (requiring $T \sim 10^4 \text{ K}$ ) is bright at the inner filaments.
- Standard shock models predict the opposite ionization gradient; CR heating provides a natural explanation for this inversion.

C. Shaping the Local CR Proton Spectrum

Local Suppression Scenario: Assuming the Solar System resides in the Local Bubble with a suppressed diffusion coefficient ( $D_0 = 10^{26} \text{ cm}^2 \text{ s}^{-1}$ ), the effective residence time of CRs is drastically reduced compared to the standard halo model.
Nearby Sources: The authors propose that the observed CR spectrum at Earth is dominated by a few nearby, recent supernovae (within $\sim 100\text{--}300 \text{ pc}$ ) rather than a superposition of distant sources.
Spectral Features: By modeling three distinct sources with ages of 1.2, 2.9, and 6.6 Myr at varying distances, the model successfully reproduces:
- The observed hardening of the spectrum around $\sim 600 \text{ GeV}$ .
- The softening at $\sim 10 \text{ TeV}$ .
- The data from CALET, AMS-02, and Voyager 1.
Grammage: The model estimates a grammage (column density traversed) consistent with Boron-to-Carbon ratios, provided backflow effects are considered.

4. Significance and Implications

Unified Explanation: The paper provides a unified physical framework where the same mechanism (CR-ISM interplay with suppressed diffusion) explains both the driving of Galactic winds (solving the mass-loading problem) and the anomalies in the local CR spectrum.
Revising Diffusion Models: It challenges the assumption of a uniform diffusion coefficient, suggesting that local suppression near sources is critical for understanding both ISM dynamics and CR propagation.
New Diagnostic Tools: The predicted "bump-tail" temperature structure offers a new diagnostic for identifying CR-heated regions in old SNRs via optical line ratios ([OIII]/H $\alpha$ ).
Galactic Evolution: The results support scenarios where CRs play a fundamental role in regulating the baryon cycle of the Galaxy, linking star formation, metal enrichment of the halo, and the chemical evolution of the Local Bubble.

In conclusion, Shimoda et al. demonstrate that escaping cosmic rays are not merely passive tracers but active agents that shape the thermal and dynamical state of the ISM, potentially resolving long-standing discrepancies in Galactic wind theory and local CR spectral observations.