JWST lensed quasar dark matter survey IV: Stringent warm dark matter constraints from the joint reconstruction of extended lensed arcs and quasar flux ratios

Résumé Technique : Enquête de Lentillage JWST IV – Contraintes sur la Matière Noire Tiède via une Reconstruction Conjointe

Énoncé du Problème
La nature de la matière noire (MN) demeure inconnue, le paradigme de la Matière Noire Froide (CDM) prédisant une hiérarchie de halos sans échelle de masse jusqu'à des valeurs très faibles. Les théories alternatives, telles que la Matière Noire Tiède (WDM), prédisent une suppression de la structure en dessous d'une échelle de diffusion libre spécifique, caractérisée par la masse de mode demi ($m_{hm}$). Le lentillage gravitationnel fort de quasars à quatre images offre une sonde unique de cette sous-structure en mesurant les perturbations dans les positions des images et les rapports de flux. Cependant, les analyses précédentes reposant uniquement sur les positions de sources ponctuelles et les rapports de flux ont été limitées par les dégénérescences entre le profil de masse du déflecteur principal (macromodèle) et la population de sous-structures. Les incertitudes liées au macromodèle (ex. : ellipticité, cisaillement externe, et structure angulaire) peuvent imiter ou masquer les signatures des sous-halos de matière noire, conduisant à des contraintes plus faibles sur les propriétés de la MN. De plus, le microlentillage stellaire a historiquement contaminé les mesures de rapport de flux pour les sources compactes, bien que cela ait été atténué par les récentes observations JWST de régions de poussière chaude étendues.

Méthodologie
Cet article présente un cadre de modélisation directe qui reconstruit simultanément les arcs lentillés étendus, les positions des images et les rapports de flux pour un échantillon de 28 quasars à quatre images. L'analyse utilise un cadre d'inférence bayésienne pour contraindre les hyper-paramètres décrivant la population de sous-structures de la MN.

• Intégration des Données : L'étude combine les rapports de flux MIRI de JWST (provenant de la région de poussière chaude du quasar, évitant ainsi le microlentillage stellaire) avec les données astrométriques et l'imagerie d'arcs lentillés étendus de JWST (NIRCam et MIRI) et HST.
• Modélisation Conjointe : Contrairement aux travaux précédents qui traitaient les arcs lentillés et les rapports de flux séparément, cette analyse emploie des poids d'échantillonnage par importance ($w_{img}$) pour propager les contraintes des arcs étendus sur les paramètres du macromodèle. Cela brise la covariance entre le macromodèle et la population de sous-structures, resserrant ainsi les contraintes sur les rapports de flux prédits par le modèle.
• Modélisation Physique :
  • Sous-structure : Le cadre génère des populations complètes de sous-halos, de halos sur la ligne de visée et d'amas globulaires en utilisant le logiciel pyHalo. Il incorpore un modèle d'évolution tidale (basé sur Du et al.) qui prédit les profils de densité des sous-halos en fonction des propriétés d'infall, tenant compte de la perte de masse et des changements de concentration.
  • Effets WDM : Le modèle modifie la fonction de masse des halos CDM et la relation concentration-masse pour tenir compte de la suppression par diffusion libre en WDM, paramétrée par la masse de mode demi $m_{hm}$.
  • Macromodèle : Le déflecteur principal est modélisé à l'aide d'un profil de loi de puissance elliptique (EPL) augmenté de perturbations multipolaires elliptiques (m=1, 3, 4) pour tenir compte de la structure angulaire (ex. : boxiness/diskiness) qui pourrait autrement passer pour de la sous-structure.
• Inférence : La fonction de vraisemblance est évaluée en générant des millions de réalisations (plus de 174 millions pour l'échantillon complet) de modèles de lentilles. L'analyse marginalise les paramètres de nuisance (macromodèle, lumière de la source, réalisations de sous-structures) pour inférer les distributions postérieures des hyper-paramètres, spécifiquement la masse de mode demi $m_{hm}$ et l'amplitude de la fonction de masse des sous-halos $\Sigma_{sub}$.

Contributions Clés

1. Briser les Dégénérescences : L'avancée méthodologique primaire est la modélisation simultanée des arcs lentillés étendus aux côtés des rapports de flux et des positions d'images. Cette approche réduit significativement l'incertitude du macromodèle, brisant ainsi la covariance entre l'abondance des sous-halos ($\Sigma_{sub}$) et l'échelle de diffusion libre ($m_{hm}$).
2. Amélioration de la Modélisation des Sous-halos : L'article implémente un modèle d'évolution tidale plus rigoureux physiquement qui prédit les propriétés des sous-halos au moment de l'infall, plutôt qu'au redshift de la lentille. Cela rend compte de la dépendance de l'arrachement tidal vis-à-vis de la concentration lors de l'infall, ce qui est particulièrement pertinent pour la WDM où les halos sont moins concentrés.
3. Sources de Perturbation Exhaustives : L'analyse inclut explicitement les amas globulaires et les perturbations multipolaires elliptiques comme sources de bruit de rapport de flux, garantissant que les contraintes sur la WDM ne soient pas biaisées par des structures à petite échelle non modélisées.
4. Taille d'Échantillon Importante : Ce travail analyse le plus grand échantillon de quasars à quatre images à ce jour pour l'inférence de la MN (28 systèmes), exploitant pleinement les capacités de l'enquête de quasars lentillés de JWST.

Résultats

• Contraintes sur la Matière Noire Tiède (WDM) :
  • En supposant un a priori log-uniforme sur l'amplitude de la fonction de masse des sous-halos, l'inclusion des arcs lentillés améliore la limite sur la masse de mode demi de 0,6 dex par rapport à l'utilisation seule des positions d'images et des rapports de flux.
  • Lors de l'application d'a priori informatifs sur $\Sigma_{sub}$ dérivés du modèle semi-analytique galacticus, l'analyse infère $m_{hm} < 10^{7,4} M_\odot$ (facteur de Bayes de 10:1).
  • Lors de l'application d'a priori basés sur des simulations N-corps (Symphony suite), la contrainte se resserre à $m_{hm} < 10^{7,2} M_\odot$.
  • Ces limites correspondent à des limites inférieures sur la masse d'une particule de relique thermique WDM de 6,5 keV et 7,4 keV, respectivement.
  • Les limites d'exclusion à 95 % sont $m_{hm} < 10^{6,8} M_\odot$ et $m_{hm} < 10^{6,7} M_\odot$, correspondant à des masses de particules $> 9,6$ keV et $> 10,2$ keV.
• Abondance de Sous-halos CDM :
  • En supposant la limite CDM ($m_{hm} < 10^5 M_\odot$), l'étude mesure la densité de masse de surface projetée des sous-halos à $1,7^{+2,6}_{-1,2} \times 10^7 M_\odot \text{kpc}^{-2}$ avec un intervalle de confiance de 95 %.
  • Cela correspond à une fraction de masse projetée de $3,2^{+5,1}_{-2,3}\%$.
  • Cette mesure est statistiquement cohérente avec les prédictions de galacticus mais montre une légère tension ($\sim 1,8\sigma$) avec les récentes simulations N-corps, qui prédisent une abondance de sous-halos plus faible.

Signification
L'article affirme que ces résultats représentent certaines des limites les plus fortes sur la WDM à ce jour et la mesure la plus précise de l'abondance des sous-halos autour des lentilles fortes. La signification réside dans la démonstration que la modélisation conjointe des arcs étendus et des rapports de flux isole efficacement le signal de la sous-structure de la matière noire des incertitudes liées au profil de masse du déflecteur principal. Les auteurs notent que leurs contraintes surpassent les limites précédentes de lentillage fort et sont compétitives ou plus fortes que d'autres sondes à petite échelle comme la forêt Lyman-$\alpha$ et les comptages de galaxies naines. Ce travail établit un cadre robuste pour les futures inférences de DM utilisant les grands échantillons de lentilles attendus des relevés à venir (Euclid, Rubin, Roman), soulignant le potentiel transformateur de la combinaison de l'imagerie spatiale à haute résolution avec une modélisation directe avancée.