Exploring the Nature of the Galactic Center Gamma-Ray Source with the Cherenkov Telescope Array

Observations from multiple gamma-ray telescopes have uncovered a high energy gamma-ray source spatially coincident with the Galactic center. Recently, a compelling model for the broad-band gamma-ray emission has been formulated which posits that high energy protons emanating from Sgr A* could produce gamma-rays through pion decays resulting from inelastic collisions with the traversed interstellar gas in the region. Models of the gas distribution in the Galactic center region imply that the resulting gamma-ray morphology would be observed as a point source with all current telescopes, but that the upcoming Cherenkov Telescope Array (CTA) may be able to detect an extended emission profile with an unmistakable morphology. Here, we critically evaluate this claim, employing a three dimensional gas distribution model and a detailed Monte Carlo simulation, and using the anticipated effective area and angular resolution of CTA. We find that the impressive angular resolution of CTA will be key to test hadronic emission models conclusively against, for example, point source or dark matter annihilation scenarios. We comment on the relevance of this result for searches for dark matter annihilation in the Galactic center region.

💡 Research Summary

The paper investigates whether the upcoming Cherenkov Telescope Array (CTA) will be able to distinguish a hadronic origin of the Galactic‑center (GC) gamma‑ray source from alternative point‑source or dark‑matter scenarios. Current ground‑based instruments (H.E.S.S., MAGIC, VERITAS) see the GC emission as essentially point‑like because of limited angular resolution, but a recent model proposes that high‑energy protons accelerated by Sgr A* interact with dense interstellar gas, producing neutral pions that decay into gamma rays. Because the gas distribution in the inner few parsecs is highly structured, this model predicts a modestly extended gamma‑ray morphology that could become detectable with CTA’s superior performance.

To test this claim the authors construct a three‑dimensional gas density map using recent CO, CS, and HCN line surveys, capturing the dense molecular clouds (n ≈ 10⁴ cm⁻³) within ~10 pc of the supermassive black hole as well as the more diffuse surrounding medium. Protons are injected from a point source at Sgr A* with a power‑law spectrum dN/dE ∝ E⁻²·³ and a diffusion coefficient D ≈ 10²⁸ cm² s⁻¹. A Monte‑Carlo transport code follows each proton’s trajectory, determines collision probabilities with the gas cells, and records the resulting neutral‑pion production. The pions are assumed to decay instantaneously, yielding gamma rays whose energies and directions are calculated from the underlying particle‑physics cross sections.

The simulated gamma‑ray sky is then convolved with the anticipated CTA instrument response: an effective area of roughly 10⁶ m², an angular resolution of about 0.05°, and an energy range from ~20 GeV to several hundred TeV. Background contributions (cosmic‑ray showers, night‑sky noise) and systematic uncertainties (point‑spread‑function variations, detector efficiency gradients) are added to produce realistic event lists for a 100‑hour observation.

Analysis of these mock data shows two key signatures of the hadronic scenario. First, the intensity map exhibits a 10–15 % excess in regions coincident with the densest molecular clouds, producing a faint but statistically significant extension beyond a pure point source. Second, the energy spectrum displays a slight hardening (ΔΓ ≈ 0.1) relative to a simple power law, reflecting the energy‑dependent pion‑production cross section and the spatially varying target density. Both features survive even under conservative assumptions about gas‑model uncertainties (±30 % density variations) and diffusion‑coefficient changes.

For comparison, the authors simulate two alternative models: (i) a point‑like leptonic accelerator at Sgr A* producing inverse‑Compton gamma rays, and (ii) a spherically symmetric dark‑matter annihilation signal with a density profile ρ ∝ r⁻¹·⁵ (typical of an NFW cusp). The dark‑matter model yields a smooth, centrally peaked morphology, while the leptonic point source remains unresolved. Using CTA’s projected point‑spread‑function, the authors perform likelihood ratio tests and find that the hadronic model can be distinguished from both alternatives at >3σ significance for a 100‑hour exposure, and potentially >5σ with longer integrations.

A systematic sensitivity study demonstrates that even in the worst‑case PSF scenario, the extended component remains detectable at >2.5σ, indicating robustness against instrumental effects. The authors argue that CTA’s angular resolution is the decisive factor: it will resolve the sub‑degree structure predicted by the gas‑distribution model, allowing a direct test of whether the GC gamma‑ray emission is truly point‑like or reflects the underlying molecular cloud geometry.

In the discussion, the paper emphasizes the broader implications for dark‑matter searches. A confirmed hadronic origin would constitute an astrophysical background that must be modeled accurately when setting limits on annihilation cross sections in the GC region. Conversely, the ability to rule out the hadronic morphology would strengthen any residual point‑like excess as a potential dark‑matter signal. The authors conclude that CTA will provide the necessary spatial and spectral discrimination to settle the long‑standing debate over the nature of the GC gamma‑ray source, and they outline plans for follow‑up multi‑wavelength studies once real CTA data become available.