Cosmic variance or galaxy bias? Disentangling finite-volume and galaxy formation effects in cosmological analysis

Résumé Technique : Découpler les effets de volume fini et de formation des galaxies dans l'analyse cosmologique

Énoncé du Problème
Les relevés actuels et futurs de redshifts de galaxies à grande échelle (par exemple, DESI, Euclid) visent à déterminer les paramètres cosmologiques avec une précision sans précédent. Cependant, l'obtention de résultats non biaisés nécessite un contrôle rigoureux des systématiques théoriques. Un défi majeur consiste à distinguer deux sources d'incertitude fortement couplées : la variance cosmique, provenant de la variance intrinsèque du champ aléatoire gaussien initial dans un volume fini, et le biais de galaxie, qui décrit les processus complexes non linéaires de formation des galaxies.

Dans les cadres perturbatifs standards, les paramètres de biais sont estimés en ajustant les observations au spectre de puissance de la matière dont la moyenne d'ensemble est calculée. Cette approche suppose implicitement que le volume du relevé représente la moyenne d'ensemble, négligeant ainsi la variance cosmique spécifique de la réalisation cosmique observée. Par conséquent, les paramètres de biais estimés deviennent dégénérés, confondant la physique réelle de la formation des galaxies avec les fluctuations statistiques du volume cosmique spécifique sondé. Cette dégénérescence complique des applications telles que la reconstruction des Oscillations Acoustiques de Baryons (BAO) et l'inférence de paramètres cosmologiques, où les hypothèses sur les valeurs de biais sont critiques.

Méthodologie
Les auteurs développent un nouveau formalisme basé sur la Théorie de la Perturbation Eulérienne (EPT) pour distinguer mathématiquement la variance cosmique du biais de galaxie.

1. Recadrage de la variance cosmique en tant que biais : Les auteurs proposent de traiter le champ de densité d'une réalisation arbitraire de volume fini ($\delta'$) comme un traceur « biaisé » du champ de densité de la moyenne d'ensemble ($\delta$). En introduisant un nouvel ensemble de paramètres de biais ($b_{1,dm}, b_{2,dm}$), ils expriment le champ de volume fini comme suit :
   $$\delta'(x) = b_{0,dm} + b_{1,dm}\delta(x) + \frac{1}{2}b_{2,dm}\delta^2(x) + \dots$$
   Dans ce cadre, à mesure que le volume tend vers l'infini, $b_{1,dm} \to 1$ et $b_{2,dm} \to 0$. Cette formulation permet de séparer le « biais de variance cosmique » du « biais de galaxie » physique ($b_{1,g}, b_{2,g}$).

2. Application au décalage BAO : Les auteurs appliquent ce formalisme au décalage non linéaire de l'échelle BAO. En utilisant les expressions analytiques dérivées de [22], ils décomposent le décalage total BAO ($\Delta \alpha$) en deux composantes distinctes :

   • Une composante pilotée par la variance cosmique intrinsèque du champ de matière noire (contrôlée par $b_{1,dm}, b_{2,dm}$).
   • Une composante pilotée par le biais des galaxies (contrôlée par $b_{1,g}, b_{2,g}$).
     Cette décomposition révèle que la variance cosmique seule peut induire un décalage BAO non négligeable, particulièrement dans les petits volumes, ce qui est dégénéré avec le décalage causé par la croissance gravitationnelle non linéaire.

3. Gestion de la contrainte d'intégrale : L'article traite de la contrainte d'intégrale (le biais introduit en forçant la densité moyenne d'un volume de relevé à zéro). Les auteurs démontrent que cet effet peut être réabsorbé dans le terme constant de l'expansion du biais ($b'_0$), garantissant qu'il n'introduit pas de paramètres de biais supplémentaires ni ne complique la séparation de la variance cosmique.

4. Validation Numérique : Pour valider les arguments analytiques, les auteurs utilisent la suite de simulations UNIT. Ils analysent quatre ensembles distincts de simulations N-corps (volumes de $2 \times 10^3$ particules dans $(1 h^{-1}\text{Gpc})^3$) :

   • 100 simulations standards.
   • 100 simulations à amplitude fixe (supprimant la variance d'amplitude).
   • 2 $\times$ 50 simulations appariées (supprimant la variance de phase).
   • 2 $\times$ 50 simulations fixes et appariées (supprimant les deux).
     Ils mesurent la position du pic BAO de manière indépendante du modèle (en ajustant une gaussienne à la fonction de corrélation résiduelle) pour quantifier la réduction de la variance.

Résultats Clés

• Découplage Analytique : Le nouveau formalisme sépare avec succès les contributions de la variance cosmique et du biais de galaxie au décalage non linéaire de la BAO. L'analyse confirme que la variance cosmique peut induire un décalage comparable à celui du biais de galaxie, nécessitant une attention particulière lors de l'utilisation d'approches perturbatives sur des volumes finis.
• Suppression de la Variance : Les résultats numériques des simulations UNIT démontrent que le contrôle des conditions initiales réduit considérablement l'incertitude sur la position du pic BAO :
  • Fixer (contrôler l'amplitude) réduit l'incertitude d'un facteur $\sim 1,5 - 2,1$.
  • Apparier (contrôler la phase) réduit l'incertitude d'un facteur $\sim 1,5 - 1,7$.
  • Combiné (Fixe-et-Appairé) : La combinaison de ces deux techniques produit un facteur de suppression de la variance de $\sim 3,4 - 5,4$ par rapport aux ensembles standards.
• Convergence : La position du pic BAO est jugée bien convergée au niveau du sous-pourcentage à travers tous les types de simulations, sans introduction apparente de systématiques par les techniques de fixation ou d'appariement par rapport à l'ensemble standard.

Signification et Stratégies Proposées
L'article conclut que les méthodes actuelles d'estimation des paramètres de biais incluent par inadvertance les effets dégénérés de la variance cosmique. Pour y remédier, les auteurs proposent deux stratégies pratiques pour les futures analyses cosmologiques :

1. Combinaison de Multi-Traceurs Optimale : Au lieu de viser uniquement la précision statistique, les analyses futures pourraient combiner des traceurs complémentaires (par exemple, les galaxies et les vides cosmiques) spécifiquement pour annuler la variance cosmique. Cela pourrait impliquer de redéfinir les échantillons pré-reconstruction pour tracer des régimes de densité complémentaires ou d'optimiser les poids dans la fonction de corrélation combinée pour imiter l'effet des réalisations appariées.
2. Reconstruction à Variance Supprimée : Les auteurs suggèrent d'appliquer des techniques de suppression de la variance (fixation et appariement) à la fois aux données observées (via des méthodes de reconstruction comme l'inversion inverse ou la reconstruction bayésienne directe) et aux catalogues de simulations utilisés pour l'estimation de la covariance. Cela permettrait l'analyse d'une version « à variance supprimée » de l'Univers observé et de son anti-Univers apparié. Ils notent que cela nécessite une mise en œuvre prudente pour éviter de biaiser le décalage BAO par le biais de la cosmologie adoptée dans le processus de reconstruction.

Les auteurs soulignent que bien que ces stratégies soient prometteuses, elles nécessitent des tests supplémentaires sur des catalogues de simulations détaillés (similaires à DESI) et une intégration avec les techniques d'analyse de pointe telles que la reconstruction BAO et la théorie des champs effectifs.