New insights on low-mass dark matter subhalo tidal tracks via numerical simulations

Résumé technique : Nouvelles perspectives sur les trajectoires de marée des sous-halos de matière noire de faible masse via des simulations numériques

Énoncé du problème
L'article traite de l'évolution des sous-halos de matière noire (DM) de faible masse ($\lesssim 10^6 M_\odot$) lorsqu'ils orbitent au sein d'une galaxie hôte, en se concentrant spécifiquement sur l'environnement de la Voie Lactée (MW). Bien que des études antérieures aient établi que les sous-halos perdent de la masse à cause de l'érosion par effet de marée (tidal stripping), il existe un débat persistant concernant la résilience des sous-halos possédant des noyaux denses (spécifiquement les « cusps prompts » avec des pentes $\gamma \geq 1,5$) et l'évolution précise de leurs paramètres structurels. Une lacune critique dans la littérature est l'absence d'une analyse détaillée de la manière dont ces paramètres évoluent aux péricentres orbitaux (où les forces de marée sont les plus fortes) par rapport aux apocentres, particulièrement lorsque le halo hôte inclut des composantes baryoniques réalistes (disque et bulbe) et un potentiel gravitationnel en évolution temporelle. De plus, la relation entre la perte de masse et la « trajectoire de marée » (tidal track) résultante des paramètres structurels (vitesse circulaire maximale $V_{max}$, son rayon $r_{max}$, et la concentration) nécessite une résolution plus élevée et une exploration de l'espace des paramètres plus large que ce qui était disponible précédemment.

Méthodologie
Les auteurs utilisent le code DASH (un intégrateur N-body optimisé pour l'érosion par marée), utilisant une version améliorée introduite par Aguirre-Santaella et al. (2023). Les simulations suivent un sous-halo de DM unique orbitant au sein d'un potentiel analytique, en évolution temporelle, représentant la Voie Lactée, qui comprend un halo de DM, un disque stellaire et un bulbe.

• Propriétés des sous-halos : Les sous-halos sont totalement sombres (sans baryons) avec une masse initiale de $10^6 M_\odot$. Ils sont initialisés avec un profil de densité NFW généralisé, faisant varier la pente interne ($\gamma$) pour inclure des profils NFW standards ($\gamma=1$) et des cusps prompts ($\gamma=1,5$).
• Espace des paramètres : L'étude fait varier la concentration initiale du sous-halo ($c$), le redshift d'accrétion ($z_{acc}$), l'énergie orbitale ($x_c$), la circularité ($\eta$) et l'angle d'inclinaison ($\theta$).
• Résolution : Les simulations utilisent jusqu'à $N=2^{25}$ particules pour garantir des résultats robustes, même lorsque les sous-halos perdent une masse significative (jusqu'à des fractions de masse liées $f_b < 10^{-3,5}$).
• Analyse : Les auteurs extraient $V_{max}$ et $r_{max}$ à chaque instant (snapshot), en distinguant les apocentres et les péricentres orbitaux. Ils dérivent une « concentration de vitesse » ($c_V$) basée sur ces quantités, ce qui évite de dépendre d'un rayon de viriel qui peut ne plus exister après l'érosion.

Contributions clés

1. Analyse de péricentre : Contrairement aux travaux précédents qui analysaient principalement les trajectoires de marée aux apocentres orbitaux (où les sous-halos sont plus proches de l'équilibre), cette étude dérive et compare explicitement les trajectoires de marée aux péricentres, là où le sous-halo subit un stress de marée maximal.
2. Évolution des cusps prompts : L'article étudie l'évolution par marée des sous-halos avec des cusps prompts ($\gamma=1,5$), une population prédite comme étant hautement résiliente mais souvent sous-résolue dans les simulations cosmologiques standards.
3. Potentiel baryonique en évolution temporelle : Les simulations intègrent un potentiel hôte qui évolue au fil du temps et inclut des composantes baryoniques, offrant une représentation plus réaliste de l'environnement de la Voie Lactée par rapport aux potentiels statiques ou uniquement de DM utilisés dans la littérature antérieure.
4. Trajectoires de marée à haute résolution : En atteignant une résolution de particules sans précédent, les auteurs cartographient l'évolution structurelle des sous-halos jusqu'à des fractions de masse liées extrêmement basses, fournissant de nouvelles fonctions d'ajustement pour les trajectoires de marée.

Résultats

• Évolution structurelle : $V_{max}$ et $r_{max}$ diminuent avec le temps, mais $r_{max}$ rétrécit plus significativement que $V_{max}$. Ce rétrécissement différentiel conduit à une augmentation continue de la concentration du sous-halo ($c_V$).
• Augmentation de la concentration : La concentration de vitesse augmente d'environ deux ordres de grandeur de l'incursion jusqu'à aujourd'hui. C'est environ un ordre de grandeur de plus que l'augmentation observée pour les halos de champ. L'augmentation est plus prononcée pour les trajectoires de marée de péricentre par rapport aux apocentres.
• Trajectoires de marée ($V_{max}$ vs $f_b$ et $r_{max}$ vs $f_b$) :
  • Pour les profils NFW standards, les auteurs confirment l'existence d'une trajectoire de marée où l'évolution structurelle dépend principalement de la fraction de masse liée ($f_b$).
  • Péricentres vs Apocentres : Les trajectoires de péricentre présentent une courbure plus forte pour des pertes de masse modérées et n'atteignent les mêmes valeurs que les trajectes d'apocentre qu'après un épuisement important ($\log_{10} f_b < -3$).
  • Cusps prompts : Les sous-halos avec des cusps prompts sont plus résilients ; à $\log_{10} f_b = -2,5$, $V_{max}$ chute à environ $40\,\%$ de sa valeur initiale pour les cusps prompts, contre environ $20\,\%$ pour les profils NFW.
• Facteurs de dispersion : La dispersion dans les trajectoires de marée est pilotée par le rapport péricentre-apocentre, la distance précise du péricentre, le redshift d'accrétion et la circularité orbitale. Le redshift d'accrétion est identifié comme le principal moteur de la dispersion dans l'évolution de la concentration de vitesse.
• Limites de résolution : Les auteurs notent que pour $\log_{10}(r_{max}/r_{max,i}) \lesssim -1,5$, les résultats peuvent être affectés par une relaxation corps-deux artificielle, limitant la fiabilité de la trajectoire de marée dans les régimes d'érosion les plus extrêmes.

Signification
L'article affirme que ces conclusions sont cruciales pour interpréter correctement les données des populations de satellites galactiques, les recherches de sous-halos via les lentilles gravitationnelles ou les courants stellaires, et les recherches indirectes de matière noire (par exemple, les rayons gamma issus de l'annihilation).

• Détection indirecte : La découverte que les concentrations des sous-halos augmentent de manière significative (d'environ 2 ordres de grandeur) implique que les signaux d'annihilation de ces sous-halos pourraient être substantiellement plus élevés que précédemment estimé, car la force du signal dépend du carré de la densité.
• Taux de survie : L'identification des cusps prompts comme étant hautement résilients suggère qu'une population plus importante de sous-halos sombres de faible masse pourrait survivre dans les régions internes de la Voie Lactée, plus que ce que l'on pensait auparavant, particulièrement ceux sur des orbites excentriques qui passent près du centre galactique.
• Précision de la modélisation : L'étude souligne que négliger les composantes baryoniques et l'évolution temporelle du potentiel de l'hôte, ainsi que négliger la dynamique de péricentre, conduit à une compréhension incomplète de l'évolution structurelle des sous-halos. Les trajectoires de marée dérivées fournissent des fonctions d'ajustement mises à jour pour les modèles semi-analytiques et les futures simulations.