Improved recipes for peculiar velocity power spectra using Evolution Mapping

Technische Samenvatting: Verbeterde recepten voor vermogensspectra van peculiar velocity met behulp van Evolution Mapping

Probleemstelling
Huidige en toekomstige sterrenstelselonderzoeken (bijv. DES, DESI, Euclid, LSST) vereisen nauwkeurige theoretische modellen van niet-lineaire structuurvorming om kosmologische informatie te extraheren uit roodverschuivingsruimte-distorsies (RSD). Een kritiek onderdeel van RSD-modellering is het divergentieveld van de snelheid ($\theta = -\nabla \cdot \mathbf{v}$), specifiek zijn autopowervermogen $P_{\theta\theta}(k)$ en zijn kruisvermogen met het dichtheidsveld $P_{\delta\theta}(k)$.

Bestaande aanpassingsvoorschriften, met name die van Bel et al. (2019, hierna B19), lijden onder drie primaire beperkingen:

1. Numerieke Convergentie: Vorige kalibraties maakten gebruik van simulatie-resoluties waarbij bemonsteringsartefacten (ontstaand uit het schatten van snelheidsvelden uit discrete deeltjesverdelingen) metingen onbetrouwbaar maakten op schalen zo groot als $k \sim 0.2\text{--}0.3 \, h\text{Mpc}^{-1}$.
2. Schijnbare $h$-afhankelijkheid: B19 parametriseren niet-lineaire evolutie met $\sigma_8$ (RMS-dichtheidsfluctuatie gesmeerd op $8 \, h^{-1}\text{Mpc}$). Omdat de gladmakende schaal afhangt van de Hubble-parameter $h$, introduceert dit onfysische afhankelijkheden van $h$. Bijgevolg leveren modellen met identieke lineaire fysica maar verschillende $h$-waarden divergerende voorspellingen op, wat de nauwkeurigheid ondermijnt bij extrapolatie over de ruimte van kosmologische parameters.
3. Functionele Vorm: Bestaande modellen vangen mode-koppelingseffecten of de niet-lineaire demping van Baryon Acoustic Oscillations (BAO) niet volledig op zonder extra vrije parameters.

Methodologie
De auteurs presenteren nieuwe aanpassingsfuncties die zijn gekalibreerd op een suite van zwaartekracht-only $N$-body-simulaties (Aletheia, AletheiaMass en AletheiaEmu) met behulp van het Evolution Mapping-kader. Dit kader stelt dat wanneer vormparameters ($\Theta_s$) vaststaan, de impact van evolutieparameters ($\Theta_e$) en roodverschuiving kan worden opgenomen in een enkele clustering-amplitudparameter, $\sigma_{12}$ (RMS-dichtheidsfluctuatie gesmeerd op 12 Mpc), die onafhankelijk is van $h$.

Belangrijke methodologische stappen omvatten:

• Convergentiestrategie: De auteurs maken gebruik van de AletheiaMass-suite, die multi-resolutiesimulaties bevat (variërende doosgroottes en deeltjesaantallen). Zij hanteren een "patching"-strategie waarbij metingen uit verschillende volumes aan elkaar worden genaaid. Een convergentieschaal $k_{\text{crop}}$ wordt gedefinieerd voor elke doos door hoge- en laag-resolutie runs te vergelijken; het uiteindelijke vermogensspectrum maakt gebruik van de meting uit het grootste volume dat voldoet aan $k < k_{\text{crop}}$. Dit levert geconvergeerde metingen op tot $k \simeq 0.56 \, \text{Mpc}^{-1}$.
• Snelheidsschatting: Om bemonsteringsartefacten te mitigeren, gebruiken de auteurs de MC-Voronoi-methode. Dit houdt in dat het snelheidsveld wordt herschikt op een uniforme glasachtige verdeling van punten (volumegewogen) gebaseerd op de dichtstbijzijnde $N$-body-deeltjes, gevolgd door Cloud-in-Cell (CIC)-toewijzing en deconvolutie.
• Functionele Vormen:
  • $P_{\theta\theta}(k)$: De auteurs stellen een ansatz voor die een gedempt lineair term combineert (met behulp van een ontgolveld lineair vermogensspectrum $P_{dw}$ om rekening te houden met niet-lineaire BAO-demping) en een additieve Gaussische term in $\ln k$ om mode-koppeling te vangen. De parameters zijn functies van $\sigma_{12}$.
  • $P_{\delta\theta}(k)$: Dit wordt gemodelleerd als het geometrische gemiddelde van het niet-lineaire dichtheidspectrum $P_{\delta\delta}(k)$ (voorspeld door de Aletheia-emulator) en het nieuwe $P_{\theta\theta}(k)$, vermenigvuldigd met een exponentiële dempingsterm afhankelijk van $\sigma_{12}$.
• Kalibratie: Parameters worden gefit met een Bayesiaans kader (nautilus) tegen de AletheiaMass-suite, die een breed scala aan $\sigma_{12}$ (0.2 tot 1.0) en roodverschuivingen ($z \approx -0.4$ tot 4.2) bestrijkt.

Belangrijkste bijdragen

1. Nieuwe Aanpassingsfuncties: Het artikel levert gekalibreerde formules voor $P_{\theta\theta}(k)$ en $P_{\delta\theta}(k)$ die expliciet afhankelijk zijn van $\sigma_{12}$ en fysische eenheden (Mpc$^{-1}$) gebruiken, waardoor schijnbare $h$-afhankelijkheden worden geëlimineerd.
2. Robuuste Convergentie: Door een conservatieve patching-strategie met multi-resolutiesimulaties te implementeren, bereiken de auteurs een nauwkeurigheid op procentniveau tot $k \simeq 0.56 \, \text{Mpc}^{-1}$, wat het betrouwbare bereik ten opzichte van eerdere werken aanzienlijk uitbreidt.
3. Validatie van Evolution Mapping: De studie bevestigt dat het uitdrukken van snelheidsstatistieken in termen van $\sigma_{12}$ de Evolution Mapping-relatie behoudt, waardoor ervoor wordt gezorgd dat kosmologieën met identieke lineaire fysica en $\sigma_{12}$ consistente niet-lineaire voorspellingen opleveren, ongeacht hun groeigeschiedenis of $h$-waarde.

Resultaten

• Nauwkeurigheid: De nieuwe voorschriften bereiken 1–2% nauwkeurigheid voor zowel $P_{\theta\theta}(k)$ als $P_{\delta\theta}(k)$ over een breed scala aan $\sigma_{12}$ en schalen tot $k \simeq 0.56 \, \text{Mpc}^{-1}$ binnen de kalibratieset.
• Vergelijking met B19: De nieuwe modellen presteren systematisch beter dan B19. B19-voorspellingen wijken significant af op kleine schalen ($k > 0.2 \, \text{Mpc}^{-1}$) vanwege resolutielimieten en tonen onfysische verschuivingen voor kosmologieën met verschillende $h$-waarden. De nieuwe modellen voorspellen correct snelheidsstatistieken over verschillende Hubble-parameters.
• Effecten van Groeigeschiedenis: Met behulp van de Aletheia-suite (die effecten van groeigeschiedenis isoleert door vormparameters vast te houden maar evolutieparameters te variëren), kwantificeren de auteurs afwijkingen van perfecte Evolution Mapping. Zij vinden dat voor de meeste kosmologieën van praktisch belang het negeren van verschillen in groeigeschiedenis fouten onder de 1% introduceert. Significante afwijkingen ($>1\%$) treden alleen op voor extreme, observationeel uitgesloten kosmologieën (bijv. Einstein-de Sitter of extreme $h$-waarden).
• Kruisspectrum: Het $P_{\delta\theta}(k)$-model profiteert van de Aletheia-emulator voor $P_{\delta\delta}(k)$, wat impliciet enige afhankelijkheid van groeigeschiedenis vastlegt, waardoor fouten verder worden gereduceerd.

Betekenis
Het artikel beweert dat de keuze van parametrisatie en eenheden cruciaal is voor de robuustheid van modellen voor snelheidsstatistieken. Door aan te tonen dat $h$-afhankelijke eenheden (zoals $\sigma_8$) systematische vertekeningen introduceren, bieden de auteurs sterke empirische ondersteuning voor het gebruik van fysische eenheden en $\sigma_{12}$ als primaire schalingsparameter. Deze aanpassingsfuncties bieden een robuuste, nauwkeurige beschrijving op procentniveau van vermogensspectra van snelheid, direct toepasbaar op RSD-modellering in huidige en toekomstige sterrenstelsel-roodverschuivingsonderzoeken. De auteurs merken op dat hoewel hun huidige aanpassingsfuncties voldoende zijn voor praktische toepassingen, toekomstig werk kan inhouden dat er een dedicated snelheidsemulator wordt geconstrueerd om residuele afwijkingen in groeigeschiedenis expliciet te modelleren en hogere precisie te bereiken over de volledige ruimte van kosmologische parameters.