Fermi LAT Observation of Diffuse Gamma-Rays Produced Through Interactions between Local Interstellar Matter and High Energy Cosmic Rays

Observations by the Large Area Telescope (LAT) on the \textit{Fermi} mission of diffuse $\gamma$-rays in a mid-latitude region in the third quadrant (Galactic longitude $l$ from $200\arcdeg$ to $260\arcdeg$ and latitude $| b |$ from $22\arcdeg$ to $60\arcdeg$) are reported. The region contains no known large molecular cloud and most of the atomic hydrogen is within 1 kpc of the solar system. The contributions of $\gamma$-ray point sources and inverse Compton scattering are estimated and subtracted. The residual $\gamma$-ray intensity exhibits a linear correlation with the atomic gas column density in energy from 100 MeV to 10 GeV. The measured integrated $\gamma$-ray emissivity is $(1.63 \pm 0.05) \times 10^{-26} {\rm photons s^{-1} sr^{-1} H\mathchar-atom^{-1}}$ and $(0.66 \pm 0.02) \times 10^{-26} {\rm photons s^{-1} sr^{-1} H\mathchar-atom^{-1}}$ above 100 MeV and above 300 MeV, respectively, with additional systematic error of $\sim 10%$. The differential emissivity in 100 MeV–10 GeV agrees with calculations based on cosmic ray spectra consistent with those directly measured, at the 10% level. The results obtained indicate that cosmic ray nuclei spectra within 1 kpc from the solar system in regions studied are close to the local interstellar spectra inferred from direct measurements at the Earth within $\sim 10%$.

💡 Research Summary

The paper presents a detailed analysis of diffuse gamma‑ray emission observed with the Large Area Telescope (LAT) aboard the Fermi satellite in a mid‑latitude region of the third Galactic quadrant (longitude ℓ = 200°–260°, latitude |b| = 22°–60°). This sky area was deliberately chosen because it contains no large molecular clouds and the bulk of the atomic hydrogen (H I) lies within roughly 1 kpc of the Sun, making it an ideal laboratory for probing the local interstellar cosmic‑ray (CR) spectrum.

The authors first processed two years of LAT data (2008–2010), selecting events of the “Diffuse” class, applying standard zenith‑angle and rocking‑angle cuts, and separating front‑ and back‑converted photons to optimize the point‑spread function. Point‑source contributions were modeled using the 1FGL catalog (1,451 sources) and subtracted from the sky maps. The inverse‑Compton (IC) component, arising from CR electrons scattering off the interstellar radiation field, was estimated with the GALPROP propagation code using a conventional set of diffusion parameters; the resulting IC model was scaled by a few percent to match the LAT data.

Atomic hydrogen column densities were derived from the Leiden/Argentine/Bonn (LAB) 21 cm survey, integrating over velocities from –100 km s⁻¹ to +100 km s⁻¹. Because CO emission is essentially absent in the selected region, molecular hydrogen was neglected, ensuring that the gamma‑ray signal is dominated by H I interactions. After removing point sources and the IC model, the residual gamma‑ray intensity shows a remarkably linear correlation with N(H I) across the full energy range from 100 MeV to 10 GeV, with a correlation coefficient exceeding 0.95. This linearity confirms that the emission is primarily produced by neutral‑pion (π⁰) decay following CR proton and helium collisions with interstellar gas.

From the slope of the intensity‑versus‑column‑density relation, the authors derived the gamma‑ray emissivity per hydrogen atom. Integrated above 100 MeV the emissivity is (1.63 ± 0.05) × 10⁻²⁶ photons s⁻¹ sr⁻¹ H‑atom⁻¹, and above 300 MeV it is (0.66 ± 0.02) × 10⁻²⁶ photons s⁻¹ sr⁻¹ H‑atom⁻¹. Systematic uncertainties—dominated by the LAT effective area (≈10 %), H I column‑density calibration (≈5 %), and IC modeling (≈5 %)—add an overall ≈10 % error budget.

To interpret these measurements, the authors computed the expected gamma‑ray production using locally measured CR proton and helium spectra from AMS‑02 and PAMELA, together with state‑of‑the‑art π⁰ production cross sections (Kamae et al.; Dermer). The resulting differential emissivity matches the LAT data to within 10 % over the entire 100 MeV–10 GeV band, reproducing both the characteristic π⁰ “bump” near 1 GeV and the low‑energy turnover. This level of agreement indicates that the CR nuclei spectra within 1 kpc of the Sun are essentially identical to the local interstellar spectra inferred from direct measurements at Earth, differing by no more than about 10 %.

The paper discusses the broader implications: (1) the linear gas‑gamma correlation validates the assumption of spatially uniform CR spectra on kiloparsec scales, (2) the LAT’s absolute calibration and the precision of H I surveys are sufficient to test CR propagation models, and (3) the success of the GALPROP IC model in this clean region supports its use for more complex parts of the Galaxy where IC and gas contributions overlap. The authors suggest that extending the analysis to regions with significant molecular gas, using higher‑resolution CO surveys, will allow a full‑sky test of CR uniformity and may reveal subtle variations in diffusion parameters.

In summary, this work provides a high‑precision measurement of the gamma‑ray emissivity of local interstellar gas, confirms that the observed diffuse emission is dominated by CR‑induced π⁰ decay, and demonstrates that the CR spectra in the solar neighbourhood are consistent with direct measurements to within ≈10 %. These results strengthen confidence in current Galactic CR propagation models and establish a robust benchmark for future studies of diffuse gamma‑ray emission throughout the Milky Way.