The contribution to Galactic Centre {\gamma}-ray excess from cluster-born millisecond pulsars. Constraints from direct N-body simulations

En utilisant des simulations N-corps directes pour modéliser la dynamique des amas globulaires et des pulsars millisecondes qu'ils hébergent, cette étude démontre que la population combinée de ces pulsars, issus d'amas survivants et disruptés, suffit à expliquer l'excès de rayons gamma du centre galactique, privilégiant ainsi une origine astrophysique plutôt que l'annihilation de matière noire.

L'essentiel

🌌 Le Mystère du "Surplus" de Lumière au Cœur de la Galaxie

Imaginez que vous regardez le centre de notre galaxie, la Voie Lactée, avec des lunettes spéciales capables de voir les rayons gamma (une lumière très énergétique, invisible à l'œil nu). Les astronomes ont remarqué quelque chose d'étrange : il y a beaucoup plus de lumière que prévu. C'est comme si, au milieu d'une ville calme, vous voyiez soudainement des milliers de feux d'artifice s'allumer sans raison apparente. Ce phénomène s'appelle le "Surplus Gamma du Centre Galactique".

Pendant longtemps, deux explications principales se battaient pour savoir d'où venait cette lumière :

1. L'hypothèse "Sci-Fi" : Ce serait la preuve de la Matière Noire qui s'annihile (des particules mystérieuses qui s'entrechoquent et disparaissent en libérant de l'énergie).
2. L'hypothèse "Astrophysique" : Ce serait simplement une foule immense d'étoiles mortes mais très actives, appelées pulsars millisecondes, qui brillent toutes ensemble.

🎢 L'Enquête : Une Simulation de "Toboggan Galactique"

L'équipe de chercheurs (Kalambay et ses collègues) a décidé de trancher en faveur de l'hypothèse des étoiles. Pour cela, ils ont créé un simulateur de réalité virtuelle géant (des simulations numériques) pour voir comment les étoiles se déplacent dans la galaxie.

Voici comment ils ont procédé, avec une image simple :

• Les Amas Globulaires (Les "Nurseries" d'étoiles) : Imaginez des boules de billard géantes composées de centaines de milliers d'étoiles, qui orbitent autour du centre de la galaxie comme des satellites.
• Le Centre Galactique (Le "Toboggan") : Au centre de la galaxie, il y a un trou noir supermassif et une forte gravité.
• Le Processus : Au fil des milliards d'années, ces boules d'étoiles s'approchent du centre. La gravité du centre agit comme un toboggan ou un vent fort qui arrache des étoiles de la boule.

🌟 Le Mécanisme : Comment les "Étoiles Zombies" arrivent au centre

Les chercheurs ont suivi deux types de voyages :

1. Les survivants : Ils ont suivi les boules d'étoiles qui existent encore aujourd'hui. En passant près du centre, elles perdent des étoiles.
2. Les disparus : Ils ont simulé des boules d'étoiles qui ont été totalement détruites il y a longtemps. Elles ont laissé derrière elles un "nuage de poussière" d'étoiles qui est tombé au centre.

Parmi ces étoiles arrachées, il y a des étoiles à neutrons. Ce sont des cadavres d'étoiles très denses. Certaines d'entre elles tournent très vite et émettent des faisceaux de lumière comme des phares de bateau : ce sont les pulsars.

🔍 La Découverte : Le Compteur ne ment pas

Les chercheurs ont fait le calcul :

• Ils ont compté combien de ces "phares" (pulsars) ont été déposés au centre de la galaxie par les boules d'étoiles (les amas).
• Ils ont estimé la luminosité de chacun.
• Le résultat : Si vous additionnez toute la lumière de ces milliers de petits phares, cela correspond exactement à la lumière mystérieuse que l'on observe !

C'est comme si vous cherchiez la source d'une odeur de café dans un immeuble. Au lieu de penser qu'il y a une machine à café magique cachée dans les murs (la matière noire), vous réalisez qu'il y a en fait 500 personnes qui ont chacune fait un café dans leur appartement. L'odeur totale est la même, mais l'explication est beaucoup plus simple et logique.

📊 Ce que disent les résultats

• La forme du surplus : La lumière observée a une forme spécifique (elle est plus brillante au centre et s'estompe doucement). La simulation montre que les étoiles arrachées par les amas détruits créent exactement cette forme.
• Le "Facteur 2" : Pour que le calcul colle parfaitement, les chercheurs pensent qu'il y a peut-être eu un peu plus d'amas d'étoiles détruits dans le passé que ce que l'on pensait (environ le double). Mais même avec cette petite hypothèse, le modèle tient la route.
• La conclusion : Il est très probable que le surplus de lumière vienne de ces pulsars (des étoiles mortes) et non de la matière noire.

🚀 Pourquoi c'est important ?

Cette étude nous dit que nous n'avons peut-être pas besoin de réinventer la physique pour expliquer ce phénomène. La nature est très efficace : elle utilise simplement les "débris" d'anciennes colonies d'étoiles pour créer une lueur brillante au centre de notre galaxie.

C'est une victoire pour l'astrophysique classique : parfois, la solution la plus simple (des étoiles) est la bonne, même quand le mystère semble très complexe.

En résumé : Le centre de la galaxie brille fort non pas à cause d'une magie sombre, mais parce qu'il est rempli de milliers de petits "phares" stellaires tombés du ciel, arrachés à leurs anciennes familles d'étoiles il y a des milliards d'années.

Résumé technique

Voici un résumé technique détaillé de l'article de recherche intitulé « The contribution to Galactic Centre γ-ray excess from cluster-born millisecond pulsars : Constraints from direct N-body simulations », rédigé en français.

1. Problématique et Contexte

Le centre galactique (GalC) présente un excès de rayonnement gamma dans la gamme des GeV (GCE), observé par le télescope Fermi-LAT. Ce signal, étendu spatialement et centré sur Sgr A*, excède les prédictions des modèles astrophysiques standards (interactions des rayons cosmiques avec le milieu interstellaire).
L'origine de ce GCE fait l'objet d'un débat intense entre deux hypothèses principales :

1. L'annihilation de matière noire (DM) : Un signal provenant de particules massives faiblement interactives (WIMPs) concentrées dans le halo de matière noire.
2. Une population de pulsars millisecondes (MSP) : Un ensemble non résolu de pulsars dans le bulbe galactique, dont l'émission collective imite le spectre observé.

L'hypothèse des MSP suggère qu'une partie significative de ces objets provient de l'évolution dynamique et de la destruction des amas globulaires (GC), qui déposent leurs étoiles et leurs restes compacts (étoiles à neutrons) dans le centre galactique au fil du temps.

2. Méthodologie

Les auteurs ont réalisé une étude dynamique complète pour réévaluer la contribution des MSP issus des amas globulaires au GCE, en utilisant des simulations N-corps directes de haute résolution.

• Simulations N-corps : Utilisation du code parallèle $\phi$-GPU, intégrant un schéma d'intégration Hermite d'ordre 4 et une évolution stellaire mise à jour (pertes de masse, supernovae, kicks nataux).
• Potentiel Galactique : Les amas évoluent dans un potentiel galactique variable dans le temps, dérivé de la simulation IllustrisTNG-100, pour capturer les effets de marée réalistes.
• Échantillons simulés :
  1. Amas actuels (MW GCs) : 12 amas globulaires réels observés, dont les orbites ont été intégrées vers le passé sur 8 Gyr.
  2. Amas détruits (Destr. GCs) : Un ensemble théorique d'amas (50 amas de $6 \times 10^4 M_\odot$et 4 amas massifs de$1,8 \times 10^5 M_\odot$) qui ont été totalement disruptés au cours des 5 premiers Gyr, déposant leurs débris dans le centre galactique.
• Suivi des Étoiles à Neutrons (NS) : Le code suit individuellement la formation et la dynamique de toutes les étoiles à neutrons produites, y compris celles qui sont éjectées de leurs amas parents par les kicks de supernova ou les effets de marée.
• Conversion NS $\to$ MSP : Comme les simulations ne modélisent pas l'évolution binaire nécessaire au recyclage des NS en MSP, les auteurs utilisent un facteur d'efficacité empirique ($k$). Ce facteur est calibré sur le rapport observé entre le nombre de MSP et le nombre d'étoiles à neutrons liées aux amas réels.
• Calcul du flux : Les distributions de phase (position et vitesse) des NS déposées dans le rayon central de 1 kpc sont converties en cartes de flux gamma, en supposant des luminosités moyennes par pulsar représentatives ($\langle L \rangle \sim 1-8 \times 10^{33}$ erg s$^{-1}$).

3. Contributions Clés

• Approche dynamique complète : Contrairement aux études précédentes qui utilisaient des fractions de rétention fixes ou des modèles statistiques simplifiés, cette étude suit la trajectoire individuelle de chaque étoile à neutrons, tenant compte des kicks nataux, de la ségrégation de masse et de l'évolution orbitale complexe.
• Inclusion des amas détruits : L'étude quantifie explicitement la contribution des amas globulaires qui n'existent plus aujourd'hui mais dont les débris ont peuplé le bulbe galactique, un canal souvent sous-estimé.
• Calibration empirique robuste : Détermination d'un coefficient médian $k \approx 0,0114$ (rapport MSP/NS liés) basé sur des catalogues observés et des prédictions théoriques récentes (incluant les découvertes FAST), appliqué de manière cohérente aux populations actuelles et détruites.

4. Résultats Principaux

• Distribution des NS :
  • La majorité des NS restent liées à leurs amas parents, mais une fraction significative (environ 6% dans la région centrale de 1 kpc) est déposée dans le bulbe galactique.
  • Les NS provenant des amas détruits présentent une morphologie axisymétrique avec des surdensités marquées dans le plan galactique et, de manière moins prononcée, perpendiculairement à celui-ci.
• Flux Gamma et GCE :
  • Les MSP issus des amas actuels fournissent déjà un signal gamma substantiel.
  • L'ajout de la population d'amas détruits augmente l'amplitude et la concentration centrale du signal.
  • Reproduction du GCE : En adoptant une luminosité moyenne par pulsar de $\langle L \rangle \approx 8 \times 10^{33}$ erg s$^{-1}$ et en augmentant légèrement (facteur $\sim 2$) la contribution estimée des amas détruits, le modèle reproduit très bien l'amplitude et la concentration spatiale du GCE observé.
  • L'accord est quasi-idéal après soustraction de la contribution du trou noir supermassif (Sgr A*).
• Sensibilité aux paramètres : Pour des luminosités plus faibles ($1-4 \times 10^{33}$erg s$^{-1}$), le modèle sous-produit le flux observé, indiquant des incertitudes sur la population d'amas primordiaux, leur masse initiale ou leur taux de disruption.

5. Signification et Conclusion

• Faveur pour l'origine astrophysique : Les résultats soutiennent fortement l'hypothèse d'une origine par les MSP pour le GCE, plutôt que l'annihilation de matière noire. La contribution combinée des MSP issus des amas survivants et détruits reproduit naturellement les propriétés observées du signal (amplitude et morphologie) sans nécessiter de physique exotique.
• Morphologie prédictive : La présence de surdensités dans le plan galactique et perpendiculairement à celui-ci, prédite par les simulations, offre des signatures observables testables par de futures campagnes d'observation (radio et gamma).
• Perspectives futures : L'étude suggère que des recherches radio profondes (avec FAST, SKA) ciblant ces surdensités, ainsi que des mesures spectroscopiques améliorées, permettront de contraindre davantage la démographie des MSP et de valider définitivement ce scénario.

En résumé, cette étude démontre que l'évolution dynamique des amas globulaires, y compris ceux qui ont disparu, suffit à expliquer le surplus de rayons gamma du centre galactique via une population de pulsars millisecondes, offrant une alternative astrophysique robuste à la matière noire.