Cosmic Shear in Effective Field Theory at Two-Loop Order: Revisiting $S_8$ in Dark Energy Survey Data

Probleemstelling

Cosmische schuif (cosmic shear), de coherente vervorming van de schijnbare vormen van verre sterrenstelsels door gravitationele lensing, is een van de krachtigste methoden om kosmologische parameters te meten, met name de dichtheid van materie ($\Omega_m$) en de amplitude van massavariaties ($\sigma_8$). De combinatie hiervan wordt vaak uitgedrukt als $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3}$.

Er bestaat echter een bekende spanning (de "$S_8$-tension") tussen metingen uit zwakke lensing (zoals de Dark Energy Survey, DES) en voorspellingen gebaseerd op de kosmische microgolfachtergrondstraling (CMB) binnen het standaard $\Lambda$CDM-model. DES-metingen geven doorgaans lagere $S_8$-waarden dan de Planck-CMB-data.

De beperkende factor bij het analyseren van kosmische schuifdata is het begrip van niet-lineaire effecten en baryonische astrofysica (zoals feedback van sterrenvorming en zwarte gaten). Traditionele analyses vertrouwen op N-body-simulaties en empirische fitting-functies (zoals Halofit of HMCode). Om systematische fouten door baryonische feedback te vermijden, worden deze analyses vaak beperkt tot grote schalen (lage golflengten, $k \lesssim 1 \, h\text{Mpc}^{-1}$), waardoor een groot deel van de beschikbare data onbenut blijft. De vraag is of de resterende schalen die robuust zijn tegenover baryonische onzekerheden, voldoende zijn om nauwkeurige beperkingen te stellen, en of deze schalen theoretisch goed beschreven kunnen worden zonder zware simulaties.

Methodologie

De auteurs introduceren een consistente toepassing van de Effectieve Veldtheorie voor de Grote Schaalstructuur (EFT) op twee-lus niveau (two-loop order) voor kosmische schuifdata. Dit is een analytische benadering die de standaard kosmologische perturbatietheorie uitbreidt naar het niet-lineaire regime.

De kernpunten van de methodologie zijn:

1. Two-Loop EFT voor Materiekracht:
   Het spectrum van de materiedichtheid ($P_m$) wordt berekend tot op twee-lus orde. De formule omvat standaard Euleriaanse perturbatietheorie (SPT) termen en hogere-afgeleide correcties (counterterms) die de invloed van niet-perturbatieve kleine schalen renormaliseren:
   $$P_m(k, z) = P^{\text{lin}}_m + P^{\text{1-loop, SPT}}_m + P^{\text{2-loop, SPT}}_m + \text{counterterm-termen}$$
   De countertermen worden genormaliseerd met een niet-lineaire schaal $k_{nl}(z)$ en gemodelleerd met een spline-functie om de tijdafhankelijkheid (roodverschuiving) flexibel te beschrijven zonder sterke aannames.

2. Lensings-countertermen (LCTs):
   Een cruciale innovatie is de behandeling van de projectie-integraal van het 3D-spectrum naar het 2D-lensing-spectrum. Zelfs als het 3D-spectrum wordt ge-renormaliseerd, blijft de lensing-integraal gevoelig voor ultraviolette (UV) modi (kleine schalen) vanwege de eindige breedte van de lensingskern. De auteurs gebruiken een Taylor-expansie van de lensingskern in de UV-limiet om een reeks vrije EFT-coëfficiënten (lensing counterterms) in te voeren. Deze marginaliseren over de onzekerheid van de kleine schalen zonder de data te hoeven afsnijden op grote schalen.

3. Intrinsieke Uitlijningen (IA):
   De intrinsieke uitlijning van sterrenstelsels (waarbij stelsels niet willekeurig georiënteerd zijn, maar door het lokale zwaartekrachtveld worden uitgelijnd) wordt gemodelleerd met behulp van Lagrangiaanse EFT voor spin-2 operatoren, wat een consistent kader biedt voor het scheiden van lensing- en IA-signalen.

4. Emulatie en Validatie:
   Omdat twee-lus berekeningen computatieel zwaar zijn, hebben de auteurs een emulator ontwikkeld die is getraind op Monte Carlo-berekeningen van de lus-integrals. Ze valideren hun pijplijn op "mock" data gegenereerd door de Aemulus-$\nu$ emulator (een niet-lineaire simulatie-emulator) en vinden dat hun model onbevooroordeelde kosmologische beperkingen oplevert.

Belangrijkste Bijdragen

• Eerste consistente analyse: Dit is de eerste analyse van publieke kosmische schuifdata (DES-Y3) die consistent gebruikmaakt van EFT op twee-lus niveau.
• Renormalisatie van kleine schalen: Door het gebruik van lensing-countertermen kunnen de auteurs data op kleinere schalen (hoger $\ell$) gebruiken dan traditionele methoden, zonder dat de analyse gevoelig wordt voor onbekende baryonische fysica.
• Model-onafhankelijkheid: De methode is niet afhankelijk van specifieke baryonische feedback-modellen; de onzekerheid wordt gemarginaliseerd via de EFT-countertermen.
• Toepasbaarheid op Beyond-$\Lambda$CDM: Omdat EFT analytisch is, kan het eenvoudig worden uitgebreid naar modellen buiten het standaardmodel (bijv. dynamische donkere energie), waar simulatie-benaderingen vaak vastlopen door gebrek aan dekking in de parameter ruimte.

Resultaten

De auteurs passen hun model toe op de DES-Y3 (Dark Energy Survey Year 3) data:

• $S_8$ Beperkingen binnen $\Lambda$CDM:
  Ze vinden een waarde van $S_8 = 0.783^{+0.038}_{-0.031}$. Dit is vergelijkbaar met de officiële DES-Y3 resultaten die op empirische modellen gebaseerd zijn ($S_8 \approx 0.793$), maar bereikt dit met een fundamenteel andere, analytische theorie.
• Oplossing van de $S_8$-spanning:
  Wanneer de DES-data wordt gecombineerd met externe priors (BAO van DESI, CMB van Planck, en supernova-data), vinden ze $S_8 = 0.824 \pm 0.029$.
  Dit resultaat is consistent met de Planck-CMB voorspelling ($S_8 \approx 0.832$). De auteurs concluderen dat er geen $S_8$-spanning is in de DES-Y3 zwakke lensing-data, mits deze correct wordt geanalyseerd met een robuuste theorie die de groei van structuur goed beschrijft. De eerdere spanning lijkt te komen uit de beperkingen van de gebruikte schaalcut-offs of de modelafhankelijkheid in eerdere analyses.
• Dynamische Donkere Energie:
  In een $w_0w_a$CDM-model (waarbij donkere energie dynamisch is) vinden ze eveneens consistente resultaten ($S_8 = 0.827 \pm 0.010$ wanneer gecombineerd met BAO en SNe), wat aangeeft dat de conclusie robuust is ongeacht het onderliggende kosmologische model.

Significantie en Toekomstperspectief

De studie toont aan dat een aanzienlijk deel van de data die in standaard schuifanalyses wordt gebruikt, perturbatief is en nauwkeurig beschreven kan worden met twee-lus EFT. Dit heeft grote gevolgen:

1. Efficiëntie: Het elimineert de noodzaak voor zware N-body-simulaties voor elke nieuwe kosmologische hypothese.
2. Nieuwe Fysica: Het opent de deur voor het testen van nieuwe fysica (zoals ultralichte axionen, vroege donkere energie) met lopende en toekomstige surveys zoals Euclid, Roman en LSST, aangezien de theorie flexibel uitbreidbaar is.
3. Robuustheid: Het bewijst dat de $S_8$-spanning mogelijk een artefact is van de analysemethoden (schaalcut-offs en baryonische modellering) in plaats van een teken van nieuwe fysica.

Kortom, deze paper markeert een paradigmaverschuiving van empirische, simulatie-afhankelijke analyses naar een fundamentele, analytische benadering voor het begrijpen van de groei van kosmische structuren.