Euclid: Optimising tomographic redshift binning for 3$\times$2pt power spectrum constraints on dark energy

Résumé technique : Optimisation du découpage tomographique en tranches de décalage vers le rouge pour les contraintes du spectre de puissance 3×2pt d'Euclid DR1

Énoncé du problème
L'analyse tomographique du signal 3×2pt (combinant le cisaillement cosmique, le regroupement angulaire des galaxies et le lentillage gravitationnel de type galaxy-galaxy) nécessite de diviser l'échantillon de galaxies observé en tranches de décalage vers le rouge (redshift). Bien que l'augmentation du nombre de tranches permette généralement de capturer plus d'informations cosmologiques, elle augmente également les coûts de calcul et le bruit de tir de Poisson. Des études antérieures examinant les stratégies de découpage « optimales » pour les relevés de l'Étape IV comme Euclid ont produit des conclusions incohérentes concernant le nombre idéal de tranches et la méthode optimale pour définir les limites des tranches (par exemple, tranches équipopulées, largeur de redshift égale ou distance comobile égale). Cette publication vise à résoudre ces incohérences en présentant une méthode basée sur des simulations pour déterminer la stratégie de découpage tomographique optimale spécifiquement pour la première version des données d'Euclid (DR1), en se concentrant sur les contraintes des paramètres de l'équation d'état de l'énergie noire $(w_0, w_a)$.

Méthodologie
Les auteurs ont développé un pipeline de simulation (SWEPT) pour générer des catalogues de galaxies fictifs de type Étape IV imitant une observation réaliste d'Euclid DR1. Le flux de travail comprend :

1. Génération de la population de galaxies : Échantillonnage d'une population de galaxies à partir d'une distribution de décalage vers le rouge $p(z)$ prédite pour Euclid et d'une densité numérique spatiale.
2. Calcul du signal théorique : Calcul des spectres de puissance 3×2pt théoriques ($C^{\gamma\gamma}_{ij}$, $C^{\delta_g\delta_g}_{ij}$, $C^{\delta_g\gamma}_{ij}$) en utilisant CosmoSIS, CAMB et halofit sous une cosmologie fiducielle plate $w_0w_a$CDM ($\Omega_m=0,3, H_0=70, w_0=-1, w_a=0$).
3. Génération de cartes : Création de cartes HEALPix 2D des champs corrélés de cisaillement cosmique et de densité de galaxies sur la sphère céleste en utilisant FLASK, approximant l'information cosmologique 3D complète.
4. Construction du catalogue fictif : Échantillonnage de Poisson des champs de densité pour placer les galaxies et assignation des observables de lentillage faible (convergence et cisaillement) à partir des cartes de champs aux positions des galaxies.
5. Injection de bruit : Introduction d'incertitudes observationnelles réalistes, incluant des erreurs de décalage vers le rouge photométrique (photo-z) distribuées selon une loi gaussienne ($\sigma_z = \sigma_{phot}(1+z)$ avec $\sigma_{phot}=0,05$) et un bruit de forme gaussien ($\sigma(\epsilon_{int})=0,3$). Un sous-ensemble de simulations a également inclus 5 % d'anomalies (outliers) catastrophiques de photo-z.
6. Analyse : Mesure des spectres de puissance Pseudo-$C_\ell$ pour diverses configurations tomographiques (variant le nombre de tranches et les définitions de limites) en utilisant NaMaster. Les contraintes sur $(w_0, w_a)$ ont été dérivées via une analyse de vraisemblance gaussienne sur grille à l'aide d'une matrice de covariance numérique dérivée de 3 000 réalisations, corrigée par le facteur de Hartlap.

Résultats clés
L'étude a évalué trois stratégies de découpage : des tranches équipopulées, des tranches espacées de manière égale en décalage vers le rouge, et des tranches espacées de manière égale en distance comobile fiduciale.

• Stratégie de découtage optimale :

  • Plein 3×2pt et regroupement angulaire : Pour le signal complet 3×2pt et la composante de regroupement angulaire seule, les tranches équipopulées produisent les contraintes les plus serrées sur $(w_0, w_a)$. Cela est vrai tant pour les simulations sans bruit que pour les simulations réalistes avec bruit.
  • Cisaillement cosmique uniquement : La stratégie optimale dépend du modèle de bruit. Dans les simulations sans bruit, les tranches à largeur de redshift égale fournissent les meilleures contraintes. Cependant, dans les simulations réalistes incluant le bruit de photo-z et de forme, les tranches espacées de manière égale en distance comobile fiduciale produisent les meilleures contraintes, bien que la différence entre les choix de découpage soit minimale (< 5 %).
  • Dominance du regroupement : La préférence pour les tranches équipopulées dans l'analyse 3×2pt complète est principalement dictée par la composante de regroupement angulaire, qui domine le pouvoir de contrainte du signal conjoint.

• Nombre de tranches et saturation de l'information :

  • Le gain d'information sature au-delà de six à huit tranches pour toutes les composantes du signal (cisaillement, regroupement, et 3×2pt complet).
  • Augmenter le nombre de tranches de 7 à 10 n'apporte que des gains marginaux (~1–2 % d'amélioration).
  • Étendre l'analyse à 13 tranches (une configuration proposée dans d'autres publications) ne produit qu'une amélioration supplémentaire d'environ 5 % de la précision par rapport à 7 tranches. Les auteurs notent que les gains marginaux de tranches supplémentaires sont probablement limités par l'exigence croissante de précision pour l'estimation de la matrice de covariance et la présence d'autres systématiques dans les données réelles.

• Impact des erreurs de photo-z catastrophiques :

  • L'introduction de 5 % d'anomalies (outliers) de photo-z sans atténuation entraîne des biais sévères.
  • Pour les analyses avec 5 ou 10 tranches, les contraintes résultantes sur $(w_0, w_a)$ sont incohérentes avec la cosmologie $\Lambda$CDM fiduciale à $\gtrsim 3\sigma$, favorisant un modèle d'énergie noire en évolution temporelle ($w_0 > -1, w_a < 0$).
  • Les distributions postérieures présentent des caractéristiques nettes et non elliptiques, indiquant que les erreurs catastrophiques non atténuées induisent des biais complexes et non triviaux dans les spectres de puissance.

Signification et affirmations
L'article affirme fournir un cadre robuste basé sur des simulations qui imite une véritable observation de relevé, offrant une évaluation plus directe des stratégies de découpage tomographique que les études théoriques précédentes. Les principales conclusions sont :

1. Les tranches équipopulées sont la stratégie recommandée pour les analyses Euclid DR1 impliquant le signal complet 3×2pt ou le regroupement angulaire, car elles fournissent les contraintes les plus robustes.
2. La composante de cisaillement cosmique est relativement insensible à la méthodologie de découpage spécifique (différences < 5 % dans les scénarios avec bruit), suggérant que le choix du découpage pour le cisaillement seul est moins critique que pour le regroupement.
3. Il existe un rendement décroissant dans l'augmentation du nombre de tranches au-delà de sept ou huit ; les gains statistiques marginaux sont probablement surcompensés par les coûts de calcul et la difficulté accrue de modéliser les systématiques et les matrices de covariance.
4. Les erreurs de photo-z catastrophiques constituent une menace critique pour l'inférence cosmologique. Si elles ne sont pas atténuées, elles peuvent induire des biais qui imitent une déviation de $\Lambda$CDM avec une signification élevée, nécessitant des stratégies de modélisation ou d'atténuation rigoureuses dans le pipeline d'analyse.

Les auteurs soulignent que, bien que leurs résultats suggèrent que les tranches équipopulées sont optimales pour les caractéristiques spécifiques d'Euclid DR1, l'insensibilité relative du cisaillement cosmique aux choix de découpage et la saturation du gain d'information suggèrent que la stratégie de découpage exacte peut être moins critique que la modélisation précise des systématiques et de la matrice de covariance.