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

Résumé technique : Contraintes par les rayons gamma sur les désintégrations d'axions de type lourds issues des spectres de blazars de Fermi-LAT et H.E.S.S.

Énoncé du problème
Les rayons gamma de très haute énergie (VHE ; $E_\gamma \ge 100$ GeV) provenant de sources extragalactiques, telles que les blazars, sont atténués lors de leur propagation à travers l'espace intergalactique en raison de la production de paires ($\gamma\gamma \to e^+e^-$) avec le fond diffus extragalactique (EBL - Extragalactic Background Light). Bien que les modèles d'EBL standards, basés sur l'évolution des galaxies et les relevés profonds, décrivent avec succès cette atténuation, les anomalies potentielles de la transparence de l'Univers ont motivé l'exploration d'une physique au-delà du Modèle Standard (BSM). Plus précisément, cette étude examine un scénario impliquant des particules de type axion (ALP) lourdes avec des masses $m_a \sim 2,5\text{--}20$ eV. Contrairement aux ALPs plus légères où les oscillations photon-ALP pourraient augmenter la transparence, les ALPs lourdes dans cette gamme de masse sont supposées se désintégrer en deux photons ($a \to \gamma\gamma$) sur des échelles de temps cosmologiques. Si les ALPs constituent une fraction significative de la matière noire, leur désintégration injecterait des photons supplémentaires dans le fond diffus cosmique infrarouge/optique, augmentant ainsi l'intensité de l'EBL et augmentant l'épaisseur optique des rayons gamma. Cet article cherche à contraindre la masse ($m_a$) et le couplage photon-ALP ($g_{a\gamma}$) de ces ALPs lourdes en testant si les spectres de blazars observés sont cohérents avec un modèle d'EBL augmenté par les contributions de désintégration des ALPs.

Méthodologie
L'analyse utilise un ensemble de données conjoint comprenant les observations du Fermi Large Area Telescope (Fermi-LAT) et du High Energy Stereoscopic System (H.E.S.S.).

• Sélection des données : L'échantillon est composé de 11 blazars sélectionnés sur la base de leur détection aux énergies VHE et à des décalages vers le rouge (redshifts) où l'absorption par l'EBL est significative ($\tau > 2$). Le jeu de données comprend 14 années de données Fermi-LAT (2008–2022) et des observations H.E.S.S.
• Segmentation temporelle : Pour tenir compte de la variabilité des sources, les auteurs ont employé l'algorithme des blocs bayésiens (Bayesian blocks) pour segmenter les courbes de lumière. Une analyse spectrale a été effectuée sur ces segments pour s'assurer que l'évolution spectrale intrinsèque ne biaisait pas les résultats. Pour les sources variables, les données contemporaines Fermi-LAT et H.E.S.S. ont été appariées ; pour les segments antérieurs à Fermi, seules les données H.E.S.S. ont été utilisées.
• Modélisation spectrale : Une analyse conjointe de maximum de vraisemblance a été réalisée à l'aide du package Gammapy. Les spectres intrinsèques des sources ont été modélisés à l'aide de fonctions de loi de puissance ou de log-parabole. L'analyse a inclus un paramètre de nuisance pour tenir compte d'un décalage cohérent de l'échelle d'énergie dans les données H.E.S.S. ($\sim 10\text{--}15\text{ \%}$ d'incertitude).
• Formalisme de désintégration des ALPs : La contribution des ALPs en désintégration à l'EBL a été modélisée en utilisant le formalisme de [54–56]. L'intensité spécifique de l'EBL induit par les ALPs dépend de la masse de l'ALP, du couplage et de la densité ($\Omega_a$). L'épaisseur optique totale a été calculée comme la somme de la contribution standard de l'EBL (en utilisant le modèle de Domínguez et al. [9] comme base, avec des comparaisons avec d'autres modèles) et de l'atténuation supplémentaire induite par la désintégration des ALPs.
• Cadre statistique : Les auteurs ont construit un rapport de vraisemblance profilé, $\lambda(m_a, g_{a\gamma})$, sommant les statistiques C à travers toutes les sources et tous les segments temporels. Les limites d'exclusion ont été dérivées au niveau de confiance de 95 % (C.L.) sous l'hypothèse que les ALPs constituent l'intégralité de la densité de matière noire ($\Omega_a = \Omega_{DM}$).

Principales contributions

• Première analyse conjointe : Ce travail représente la première analyse des signatures de désintégration d'ALPs lourdes utilisant un cadre de vraisemblance unifié combinant les données de Fermi-LAT et de H.E.S.S., maximisant la robustesse statistique sur une large gamme d'énergies (MeV à TeV).
• Couverture de l'espace des paramètres : L'étude balaie la gamme de masse des ALPs $m_a \approx 2,5\text{--}20$ eV. Cette gamme est physiquement motivée car les photons de désintégration dans cette fenêtre de masse tombent dans le domaine optique et proche infrarouge, ce qui affecte fortement l'atténuation des rayons gamma TeV via la production de paires.
• Gestion des systématiques : L'analyse intègre explicitement les incertitudes d'échelle d'énergie de H.E.S.S. et teste la robustesse des résultats par rapport à cinq modèles d'EBL différents (Domínguez, Franceschini, Finke, Saldana-Lopez, et Porras-Bedmar & Meyer).

Résultats

• Limites d'exclusion : L'analyse combinée produit des limites d'exclusion à 95 % C.L. sur le couplage photon-ALP $g_{a\gamma}$. Les contraintes les plus fortes sont trouvées aux masses plus élevées où la durée de vie de désintégration est plus courte, entraînant une émissivité de photons plus élevée.
  • À $m_a \approx 15$ eV, la limite atteint $g_{a\gamma} \sim 7 \times 10^{-12}$ GeV$^{-1}$.
  • À $m_a \approx 20$ eV, la limite est d'environ $g_{a\gamma} \sim 5,4 \times 10^{-12}$ GeV$^{-1}$.
• Cohérence avec l'hypothèse nulle : Les données ne montrent aucune préférence statistique significative pour une composante induite par les ALPs par rapport au modèle d'EBL standard. La grande majorité des segments de sources individuelles et l'ensemble des données combinées sont cohérents avec l'hypothèse nulle (pas de désintégration d'ALP).
• Robustesse du modèle : Les contraintes dérivées sont robustes face aux incertitudes de modélisation de l'EBL, les variations entre les différents modèles d'EBL restant dans un facteur d'environ 1,5.
• Comparaison avec les revendications précédentes : Les résultats défavorisent l'espace des paramètres suggéré par des travaux précédents (par exemple, [55, 67]) qui proposaient que les désintégrations d'ALPs pourraient expliquer des excès dans le fond optique cosmique ou les spectres de l'EBL. Plus précisément, la région suggérée par [55] ($m_a \sim 8\text{--}20$ eV, $g_{a\gamma} \sim (3\text{--}6) \times 10^{-11}$ GeV$^{-1}$) est exclue sous l'hypothèse $\Omega_a = \Omega_{DM}$.

Signification et affirmations
L'article affirme fournir les contraintes les plus fortes à ce jour sur les ALPs lourdes dans la gamme $m_a \sim 15\text{--}20$ eV basées sur les observations de rayons gamma. En fermant une région de l'espace des paramètres auparavant non contrainte, ces résultats offrent un test complémentaire aux recherches spectroscopiques de raies de désintégration étroites et aux expériences de laboratoire. Les auteurs soulignent que leur méthode repose sur la contribution intégrée de la désintégration à l'EBL et la modification résultante de l'épaisseur optique des rayons gamma, ce qui la rend indépendante des hypothèses concernant les facteurs D spécifiques des sources ou les fonds de raies étroites. Bien que l'analyse ne trouve aucune preuve de désintégration d'ALPs lourdes, elle établit des limites supérieures strictes qui contraignent la viabilité de telles particules en tant que composante dominante de la matière noire capable d'altérer significativement la transparence des rayons gamma de l'Univers. Ce travail démontre l'efficacité des analyses de rayons gamma multi-instruments pour sonder la physique au-delà du Modèle Standard et établit une base pour de futures recherches plus sensibles avec les instruments de prochaine génération comme le Cherenkov Telescope Array (CTAO).