Detection of Gamma Rays From a Starburst Galaxy

Starburst galaxies exhibit in their central regions a highly increased rate of supernovae, the remnants of which are thought to accelerate energetic cosmic rays up to energies of ~ 10^15 eV. We report the detection of gamma rays – tracers of such cosmic rays – from the starburst galaxy NGC 253 using the H.E.S.S. array of imaging atmospheric Cherenkov telescopes. The gamma-ray flux above 220 GeV is F = (5.5 +/- 1.0stat +/- 2.8sys) x 10^-13 ph. s-1 cm-2, implying a cosmic-ray density about three orders of magnitude larger than that in the center of the Milky Way. The fraction of cosmic-ray energy channeled into gamma rays in this starburst environment is 5 times larger than that in our Galaxy.

💡 Research Summary

The paper reports the first detection of very‑high‑energy (VHE) gamma‑ray emission from the nearby starburst galaxy NGC 253 using the H.E.S.S. (High Energy Stereoscopic System) array of imaging atmospheric Cherenkov telescopes. Starburst galaxies are characterized by intense, centrally concentrated star formation and a correspondingly high supernova rate. Supernova remnants are widely believed to accelerate cosmic‑ray (CR) particles up to the “knee” of the CR spectrum (~10¹⁵ eV). Because CRs are charged, they are deflected by magnetic fields and cannot be traced directly back to their sources; however, when CRs interact with interstellar gas or radiation fields they produce neutral pions (π⁰) that decay into gamma rays. Consequently, gamma‑ray observations provide a unique probe of the CR population in distant astrophysical environments.

Observations and Data Analysis
The H.E.S.S. instrument, located in the Southern Hemisphere, consists of four 12‑meter telescopes that record the Cherenkov light from air showers initiated by gamma rays in the 100 GeV–10 TeV range. The NGC 253 field was observed for a total live time of roughly 30 hours under good atmospheric conditions. Standard H.E.S.S. analysis pipelines—both the “reflected‑region” background method and a model‑based likelihood approach—were applied to extract a statistically significant excess of events from the galaxy’s central region. Systematic uncertainties were carefully evaluated by varying the atmospheric model, telescope optical efficiency, and analysis cuts.

Results
A gamma‑ray flux above 220 GeV was measured as

F(>220 GeV) = (5.5 ± 1.0 (stat) ± 2.8 (sys)) × 10⁻¹³ ph s⁻¹ cm⁻².

The spectrum is compatible with a simple power‑law (photon index ≈ 2.2), typical for hadronic CR interactions. Using the known distance to NGC 253 (≈ 3.5 Mpc) and assuming isotropic emission, the corresponding gamma‑ray luminosity is ≈ 1 × 10³⁹ erg s⁻¹. By adopting standard hadronic interaction models, the authors infer a CR energy density in the starburst core of roughly 10³ eV cm⁻³, i.e., three orders of magnitude higher than the CR density measured near the Galactic Center (~10 eV cm⁻³). Moreover, the fraction of CR power that is converted into gamma‑ray emission (the “calorimetric efficiency”) is about five times larger than in the Milky Way, indicating that the dense gas and intense radiation fields in NGC 253’s starburst region force CRs to lose energy more rapidly via pion production.

Interpretation and Implications
These findings provide direct observational support for the long‑standing hypothesis that starburst nuclei act as efficient CR factories. The elevated CR density is consistent with the high supernova rate (≈ 0.1–0.3 yr⁻¹) inferred from infrared and radio observations, which is an order of magnitude larger than in a typical spiral galaxy. The enhanced calorimetric efficiency suggests that NGC 253 is approaching the “CR calorimeter” regime, where most of the injected CR energy is dissipated locally rather than escaping into the intergalactic medium. This has important consequences for galaxy evolution: CR pressure can drive outflows, heat the interstellar medium, and influence the star‑formation feedback loop.

Future Prospects
The authors emphasize that the current measurement is limited by statistical uncertainties and the relatively narrow energy band accessible to H.E.S.S. The upcoming Cherenkov Telescope Array (CTA), with an order‑of‑magnitude improvement in sensitivity and a broader energy coverage (≈ 20 GeV–100 TeV), will be able to resolve the spatial distribution of the gamma‑ray emission, test for spectral curvature, and possibly detect the transition from hadronic to leptonic processes. Complementary observations at other wavelengths—radio synchrotron (tracing CR electrons), infrared (tracing dust‑heated radiation fields), and X‑ray (tracing hot gas and supernova remnants)—will enable a multi‑messenger model of CR injection, propagation, and loss in starburst environments.

Conclusion
The detection of VHE gamma rays from NGC 253 confirms that the intense star‑forming core of a starburst galaxy hosts a CR population far more energetic than that of the Milky Way. The measured gamma‑ray flux translates into a CR energy density ~1000 times higher than in our Galactic Center and a calorimetric efficiency ~5 × greater. These results cement starburst nuclei as key contributors to the high‑energy particle budget of the Universe and open a new observational window on the interplay between star formation, supernova feedback, and cosmic‑ray physics.