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

Technische Zusammenfassung: Testen von Mehr-Felder-Inflation mit dem Galaxien-Leistungsspektrum und der Bispektrum

Problemstellung
Primordiale Nicht-Gaußsche Verteilung (PNG), insbesondere vom lokalen Typ, dient als kritische Sonde zur Unterscheidung zwischen Ein-Felder- und Mehr-Felder-Modellen der Inflation. Während der kosmische Mikrowellenhintergrund (CMB) enge Grenzen für lokale PNG gesetzt hat (nahezu kosmische-Varianz-begrenzt), versprechen Galaxien-Übersichten der Stufe IV wie Euclid, diese Grenzen zu übertreffen, indem sie 2-Punkt- (Leistungsspektrum) und 3-Punkt- (Bispektrum) Klusterungsstatistiken über riesige kosmologische Volumina messen. Die Extraktion robuster Einschränkungen für den PNG-Parameter $f_{\text{NL}}$ aus der Galaxienklusterung erfordert jedoch ein rigoroses gemeinsames Analyse-Framework, das Rotverschiebungsraum-Verzerrungen (RSD), nichtlineare gravitative Evolution und Galaxien-Bias berücksichtigt. Eine Haupt Herausforderung besteht in der Entartung zwischen $f_{\text{NL}}$ und dem PNG-induzierten Bias-Parameter $b_\phi$, die ohne informative Priors oder eine umfassende Behandlung theoretischer Systematiken schwer zu brechen ist.

Methodik
Die Autoren validieren ein theoretisches Modellierungs- und Inferenz-Pipeline, das für die spektroskopische H$\alpha$-Galaxienstichprobe von Euclid entwickelt wurde. Die Arbeit nutzt:

• Simulationen: Eine Suite von Abacus-PNG N-Körper-Simulationen mit lokalen PNG-Anfangsbedingungen ($f_{\text{NL}} = \{0, 30\}$) und entsprechenden gaußschen Gegenstücken. Diese werden mit Galaxien besetzt, wobei eine Halo-Besetzungsverteilung (HOD) verwendet wird, die auf die Reproduktion der Klustereigenschaften der Euclid Flagship 2 (FS2)-Simulation abgestimmt ist. Ein Abundanz-Matching (AM)-Verfahren wird angewendet, um Unterschiede in den Halo-Massen-Definitionen zwischen den FS2- und Abacus-Simulationen auszugleichen.
• Theoretische Modellierung: Die Analyse verwendet die Effektivfeldtheorie der großskaligen Struktur (EFTofLSS). Das Leistungsspektrum im Rotverschiebungsraum wird im Ein-Schleifen-Ordnung modelliert, während das Bispektrum im Baum-Level modelliert wird. Das Framework umfasst:
  • Eine renormierte Bias-Expansion, die um PNG-Operatoren erweitert wurde.
  • Rotverschiebungsraum-Verzerrungen, die über eine Geschwindigkeits-Moment-Expansion modelliert werden.
  • Gegenterme zur Absorption ultravioletter Sensitivitäten und zur Modellierung von Geschwindigkeitseffekten auf kleinen Skalen (Finger of God), wobei speziell ein Gegenterm übernommen wird, der vom Geschwindigkeitsdifferenz-Generator (VDG)-Framework für das Bispektrum motiviert ist.
  • Infrarot-Resummation und Alcock–Paczynski-Effekte.
• Inferenz-Pipeline: Ein neuer Code, CosmoSIS-gClust, wird eingeführt, um gemeinsame Likelihood-Analysen der Leistungsspektrum-Multipole ($P_0, P_2, P_4$) und der Bispektrum-Multipole ($B_0, B_2$) durchzuführen. Die Analyse marginalisiert über Störparameter (Gegenterme, Schussrauschen) und testet verschiedene Priors für den PNG-Bias-Parameter $b_\phi$.
• Metriken: Die Gültigkeit des Modells wird mit der Bias-Figur (FoB) zur Überprüfung auf systematische Fehler und mit der Merit-Figur (FoM) zur Quantifizierung des Informationsgehalts bewertet. Skalenabschnitte ($k_{\text{max}}$) werden optimiert, um eine unverzerrte Wiederherstellung kosmologischer Parameter sicherzustellen.

Hauptbeiträge

1. Pipeline-Validierung: Das Papier liefert eine gründliche Validierung einer gemeinsamen Analyse-Pipeline für Leistungsspektrum und Bispektrum für lokale PNG im Rotverschiebungsraum, die speziell für Euclid-ähnliche Daten zugeschnitten ist.
2. Bias-Prior-Strategie: Die Autoren zeigen, dass ein „prior-agnostischer" Ansatz (unter Verwendung breiter flacher Priors für $b_\phi$) zu verzerrten oder unbeschränkten $f_{\text{NL}}$-Schätzungen führt, bedingt durch die starke Entartung mit $b_\phi$. Sie schlagen einen physikalisch motivierten Gaußschen Prior für den Parameter $p_{b_\phi}$ vor (wobei $b_\phi = 2\delta_c(b_1 - p_{b_\phi})$ gilt), der sowohl Vorhersagen der universellen Halo-Massenfunktion (UHMF) als auch des Galaxienbildungsmodells (GFM) einschließt. Diese Wahl liefert unverzerrte Einschränkungen unter Berücksichtigung theoretischer Unsicherheiten.
3. Optimierung der Skalenabschnitte: Die Studie bestimmt optimale Skalenabschnitte ($k_{\text{max}}$) für beide Statistiken. Für das Leistungsspektrum bleiben Ein-Schleifen-Modelle bis $k \approx 0,30 \, h \text{ Mpc}^{-1}$ unverzerrt. Für das Bispektrum bleiben Baum-Level-Modelle bis $k \approx 0,12 \, h \text{ Mpc}^{-1}$ unverzerrt (bei höheren Rotverschiebungen niedriger), wobei darüber hinaus das Fehlen von Schleifenkorrekturen zu Verzerrungen führt.
4. Null-Test: Ein Null-Test an Simulationen mit gaußschen Anfangsbedingungen ($f_{\text{NL}}=0$) bestätigt, dass die Pipeline PNG nicht fälschlicherweise detektiert oder die Standard-$\Lambda$CDM-Parameter verzerrt.

Ergebnisse

• Detektionsbedeutung: In einem „prior-agnostischen" Setup ist eine hochsignifikante Detektion des dominanten PNG-Terms ($f_{\text{NL}} b_\phi$) in einzelnen Rotverschiebungsbins schwierig. Mit dem vorgeschlagenen informativen Prior für $b_\phi$ gewinnt die gemeinsame Analyse jedoch den fiduziellen Wert $f_{\text{NL}} = 30$ mit weniger als $1\sigma$ Verzerrung zurück.
• Informationsgewinn:
  • Bei fester Rotverschiebung und festen Skalenabschnitten reduziert das Bispektrum allein den Fehler für $f_{\text{NL}}$ ($\sigma(f_{\text{NL}})$) um etwa 29–46 % im Vergleich zum Leistungsspektrum allein.
  • Die gemeinsame Leistungsspektrum-Bispektrum-Analyse schränkt die Grenzen weiter um zusätzliche 8–13 % im Vergleich zum Bispektrum allein ein.
  • Der kumulative Gewinn an Präzision von $z=0,8$ bis $z=1,7$ beträgt für den gemeinsamen Fall einen Faktor von $\sim 2,3$.
• Bedeutung der Multipole: Das Bispektrum-Quadrupol ($B_2$) wird als wesentlicher Beitrag identifiziert, der Entartungen zwischen $f_{\text{NL}}$, $\omega_c$ und Bias-Parametern aufbricht, die Analysen nur mit dem Monopol begrenzen.
• Rotverschiebungsabhängigkeit: Die Einschränkungen verbessern sich signifikant mit der Rotverschiebung. Die stärksten Ergebnisse werden bei $z=1,7$ erzielt, was eine $1,9\sigma$-Detektion für $f_{\text{NL}} b_\phi$ (prior-agnostisch) und eine $2,35\sigma$-Detektion für $f_{\text{NL}}$ (prior-basiert) ergibt.
• Empfindlichkeit auf großen Skalen: Die Analyse unterstreicht die kritische Bedeutung von Moden auf großen Skalen; das Ausschließen der ersten paar $k$-Bins verschlechtert die Einschränkungen für $f_{\text{NL}}$ und führt zu erheblichen Verzerrungen, insbesondere für das Bispektrum.

Bedeutung und Behauptungen
Die Autoren stellen fest, dass diese Arbeit als Vorbereitung für die Analyse der spektroskopischen Euclid-Stichprobe dient und sich auf die Validierung theoretischer Systematiken konzentriert, anstatt endgültige Vorhersagen für die gesamte Übersicht zu liefern. Das Papier behauptet, dass gemeinsame Analysen von Leistungsspektrum und Bispektrum „gute Aussichten für das Testen von Mehr-Felder-Inflation" bieten, vorausgesetzt, eine End-zu-End-Validierung mit beobachtungsbezogenen Systematiken und realistischen Übersichtsgeometrien wird durchgeführt. Die Ergebnisse bestätigen, dass Übersichten der Stufe IV mit robuster Modellierung und angemessenen Priors für PNG-Biases lokale PNG effektiv einschränken können und möglicherweise aktuelle CMB-Grenzen übertreffen. Die Autoren weisen ausdrücklich darauf hin, dass die berichteten Einschränkungen konservativ sind, da sie eine Teilmenge des Übersichtsvolumens nutzen und noch keine vollständigen beobachtungsbezogenen Systematiken oder Kreuzkorrelationen mit anderen Tracern einbeziehen.