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

Technische Zusammenfassung: Optimierung des tomografischen Redshift-Binning für Euclid DR1 3×2pt Power-Spektrum-Constraints

Problemstellung
Die tomografische Analyse des 3×2pt-Signals (die Kombination aus Cosmic Shear, Winkelsystematik der Galaxienverteilung und Galaxy-Galaxy Lensing) erfordert die Unterteilung der beobachteten Galaxienprobe in Redshift-Bins (Rotverschiebungsschub). Während eine Erhöhung der Anzahl der Bins im Allgemeinen mehr kosmologische Informationen erfasst, erhöht dies auch die Rechenkosten und das Poisson-Schrotrauschen (Shot Noise). Vorherige Studien, die die „optimale“ Binning-Strategie für Stage-IV-Durchmusterungen wie Euclid untersuchten, lieferten inkonsistente Ergebnisse hinsichtlich der idealen Anzahl an Bins und der optimalen Methode zur Definition der Bin-Grenzen (z. B. äquipopuliert, gleiche Rotverschiebungsbreite oder gleicher Eigenabstand/Comoving Distance). Diese Arbeit zielt darauf ab, diese Inkonsistenzen durch eine simulationsbasierte Methode zu lösen, um die optimale tomografische Binning-Strategie spezifisch für die erste Euclid-Datenveröffentlichung (DR1) zu bestimmen, mit Fokus auf die Constraints für die Parameter der Dunklen Energie $(w_0, w_a)$.

Methodik
Die Autoren entwickelten eine Simulations-Pipeline (SWEPT), um Mock-Galaxienkataloge zu generieren, die eine realistische Euclid-DR1-Beobachtung nachahmen. Der Arbeitsablauf umfasst:

1. Generierung der Galaxienpopulation: Stichprobenziehung einer Galaxienpopulation aus einer vorhergesagten Euclid-ähnlichen Rotverschiebungsverteilung $p(z)$ und räumlichen Nummerndichte.
2. Berechnung des theoretischen Signals: Berechnung der theoretischen 3×2pt-Power-Spektren ($C^{\gamma\gamma}_{ij}$, $C^{\delta_g\delta_g}_{ij}$, $C^{\delta_g\gamma}_{ij}$) unter Verwendung von CosmoSIS, CAMB und halofit unter einer fiduziellen flachen $w_0w_a$CDM-Kosmologie ($\Omega_m=0.3, H_0=70, w_0=-1, w_a=0$).
3. Map-Generierung: Erstellung von 2D-HEALPix-Maps der korrelierten Cosmic-Shear- und Galaxendichtefelder am Himmel mittels FLASK, welche die volle 3D-kosmologische Information approximieren.
4. Konstruktion des Mock-Katalogs: Poisson-Sampling der Dichtefelder zur Platzierung von Galaxien und Zuweisung von schwachen Lins-Observablen (Konvergenz und Shear) aus den Feld-Maps an den Positionen der Galaxien.
5. Rauschinjektion: Einführung realistischer Beobachtungsunsicherheiten, einschließlich Gauß-verteilter Photometrischer-Rotverschiebungsfehler (photo-z errors) ($\sigma_z = \sigma_{phot}(1+z)$ mit $\sigma_{phot}=0.05$) und Gauß-verteiltem Shape-Rauschen ($\sigma(\epsilon_{int})=0.3$). Eine Teilmenge der Simulationen beinhaltete zudem 5 % katastrophale photo-z Outlier.
6. Analyse: Messung der Pseudo-$C_\ell$-Power-Spektren für verschiedene tomografische Konfigurationen (variable Bin-Anzahl und Grenzdefinitionen) mittels NaMaster. Die Constraints auf $(w_0, w_a)$ wurden über eine gitterbasierte Gaußsche Likelihood-Analyse mit einer numerischen Kovarianzmatrix abgeleitet, die aus 3.000 Realisierungen gewonnen und mittels des Hartlap-Faktors korrigiert wurde.

Kernergebnisse
Die Studie evaluierte drei Binning-Strategien: äquipopulierte Bins, in Rotverschiebung gleichmäßig beabstandete Bins und in den fiduziellen Eigenabstand (comoving distance) gleichmäßig beabstandete Bins.

• Optimale Binning-Strategie:

  • Volles 3×2pt und Winkel-Clustering: Sowohl für das volle 3×2pt-Signal als auch für die Komponente des Winkel-Clusterings liefern äquipopulierte Bins die engsten Constraints auf $(w_0, w_a)$. Dies gilt sowohl für rauschfreie als auch für realistische verrauschte Simulationen.
  • Nur Cosmic Shear: Die optimale Strategie hängt vom Rauschmodell ab. In rauschfreien Simulationen liefern gleichmäßige Rotverschiebungsbreiten die besten Constraints. In realistischen Simulationen, die photo-z- und Shape-Rauschen enthalten, liefern jedoch im fiduziellen Eigenabstand gleichmäßig beabstandete Bins die besten Constraints, wobei der Unterschied zwischen den Binning-Wahlen minimal ist (< 5 %).
  • Dominanz des Clusterings: Die Bevorzugung von äquipopulierten Bins in der vollen 3×2pt-Analyse wird primär durch die Winkel-Clustering-Komponente getrieben, welche die Constraint-Stärke des kombinierten Signals dominiert.

• Anzahl der Bins und Informationssättigung:

  • Der Informationsgewinn sättigt jenseits von sechs bis acht Bins für alle Signal-Komponenten (Shear, Clustering und volles 3×2pt).
  • Eine Erhöhung der Bin-Anzahl von 7 auf 10 liefert nur marginale Gewinne (~1–2 % Verbesserung).
  • Eine Erweiterung der Analyse auf 13 Bins (eine in anderer Literatur vorgeschlagene Konfiguration) führt im Vergleich zu 7 Bins nur zu einer weiteren Verbesserung der Präzision um etwa 5 %. Die Autoren merken an, dass die marginalen Gewinne durch zusätzliche Bins wahrscheinlich durch die steigenden Anforderungen an die Genauigkeit der Kovarianzmatrix-Schätzung und das Vorhandensein anderer Systematiken in realen Daten begrenzt sind.

• Auswirkung katastrophaler Photo-z-Fehler:

  • Die Einführung von 5 % katastrophalen photo-z Outliern ohne Mitigierung führt zu schweren Bias (Verzerrungen).
  • Für Analysen mit 5 oder 10 Bins sind die resultierenden Constraints auf $(w_0, w_a)$ inkonsistent mit der fiduziellen $\Lambda$CDM-Kosmologie bei $\gtrsim 3\sigma$ und favorisieren ein zeitentwickelndes Dunkle-Energie-Modell ($w_0 > -1, w_a < 0$).
  • Die Posterior-Verteilungen zeigen scharfe, nicht-elliptische Merkmale, was darauf hindeutet, dass unmitigierte katastrophale Fehler komplexe, nicht-triviale Biases in die Power-Spektren induzieren.

Bedeutung und Ansprüche
Das Paper beansprucht, einen robusten, simulationsbasierten Rahmen bereitzustellen, der eine reale Beobachtung eines Surveys nachahmt und somit eine direktere Bewertung der tomografischen Binning-Strategien ermöglicht als bisherige theoretische Studien. Die primären Schlussfolgerungen sind:

1. Äquipopulierte Bins sind die empfohlene Strategie für Euclid DR1-Analysen, die das volle 3×2pt-Signal oder das Winkel-Clustering beinhalten, da sie die robustesten Constraints liefern.
2. Die Cosmic-Shear-Komponente ist relativ unempfindlich gegenüber der spezifischen Binning-Methodik (Unterschiede < 5 % in verrauschten Szenarien), was darauf hindeutet, dass die Wahl des Binning für Shear allein weniger kritisch ist als für das Clustering.
3. Es gibt einen abnehmenden Grenznutzen bei der Erhöhung der Bin-Anzahl über sieben oder acht hinaus; die marginalen statistischen Gewinne werden wahrscheinlich durch die Rechenkosten und die erschwerte genaue Modellierung von Systematiken und Kovarianzmatrizen übertroffen.
4. Katastrophale photo-z-Fehler stellen eine kritische Bedrohung für die kosmologische Inferenz dar. Wenn sie unmitigiert bleiben, können sie Biases induzieren, die eine Abweichung von $\Lambda$CDM mit hoher Signifikanz vortäuschen, was eine rigorose Modellierung oder Mitigierungsstrategien im Analyse-Pipeline erforderlich macht.

Die Autoren betonen, dass obwohl ihre Ergebnisse nahelegen, dass äquipopulierte Bins für die spezifischen Charakteristika von Euclid DR1 optimal sind, die relative Unempfindlichkeit des Cosmic Shear gegenüber Binning-Wahlen und die Sättigung des Informationsgewinns darauf hindeuten, dass die exakte Binning-Strategie weniger kritisch sein kann als die genaue Modellierung von Systematiken und der Kovarianzmatrix.