Dynamical Systematics for Time Delay Lenses and the Impact on the Hubble Constant

Résumé technique : Systématique dynamique pour les lentilles à retard de temps et l'impact sur la constante de Hubble

Énoncé du problème
Les lentilles gravitationnelles à retard de temps offrent une méthode pour mesurer la constante de Hubble ($H_0$) indépendamment de l'échelle de distance cosmique. Cependant, ces mesures souffrent d'une dégénérescence fondamentale où $H_0 \propto 1 - \kappa_E$, avec $\kappa_E$ étant la convergence près du rayon d'Einstein. Les données de lentille contraignent le rayon d'Einstein ($R_E$) et une quantité adimensionnelle $\xi$ liée à la pente du profil de masse, mais elles ne peuvent pas déterminer $\kappa_E$ sans hypothèses supplémentaires. Des mesures de dynamique stellaire, spécifiquement les dispersions de vitesse centrales ($\sigma$), sont utilisées pour lever cette dégénérescence. La précision de $H_0$ dérivée de cette méthode est limitée par la précision des mesures cinématiques ; atteindre une précision de 2 % pour $H_0$ nécessite que les incertitudes systématiques de la dispersion de vitesse au carré ($\Delta\sigma^2/\sigma^2$) soient inférieures à 2 %. Cet article étudie la plage d'incertitudes systématiques potentielles affectant huit lentilles à retard de temps utilisées dans les analyses cosmologiques, arguant que les estimations actuelles de $H_0$ pourraient être biaisées par des systématiques non prises en compte dans la modélisation dynamique.

Méthodologie
Les auteurs analysent huit lentilles à retard de temps (DES J0408−5354, HE 0435−1223, PG 1115+080, RX J1131−1231, SDSS J1206+4332, B1608+656, WFI 2033−4723, et WGD 2038−4008) en utilisant un modèle dynamique à symétrie sphérique. L'étude emploie une distribution stellaire de Hernquist enchâssée dans un modèle de masse de Sphère Isotherme Singulière (SIS) pour isoler les effets systématiques. L'analyse se concentre sur quatre catégories principales de systématiques :

1. Modélisation de la Fonction d'Étalement du Point (PSF) : Les auteurs testent la sensibilité des mesures de dispersion de vitesse aux erreurs de la largeur à mi-hauteur (FWHM) de la PSF, au décentrage de l'ouverture, et à l'hypothèse d'une PSF gaussienne par rapport à des profils de Moffat plus réalistes.
2. Dispersion mesurée vs réelle : L'étude distingue la dispersion de vitesse mesurée ($\sigma_*$, dérivée de l'ajustement d'une gaussienne aux raies spectrales) de la vitesse efficace quadratique moyenne ($v_{rms}$) requise pour les équations de Jeans. Les auteurs quantifient le biais introduit par les distributions de vitesse non gaussiennes, caractérisées par le moment de Gauss-Hermite $h_4$.
3. Modèles d'anisotropie : Les auteurs comparent des modèles d'anisotropie constante, le modèle d'Osipkov-Merritt (O-M), et des profils alternatifs (Mamon & Lokas ; Cuddeford) pour déterminer comment le choix de l'anisotropie orbitale affecte la dispersion de vitesse inférée et $H_0$.
4. Modèles photométriques et de populations stellaires : L'analyse examine l'impact de l'utilisation de différents profils de densité stellaire (Sérsic vs Hernquist vs Jaffe), du choix de l'axe du rayon effectif (majeur vs intermédiaire), et des complications découlant des gradients de couleur et des populations stellaires multiples dans les galaxies de type précoce.

De plus, l'article traite du traitement statistique de l'échantillon de lentilles, arguant que les galaxies de type précoce constituent une population homogène en termes de dynamique et de populations stellaires, ce qui implique que certaines erreurs systématiques devraient être marginalisées sur l'ensemble de l'échantillon plutôt que traitées comme des erreurs aléatoires indépendantes pour chaque lentille.

Contributions clés et résultats
Les auteurs constatent que les incertitudes systématiques dans la modélisation dynamique peuvent facilement dépasser le seuil de 2 % requis pour des mesures compétitives de $H_0$. Les conclusions spécifiques incluent :

• Systématiques de la PSF : Pour les observations au sol, les biais dans la largeur et la forme de la PSF peuvent induire des changements fractionnaires de $\sigma^2$ de 1 à 5 %. Sous-estimer la largeur de la PSF ou supposer un profil gaussien alors que la vraie PSF possède des ailes étendues (Moffat) conduit à sous-estimer les dispersions de vitesse, poussant les modèles vers des distributions de masse moins concentrées et des valeurs de $H_0$ plus faibles. Le décentrage de l'ouverture a un effet moindre (typiquement < 0,3 %) mais biaise également les résultats vers le bas.
• Définition de la dispersion : La différence entre la $\sigma_*$ mesurée et la $v_{rms}$ requise est significative. Pour les modèles isotropes, cette différence varie de ~1,8 % à 5,6 % à travers l'échantillon. Si les galaxies de type précoce sont modérément anisotropes radialement (comme on peut s'y attendre), la vraie $v_{rms}$ est plus grande que la $\sigma_*$ mesurée, ce qui conduit à des sous-estimations de $H_0$ si elle n'est pas corrigée.
• Modèles d'anisotropie : Le choix du modèle d'anisotropie a un impact profond. Bien que les modèles d'anisotropie constante et les modèles O-M avec de grands rayons de rupture se situent dans une enveloppe étroite, les modèles permettant une anisotropie radiale modérée dans les régions centrales (par exemple, les modèles de Cuddeford avec $\beta_0 > 0$) peuvent augmenter la dispersion de vitesse inférée de 2 à 18 %. Cela déplace la distribution de masse pour qu'elle soit moins concentrée, pouvant biaiser $H_0$ vers le haut si elle n'est pas correctement contrainte. De plus, les modèles produisant des valeurs de $h_4$ plus élevées (cohérentes avec les observations des galaxies de type précoce) nécessitent des corrections plus importantes entre $\sigma_*$ et $v_{rms}$.
• Profils photométriques : Les décalages entre le profil de masse stellaire supposé et la distribution réelle (par exemple, utiliser un profil de Sérsic avec un indice incorrect) peuvent causer des variations de $\Delta\sigma^2/\sigma^2$ de 1 à 10 %. De plus, la présence de gradients de couleur implique que la densité numérique stellaire $n$ dans l'équation de Jeans devrait être pondérée par la largeur équivalente des raies spectrales plutôt que par le profil de lumière à large bande, un facteur actuellement non pris en compte dans les analyses standards.
• Homogénéité de la population : L'article soutient que traiter les paramètres tels que le rayon d'anisotropie ou le choix du modèle stellaire comme des variables aléatoires indépendantes pour chaque lentille sous-estime l'incertitude totale. Puisque les galaxies de type précoce sont une population homogène, les erreurs systématiques dans ces paramètres sont corrélées à travers l'échantillon. La marginalisation de ces paramètres au niveau de la population (plutôt que par lentille) empêche la réduction artificielle de l'incertitude par un facteur $N^{-1/2}$.

Signification et revendications
L'article affirme que la précision actuelle des estimations de $H_0$ issues des lentilles à retard de temps est probablement limitée par les incertitudes systématiques de la modélisation dynamique plutôt que par le bruit statistique. Les auteurs démontrent que des biais de $\sigma^2$ de 1 à 18 % sont possibles selon les choix de modélisation, ce qui est comparable ou supérieur à la précision de 2 % requise pour $H_0$.

Les auteurs soulignent que :

1. Le calibrage absolu importe : Contrairement aux études de dynamique stellaire pure qui reposent sur des rapports, le lentillage dépend du calibrage absolu des dispersions de vitesse, ce qui le rend très sensible aux erreurs systématiques de modélisation de la PSF et de définition de la dispersion.
2. Dépendances du modèle : Le choix du modèle d'anisotropie et le traitement de la fonction de distribution des vitesses (spécifiquement la correction $h_4$) sont critiques. Les modèles actuels échouent souvent à reproduire les valeurs de $h_4$ observées dans les galaxies de type précoce, suggérant que les corrections actuelles sont des sous-estimations.
3. Traitement statistique : L'hypothèse d'une population de lentilles hétérogène conduit à une sous-estimation des incertitudes finales. Un modèle de population homogène nécessite de marginaliser les paramètres systématiques sur l'ensemble de l'échantillon, ce qui préserve la contribution de l'erreur systématique quel que soit la taille de l'échantillon.

L'article conclut que l'obtention d'une mesure de $H_0$ à 2 % nécessite de contrôler ces incertitudes systématiques, particulièrement le calibrage absolu des dispersions de vitesse et la modélisation correcte de l'anisotropie orbitale et des populations stellaires. Bien que la cinématique résolue spatialement (données 2D) offre une voie pour briser les dégénérescences, les biais systématiques fondamentaux concernant la définition de la dispersion de vitesse et l'homogénéité de la population doivent toujours être adressés.