The Sensitivity of Substructure Lensing to SIDM Core-collapse Model Variation

Résumé technique : La sensibilité du lentillage de sous-structures aux variations du modèle d'effondrement du cœur de la matière noire auto-interagissante

Énoncé du problème
Le lentillage gravitationnel fort est devenu une sonde puissante pour la sous-structure de la matière noire, offrant un test spécifique pour les modèles de matière noire auto-interagissante (SIDM). Dans les scénarios SIDM avec des collisions élastiques, le transfert de chaleur conduit à la formation de cœurs isothermes, suivis d'une « catastrophe gravothermique » (effondrement du cœur) qui produit des cœurs centraux denses. Ces halos compacts sont prédits pour laisser des empreintes distinctes sur les observables de lentillage. Cependant, la prédiction précise de ces signatures est entravée par la difficulté de modéliser la chronologie de l'effondrement du cœur et l'état final des sous-halos. Les approches actuelles reposent souvent sur des hypothèses simplificatrices, telles que le traitement de l'effondrement du cœur comme une transition instantanée ou l'utilisation d'assignations probabilistes pour les états d'effondrement. Les auteurs identifient un décalage entre les pipelines d'inférence par lentillage existants et les modèles physiques sous-jacents, notant que de petites variations dans la modélisation des sous-halos en cours d'effondrement du cœur peuvent entraîner des différences significatives dans la puissance du lentillage. Sans quantifier ces incertitudes de modélisation, la capacité des observations de lentillage à contraindre la nature de la matière noire est limitée.

Méthodologie
Les auteurs quantifient ce défi de modélisation en testant diverses modifications des profils de densité des halos en effondrement du cœur SIDM dans le contexte du lentillage de sous-structures. Leur méthodologie comprend :

1. Simulation et calibration : Ils utilisent des données de simulation $N$-corps (générées avec le code Arepo et le module SIDM de la référence [20]) pour calibrer un traitement approximatif des effets de marée sur l'effondrement gravothermique du cœur. Ils dérivent un « facteur d'accélération de l'effondrement » pour tenir compte de la différence dans les échelles de temps d'effondrement entre les halos isolés et les sous-halos évoluant dans un potentiel hôte.
2. Réalisations de lentillage : En utilisant le package open-source pyHalo, ils génèrent 200 réalisations aléatoires de systèmes de lentilles à l'échelle galactique (à $z_{lentille}=0,5$) contenant des sous-halos de déflecteur principal et des halos de ligne de visée (LOS). Ils remplacent les profils de matière noire froide (CDM) standard par des profils SIDM basés sur une solution fluide gravothermique calibrée sur des simulations $N$-corps.
3. Variations de modèle : L'étude teste plusieurs choix de modélisation :
   • Continu vs Instantané : Comparaison des profils de densité évoluant de manière fluide avec des modèles binaires où les halos sont soit totalement non effondrés, soit totalement effondrés.
   • Seuils d'effondrement : Variation du moment où le profil de densité est « gelé » (arrêt) ou « boosté » vers l'état final, testant différents états finaux de l'effondrement.
   • Probabiliste vs Déterministe : Comparaison du suivi de l'évolution individuelle des halos avec l'assignation aléatoire d'états d'effondrement basée sur des tranches de masse.
   • Profils initiaux : Test d'halos non effondrés modélisés soit comme des TNFW (Truncated NFW) standard, soit comme ayant atteint la phase de cœur maximal.
4. Analyse statistique : Ils calculent le potentiel de lentillage gravitationnel multiplan effectif, en déterminant la divergence ($\kappa_{div}$) et le rotationnel ($\kappa_{curl}$) du champ de déflexion effectif. La statistique principale est la fonction de corrélation à deux points de ces champs, décomposée en moments multipolaires (monopôle et quadrupôle) à l'aide de bases orthonormées pour dissocier les effets de la ligne de visée de ceux du lentille principal.

Contributions clés

• Quantification de l'incertitude de modélisation : L'article fournit une quantification systématique de la manière dont des hypothèses simplificatrices spécifiques dans la modélisation de l'effondrement du cœur SIDM affectent les statistiques de lentillage.
• Calibration des marées : Il introduit une méthode probabiliste pour tenir compte de l'accélération/décélération par marée de l'effondrement du cœur dans les sous-halos, calibrée par rapport aux simulations $N$-corps.
• Comparaison binaire vs continue : Il compare rigoureusement l'approximation de « l'effondrement binaire » (courante dans la littérature) avec un modèle gravothermique continu et évolutif dans le temps, mettant en évidence les échecs de cette approximation.
• Sensibilité du profil : L'étude démontre que le choix du profil de densité initial (TNFW vs à cœur) et la méthode d'assignation des états d'effondrement modifient considérablement le signal de lentillage prédit.

Résultats

• Sensibilité à la densité vs au temps : La fonction de corrélation à deux points (spécifiquement les moments multipolaires) est très sensible à la densité centrale des halos effondrés mais relativement insensible à l'échelle de temps exacte de l'effondrement. Les variations du temps d'effondrement (par exemple, arrêt de l'évolution à 20 % du temps total) ne produisent pas de résultats significativement différents par rapport à l'évolution complète, à condition que le profil de densité final soit précis.
• Suppression par les modèles binaires : Les modèles d'effondrement du cœur binaire suppriment systématiquement la puissance à petite échelle dans les statistiques de lentillage par rapport au modèle d'évolution temporelle complète. Cette suppression est la plus sévère lorsque :
  • Les halos non effondrés sont modélisés avec des profils initiaux à cœur plutôt qu'avec des profils TNFW.
  • Les états d'effondrement sont assignés de manière aléatoire plutôt que basés sur l'évolution individuelle des halos.
• Dépendance à l'échelle : Les différences entre les modèles sont les plus discernables aux petites échelles de longueur (inférieures à 0,1 seconde d'arc). Le moment quadrupôle du champ de divergence ($\xi_{div,2}$) reste sensible aux variations de modèle à des échelles plus grandes que le monopôle ou les moments de rotationnel.
• Corrélation entre l'effondrement et la puissance de lentillage : L'étude constate que les halos ayant un effet de lentillage plus prononcé (masse/concentration plus élevée) tendent à s'effondrer plus rapidement. L'aléatorisation des états d'effondrement supprime l'amélioration SIDM car elle découple l'état d'effondrement des propriétés de l'halo qui pilotent le signal de lentillage.
• Validation par $N$-corps : Le modèle approximatif d'accélération par marée des auteurs produit des statistiques de lentillage cohérentes avec les simulations $N$-corps dans un écart-type, validant son utilisation comme approximation de premier ordre pour cette étude.

Signification et affirmations
L'article affirme que la fonction de corrélation à deux points est une sonde sensible des simplifications spécifiques apportées dans la modélisation de l'effondrement du cœur SIDM. Les auteurs affirment que :

1. Robustesse de l'état final : Des incertitudes systématiques modestes dans le timing de l'effondrement du cœur peuvent ne pas affecter significativement les prédictions de lentillage, à condition que le profil de densité final de l'halo effondré soit modélisé avec précision.
2. Limites des approximations binaires : Bien que les modèles d'effondrement binaire soient efficaces sur le plan computationnel, ils peuvent introduire des biais significatifs en supprimant la puissance à petite échelle, en particulier s'ils supposent des conditions initiales à cœur ou aléatorisent les états d'effondrement.
3. Nécessité d'une motivation physique : Les futures analyses de lentillage visant à tester le SIDM doivent faire des choix intentionnels et physiquement motivés concernant les profils de sous-structures et de halos de ligne de visée. Les auteurs recommandent que de tels travaux incluent des tests analogues de leurs choix de profils et comparent leurs profils d'angle de déflexion avec des modèles d'évolution complète physiquement motivés.
4. Dépendance observationnelle : La détectabilité de ces différences de modélisation dépend fortement de la sonde observationnelle et des échelles analysées. L'article ne prétend pas que les données actuelles peuvent exclure définitivement des modèles SIDM spécifiques, mais plutôt que l'interprétation de telles données repose fortement sur la robustesse des modèles de sous-structure sous-jacents.

Les auteurs concluent que l'affinement des modèles d'effondrement du cœur SIDM dans des environnements réalistes est essentiel, car les prédictions testables peuvent fortement dépendre de l'état final des halos effondrés. Ils suggèrent que leur cadre peut être utilisé pour intégrer ces modèles d'évolution complète dans les prédictions des anomalies de rapports de flux et d'autres observables de lentillage dans des travaux futurs.