Improved Cosmic-Ray Injection Models and the Galactic Center Gamma-Ray Excess

TL;DR

This study improves cosmic-ray injection models near the Galactic Center, revealing their impact on gamma-ray excess interpretations and highlighting the importance of outflows and diffuse emission modeling in understanding the GC gamma-ray data.

Contribution

It introduces a non-linear cosmic-ray injection model based on 3D H2 density, significantly affecting the interpretation of the GC gamma-ray excess.

Findings

GC excess intensity decreases by a factor of 2 with larger sideband normalization.

Adding an excess template suppresses diffuse emission in smaller regions.

Fermi-LAT data favor outflows of several hundred km/s from the GC.

Abstract

Fermi-LAT observations of the Galactic Center (GC) have revealed a spherically- symmetric excess of GeV gamma rays extending to at least 10 deg from the dynamical center of the Galaxy. A critical uncertainty in extracting the intensity, spectrum, and morphology of this excess concerns the accuracy of astrophysical diffuse gamma-ray emission models near the GC. Recently, it has been noted that many diffuse emission models utilize a cosmic-ray injection rate far below that predicted based on the observed star formation rate in the Central Molecular Zone. In this study, we add a cosmic-ray injection component which non-linearly traces the Galactic H2 density determined in three-dimensions, and find that the associated gamma-ray emission is degenerate with many properties of the GC gamma-ray excess. In models that utilize a large sideband (40x40 deg surrounding the GC) to normalize the best-fitting diffuse emission models, the intensity of the GC excess decreases by a factor of 2, and the morphology of the excess becomes less peaked and less spherically symmetric. In models which utilize a smaller region of interest (15x15 deg) the addition of an excess template instead suppresses the intensity of the best-fit astrophysical diffuse emission, and the GC excess is rather resilient to changes in the details of the astrophysical diffuse modeling. In both analyses, the addition of a GC excess template still provides a statistically significant improvement to the overall fit to the gamma-ray data. We also implement advective winds at the GC, and find that the Fermi-LAT data strongly prefer outflows of order several hundred km/s, whose role is to efficiently advect low-energy cosmic rays from the Galactic center. Finally, we perform numerous tests of our models, and conclude that they significantly improve our understanding of multi-wavelength non-thermal emission from the GC.