Euclid: Exploring observational systematics in cluster cosmology -- a comprehensive analysis of cluster counts and clustering

Résumé technique : Euclid : Exploration des effets systématiques en cosmologie des amas

Énoncé du problème
Les amas de galaxies servent de sondes cosmologiques critiques pour contraindre le paramètre de densité de matière ($\Omega_m$) et l'amplitude des fluctuations de densité ($\sigma_8$). Bien que les analyses traditionnelles reposent sur les comptages d'amas, la distribution spatiale des amas (leur agrégation) fournit des informations complémentaires pour briser les dégénérescences dans les paramètres cosmologiques et de relation masse-observable. Cependant, la future mission Euclid, qui vise à cartographier $\sim 10^5$ amas jusqu'à $z=2$, exige une précision inférieure au pourcent. Un défi majeur réside dans la caractérisation de l'impact des effets systématiques observationnels et de modélisation — spécifiquement les incertitudes sur les redshifts photométriques (photo-$z$) et les distorsions de l'espace des redshifts (RSD) — sur les comptages d'amas et les statistiques d'agrégation, ainsi que sur leurs matrices de covariance. De plus, la nécessité de combiner ces sondes pour maximiser le pouvoir de contrainte exige une compréhension rigoureuse de leurs corrélations croisées et de la dépendance des matrices de covariance par rapport à la cosmologie.

Méthodologie
Les auteurs utilisent un cadre basé sur 1 000 catalogues de cônes de lumière de type Euclid, générés à l'aide du code PINOCCHIO (employant la théorie des perturbations lagrangiennes d'ordre trois et l'effondrement ellipsoïdal). Ces simulations couvrent une surface céleste de 10 313 deg$^2$ et des redshifts de $z=0$ à $2.5$, contenant environ $10^5$ halos par catalogue avec $M_{vir} \ge 10^{13.5} h^{-1} M_{\odot}$.

L'étude emploie les étapes méthodologiques suivantes :

1. Modélisation théorique : Les auteurs modélisent les comptages d'amas et le monopôle de la fonction de corrélation à deux points (2PCF), en intégrant les effets des erreurs de photo-$z$ (modélisées comme une dispersion gaussienne), des RSD (effet Kaiser), de la resommation infrarouge (IR) pour l'amortissement des oscillations acoustiques de baryons (BAO) et des distorsions géométriques (GD). Ils modélisent également la masse moyenne par lentille gravitationnelle faible comme sonde complémentaire pour briser les dégénérescences dans les relations d'échelle.
2. Analyse de covariance : Les covariances numériques sont calculées directement à partir des 1 000 simulations de type « mock ». Elles sont comparées à des modèles semi-analytiques (adaptés de Fumagalli et al. 2021 et Euclid Collaboration : Fumagalli et al. 2024) pour valider le traitement du bruit de Poisson, de la variance d'échantillon et des corrélations croisées.
3. Analyse de vraisemblance : Une inférence bayésienne est réalisée à l'aide de PolyChord pour explorer les distributions a posteriori de $\Omega_m$, $\sigma_8$ et de quatre paramètres de relation masse-observable. L'analyse évalue l'Indice de Mérite (FoM) et l'Indice d'Incohérence (IoI) à travers diverses combinaisons de sondes (comptages, agrégation, lentille gravitationnelle faible) et configurations systématiques.
4. Tests systématiques : L'impact de la négligence de systématiques spécifiques (par exemple, RSD, erreurs de photo-$z$) dans la modélisation des statistiques et de leurs matrices de covariance est quantifié en comparant les résultats à une référence de « modèle complet ».

Contributions et résultats clés

• Indépendance des sondes : L'analyse confirme que les comptages d'amas et l'agrégation des amas sont des observables effectivement non corrélés. Le coefficient de corrélation de Pearson entre leurs résidus de vraisemblance logarithmique est $\rho = -0.011 \pm 0.031$, compatible avec zéro. Cette indépendance découle du fait que les comptages mesurent la densité numérique comobile moyenne, tandis que l'agrégation mesure les fluctuations autour de cette moyenne.
• Pouvoir de contrainte de la combinaison : La combinaison de l'agrégation d'amas avec les comptages améliore considérablement les contraintes cosmologiques. Le FoM augmente de plus de 300 % par rapport aux comptages seuls. La combinaison des trois sondes (comptages, agrégation et masse par lentille gravitationnelle faible) produit les contraintes les plus serrées et les plus précises, l'analyse multi-sondes étant moins sensible aux effets de volume a priori que les combinaisons à deux sondes.
• Impact des incertitudes de photo-$z$ : Les erreurs de redshift photométrique constituent l'effet systématique dominant. Une incertitude optimiste ($\sigma_{z0} = 0.005$) élargit les distributions a posteriori de moins de 10 %, tandis qu'une incertitude plus réaliste ($\sigma_{z0} = 0.01$) les élargit de 20 à 30 %. Les erreurs de photo-$z$ réduisent également significativement la variance d'échantillon dans la covariance des comptages.
• Impact des distorsions secondaires : Bien que les RSD, la resommation IR et les distorsions géométriques aient un impact plus faible sur le FoM (5 % pour l'agrégation seule, ~10 % pour la combinaison complète) comparé aux erreurs de photo-$z$, la négligence des RSD introduit un biais significatif (décalage de $\sim 2\sigma$) dans les distributions a posteriori de $\Omega_m$–$\sigma_8$.
• Dépendance à l'échelle : Les données d'agrégation aux échelles inférieures à $60 h^{-1}$ Mpc fournissent un pouvoir de contrainte supplémentaire, tandis que les échelles supérieures à l'échelle des oscillations acoustiques ($>130 h^{-1}$ Mpc) ajoutent une information négligeable sur $\Omega_m$ et $\sigma_8$.
• Modélisation de la covariance : L'étude démontre que la négligence des effets RSD et photo-$z$ dans la covariance des comptages conduit à une sous-estimation de 10 à 20 % de l'amplitude a posteriori. Les auteurs proposent un modèle de spectre de puissance efficace (sommant les termes monopôle et quadrupôle) qui capture avec précision la covariance d'échantillon sans le coût computationnel d'une modélisation anisotrope complète.
• Covariance dépendante de la cosmologie : L'utilisation d'une matrice de covariance dépendante de la cosmologie améliore les contraintes d'agrégation d'environ 40 %. Cependant, lorsque les comptages sont combinés à l'agrégation, ce gain supplémentaire disparaît, suggérant que les comptages encapsulent déjà l'information encodée dans la dépendance cosmologique de la covariance de l'agrégation.

Signification
L'article établit un cadre robuste pour l'analyse conjointe des comptages et de l'agrégation d'amas pour la mission Euclid. Il met en évidence que, si la combinaison de ces sondes améliore drastiquement les contraintes cosmologiques, elle accroît simultanément la sensibilité aux effets systématiques. Le travail souligne la nécessité de modéliser avec précision les incertitudes de photo-$z$ et les RSD non seulement dans les statistiques elles-mêmes, mais également au sein des matrices de covariance, afin d'éviter des contraintes biaisées. En validant l'indépendance des comptages et de l'agrégation et en démontrant l'efficacité de l'approche combinée, l'étude fournit des prescriptions essentielles pour maximiser le retour scientifique de la cosmologie des amas d'Euclid, tout en mettant en garde que des contraintes plus étroites issues d'analyses combinées nécessitent un contrôle systématique encore plus rigoureux que les analyses à sonde unique.