🔭 astrophysics

The Diffuse Gamma-ray Sky of a Milky Way Analogue: Local Diversity and Global Constraints

Using the Rhea suite of CR-MHD simulations, this study demonstrates that a Milky Way analogue naturally reproduces observed diffuse gamma-ray luminosities and spectral properties, revealing that while gas density fluctuations primarily dictate emission morphology dependent on observer location, cosmic ray transport physics governs the spectral characteristics and angular power spectra.

Original authors: Karin Kjellgren, Philipp Girichidis, Maria Werhahn, Ralf S. Klessen, Christoph Pfrommer, Juan Soler, Brian Reville, Jim Hinton, Patrick Hennebelle, Noé Brucy, Simon C. O. Glover

Published 2026-03-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Karin Kjellgren, Philipp Girichidis, Maria Werhahn, Ralf S. Klessen, Christoph Pfrommer, Juan Soler, Brian Reville, Jim Hinton, Patrick Hennebelle, Noé Brucy, Simon C. O. Glover

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 the Milky Way galaxy not as a static picture, but as a bustling, glowing city at night. Now, imagine that the "lights" we see aren't from streetlamps or houses, but from invisible, high-speed particles called Cosmic Rays crashing into the air (gas) around them. When these particles hit the gas, they create a flash of light called a gamma-ray.

This paper is like a team of astronomers building a super-advanced, 3D video game simulation of a galaxy that looks just like ours. They wanted to answer a big question: If we were standing in different neighborhoods of this galaxy, what would the "gamma-ray sky" look like?

Here is the story of their discovery, broken down into simple concepts:

1. The Invisible City and the "Flashbulbs"

Think of the galaxy as a city filled with gas (the air) and stars. Cosmic rays are like invisible, super-fast bullets flying everywhere.

The Collision: When a cosmic ray bullet hits a gas molecule, it creates a tiny explosion that releases a gamma-ray photon (a flash of light).
The Problem: We live inside this galaxy, specifically in a quiet neighborhood called the "Local Bubble" (a giant, empty cavity created by old supernova explosions). Because we are stuck in one spot, we can only see the sky from our backyard. We don't know if the view looks different if you stand on the other side of the city.

2. The Simulation: A "Time-Travel" Camera

The researchers used a powerful computer simulation (called the Rhea suite) to build a fake galaxy that behaves exactly like the real one.

The Magic Trick: Instead of just looking at the galaxy from the outside, they placed "virtual observers" in 18 different locations inside the simulation. Some were in crowded, dense gas clouds; others were in quiet, empty bubbles like ours.
The Result: They generated a "gamma-ray map" for each observer, showing exactly what the sky would look like from that specific spot.

3. The Big Surprise: Your Neighborhood Matters Most

The most exciting finding is that where you stand changes everything.

The Analogy: Imagine looking at a city skyline. If you stand in a park, you see a clear view of the tall skyscrapers in the distance. If you stand in a dense forest, you only see the trees right next to you, and the city is hidden.
The Discovery: The researchers found that the "gamma-ray sky" looks completely different depending on your location.
- Near the Galactic Center: The sky is bright and chaotic.
- In a "Local Bubble" (like ours): The sky is dominated by the immediate surroundings. The features you see (filaments and loops of light) are mostly caused by gas right next to you (within 2,000 light-years), not by the whole galaxy.
- The Takeaway: The "shape" of the gamma-ray sky is mostly determined by the local gas density, not by the total amount of cosmic rays in the galaxy.

4. The "Fog" vs. The "Lamps"

The paper also investigated what creates the patterns in the sky.

Cosmic Rays (The Fog): These particles are like a thick, smooth fog that spreads out evenly. They don't have sharp edges.
Gas (The Lamps): The gas is clumpy, like streetlamps or buildings.
The Conclusion: The gamma-ray sky looks like the gas, not the cosmic rays. The "fog" of cosmic rays is everywhere, but the "flash" only happens where the "lamps" (gas) are. So, the patterns we see in the sky are actually a map of the gas clouds, not a map of the cosmic rays themselves.

5. Matching the Real World

Finally, they compared their fake galaxy to the real data we have from telescopes like Fermi-LAT.

The Score: Their simulation matched the real universe surprisingly well!
The Secret Sauce: To get the match right, they had to tune how fast the cosmic rays move. They found that cosmic rays move faster at higher energies (a specific mathematical rule called $\delta = 0.5$ ). If they used a slower rule, the simulation didn't match reality.

Why Does This Matter?

This paper is like a guidebook for understanding our place in the universe.

It explains our view: It tells us that the gamma-ray sky we see is heavily influenced by our local neighborhood (the Local Bubble). We can't just look at the sky and assume we see the whole galaxy; we have to account for the "fog" right outside our window.
It validates our physics: The fact that a computer simulation, built from first principles (without cheating or tweaking numbers to fit the data), produced a sky that looks like the real one proves that our understanding of how cosmic rays and gas interact is correct.

In a nutshell: The universe is a giant, glowing city. The light we see isn't just a uniform glow; it's a complex pattern shaped by the local "streets" and "buildings" (gas clouds) around us. To understand the whole city, we first have to understand our own backyard.

1. Problem Statement

Interpreting the diffuse gamma-ray emission of the Milky Way is challenging because the observed signal depends on a complex interplay between the three-dimensional distribution of cosmic rays (CRs), the interstellar medium (ISM) gas distribution, and the specific location of the observer within the Galaxy. While instruments like Fermi-LAT provide high-resolution all-sky maps, deriving global properties is difficult due to our position inside the Local Bubble and the inability to view the Galaxy from an external perspective. Previous studies often relied on static models or simplified CR propagation assumptions. This work aims to bridge the gap by using fully self-consistent, time-evolving magnetohydrodynamic (MHD) simulations that include CR feedback to investigate how local environmental diversity and global CR transport physics shape the gamma-ray sky.

2. Methodology

Simulation Setup (Rhea Suite):

Code: The authors used the moving-mesh code AREPO to simulate isolated, Milky Way-analogue galaxies.
Physics: The simulations include a self-consistent treatment of star formation, supernova (SN) feedback, magnetic fields, and CRs.
- CRs: Treated in a "gray" approximation, evolving the total integrated CR proton energy density. CRs are injected at 10% of SN explosion energy and transported via advection and diffusion along magnetic field lines.
- Magnetic Fields: Initialized as toroidal fields (3 $\mu$ G in the fiducial CRMHD run).
- Chemistry: Uses the NL97 non-equilibrium chemical network to track atomic, ionized, and molecular hydrogen.
Resolution: Gas mass resolution of $3000 M_\odot$ , with volume refinement in the disk limiting cell sizes to $(100 \text{ pc})^3$ . Simulations evolved for 2 Gyr.

Gamma-Ray Emission Modeling:

Post-Processing: Gamma-ray emission was calculated using crayon+, a code that computes steady-state CR spectra and gamma-ray emissivity in every computational cell.
Mechanism: The study focuses exclusively on neutral pion decay ( $\pi^0 \to \gamma\gamma$ ) resulting from CR proton interactions with thermal gas. Leptonic processes (Inverse Compton and bremsstrahlung) were excluded, as pion decay dominates the 0.1–100 GeV range.
CR Transport Parameters: The diffusion coefficient was modeled as $D(E) = D_0 (E/E_0)^\delta$ . The fiducial value used was $\delta = 0.5$ , with a comparison run using $\delta = 0.3$ .
Observer Perspective: The authors generated all-sky maps (Mollweide projections) and angular power spectra (APS) for multiple observer positions located within low-density bubbles (mimicking the Local Bubble) at various galactocentric radii.

3. Key Contributions

Self-Consistent Framework: This is one of the first studies to generate gamma-ray sky maps from a galaxy simulation where CRs dynamically drive outflows and shape the ISM, rather than using static gas distributions.
Local Diversity Analysis: The paper systematically quantifies how the gamma-ray sky morphology changes for observers at different locations, specifically within Local-Bubble-like cavities.
Decoupling Gas and CR Structure: By comparing the Angular Power Spectra (APS) of the gamma-ray emission, gas column density, and CR energy density, the authors demonstrate which physical component primarily drives the observed structure.
Validation of Transport Physics: The study tests different CR diffusion scaling indices ( $\delta$ ) against observational data to determine which transport model best reproduces the Milky Way's spectral and structural properties.

4. Key Results

Global Luminosity and Spectra:

The simulations naturally reproduce Milky Way-like gamma-ray luminosities and spectral slopes without fine-tuning source distributions.
The total luminosity is relatively stable over time, but the spectral shape is highly sensitive to the diffusion coefficient scaling.
$\delta = 0.5$ provides the best match to Fermi-LAT observations and GALPROP models, accurately reproducing the high-energy slope.
$\delta = 0.3$ significantly overestimates luminosity at energies $>1$ GeV.

Morphology and Local Environment:

Observer Dependence: The detailed morphology of the gamma-ray sky varies drastically depending on the observer's location. Even observers in nearby bubbles (separated by a few hundred parsecs) see different filamentary structures and loops.
Local Dominance: At high galactic latitudes ( $|b| > 20^\circ$ ), the gamma-ray emission is dominated by gas within $\le 2$ kpc of the observer. The local environment (e.g., nearby superbubbles, spurs) dictates the sky structure in these regions.
Gas vs. CR Structure: The Angular Power Spectrum (APS) of the gamma-ray emission closely follows the gas column density rather than the CR energy density. Because CRs diffuse and smooth out gradients, the gas distribution (which is clumpy) is the primary driver of the emission morphology.

Comparison with Observations:

All-Sky Spectrum: The simulated hadronic component matches the Fermi-LAT data well, lying slightly below the total observed flux (as expected, given the exclusion of leptonic components).
Regional Spectra: The model reproduces spectra for specific sky regions (e.g., inner Galaxy, high latitudes) observed by Fermi-LAT, LHAASO, and Tibet AS $\gamma$ .
APS Comparison:
- On large angular scales, the simulation matches observations well.
- On small scales (high multipoles), the simulation underestimates power compared to observations. This is attributed to limited simulation resolution (cell sizes $\sim 75-100$ pc) and the steady-state approximation, which smooths out small-scale spectral variations.
- Masking the galactic midplane ( $|b| > 5^\circ$ ) improves the agreement between the simulation and observations, particularly for the $\delta = 0.3$ model, suggesting that CR transport physics may vary with latitude.

5. Significance and Conclusions

Validation of Simulations: The Rhea simulations demonstrate that fully self-consistent CR-MHD galaxy models can naturally reproduce the global gamma-ray properties of the Milky Way, validating the use of such simulations for interpreting high-energy astrophysics.
Role of Local Environment: The study confirms that the "local" view of the gamma-ray sky is not representative of the global galaxy. The complex distribution of gas in the immediate vicinity of the observer is the dominant factor in shaping the morphology of the emission, especially at high latitudes.
CR Transport Constraints: The results strongly favor a diffusion coefficient energy scaling of $\delta \approx 0.5$ (consistent with AMS-02 secondary-to-primary ratios) over $\delta = 0.3$ for reproducing the observed spectral slope and luminosity.
Future Directions: The paper highlights the need for higher-resolution simulations to capture small-scale structure and suggests that future work should include leptonic emission (Inverse Compton) to fully explain the high-latitude and high-energy gamma-ray sky.

In summary, this work establishes a robust, physics-based framework for interpreting diffuse gamma-ray observations, emphasizing that understanding the observer's local surroundings is as critical as understanding global CR transport parameters.