Large-scale Galactic diffuse gamma rays observed with the Fermi Gamma-Ray Space Telescope

The LAT instrument on the Fermi Gamma-Ray Space Telescope is performing an all-sky survey from 20 MeV to 300 GeV with unprecedented statistics and angular resolution. This is providing a wealth of new information on the non-thermal emission from the Galactic interstellar medium with implications for cosmic rays and Galactic structure. First results at intermediate latitudes have already shown good agreement with predictions based on direct measurements of cosmic rays, suggesting that at least the local (within about 1 kpc from the Sun) gamma-ray emission is understood. We will present the first spectra from regions over the sky using the LAT data, and profiles for selected energies. The aim here is to evaluate the agreement with the models and assess what we can expect to learn as this analysis matures.

💡 Research Summary

The paper presents the first large‑scale analysis of diffuse Galactic gamma‑ray emission using data from the Large Area Telescope (LAT) aboard the Fermi Gamma‑Ray Space Telescope. LAT continuously surveys the entire sky in the 20 MeV–300 GeV energy range with unprecedented photon statistics and an angular resolution of roughly 0.1°, enabling a detailed study of the non‑thermal radiation that pervades the Milky Way’s interstellar medium.

The authors first describe the data selection and processing pipeline. Events collected over the first six months of the mission are filtered to retain only the cleanest photon class, and contaminations from Earth albedo, solar flares, and instrumental backgrounds are removed. The resulting photon list is binned in energy (four broad bands: 100 MeV–1 GeV, 1–10 GeV, 10–100 GeV, and 100–300 GeV) and mapped onto Galactic coordinates. For each band, intensity maps are generated and then projected onto latitude and longitude profiles, focusing especially on intermediate latitudes (|b|≈10°–20°) where the line‑of‑sight integration is dominated by local interstellar material.

To interpret the observations, the authors employ the GALPROP code, a widely used numerical model of cosmic‑ray (CR) propagation and associated gamma‑ray production. GALPROP takes as inputs the locally measured CR spectra (protons, helium, electrons, etc.), three‑dimensional distributions of atomic (HI) and molecular (H₂) gas, and the interstellar radiation field (IR, optical, CMB). It then computes the gamma‑ray emissivity from three processes: neutral‑pion decay (π⁰ → 2γ) following CR‑gas collisions, bremsstrahlung from CR electrons in the gas, and inverse‑Compton scattering of CR electrons on the radiation field. The model predictions are convolved with the LAT instrument response and compared directly to the observed maps.

The comparison yields several key findings. At intermediate latitudes and in the 100 MeV–10 GeV range, the LAT data agree with the GALPROP predictions to within a few percent. This close match indicates that the locally measured CR spectrum, when propagated a distance of roughly 1 kpc from the Sun, reproduces the observed gamma‑ray intensity, confirming that the basic ingredients of the model (CR source spectra, diffusion coefficient, gas maps) are adequate for the solar neighbourhood.

In the higher‑energy bands (>10 GeV) the observed intensity exceeds the model by about 10–20 %. The authors discuss three plausible contributors to this excess: (1) a spatial variation of the diffusion coefficient (δ) that is not captured by the globally uniform value used in the baseline model; (2) uncertainties in the gas column density, especially in the inner Galaxy where CO‑to‑H₂ conversion factors are poorly constrained; and (3) an additional, unresolved population of gamma‑ray sources (e.g., faint supernova‑remnant shells, pulsar wind nebulae, or dark matter subhalos) that would add a hard component to the diffuse emission.

Latitude profiles reveal a steep rise in intensity toward the Galactic plane at low energies, reflecting the dominance of π⁰ decay in dense gas. At higher energies the profiles flatten, consistent with a larger relative contribution from inverse‑Compton scattering, which depends on the broadly distributed interstellar radiation field rather than on gas density. Longitude profiles show enhancements coincident with known spiral‑arm tangents, indicating that the gamma‑ray emissivity traces the large‑scale distribution of interstellar material.

The paper’s broader scientific implications are twofold. First, the excellent agreement at low and intermediate energies validates the use of locally measured CR spectra as a baseline for Galactic diffuse emission models, providing confidence for background estimates in searches for exotic signals (e.g., dark matter annihilation). Second, the modest high‑energy discrepancy highlights the need for more sophisticated propagation scenarios, such as spatially varying diffusion, anisotropic transport, or localized re‑acceleration, and underscores the importance of improving ancillary data sets (high‑resolution HI, CO, and dust maps, as well as refined interstellar radiation field models).

Looking ahead, the authors outline a roadmap for future work. Continued LAT observations will increase photon statistics, reducing statistical uncertainties and allowing finer spatial binning. Integration of newer gas surveys (e.g., HI4PI, the latest CO surveys) and Planck dust‑derived column densities will sharpen the model inputs. Multi‑wavelength synergy—combining radio synchrotron data, X‑ray observations of supernova remnants, and TeV gamma‑ray measurements from ground‑based arrays—will help disentangle the contributions of distinct source classes. Ultimately, a three‑dimensional reconstruction of the Galactic CR density and spectrum will become feasible, opening the door to precise studies of the Milky Way’s structure, star‑formation history, and the mechanisms that accelerate cosmic rays throughout the Galaxy.