Euclid preparation: Testing multi-field inflation with galaxy power spectrum and bispectrum

Résumé technique : Tester l'inflation multifield avec le spectre de puissance et le bispectre des galaxies

Énoncé du problème
La non-gaussianité primordiale (PNG), en particulier de type local, sert de sonde critique pour distinguer les modèles d'inflation à champ unique des modèles à champs multiples. Bien que le fond diffus cosmologique (CMB) ait imposé des contraintes strictes sur la PNG locale (presque limitées par la variance cosmique), les relevés de galaxies de phase IV comme Euclid promettent de dépasser ces limites en mesurant les statistiques de regroupement à deux points (spectre de puissance) et à trois points (bispectre) sur de vastes volumes cosmologiques. Cependant, l'extraction de contraintes robustes sur le paramètre PNG $f_{\text{NL}}$ à partir du regroupement des galaxies nécessite un cadre d'analyse conjointe rigoureux qui prenne en compte les distorsions de l'espace des vitesses (RSD), l'évolution gravitationnelle non linéaire et le biais des galaxies. Un défi majeur réside dans la dégénérescence entre $f_{\text{NL}}$ et le paramètre de biais induit par la PNG, $b_\phi$, qui est difficile à briser sans des a priori informatifs ou un traitement complet des systématiques théoriques.

Méthodologie
Les auteurs valident un pipeline de modélisation théorique et d'inférence conçu pour l'échantillon de galaxies H$\alpha$ spectroscopiques d'Euclid. Le travail utilise :

• Simulations : Une suite de simulations N-corps Abacus-PNG avec des conditions initiales de type PNG local ($f_{\text{NL}} = \{0, 30\}$) et leurs équivalents gaussiens. Ces simulations sont peuplées de galaxies à l'aide d'une distribution d'occupation des halos (HOD) calibrée pour reproduire les propriétés de regroupement de la simulation Euclid Flagship 2 (FS2). Une procédure d'appariement d'abondance (AM) est appliquée pour concilier les différences dans les définitions de la masse des halos entre les simulations FS2 et Abacus.
• Modélisation théorique : L'analyse emploie la Théorie des Champs Effective de la Structure à Grande Échelle (EFTofLSS). Le spectre de puissance dans l'espace des vitesses est modélisé à l'ordre une boucle, tandis que le bispectre est modélisé au niveau arbre. Le cadre inclut :
  • Une expansion de biais renormalisée étendue pour inclure les opérateurs PNG.
  • Les distorsions de l'espace des vitesses modélisées via un développement en moments de vitesse.
  • Des termes de contre-réaction pour absorber les sensibilités ultraviolettes et modéliser les effets de vitesse à petite échelle (doigts de Dieu), en adoptant spécifiquement un terme de contre-réaction motivé par le cadre du générateur de différence de vitesse (VDG) pour le bispectre.
  • Une resommation infrarouge et les effets Alcock–Paczynski.
• Pipeline d'inférence : Un nouveau code, CosmoSIS-gClust, est introduit pour effectuer des analyses de vraisemblance conjointe des multipôles du spectre de puissance ($P_0, P_2, P_4$) et des multipôles du bispectre ($B_0, B_2$). L'analyse marginalise sur les paramètres de nuisance (termes de contre-réaction, bruit de Poisson) et teste divers a priori sur le paramètre de biais PNG $b_\phi$.
• Métriques : La validité du modèle est évaluée à l'aide de la Figure de Biais (FoB) pour vérifier les erreurs systématiques et de la Figure de Mérite (FoM) pour quantifier le contenu informationnel. Les coupures d'échelle ($k_{\text{max}}$) sont optimisées pour assurer une récupération non biaisée des paramètres cosmologiques.

Contributions clés

1. Validation du pipeline : L'article fournit une validation approfondie d'un pipeline d'analyse conjointe du spectre de puissance et du bispectre pour la PNG locale dans l'espace des vitesses, spécifiquement adapté aux données de type Euclid.
2. Stratégie d'a priori sur le biais : Les auteurs démontrent qu'une approche « agnostique vis-à-vis des a priori » (utilisant des a priori plats larges sur $b_\phi$) conduit à des estimations de $f_{\text{NL}}$ biaisées ou non contraintes en raison de la forte dégénérescence avec $b_\phi$. Ils proposent et valident un a priori gaussien physiquement motivé sur le paramètre $p_{b_\phi}$ (où $b_\phi = 2\delta_c(b_1 - p_{b_\phi})$) qui englobe à la fois les prédictions de la Fonction de Masse Universelle des Halos (UHMF) et du Modèle de Formation des Galaxies (GFM). Ce choix produit des contraintes non biaisées tout en tenant compte de l'incertitude théorique.
3. Optimisation de la coupure d'échelle : L'étude détermine les coupures à petite échelle optimales ($k_{\text{max}}$) pour les deux statistiques. Pour le spectre de puissance, les modèles à une boucle restent non biaisés jusqu'à $k \approx 0,30 \, h \text{ Mpc}^{-1}$. Pour le bispectre, les modèles au niveau arbre restent non biaisés jusqu'à $k \approx 0,12 \, h \text{ Mpc}^{-1}$ (plus faible à des décalages vers le rouge plus élevés), au-delà de quoi l'absence de corrections de boucle introduit un biais.
4. Test nul : Un test nul sur les simulations de conditions initiales gaussiennes ($f_{\text{NL}}=0$) confirme que le pipeline ne détecte pas artificiellement la PNG ni ne biaise les paramètres standards du $\Lambda$CDM.

Résultats

• Significativité de la détection : Dans une configuration « agnostique vis-à-vis des a priori », la détection à haute significativité du terme PNG dominant ($f_{\text{NL}} b_\phi$) dans des tranches de décalage vers le rouge uniques est difficile. Cependant, avec l'a priori informatif proposé sur $b_\phi$, l'analyse conjointe récupère la valeur fiduciale $f_{\text{NL}} = 30$ avec un biais inférieur à $1\sigma$.
• Gain informationnel :
  • À décalage vers le rouge et coupures d'échelle fixes, le bispectre seul réduit l'erreur sur $f_{\text{NL}}$ ($\sigma(f_{\text{NL}})$) d'environ 29 à 46 % par rapport au spectre de puissance seul.
  • L'analyse conjointe spectre de puissance-bispectre resserre davantage les contraintes d'un supplément de 8 à 13 % par rapport au bispectre seul.
  • Le gain cumulatif de précision de $z=0,8$ à $z=1,7$ est d'un facteur d'environ 2,3 pour le cas conjoint.
• Importance des multipôles : Le quadrupôle du bispectre ($B_2$) est identifié comme un contributeur clé, brisant les dégénérescences entre $f_{\text{NL}}$, $\omega_c$ et les paramètres de biais qui limitent les analyses basées uniquement sur le monopôle.
• Dépendance au décalage vers le rouge : Les contraintes s'améliorent considérablement avec le décalage vers le rouge. Les résultats les plus forts sont obtenus à $z=1,7$, produisant une détection à $1,9\sigma$ pour $f_{\text{NL}} b_\phi$ (agnostique vis-à-vis des a priori) et une détection à $2,35\sigma$ pour $f_{\text{NL}}$ (basée sur les a priori).
• Sensibilité aux grandes échelles : L'analyse souligne l'importance critique des modes à grande échelle ; exclure les premières tranches $k$ dégrade les contraintes sur $f_{\text{NL}}$ et introduit des biais significatifs, en particulier pour le bispectre.

Signification et affirmations
Les auteurs déclarent que ce travail sert de préparation à l'analyse de l'échantillon spectroscopique d'Euclid, en se concentrant sur la validation des systématiques théoriques plutôt que sur la fourniture de prévisions finales pour l'ensemble du relevé. L'article affirme que les analyses conjointes du spectre de puissance et du bispectre offrent de « fortes perspectives pour tester l'inflation multifield », à condition qu'une validation de bout en bout avec les systématiques observationnelles et les géométries de relevé réalistes soit effectuée. Les résultats confirment que, grâce à une modélisation robuste et à des a priori appropriés sur les biais PNG, les relevés de phase IV peuvent contraindre efficacement la PNG locale, potentiellement en dépassant les limites actuelles du CMB. Les auteurs notent explicitement que les contraintes rapportées sont conservatrices, car elles utilisent un sous-ensemble du volume du relevé et n'intègrent pas encore l'ensemble des systématiques observationnelles ou les corrélations croisées avec d'autres traceurs.