Stellar-like Galactic center excess challenges particle dark matter

Résumé Technique : Les défis de l'excès galactique central de type stellaire pour la matière noire particulaire

Énoncé du Problème
Le Centre Galactique (CG) est une cible primaire pour la détection indirecte de la matière noire (DM) via les rayons gamma, car il est prédit comme abritant la densité de DM la plus élevée. Cependant, la région est complexifiée par l'« Excès du Centre Galactique » (GCE), une composante inattendue de rayons gamma culminant à des énergies de l'ordre du GeV, découverte par le télescope Fermi Large Area Telescope (LAT). Bien qu'initialement interprété comme une signature de l'annihilation de la DM de relique thermique, des analyses récentes suggèrent que la morphologie du GCE est mieux décrite par une population de pulsars millisecondes (MSP) non résolus dans le renflement stellaire plutôt que par un halo de DM sphériquement symétrique. Un défi majeur pour contraindre la DM dans cette région est le « mauvais modelage » (mismodeling) des fonds d'émission diffuse galactique, ce qui peut biaiser les méthodes statistiques utilisées pour distinguer les sources ponctuelles (comme les MSP) des signaux diffus (comme la DM). Des études antérieures ont souvent écarté le CG des recherches de DM en raison de ces incertitudes astrophysiques.

Méthodologie
Cet article réévalue le CG comme cible pour les recherches de DM en employant un cadre hybride entièrement auto-cohérent qui combine deux techniques statistiques avancées : le ajustement de modèles adaptatifs (utilisant le code skyFACT) et les statistiques de comptage de photons (utilisant la méthode 1pPDF).

1. Approche Hybride :

   • skyFACT (Optimisation) : Les auteurs utilisent d'abord skyFACT pour optimiser le modèle de la partie interne de la Galaxie. Cette étape minimise le mauvais modelage de l'émission diffuse galactique en introduisant des paramètres de nuisance qui permettent une modulation spectrale et spatiale des modèles (templates). Crucialement, les auteurs ajustent les données avec un modèle mixte comprenant à la fois un signal de DM et une composante de renflement stellaire simultanément. Cela garantit que les modèles de fond diffus « optimisés » résultants sont sur mesure pour le profil de densité, la masse et le canal d'annihilation de DM spécifiques testés.
   • 1pPDF (Dérivation des Contraintes) : Les modèles optimisés de skyFACT sont ensuite injectés dans l'analyse 1pPDF. Cette méthode analyse la distribution de probabilité du nombre de photons dans chaque pixel pour démêler les contributions des fonds isotropes, de l'émission diffuse, des sources résolues et des sources ponctuelles non résolues (sous le seuil de détection). En incluant les modèles de renflement stellaire et de DM comme composantes lisses aux côtés d'une distribution de comptage de sources paramétrique ($dN/dS$), la méthode quantifie l'espace restant pour un signal de DM tout en tenant compte des MSP sous le seuil de détection.

2. Détails de la Modélisation :

   • Modèles de DM : L'étude teste des modèles génériques de WIMP avec des masses allant de 10 GeV à 1 TeV, en se concentrant sur les canaux d'annihilation hadroniques ($b\bar{b}$) et leptoniques ($\tau^+\tau^-$).
   • Profils de Densité : Quatre profils de densité de DM sont testés : un NFW standard ($\gamma=1$), un NFW contracté ($\gamma=1.26$), un profil d'Einasto et un profil de Burkert à cœur (cored).
   • Renflement Stellaire : La composante stellaire est modélisée à l'aide d'un modèle de renflement boxy de type « Coleman » combiné à un renflement nucléaire, mis à jour par rapport aux travaux précédents pour refléter la morphologie préférée pour un GCE purement stellaire.
   • Données : L'analyse utilise 639 semaines de données Fermi-LAT (événements P8R3 ULTRACLEANVETO). L'ajustement skyFACT couvre une région de $40^\circ \times 40^\circ$, tandis que la dérivation des contraintes 1pPDF se concentre sur une région de la partie interne de la Galaxie de $20^\circ \times 20^\circ$ (excluant l'intérieur de $\pm 1^\circ$) dans la fenêtre d'énergie 1,6–5,9 GeV, là où le GCE culmine.

Contributions Clés

• Ajustement Simultané : Contrairement aux analyses précédentes qui traitaient souvent le GCE soit comme de la DM pure, soit comme un composant stellaire pur, ce travail ajuste un modèle mixte où les composantes de DM et de renflement stellaire sont autorisées à contribuer simultanément. Cela traite directement la dégénérescence entre les deux hypothèses.
• Gestion des Sources Sous le Seuil : La méthodologie tient explicitement compte des sources ponctuelles non résolues (MSPs sous le seuil) via la distribution de comptage de sources 1pPDF, empêchant qu'elles ne soient attribuées à tort à un signal de DM ou qu'elles ne laissent des résidus structurés.
• Validation de la Robustesse : Les auteurs effectuent des simulations étendues pour valider le pipeline. Ils démontrent que l'ajustement skyFACT réduit le mauvais modelage résiduel au niveau du bruit de Poisson de $2\sigma$ et que les limites supérieures dérivées sont compatibles avec les tests d'hypothèse nulle (c'est-à-dire qu'aucun signal de DM spurious n'est récupéré à partir de données simulées purement stellaires).

Résultats

• Contraintes de DM : L'analyse produit des limites supérieures strictes au niveau de confiance (C.L.) de 95 % sur la section efficace d'annihilation de la DM ($\langle \sigma v \rangle$).
  • Pour le profil NFW contracté (NFW126), les limites sont les plus strictes, excluant la section efficace de la relique thermique pour des masses $\lesssim 300$ GeV (canal hadronique) et $\lesssim 80$ GeV (canal leptonique).
  • Pour les profils NFW standard et Einasto, les limites sont légèrement plus faibles mais restent compétitives.
  • Pour le profil Burkert (à cœur), les limites sont nettement plus faibles (plusieurs ordres de grandeur plus élevées) en raison de la densité de DM prédite plus faible dans la région centrale.
• Préférence de Signal : Dans les ajustements skyFACT, la composante du renflement stellaire est trouvée comme étant une source d'émission significativement plus élevée que toute composante potentielle de DM (supprimée par des ordres de grandeur). Très peu de modèles de DM montrent une évidence statistique ($>3\sigma$) pour être un ajout nécessaire à un GCE purement stellaire.
• Comparaison avec d'autres Messagers :
  • Pour le canal $b\bar{b}$ et des masses allant jusqu'à 500 GeV, les limites de rayons gamma du CG (en supposant NFW126) sont environ trois fois plus sensibles que les contraintes combinées des galaxies naines sphéroïdales.
  • Les résultats sont compétitifs avec les contraintes d'antiprotons d'AMS-02 pour les masses inférieures à 200 GeV, étant jusqu'à un ordre de grandeur plus contraignants à 10–100 GeV.

Signification et Revendications
L'article affirme démontrer que le CG reste une cible puissante pour les recherches indirectes de DM, à condition que l'analyse prenne rigoureusement en compte la dégénérescence entre la DM et une population de renflement stellaire, ainsi que le mauvais modelage des fonds diffus. En combinant l'ajustement de modèles adaptatifs avec les statistiques de comptage de photons, les auteurs montrent qu'il est possible d'extraire des contraintes de DM robustes même en présence d'un GCE dominant d'origine stellaire.

Les auteurs notent modestement que leurs contraintes sont dérivées d'une fenêtre d'énergie étroite (1,6–5,9 GeV) pour garantir des résultats conservateurs, ce qui limite la sensibilité aux particules de DM très lourdes ou à des caractéristiques spectrales spécifiques (comme la raie Higgsino) qui pourraient être mieux sondées à des énergies plus élevées. Ils concluent que bien que leurs résultats remettent en question l'existence de la DM de relique thermique dans le CG pour des masses inférieures à ~300 GeV (selon le profil du halo), la méthodologie est suffisamment générale pour être adaptée à des profils de DM non sphériques ou à des formes spectrales différentes dans des travaux futurs. L'étude renforce l'interprétation selon laquelle le GCE est probablement d'origine stellaire, laissant peu de place pour une contribution significative de la DM dans l'espace des paramètres testé.