Gamma-Ray Constraints on Heavy Axion-Like-Particle Decays from Fermi-LAT and H.E.S.S. Blazar Spectra

Technical Summary: Gamma-Ray Constraints on Heavy Axion-Like-Particle Decays from Fermi-LAT and H.E.S.S. Blazar Spectra

Problem Statement
Very-high-energy (VHE; $E_\gamma \ge 100$ GeV) gamma rays from extragalactic sources, such as blazars, are attenuated during propagation through intergalactic space due to pair production ($\gamma\gamma \to e^+e^-$) with photons of the Extragalactic Background Light (EBL). While standard EBL models based on galaxy evolution and deep surveys successfully describe this attenuation, potential anomalies in the transparency of the Universe have motivated the exploration of Beyond-the-Standard-Model (BSM) physics. Specifically, this study investigates a scenario involving heavy Axion-Like Particles (ALPs) with masses $m_a \sim 2.5\text{--}20$ eV. Unlike lighter ALPs where photon-ALP oscillations might increase transparency, heavy ALPs in this mass range are hypothesized to decay into two photons ($a \to \gamma\gamma$) over cosmological timescales. If ALPs constitute a significant fraction of dark matter, their decay would inject additional photons into the cosmic infrared/optical background, thereby enhancing the EBL intensity and increasing the gamma-ray optical depth. This paper seeks to constrain the mass ($m_a$) and photon coupling ($g_{a\gamma}$) of such heavy ALPs by testing whether the observed blazar spectra are consistent with an EBL model augmented by ALP decay contributions.

Methodology
The analysis utilizes a joint dataset comprising observations from the Fermi Large Area Telescope (Fermi-LAT) and the High Energy Stereoscopic System (H.E.S.S.).

• Data Selection: The sample consists of 11 blazars selected based on their detection at VHE energies and redshifts where EBL absorption is significant ($\tau > 2$). The dataset includes 14 years of Fermi-LAT data (2008–2022) and H.E.S.S. observations.
• Temporal Segmentation: To account for source variability, the authors employed the Bayesian blocks algorithm to segment light curves. Spectral analysis was performed on these segments to ensure that intrinsic spectral evolution did not bias the results. For variable sources, contemporaneous Fermi-LAT and H.E.S.S. data were matched; for pre-Fermi segments, only H.E.S.S. data were used.
• Spectral Modeling: A joint maximum-likelihood analysis was performed using the Gammapy package. The intrinsic source spectra were modeled using power-law or log-parabola functions. The analysis included a nuisance parameter to account for a coherent energy-scale shift in H.E.S.S. data ($\sim 10\text{--}15\%$ uncertainty).
• ALP Decay Formalism: The contribution of decaying ALPs to the EBL was modeled using the formalism from [54–56]. The specific intensity of the ALP-induced EBL depends on the ALP mass, coupling, and density ($\Omega_a$). The total optical depth was calculated as the sum of the standard EBL contribution (using the Domínguez et al. [9] model as a baseline, with comparisons to other models) and the additional attenuation induced by ALP decay photons.
• Statistical Framework: The authors constructed a profile likelihood ratio, $\lambda(m_a, g_{a\gamma})$, summing the C-statistics across all sources and time segments. Exclusion limits were derived at the 95% confidence level (C.L.) under the assumption that ALPs constitute the entire dark matter density ($\Omega_a = \Omega_{DM}$).

Key Contributions

• First Joint Analysis: This work represents the first analysis of heavy ALP decay signatures using a unified likelihood framework combining Fermi-LAT and H.E.S.S. data, maximizing statistical robustness across a wide energy range (MeV to TeV).
• Parameter Space Coverage: The study scans the ALP mass range of $m_a \approx 2.5\text{--}20$ eV. This range is physically motivated because decay photons in this mass window fall in the optical to near-infrared band, which most strongly affects TeV gamma-ray attenuation via pair production.
• Systematic Handling: The analysis explicitly incorporates H.E.S.S. energy-scale uncertainties and tests the robustness of results against five different EBL models (Domínguez, Franceschini, Finke, Saldana-Lopez, and Porras-Bedmar & Meyer).

Results

• Exclusion Limits: The combined analysis yields 95% C.L. exclusion limits on the photon-ALP coupling $g_{a\gamma}$. The strongest constraints are found at higher masses where the decay lifetime is shorter, leading to higher photon emissivity.
  • At $m_a \approx 15$ eV, the limit reaches $g_{a\gamma} \sim 7 \times 10^{-12}$ GeV$^{-1}$.
  • At $m_a \approx 20$ eV, the limit is approximately $g_{a\gamma} \sim 5.4 \times 10^{-12}$ GeV$^{-1}$.
• Null Hypothesis Consistency: The data show no significant statistical preference for an ALP-induced component over the standard EBL model. The vast majority of individual source segments and the combined dataset are consistent with the null hypothesis (no ALP decay).
• Model Robustness: The derived constraints are robust against uncertainties in EBL modeling, with variations across different EBL templates remaining within a factor of roughly 1.5.
• Comparison with Previous Claims: The results disfavor the parameter space suggested in previous works (e.g., [55, 67]) that proposed ALP decays could explain excesses in the cosmic optical background or EBL spectra. Specifically, the region suggested by [55] ($m_a \sim 8\text{--}20$ eV, $g_{a\gamma} \sim (3\text{--}6) \times 10^{-11}$ GeV$^{-1}$) is excluded under the assumption $\Omega_a = \Omega_{DM}$.

Significance and Claims
The paper claims to provide the strongest constraints to date on heavy ALPs in the $m_a \sim 15\text{--}20$ eV range based on gamma-ray observations. By closing a previously unconstrained region of parameter space, these results offer a complementary test to spectroscopic searches for narrow decay lines and laboratory experiments. The authors emphasize that their method relies on the integrated decay contribution to the EBL and the resulting modification of gamma-ray optical depth, making it independent of assumptions regarding specific source D-factors or narrow-line backgrounds. While the analysis does not find evidence for heavy ALP decay, it establishes stringent upper limits that constrain the viability of such particles as a dominant dark matter component capable of significantly altering the gamma-ray transparency of the Universe. The work demonstrates the efficacy of multi-instrument gamma-ray analyses in probing BSM physics and sets a baseline for future, more sensitive searches with next-generation instruments like the Cherenkov Telescope Array (CTAO).