Measurement of Substructure from the Kinematics of the GD-1 Stellar Stream

Résumé technique : Mesure de la sous-structure à partir de la cinématique du courant stellaire GD-1

Énoncé du problème
Dans le paradigme de la Matière Noire Froide (CDM), les galaxies devraient abriter une vaste population de sous-halos de matière noire, allant des échelles galactiques jusqu'aux échelles de la masse de la Terre. Cependant, les sous-halos situés en dessous du seuil de formation des galaxies ($\approx 10^8 M_\odot$) sont censés être dépourvus de matière baryonique, ce qui les rend invisibles à l'observation directe. Par conséquent, tester le paradigme CDM sur de petites échelles nécessite des traceurs gravitationnels indirects. Les courants stellaires, qui sont des structures cinématiquement froides formées par la disruption tidale de amas globulaires, sont sensibles aux perturbations causées par des sous-halos de faible masse. Bien que des études précédentes aient utilisé les lacunes morphologiques et les excroissances (spurs) dans les courants comme GD-1 pour inférer les interactions avec les sous-halos, ce travail répond à la nécessité de contraindre la dispersion de vitesse intrinsèque du courant. Une dispersion de vitesse mesurée significativement plus élevée que celle prédite par des potentiels lisses de la Voie Lactée indiquerait un chauffage piloté par la sous-structure de la matière noire.

Méthodologie
Les auteurs analysent le courant stellaire GD-1 en utilisant une combinaison de données observationnelles et de modélisation par théorie des perturbations :

1. Collecte de données et mesure de la dispersion de vitesse :

   • L'étude utilise 160 étoiles membres uniques identifiées à partir du catalogue de GD-1 piloté par les données (Starkman et al. 2025) avec une probabilité d'appartenance $>75\%$.
   • Les vitesses radiales ($v_r$) sont compilées à partir de quatre relevés spectroscopiques : DESI (première version des données), SDSS DR9, LAMOST DR8 et MMT Hectochelle.
   • Pour limiter la contamination par les étoiles binaires, les auteurs utilisent des observations répétées de 34 étoiles à travers différents instruments.
   • Un modèle bayésien hiérarchique est employé pour séparer les erreurs d'observation, le saut de vitesse induit par les binaires ($\sigma_{binary}$) et la dispersion de vitesse intrinsèque du courant ($\sigma_{\Delta v_r}$). Le courant est segmenté en trois segments le long de sa longueur ($\phi_1 \in [-60, -40], [-40, -20], [-20, 0]$ degrés) afin de mesurer la dispersion locale.

2. Modélisation dynamique :

   • Les auteurs utilisent le cadre streamsculptor, qui emploie la théorie des perturbations pour modéliser l'effet cumulatif de centaines d'impacts de sous-halos sur le courant.
   • Des modèles de courants non perturbés sont générés à l'aide d'une approche de pulvérisation de particules (particle-spray) pour des masses de progéniteurs de $2 \times 10^4 M_\odot$ et $10^5 M_\odot$ et des âges dynamiques de 3,5, 5 et 8 Gyr, au sein d'un potentiel statique de la Voie Lactée (MilkyWayPotential2022).
   • Les sous-halos sont modélisés comme des profils de Hernquist. L'étude explore deux scénarios pour lever les dégénérescences entre la densité numérique des sous-halos et leur structure interne :
     • Modèle I : Suppose un scénario de Matière Noire Tiède (WDM) avec une suppression dans la fonction de masse des sous-halos (SHMF) en dessous d'une masse de mode demi (half-mode mass, $M_{hm}$) et une relation masse-taille constante ($r_s/r_{s,cdm}$).
     • Modèle II : Suppose une SHMF standard de type CDM ($M_{hm}=0$) mais permet à la relation masse-taille de varier selon une loi de puissance brisée en fonction de la masse, testant ainsi une concentration dépendante de la masse.

3. Inférence :

   • Les auteurs utilisent l'Évaluation Bayésienne Approchée (ABC) pour comparer les simulations de réalisations de courants aux dispersions de vitesse observées.
   • Les simulations sont générées pour diverses populations de paramètres de sous-halos (normalisation, $M_{hm}$, concentration). Les échantillons sont acceptés si la distance euclidienne entre la dispersion de vitesse simulée et observée (ainsi que les contraintes de largeur du courant) tombe dans une tolérance $\epsilon$.

Contributions clés et résultats

• Dispersion de vitesse intrinsèque : L'étude mesure la dispersion de vitesse radiale intrinsèque de GD-1 à environ $2\text{--}5$ km s$^{-1}$ sur toute sa longueur. La région présentant la dispersion la plus élevée ($\phi_1 \in [-40, -20]$ degrés) montre une valeur de $4,8^{+1,2}_{-0,7}$ km s$^{-1}$. Cela représente un écart de $4\sigma$ par rapport aux modèles de courants non perturbés formés dans un potentiel lisse, qui prédisent des dispersions de seulement $\approx 0,7\text{--}1,1$ km s$^{-1}$.
• Contraintes sur la population de sous-halos : La dispersion observée peut être expliquée par une population de sous-halos de matière noire. La fraction de la masse de la Voie Lactée contenue dans les sous-halos est déduite à $f_{sub} = 0,05^{+0,08}_{-0,03}$ (confiance de 68 %), ce qui est cohérent avec les attentes de la CDM ($\approx 7,5 \%$).
• Compacité des sous-halos : Un résultat principal est la préférence pour des sous-halos plus compacts que les prédictions standard de la CDM.
  • Dans le Modèle I, si la fraction de masse des sous-halos est de $10\%$, le rayon d'échelle est contraint à $r_s/r_{s,cdm} = 0,3^{+0,3}_{-0,2}$ (confiance de 68 %).
  • Dans le Modèle II, il existe une préférence pour des sous-halos de plus en plus compacts à faible masse ($M < 10^7 M_\odot$). Pour les sous-halos ayant une $M_{200} = 10^6 M_\odot$, la concentration $c_{-2}$ est contrainte à $>190$ avec une confiance de 68 %, ce qui est nettement supérieur à la prédiction de la CDM d'environ $43$.
• Scénarios alternatifs : Les données peuvent être expliquées soit par de nombreux impacts avec des sous-halos de faible masse, soit par un impact unique avec un sous-halo très compact ($M \gtrsim 10^8 M_\odot$). Les auteurs notent qu'un impact unique avec un sous-halo compact est suffisant pour générer la dispersion observée, ce qui est cohérent avec les interprétations précédentes de la caractéristique d'excroissance (spur) de GD-1.

Signification et affirmations
L'article affirme fournir une nouvelle contrainte indépendante sur la sous-structure de la matière noire en utilisant la cinématique des courants stellaires, complétant les contraintes précédentes dérivées de la morphologie de la densité des courants. Les résultats suggèrent une possible déviation de la relation masse-taille standard de la CDM à faibles masses de sous-halos, là où les sous-halos sont censés être purement dominés par la matière noire. Plus précisément, la préférence pour des concentrations plus élevées (ou des rayons d'échelle plus petits) que les prédictions de la CDM s'aligne avec des théories alternatives de la matière noire, telles que la Matière Noire Auto-Interagissante (SIDM) ou la Matière Noire Atomique, qui prédisent des densités centrales accrues.

Les auteurs soulignent que leurs contraintes sont empiriques et agnostiques vis-à-vis de modèles spécifiques de physique des particules, servant de référence pour tester diverses théories de la matière noire. Ils reconnaissent que des dégénérescences existent entre le nombre d'impacts de sous-halos et leur concentration (la "dégénérescence taux-concentration"), qui pourraient être levées dans des travaux futurs en modélisant conjointement la densité et la cinématique du courant. L'étude conclut que, bien que la fraction de masse des sous-halos soit cohérente avec la CDM, la structure interne de ces sous-halos pourrait nécessiter une physique au-delà du paradigme standard de la CDM, particulièrement à faibles masses.