Accurate modeling for 3$\times$2pt analyses in Roman and Rubin: a study of model approximations

Résumé Technique : Modélisation précise pour les analyses 3×2pt dans Roman et Rubin

Énoncé du Problème
La cosmologie moderne repose largement sur les analyses « 3×2pt », qui combinent le regroupement de galaxies (clustering), le cisaillement gravitationnel des galaxies (galaxy-galaxy lensing) et le cisaillement de la matière (cosmic shear) pour contraindre les paramètres cosmologiques. Bien que ce cadre ait été appliqué avec succès aux relevés de l'Étape III (ex. DES, KiDS), les relevés de l'Étape IV à venir — spécifiquement le High Latitude Imaging Survey (HLIS) du télescope spatial Nancy Grace Roman et le Legacy Survey of Space and Time (LSST) de l'Observatoire Vera C. Rubin — fourniront une puissance statistique sans précédent. Les auteurs postulent que les approximations de modélisation acceptables pour les données de l'Étape III pourraient introduire des biais significatifs dans les contraintes de l'Étape IV. Cet article étudie systématiquement trois approximations de modélisation spécifiques couramment utilisées dans les analyses 3×2pt :

1. L'application de l'approximation de Limber pour projeter les spectres de puissance 3D sur les échelles angulaires.
2. La négligence des distorsions dans l'espace des redshifts (RSD) dans le champ de surdensité des galaxies.
3. L'utilisation de modèles moins précis pour le spectre de puissance de la matière non linéaire (comparant spécifiquement HMCode à EuclidEmulator2).

Méthodologie
L'étude utilise un cadre d'analyse de vraisemblance simulé implémenté dans CoCoA (CosmoLike-Cobaya Architecture), un pipeline développé par le Roman HLIS Cosmology PIT. Pour garantir la robustesse, les auteurs ont d'abord validé CoCoA par rapport à la bibliothèque établie Core Cosmology Library (CCL), qui sert de standard pour le Rubin LSST Dark Energy Science Collaboration (DESC).

• Configurations de Relevés : L'analyse prévoit les contraintes pour un ensemble de données Roman HLIS sur 5 ans (tier moyen) et un ensemble de données Rubin LSST sur 10 ans. Les deux utilisent 10 bacs de galaxies lentilles et des nombres variables de bacs de galaxies sources (8 pour Roman, 5 pour Rubin).
• Cadre de Modélisation : L'analyse suppose une cosmologie $\Lambda$CDM plate avec cinq paramètres libres ($A_s, n_s, H_0, \Omega_b, \Omega_m$). Les paramètres de nuisance incluent les décalages de photo-z, les biais de calibration de cisaillement, l'alignement intrinsèque (modèle NLA) et le biais de galaxie linéaire.
• Coupures d'Échelle (Scale Cuts) : Pour atténuer les incertitudes liées au feedback baryonique et au biais de galaxie non linéaire, les auteurs appliquent des coupures d'échelle conservatrices. En espace réel, cela correspond à un rayon transverse minimal de 6 Mpc/h pour le cisaillement et de 8 Mpc/h pour le regroupement. En espace harmonique, cela se traduit par des plages de multipôles spécifiques ($20 < \ell < 4000$).
• Conception Expérimentale : Les auteurs définissent cinq scénarios de modèles pour isoler les effets des approximations :
  1. Référence (Baseline) : Non-Limber pour le regroupement, RSD inclus, EuclidEmulator2 pour le spectre de puissance non linéaire.
  2. Limber : Utilise l'approximation de Limber pour le regroupement.
  3. noRSD : Omet les corrections RSD.
  4. HMCode : Utilise HMCode au lieu de EuclidEmulator2.
  5. Combiné : Utilise Limber, omet les RSD et utilise HMCode.

Résultats Clés
L'étude quantifie l'impact de ces approximations sur les vecteurs de données (via $\Delta\chi^2$) et sur les contraintes cosmologiques résultantes (biais dans le plan $\Omega_m - S_8$).

1. Validation du Code : La comparaison entre CoCoA et CCL a montré un excellent accord, avec des valeurs de $\Delta\chi^2$ bien inférieures à 1 pour les réglages fiduciaux, confirmant la fiabilité du pipeline de modélisation.
2. Impacts Individuels des Approximations :
   • Approximation de Limber : Négliger le calcul non-Limber pour le regroupement de galaxies induit des biais de ~0,7–0,8$\sigma$ pour Roman et de ~1,8–2,0$\sigma$ pour Rubin. L'effet est plus prononcé en espace harmonique qu'en espace réel car les coupures d'échelle retirent moins de modes affectés à grande échelle en espace harmonique.
   • Distorsions dans l'Espace des Redshifts (RSD) : Négliger les RSD produit des biais comparables ou légèrement supérieurs à l'approximation de Limber, atteignant ~0,7–0,8$\sigma$ pour Roman et ~2,0–2,3$\sigma$ pour Rubin. Les auteurs notent que les RSD sont l'effet dominant entre les deux, car l'approximation de Limber supprime intrinsèquement les contributions des RSD.
   • Spectre de Puissance Non Linéaire : Passer de EuclidEmulator2 à HMCode entraîne des biais plus faibles (~0,6$\sigma$ pour Roman, ~1,0–1,2$\sigma$ pour Rubin). Ceci est dû au fait que les différences entre ces modèles se situent principalement sur les petites échelles, qui sont largement éliminées par les coupures d'échelle fiduciaux.
3. Effets Combinés : Lorsque les trois approximations sont appliquées simultanément, les biais s'accumulent. Pour Rubin, le biais combiné atteint ~2,8$\sigma$ (espace réel) et ~3,1$\sigma$ (espace harmonique). Pour Roman, le biais combiné est de ~1,0–1,1$\sigma$.
4. Sensibilité du Relevé : Les biais sont systématiquement plus élevés pour Rubin que pour Roman en raison de la plus grande capacité de contrainte de Rubin. Les auteurs notent que bien qu'ils aient supposé un « tier moyen » pour Roman, une configuration de « tier large » produirait probablement des contraintes encore plus serrées et potentiellement des biais plus importants si les approximations sont utilisées.

Signification et Conclusions
L'article conclut que pour les relevés de l'Étape IV, les approximations standards utilisées pour les analyses de l'Étape III sont insuffisantes. Négliger l'un des trois effets étudiés peut induire des biais cosmologiques approchant ou dépassant 1$\sigma$ pour Roman et 2$\sigma$ pour Rubin.

Les auteurs soulignent trois implications principales pour les futures analyses de l'Étape IV :

1. Algorithmes Non-Limber : Il est crucial d'optimiser les algorithmes non-Limber pour garantir qu'ils soient numériquement stables et efficaces sur le plan computationnel, car l'approximation de Limber n'est plus suffisante pour le regroupement de galaxies dans les relevés de haute précision.
2. Modélisation des RSD : L'effet RSD à l'ordre linéaire a un impact important sur les vecteurs de données et les contraintes, motivant le développement et l'inclusion de modèles de RSD plus réalistes.
3. Spectre de Puissance Non Linéaire : Le choix du modèle de spectre de puissance de la matière non linéaire est significatif, particulièrement lorsqu'il est combiné à des stratégies spécifiques de coupure d'échelle.

Cette étude sert de validation du pipeline CoCoA par rapport à CCL et fournit une base quantitative exigeant une modélisation de plus haute fidélité pour les publications de données finales des relevés Roman et Rubin. Les auteurs suggèrent que les travaux futurs devraient étendre ces conclusions pour inclure le biais de galaxie non linéaire, des modèles d'alignement intrinsèque plus complexes et le feedback baryonique, ainsi qu'étudier la combinaison des tiers moyen et large de Roman.