Improved recipes for peculiar velocity power spectra using Evolution Mapping

Technische Zusammenfassung: Verbesserte Rezepte für Leistungsspektren peculiarer Geschwindigkeiten mittels Evolution Mapping

Problemstellung
Aktuelle und kommende Galaxien-Surveys (z. B. DES, DESI, Euclid, LSST) erfordern präzise theoretische Modelle der nichtlinearen Strukturbildung, um kosmologische Informationen aus Rotverschiebungsraum-Verzerrungen (RSD) zu extrahieren. Eine kritische Komponente der RSD-Modellierung ist das Geschwindigkeitsdivergenzfeld ($\theta = -\nabla \cdot \mathbf{v}$), insbesondere sein Autoleistungsspektrum $P_{\theta\theta}(k)$ und sein Kreuzleistungsspektrum mit dem Dichtefeld $P_{\delta\theta}(k)$.

Bestehende Anpassungsvorschriften, insbesondere die von Bel et al. (2019, im Folgenden B19), leiden unter drei primären Einschränkungen:

1. Numerische Konvergenz: Frühere Kalibrierungen beruhten auf Simulationsauflösungen, bei denen Sampling-Artefakte (die aus der Schätzung von Geschwindigkeitsfeldern aus diskreten Partikelverteilungen resultieren) Messungen bei Skalen von bis zu $k \sim 0,2\text{--}0,3 \, h\text{Mpc}^{-1}$ unzuverlässig machten.
2. Spurlose $h$-Abhängigkeit: B19 parametrisieren die nichtlineare Evolution mittels $\sigma_8$ (RMS-Dichtefluktuation, geglättet auf $8 \, h^{-1}\text{Mpc}$). Da die Glättungsskala vom Hubble-Parameter $h$ abhängt, führt dies zu unphysikalischen Abhängigkeiten von $h$. Folglich liefern Modelle mit identischer linearer Physik, aber unterschiedlichen $h$-Werten divergierende Vorhersagen, was die Genauigkeit bei der Extrapolation über den kosmologischen Parameterraum hinweg beeinträchtigt.
3. Funktionale Form: Bestehende Modelle erfassen Modenkopplungseffekte oder die nichtlineare Dämpfung von Baryonischen Akustischen Oszillationen (BAO) ohne zusätzliche freie Parameter nicht vollständig.

Methodik
Die Autoren stellen neue Anpassungsfunktionen vor, die an einer Suite von rein gravitativen $N$-Körper-Simulationen (Aletheia, AletheiaMass und AletheiaEmu) unter Verwendung des Evolution Mapping-Rahmens kalibriert wurden. Dieses Rahmenwerk geht davon aus, dass, wenn Formparameter ($\Theta_s$) fixiert sind, der Einfluss von Evolutionsparametern ($\Theta_e$) und der Rotverschiebung in einen einzigen Clustering-Amplitudenparameter, $\sigma_{12}$ (RMS-Dichtefluktuation, geglättet auf 12 Mpc), absorbiert werden kann, der unabhängig von $h$ ist.

Wichtige methodische Schritte umfassen:

• Konvergenzstrategie: Die Autoren nutzen die AletheiaMass-Suite, die Multiresolution-Simulationen (mit variierenden Boxgrößen und Partikelzahlen) beinhaltet. Sie verwenden eine „Patching"-Strategie, bei der Messungen aus verschiedenen Volumina zusammengesetzt werden. Eine Konvergenzskala $k_{\text{crop}}$ wird für jede Box definiert, indem Hoch- und Niederauflösungsläufe verglichen werden; das finale Leistungsspektrum verwendet die Messung aus dem größten Volumen, das $k < k_{\text{crop}}$ erfüllt. Dies liefert konvergierte Messungen bis zu $k \simeq 0,56 \, \text{Mpc}^{-1}$.
• Geschwindigkeitsschätzung: Um Sampling-Artefakte zu mindern, verwenden die Autoren die MC-Voronoi-Methode. Dies beinhaltet eine Neuabtastung des Geschwindigkeitsfeldes auf eine gleichmäßige glasartige Punktverteilung (volumengewichtet) basierend auf dem nächsten $N$-Körper-Partikel, gefolgt von einer Cloud-in-Cell (CIC)-Zuweisung und einer Dekonvolution.
• Funktionale Formen:
  • $P_{\theta\theta}(k)$: Die Autoren schlagen einen Ansatz vor, der einen gedämpften linearen Term (unter Verwendung eines entwackelten linearen Leistungsspektrums $P_{dw}$ zur Berücksichtigung der nichtlinearen BAO-Dämpfung) und einen additiven Gaußschen Term in $\ln k$ kombiniert, um Modenkopplung zu erfassen. Die Parameter sind Funktionen von $\sigma_{12}$.
  • $P_{\delta\theta}(k)$: Dies wird als geometrisches Mittel des nichtlinearen Dichtespektrums $P_{\delta\delta}(k)$ (vorhergesagt vom Aletheia-Emulator) und des neuen $P_{\theta\theta}(k)$ modelliert, multipliziert mit einem exponentiellen Dämpfungsterm, der von $\sigma_{12}$ abhängt.
• Kalibrierung: Parameter werden unter Verwendung eines Bayesschen Rahmens (nautilus) gegen die AletheiaMass-Suite angepasst, die eine breite Palette von $\sigma_{12}$ (0,2 bis 1,0) und Rotverschiebungen ($z \approx -0,4$ bis $4,2$) abdeckt.

Hauptbeiträge

1. Neue Anpassungsfunktionen: Die Arbeit liefert kalibrierte Formeln für $P_{\theta\theta}(k)$ und $P_{\delta\theta}(k)$, die explizit von $\sigma_{12}$ abhängen und physikalische Einheiten (Mpc$^{-1}$) verwenden, wodurch spurlose $h$-Abhängigkeiten eliminiert werden.
2. Robuste Konvergenz: Durch die Implementierung einer konservativen Patching-Strategie mit Multiresolution-Simulationen erreichen die Autoren eine Genauigkeit im Prozentbereich bis zu $k \simeq 0,56 \, \text{Mpc}^{-1}$ und erweitern damit den zuverlässigen Bereich erheblich im Vergleich zu früheren Arbeiten.
3. Validierung von Evolution Mapping: Die Studie bestätigt, dass die Darstellung von Geschwindigkeitsstatistiken in Bezug auf $\sigma_{12}$ die Evolution-Mapping-Beziehung erhält, sodass Kosmologien mit identischer linearer Physik und $\sigma_{12}$ konsistente nichtlineare Vorhersagen liefern, unabhängig von ihrer Wachstumsgeschichte oder ihrem $h$-Wert.

Ergebnisse

• Genauigkeit: Die neuen Vorschriften erreichen 1–2 % Genauigkeit sowohl für $P_{\theta\theta}(k)$ als auch für $P_{\delta\theta}(k)$ über einen weiten Bereich von $\sigma_{12}$ und Skalen bis zu $k \simeq 0,56 \, \text{Mpc}^{-1}$ innerhalb des Kalibrierungssatzes.
• Vergleich mit B19: Die neuen Modelle übertreffen B19 systematisch. B19-Vorhersagen weichen bei kleinen Skalen ($k > 0,2 \, \text{Mpc}^{-1}$) aufgrund von Auflösungsgrenzen signifikant ab und zeigen unphysikalische Verschiebungen für Kosmologien mit unterschiedlichen $h$-Werten. Die neuen Modelle sagen Geschwindigkeitsstatistiken über verschiedene Hubble-Parameter hinweg korrekt voraus.
• Effekte der Wachstumsgeschichte: Unter Verwendung der Aletheia-Suite (die Effekte der Wachstumsgeschichte isoliert, indem Formparameter fixiert, aber Evolutionsparameter variiert werden), quantifizieren die Autoren Abweichungen von der perfekten Evolution-Mapping. Sie stellen fest, dass für die meisten von praktischem Interesse liegenden Kosmologien das Vernachlässigen von Unterschieden in der Wachstumsgeschichte Fehler unter 1 % einführt. Signifikante Abweichungen ($>1\%$) treten nur für extreme, durch Beobachtungen ausgeschlossene Kosmologien auf (z. B. Einstein-de Sitter oder extreme $h$-Werte).
• Kreuzspektrum: Das $P_{\delta\theta}(k)$-Modell profitiert vom Aletheia-Emulator für $P_{\delta\delta}(k)$, der implizit einige Abhängigkeiten von der Wachstumsgeschichte erfasst und die Fehler weiter reduziert.

Bedeutung
Die Arbeit behauptet, dass die Wahl der Parametrisierung und der Einheiten für die Robustheit von Modellen für Geschwindigkeitsstatistiken entscheidend ist. Indem sie nachweisen, dass $h$-abhängige Einheiten (wie $\sigma_8$) systematische Verzerrungen einführen, liefern die Autoren starke empirische Unterstützung für die Verwendung physikalischer Einheiten und $\sigma_{12}$ als primären Skalierungsparameter. Diese Anpassungsfunktionen bieten eine robuste, auf Prozentniveau genaue Beschreibung von Geschwindigkeitsleistungsspektren, die direkt auf RSD-Modellierung in aktuellen und zukünftigen Galaxien-Rotverschiebungssurveys anwendbar ist. Die Autoren weisen darauf hin, dass ihre aktuellen Anpassungsfunktionen zwar für praktische Anwendungen ausreichend sind, zukünftige Arbeiten jedoch die Konstruktion eines dedizierten Geschwindigkeits-Emulators beinhalten könnten, um verbleibende Abweichungen der Wachstumsgeschichte explizit zu modellieren und eine höhere Präzision über den gesamten kosmologischen Parameterraum zu erreichen.