Accurate Galaxy Cluster Shear and Mass Calibration for LSST with AnaCal

Résumé technique : Calibration précise du cisaillement et de la masse des amas de galaxies pour LSST avec AnaCal

Énoncé du problème
La cosmologie des amas de galaxies repose sur l'établissement d'une relation précise entre les observables des amas et leurs masses, un processus connu sous le nom de calibration de masse. Le cisaillement gravitationnel faible est l'outil standard pour cette calibration car il sonde la masse totale (incluant la matière noire) sans dépendre d'hypothèses d'équilibre dynamique ou thermique. Cependant, un défi significatif surgit dans l'environnement des amas : le signal de cisaillement à proximité des centres d'amas peut atteindre le régime non linéaire (par exemple, $|g| \sim 0,15$), où de nombreux estimateurs de cisaillement reposent sur des hypothèses perturbatives qui nécessitent une validation explicite. De plus, le lentillage des amas introduit des complexités uniques, absentes des mesures de cisaillement cosmique standard, notamment :

• Régime de cisaillement élevé : Le cisaillement tangentiel dans les régions internes ($\sim 0,2$ Mpc) est plusieurs fois plus important que les valeurs typiquement utilisées pour les tests de calibration basés sur la perturbation.
• Déflexion inhomogène : Le champ de lentillage de l'amas déplace les positions apparentes des galaxies d'arrière-plan, modifiant leurs configurations de chevauchement et de mélange (blending), ce qui impacte la détection et la mesure de forme.
• Dépendance radiale : Le champ de convergence varie radialement, ce qui peut induire une réponse de cisaillement moyenne décroissante radialement.
• Effets d'ordre supérieur : Près du centre de l'amas, des distorsions d'ordre supérieur (par exemple, la flexion) et la lumière anisotrope provenant des halos intra-amas peuvent introduire des signaux parasites.

Bien que la collaboration LSST Dark Energy Science Collaboration (DESC) ait sélectionné AnaCal et la métadétection comme méthodes fiduciels d'estimation du cisaillement, leur performance n'avait pas été rigoureusement testée sous les conditions spécifiques de cisaillement élevé et de déflexion inhomogène caractéristiques des lentilles d'amas.

Méthodologie
Les auteurs utilisent des simulations d'images représentatives des données de 10 ans de LSST pour tester l'estimateur de cisaillement AnaCal. La méthodologie consiste en :

1. Cadre de simulation : En utilisant descwl-shear-sims et xlens, l'équipe a généré des simulations d'images de lentillage d'amas réalistes. Ils ont modélisé 9 masses de halo allant de $10^{14} M_\odot$ à $10^{15} M_\odot$ comme des halos NFW à une décalage vers le rouge de la lentille $z_l = 0,25$, avec des galaxies sources à $z_s = 1,0$.

   • Les simulations comprenaient 2178 réalisations par masse, couvrant environ 1700 deg$^2$ au total.
   • Pour isoler les effets d'estimation du cisaillement, les simulations ont exclu les galaxies membres de l'amas, la lumière intra-amas et les erreurs de calibration de redshift, se concentrant strictement sur l'impact du champ de déflexion à haut cisaillement et inhomogène.
   • Une PSF de type Moffat circulaire constante ($FWHM = 0,8''$) a été utilisée pour représenter la qualité d'image attendue de l'Observatoire Rubin.

2. Implémentation d'AnaCal : L'étude emploie le cadre AnaCal, qui dérive analytiquement la réponse de cisaillement de l'ensemble du pipeline de mesure (observables, détection et sélection) en utilisant l'auto-différentiation et une base de shapelets.

   • Observables : Les images sont projetées sur des noyaux de détection et des fonctions de base de shapelets polaires.
   • Calcul de la réponse : La réponse de cisaillement est calculée en différenciant les observables par rapport à des perturbations de cisaillement contrôlées.
   • Correction du bruit : Un schéma de « renoising » est appliqué pour corriger le biais de bruit en ajoutant des réalisations de bruit synthétique aux images.
   • Pondération : Des poids de sélection et de détection sont appliqués en fonction du rapport signal sur bruit (SNR) et de la taille de la galaxie, conçus pour être invariants sous des rotations de 90 degrés.

3. Analyse et ajustement :

   • Estimation du cisaillement : Les ellipticités tangentielles et croisées sont calculées, et la matrice de réponse est pivotée dans le référentiel tangentiel. Crucialement, les auteurs utilisent la réponse de cisaillement dans chaque classe angulaire séparément plutôt que de moyenner sur toutes les classes afin d'éviter le biais dû à la dépendance radiale du champ de convergence.
   • Ajustement de la masse : Une approche de modélisation directe ajuste le profil de cisaillement réduit tangentiel observé ($g_T$) à des modèles NFW en utilisant une fonction de vraisemblance qui incorpore une matrice de covariance analytique dérivée des prévisions de 10 ans de LSST.
   • Quantification du biais : Un biais multiplicatif de cisaillement ($m$) est ajusté selon la forme $m = a + b(g_{true})^2$, testant les biais linéaires et non linéaires du troisième ordre.

Résultats clés

• Récupération du cisaillement : AnaCal récupère le cisaillement d'entrée avec un biais minimal, même à des niveaux de cisaillement modérément élevés ($|g| \sim 0,15$).
• Tendance de la réponse radiale : L'étude confirme une réponse de cisaillement moyenne décroissante radialement, pilotée par la dépendance radiale du champ de convergence. Les auteurs démontrent que si cet effet n'est pas modélisé (c'est-à-dire si une réponse moyenne globale est utilisée), cela biaise l'inférence du cisaillement.
• Biais du troisième ordre : Un biais d'estimation du cisaillement positif est détecté au troisième ordre du cisaillement réduit près du centre de l'amas.
• Précision de la calibration de la masse : Malgré la présence de régimes de cisaillement élevé et de biais du troisième ordre, l'impact sur l'estimation finale de la masse est négligeable. Pour les masses d'amas de lentille dans la plage $[10^{14} M_\odot, 10^{15} M_\odot]$ (en adoptant une coupure d'échelle de $\sim 0,2$ Mpc à $z=0,25$), le biais moyen de masse d'amas résultant est de $0,24 \pm 0,26\%$ dans des conditions idéales. Ce faible biais est attribué au fait qu'une petite fraction des galaxies se trouve dans le régime de cisaillement élevé, et que ces mesures sont pondérées à la baisse par la matrice de covariance.

Signification et affirmations
L'article affirme que ces résultats démontrent que AnaCal est un outil robuste pour une calibration précise de la masse des amas à l'ère de LSST. Plus précisément :

• Il valide que l'estimation du cisaillement par perturbation (via AnaCal) reste précise même dans le régime de cisaillement élevé typique des centres d'amas, à condition que la dépendance radiale de la réponse soit correctement modélisée.
• Il confirme que les non-linéarités et les effets de cisaillement élevé ne dégradent pas de manière significative la précision des estimations de masse pour les échantillons d'amas cosmologiquement pertinents.
• Ce travail constitue une étape de validation critique pour le LSST DESC, garantissant que les méthodes d'estimation du cisaillement sélectionnées peuvent répondre aux exigences strictes de la cosmologie des amas sans introduire d'erreurs systématiques significatives.

Les auteurs notent explicitement que leurs simulations n'ont pas inclus les effets des membres de l'amas, de la lumière intra-amas ou des distorsions d'ordre supérieur comme la flexion, isolant ainsi la performance de l'estimateur de cisaillement face au principal défi que constituent le cisaillement élevé et les champs de déflexion inhomogènes.