Constraining the slow-diffusion zone size and electron injection spectral index for the Geminga pulsar halo

1. Problématique

L'étude des halos de pulsars en rayons gamma, en particulier celui de Geminga, est cruciale pour comprendre la propagation des rayons cosmiques dans le milieu interstellaire (MIS) et résoudre l'énigme de l'excès de positrons observé par les expériences spatiales (comme AMS-02).

• Contexte : Les mesures précédentes du halo de Geminga par l'expérience HAWC (High-Altitude Water Cherenkov) ont permis d'estimer un coefficient de diffusion très faible (deux ordres de grandeur inférieur à la moyenne galactique), suggérant une « zone de diffusion lente » autour du pulsar.
• Limites des travaux antérieurs : Les études précédentes se concentraient principalement sur la mesure du coefficient de diffusion, sans contraindre de manière rigoureuse la taille de cette zone de diffusion lente ($r^*$) ni la forme spectrale d'injection des électrons ($p$) provenant de la nébuleuse de vent de pulsar (PWN). De plus, la relation entre ces paramètres et l'excès de positrons restait incertaine.
• Objectif : Utiliser les nouvelles données de morphologie spatiale et de spectre en énergie de HAWC pour contraindre simultanément la taille de la zone de diffusion lente et l'indice spectral d'injection des électrons, afin de tester la viabilité de Geminga comme source principale de l'excès de positrons.

2. Méthodologie

L'auteur utilise un modèle de diffusion à deux zones pour décrire la propagation des électrons et positrons :

1. Zone interne (lente) : Une sphère de rayon $r^*$ autour du pulsar où le coefficient de diffusion est $D_1(E) = D_{100}(E/100 \text{ TeV})^\delta$.
2. Zone externe (rapide) : Au-delà de $r^*$, le coefficient de diffusion reprend la valeur du MIS standard ($D_{ISM}$).

Équations et paramètres :

• L'équation de transport (diffusion-pertes) est résolue numériquement par la méthode des volumes finis.
• Le spectre d'injection est modélisé par une loi de puissance avec une coupure super-exponentielle : $q(E) \propto E^{-p} \exp[-(E/E_c)^2]$.
• Les données utilisées sont la morphologie bidimensionnelle (profil de brillance de surface, SBP) et le spectre en énergie (1–100 TeV) fournis par HAWC.
• Stratégie d'ajustement (Fitting) :
  • La taille de la zone lente ($r^*$) est contrainte par l'ajustement du profil de brillance de surface (SBP) via un test du $\chi^2$.
  • L'indice spectral d'injection ($p$) est contraint par l'ajustement du spectre en rayons gamma, en utilisant un test de rapport de vraisemblance pour déterminer les limites de confiance (2$\sigma$ et 3$\sigma$).
  • Différents scénarios pour l'indice de dépendance énergétique de la diffusion ($\delta = 0, 1/3, 1/2$) sont explorés.

3. Résultats Clés

A. Contrainte sur la taille de la zone de diffusion lente ($r^*$)

• Limite inférieure : L'ajustement du profil de brillance de surface (SBP) montre que les modèles avec $r^* \lesssim 30$ pc sont exclus avec un niveau de confiance supérieur à 95 %. Une zone trop petite entraîne une chute trop rapide du flux aux grandes distances angulaires, incompatible avec les données HAWC.
• Limite supérieure : En combinant l'ajustement du spectre gamma et la contrainte physique que l'indice d'injection $p$ ne doit pas être inférieur à 1 (basé sur les simulations de chocs et les observations de PWN), l'auteur déduit que $r^*$ ne peut dépasser 70 pc. Au-delà de cette valeur, l'indice $p$ requis pour ajuster le spectre tomberait en dessous de 1, ce qui est physiquement improbable.
• Conclusion sur $r^*$ : La taille de la zone de diffusion lente autour de Geminga est contrainte à l'intervalle 30 – 70 pc.

B. Contrainte sur l'indice spectral d'injection ($p$)

• En fixant la taille de la zone lente à sa limite inférieure ($r^* = 30$ pc) et en supposant une diffusion indépendante de l'énergie ($\delta = 0$), l'auteur obtient la valeur maximale permise pour l'indice spectral $p$.
• La limite supérieure à 3$\sigma$ est $p \le 2.17$.
• Pour des valeurs plus grandes de $\delta$ (ex: $\delta = 1/2$), la contrainte sur $p$ devient plus stricte (limite à 1.65).
• L'énergie de coupure du spectre d'injection est déterminée à environ 100 TeV.

C. Implications pour l'excès de positrons

• L'auteur projette le flux de positrons de Geminga vers la Terre en utilisant les paramètres contraints ($r^*=30$ pc, $p=2.17$, $\delta=0$).
• Résultat : Dans ce scénario « optimiste », le spectre de positrons prédit correspond remarquablement bien aux mesures d'AMS-02 dans la gamme 50 – 500 GeV, tant en forme spectrale qu'en flux absolu.
• Cela suggère que Geminga pourrait être une source majeure (voire dominante) de l'excès de positrons observé, à condition que la diffusion soit lente et que le spectre d'injection soit suffisamment dur (mais compatible avec les limites de 3$\sigma$).

4. Contributions et Significativité

• Précision accrue : Cette étude dépasse la simple mesure du coefficient de diffusion pour fournir des contraintes rigoureuses sur la géométrie de la zone de diffusion lente et la physique d'injection des particules, grâce aux données améliorées de HAWC.
• Résolution de dégénérescence : L'article démontre la corrélation inverse entre la taille de la zone lente ($r^*$) et l'indice spectral ($p$). Une zone plus petite permet un spectre d'injection plus dur, ce qui est crucial pour expliquer l'excès de positrons sans violer les contraintes des rayons gamma.
• Validation du modèle à deux zones : Les résultats soutiennent fortement le modèle de diffusion à deux zones, excluant les scénarios avec une zone de diffusion très petite (< 30 pc) qui étaient parfois envisagés.
• Impact sur la physique des rayons cosmiques : En établissant que Geminga peut expliquer l'excès de positrons dans le cadre des contraintes observationnelles actuelles, l'article renforce l'hypothèse selon laquelle les pulsars proches et d'âge moyen sont les sources de cet excès, plutôt que de nouvelles physiques (comme l'annihilation de matière noire).
• Perspectives futures : L'auteur note que la dégénérescence entre $\delta$, $r^*$ et $p$ pourrait être levée par les futures mesures de morphologie dépendante de l'énergie fournies par l'observatoire LHAASO.

En résumé, ce travail établit un cadre cohérent reliant la morphologie du halo de Geminga, son spectre gamma et le flux de positrons sur Terre, validant le rôle de Geminga comme source principale de l'excès de positrons sous des conditions physiques spécifiques mais plausibles.