Towards an application of fourth-order shear statistics I. The information content of $\langle M_\mathrm{ap}^4 \rangle$

Probleemstelling

Cosmische schuif (cosmic shear) is een krachtige tool om de grote schaalstructuur van het heelal en kosmologische parameters te bestuderen. Traditioneel zijn tweede-orde statistieken (zoals de schuif-correlatiefunctie) de standaard, omdat ze de Gaussische informatie in het veld vastleggen. Echter, naarmate de waarnemingsprecisie toeneemt (bijv. bij Stage III en IV surveys zoals DES-Y3 en Euclid), wordt de informatie in niet-Gaussische fluctuaties steeds belangrijker.

Hoewel derde-orde statistieken al succesvol worden toegepast, is de toepassing van vierde-orde schuifstatistieken (vierpunt-correlatiefuncties of 4PCF) nog beperkt door theoretische complexiteit en computatievere eisen. Dit artikel richt zich op het ontwikkelen van een theoretisch raamwerk en een numerieke pijplijn om vierde-orde aperture-statistieken (⟨𝑀⁴ₐₚ⟩) te berekenen en te evalueren, met als doel te bepalen welke extra informatie deze statistieken bieden ten opzichte van tweede- en derde-orde statistieken.

Methodologie

De auteurs hanteren een aanpak die theoretische afleidingen combineert met numerieke integratie en forecast-analyses:

1. Theoretisch Raamwerk:

   • De relatie wordt gelegd tussen de $n$-de orde aperture-statistieken (⟨𝑀ⁿₐₚ⟩), de $n$-punt correlatiefuncties (nPCF) van het schuifveld, en de $n$-punt polyspectra (zoals het trispectrum voor $n=4$).
   • Specifiek voor de vierde orde worden de filterfuncties afgeleid die nodig zijn om de 4PCF te integreren naar aperture-statistieken. Dit gebeurt via een projectie op een "×-projectie" (adaptief aan Porth et al. 2024) om coördinaat-singulariteiten te vermijden.
   • De auteurs onderscheiden tussen de verbonden (connected) en ongebonden (disconnected) delen van de statistieken. Voor een Gaussisch willekeurig veld (GRF) bestaat de vierde-orde statistiek uitsluitend uit ongebonden termen die afgeleid kunnen worden uit tweede-orde statistieken.

2. Numerieke Integratiepijplijn:

   • Een numerieke integratie wordt uitgevoerd over de 4PCF van een model voor een Gaussisch willekeurig schuifveld (GRF).
   • De integratie wordt uitgevoerd in een 5-dimensionale ruimte (naast de oorspronkelijke 8 dimensies) door gebruik te maken van symmetrieën en het integreren over de positie van het zwaartepunt van het viereck (CMV) en de oriëntatie.
   • Er wordt aandacht besteed aan de binningsgeometrie (binning accuracy). De auteurs testen verschillende radiale en hoekige binningsstrategieën (logaritmisch gespatieerd) om numerieke fouten en E/B-mode mixing te minimaliseren.
   • De nauwkeurigheid wordt getest door de berekende vierde-orde statistieken te vergelijken met de theoretische verwachtingen voor een GRF (waarbij ⟨𝑀⁴ₐₚ⟩ = 3⟨𝑀²ₐₚ⟩²).

3. Fisher Forecast:

   • Een Fisher-informatiematrix-analyse wordt uitgevoerd om het constrainerend vermogen te voorspellen voor een opstelling vergelijkbaar met DES-Y3.
   • De observabelen zijn gelijke-schaal aperture-statistieken voor de tweede, derde en vierde orde.
   • Een nieuwe, robuuste methode wordt ontwikkeld om afgeleiden van de observabelen naar kosmologische parameters ($\Omega_m, S_8, h, w_0$) te schatten vanuit een simulatiesuite (cosmo-SLICS) met willekeurig verdeelde kosmologische sets, in plaats van gebruik te maken van onstabiele eindige-differentie-methoden.

Belangrijkste Bijdragen

• Theoretische Afleiding: De eerste expliciete afleiding van de filterfuncties voor vierde-orde aperture-statistieken in termen van de 4PCF en het trispectrum, inclusief de decompositie in E- en B-modi.
• Numerieke Pijplijn: De ontwikkeling van een geoptimaliseerde integratiepijplijn die numerieke bias (voornamelijk door eindige integratiegebieden en binningsfouten) tot een minimum beperkt.
• Nieuwe Schattingstechniek: Een methode om afgeleiden voor Fisher-forecasts te berekenen op basis van simulaties met onregelmatige sampling (Latin Hypercube), wat essentieel is voor het omgaan met ruis in simulaties.
• Validatie: Een grondige validatie van de methode op een GRF, waarbij wordt aangetoond dat de berekende statistieken binnen de verwachte foutmarges vallen.

Resultaten

1. Numerieke Nauwkeurigheid:

   • Door numerieke effecten te mitigeren (zoals E/B-mode mixing door onvoldoende integratiebereik en binningsfouten), bereiken de auteurs een precisie van twee procent op de vierde-orde aperture-statistieken voor een GRF.
   • Deze precisie ligt ver onder het ruisbudget van Stage IV surveys, wat aangeeft dat de methode robuust is voor toekomstige toepassingen.
   • Het gebruik van "hypercubische" bins en het integreren van de 4PCF binnen bredere radiale bins (in plaats van op vaste punten) vermindert de dispersie aanzienlijk.

2. Fisher Forecast (DES-Y3):

   • Voor een opstelling vergelijkbaar met DES-Y3 (niet-tomografisch, gelijke schaal) wordt geen significante verbetering gevonden in het constrainerend vermogen voor kosmologische parameters wanneer vierde-orde statistieken worden toegevoegd aan een gezamenlijke analyse van tweede- en derde-orde statistieken.
   • De beperkte winst wordt toegeschreven aan een aantal aannames:
     • De aanneming van een Gaussische likelihood (wat niet ideaal is voor hogere-orde statistieken).
     • Het gebruik van alleen gelijke-schaal statistieken (meerdere schalen zouden meer informatie kunnen opleveren).
     • Het marginaliseren over parameters zoals $h$ en $w_0$ zonder tomografie.
     • De sterke correlatie tussen schalen en de snelle afname van het signaal op kleine schalen (waar de vierde-orde informatie het grootst is, maar waar simulaties onnauwkeurig kunnen zijn).

Betekenis en Conclusie

Dit werk vormt de eerste stap in een serie van twee artikelen (het tweede artikel, P25, presenteert de daadwerkelijke detectie in DES-Y3 data).

• Voor Stage III surveys (zoals DES): De toegevoegde waarde van vierde-orde statistieken voor het verbeteren van kosmologische constraints lijkt beperkt, tenzij er een grote set simulaties of een zeer nauwkeurig theoretisch model beschikbaar is om de fouten te beheersen.
• Voor Stage IV surveys (zoals Euclid): De methode is cruciaal voor toekomstige surveys met een extreem hoog signaal-ruisverhouding. Hoewel de vierde-orde statistieken mogelijk niet de kosmologische parameters zelf sterk verbeteren, kunnen ze dienen als een krachtige consistentiecheck en een test voor systematische fouten (zoals intrinsieke uitlijning of baryonische effecten). Als ze degenereren met lagere-orde statistieken, kunnen ze helpen bij het kalibreren van storende parameters (nuisance parameters).

Samenvattend biedt dit artikel het noodzakelijke theoretische en numerieke fundament om vierde-orde schuifstatistieken toe te passen in de moderne kosmologie, en bevestigt het dat deze statistieken, hoewel ze nu nog beperkt zijn door datakwaliteit en modelcomplexiteit, essentieel zullen zijn voor de precisiecosmologie van de toekomst.