Euclid preparation. Galaxy power spectrum modelling in redshift space

Voor de voorbereiding op de Euclid-missie toont dit onderzoek aan dat het ${\rm VDG_ \infty}$-model en de BACCO-emulator de effecten van roodverschuivingsruimtestructuur op de galaxiekrachtspectrum beter modelleren dan het EFT-model, waardoor nauwkeurigere kosmologische parameterschattingen mogelijk zijn tot kleinere schalen.

De kern

De Euclid-missie: Een kosmische zoektocht en de perfecte kaart

Stel je voor dat het heelal een gigantische, donkere oceaan is. De sterrenstelsels zijn als schitterende eilanden die verspreid liggen in deze duisternis. De Europese ruimtevaartorganisatie ESA heeft een nieuwe missie, Euclid, die als een superkrachtige verrekijker door deze oceaan vaart. Het doel? Een driedimensionale kaart maken van miljarden van deze eilanden om te begrijpen hoe het heelal is opgebouwd en hoe het zich gedraagt.

Maar hier zit een addertje onder het gras: als we naar deze sterrenstelsels kijken, zien we ze niet precies waar ze echt zijn. Ze lijken een beetje verschoven, alsof ze op een vervormde spiegel worden weerspiegeld. Dit komt door twee dingen:

1. Het heelal breidt zich uit (ze bewegen weg van ons).
2. Ze hebben hun eigen "eigenbeweging" door de zwaartekracht van andere sterrenstelsels (net als auto's die in een file versnellen of vertragen).

Deze vervorming noemen wetenschappers roodverschuiving-distorsies. Om de echte kaart te maken, moeten we een wiskundig model hebben dat deze vervorming perfect kan "ontwarren".

In dit artikel kijken drie verschillende teams (of modellen) naar hoe ze deze wiskundige puzzel het beste kunnen oplossen. Ze testen hun modellen op een virtuele simulatie van het heelal (een soort "videospel" van het heelal dat ze zelf hebben gemaakt) om te zien wie de beste kaart tekent.

De drie kampioenen in de ring

De auteurs vergelijken drie verschillende methoden om de vervorming te berekenen:

1. De EFT (Effective Field Theory):

   • De analogie: Stel je voor dat je probeert de beweging van een drukke menigte te voorspellen door alleen te kijken naar hoe mensen zich gedragen als ze rustig staan, en dan kleine correcties toe te voegen voor mensen die rennen.
   • Hoe het werkt: Dit is een heel strakke, theoretische formule. Het werkt fantastisch op grote schaal (waar de menigte rustig is), maar begint te struikelen als je te dicht bij de chaos komt (op kleine schaal). Het model probeert de chaos te beschrijven met "correctie-termen", maar op een gegeven moment is de formule te simpel voor de complexiteit.

2. De VDG∞ (Velocity Difference Generator):

   • De analogie: In plaats van alleen te kijken naar de rustige menigte, kijkt dit model ook naar hoe mensen zich gedragen als ze in een drukke file zitten. Het gebruikt een slimme formule die specifiek is ontworpen om die "file-gedrag" (de snelle bewegingen) nauwkeurig te beschrijven.
   • Hoe het werkt: Dit model is een hybride tussen pure theorie en slimme observatie. Het beschrijft de kleine, chaotische bewegingen veel beter dan de EFT. Het kan dus verder "inzoomen" op de details zonder dat de formule in elkaar stort.

3. De BACCO Emulator:

   • De analogie: Stel je voor dat je in plaats van een formule te schrijven, een supercomputer hebt die miljoenen keren een virtueel heelal heeft nagebootst. Als je een vraag hebt, kijkt de computer in zijn enorme database van eerdere simulaties en zegt: "Ah, dit lijkt op situatie X, hier is het antwoord."
   • Hoe het werkt: Dit is een hybride model. Het gebruikt de theorie als basis, maar vult de gaten op met data uit superkrachtige computersimulaties. Het is als een GPS die niet alleen op kaarten kijkt, maar ook live verkeersinformatie gebruikt.

Wat hebben ze ontdekt?

De onderzoekers hebben hun modellen getest op vier verschillende tijdstippen in de geschiedenis van het heelal (van 0,9 tot 1,8 miljard jaar geleden). Ze keken naar drie dingen:

• Hoe nauwkeurig is de kaart? (Komen de coördinaten overeen met de werkelijkheid?)
• Hoe scherp is de focus? (Kunnen we de details van het heelal goed zien?)
• Hoe goed past het model? (Ziet het er logisch uit?)

De resultaten:

• De EFT (De theorie-purist): Werkt goed zolang je niet te dicht bij de chaos komt. Zodra ze te ver inzoomen (op schalen kleiner dan een bepaalde grens), begint het model fouten te maken. Het wordt "vooroordeelsvol" en geeft een verkeerde kaart. Het is alsof je probeert een foto van een snel bewegende auto te maken met een trage camera; het resultaat wordt wazig.
• De VDG∞ (De slimme formule): Deze wint het van de EFT. Door de "file-gedrag" beter te begrijpen, kan het model veel verder inzoomen zonder fouten. Het levert een scherpere kaart op en geeft betere antwoorden over de snelheid van het heelal en de hoeveelheid donkere materie.
• De BACCO Emulator (De database-kijker): Deze is ook heel goed, vooral op grote schaal. Het is echter een beetje "verzadigd": op een zeker punt haalt het niet meer veel extra informatie uit het inzoomen. Het is alsof je al zo veel details hebt dat het toevoegen van nog meer pixels niet meer helpt.

De grote les voor de toekomst

Het belangrijkste nieuws voor de Euclid-missie is dit:
Als we de VDG∞-methode gebruiken, kunnen we veel meer informatie uit de data halen dan met de oude methoden. We kunnen "dieper" inzoomen op het heelal zonder dat de resultaten onbetrouwbaar worden.

• BACCO is geweldig, maar vraagt veel rekenkracht (het is zwaar om te draaien).
• EFT is licht en snel, maar niet nauwkeurig genoeg voor de diepste details.
• VDG∞ is de gouden middenweg: het is nauwkeurig, snel genoeg, en haalt het meeste uit de data.

Conclusie in één zin:
Om de perfecte kaart van ons heelal te maken met de Euclid-ruimtesonde, moeten we de "slimme formule" (VDG∞) gebruiken die de chaos van de sterrenstelsels beter begrijpt dan de oude theorieën, zodat we de geheimen van donkere energie en de uitdijing van het heelal eindelijk kunnen ontcijferen.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Galaxy power spectrum modelling in redshift space" van de Euclid Collaboration, geschreven in het Nederlands.

Probleemstelling

De analyse van grote schaalstructuur (LSS) surveys, zoals de komende Euclid-missie, is cruciaal voor het begrijpen van de geometrie en de groei van structuren in het universum. Een van de belangrijkste bronnen van kosmologische informatie is de galaxiekracht-spectrum (power spectrum) in roodshift-ruimte. Echter, waarnemingen in roodshift-ruimte zijn vervormd door de eigenbewegingen (peculiare snelheden) van sterrenstelsels, een fenomeen bekend als Roodshift-Space Distortions (RSD).

De uitdaging ligt in het nauwkeurig modelleren van deze RSD-effecten, vooral op kleine schalen (hoge golfgetallen $k$), waar niet-lineaire fysica en de overgang van de lineaire naar de niet-lineaire regime de theorie complex maken. Bestaande modellen moeten een balans vinden tussen het benutten van maximale informatie (door tot kleine schalen te gaan) en het vermijden van systematische fouten (bias) die ontstaan door onnauwkeurige modellen. Voor de Euclid-missie is het essentieel om te bepalen welk model de beste prestaties levert voor de analyse van emissielijn-galaxies (H$\alpha$) in het roodshift-bereik $0.9 \le z \le 1.8$.

Methodologie

De auteurs vergelijken drie verschillende modelleringbenaderingen voor het anisotrope galaxiekracht-spectrum, allemaal gebaseerd op een één-lus perturbatietheorie (one-loop perturbation theory) expansie, maar verschillend in hun behandeling van de mapping van reële naar roodshift-ruimte:

1. EFT (Effective Field Theory): Een volledig perturbatieve aanpak waarbij de exponentiële prefactor in de roodshift-ruimte mapping wordt uitgebreid. Kleine schaal-effecten worden gemodelleerd via een reeks "counterterms" die de ultraviolette (UV) onvolkomenheden en Finger-of-God (FoG) onderdrukking absorberen.
2. VDG$\infty$ (Velocity Difference Generator at infinity): Een model dat de bijdrage van geventileerde snelheden beschrijft via een analytische dempingsfunctie (afgeleid van de resummation van niet-lineaire termen). Dit model behandelt de real-naar-roodshift mapping anders dan EFT, met een focus op de fysica van de snelheidsverschilverdeling.
3. BACCO (Hybride Emulator): Een hybride aanpak die een Lagrangiaanse bias-expansie combineert met componenten die zijn geëmulleerd op basis van hoge-resolutie N-body simulaties. Dit model gebruikt daadwerkelijke simulatie-data om de complexe niet-lineaire effecten te vangen.

Data en Validatie:

• De modellen worden getest op mock catalogi van H$\alpha$-emitterende galaxies, gegenereerd uit de Flagship I N-body simulatie (een simulatie met twee triljoen deeltjes).
• Er worden vier roodshift-bins gebruikt: $z = 0.9, 1.2, 1.5, 1.8$.
• De analyse omvat de Legendre-multipolen ($\ell=0, 2, 4$) van het kracht-spectrum.
• De prestaties worden beoordeeld aan de hand van drie metrieken:
  • FoB (Figure of Bias): Meet de nauwkeurigheid van het herstel van de ingevoerde kosmologische parameters.
  • FoM (Figure of Merit): Meet de statistische kracht (precisie) van de parameterbeperkingen.
  • p-waarde: Beoordeelt de goodness-of-fit.

De analyse varieert de maximale golfmodus ($k_{max}$) om de geldigheidsbereik van elk model te bepalen.

Belangrijkste Resultaten

1. Geldigheidsbereik en Bias:

   • Het EFT-model is robuust op grote schalen, maar vertoont significante bias in de kosmologische parameters voor $k_{max} \gtrsim 0.25 \, h \, \text{Mpc}^{-1}$. De FoB overschrijdt de 2$\sigma$ drempel, wat aangeeft dat het model te vroeg "breekt" op niet-lineaire schalen.
   • Zowel VDG$\infty$ als BACCO presteren aanzienlijk beter en blijven binnen de 2$\sigma$ grenzen tot $k_{max} \lesssim 0.35 \, h \, \text{Mpc}^{-1}$ (en zelfs hoger voor VDG$\infty$).
   • Het EFT-model toont een toenemende bias, vooral bij lagere roodshifts ($z=0.9, 1.2$), terwijl VDG$\infty$ en BACCO stabiel blijven.

2. Beperkende Kracht (Constraining Power):

   • BACCO levert de beste beperkingen op schalen $k_{max} \lesssim 0.20 \, h \, \text{Mpc}^{-1}$, maar vertoont een verzadiging (plateau) in de informatie-inwinning bij $k_{max} \gtrsim 0.30 \, h \, \text{Mpc}^{-1}$, vooral bij hogere roodshifts.
   • VDG$\infty$ blijft de beperkende kracht verbeteren tot $k_{max} \gtrsim 0.30 \, h \, \text{Mpc}^{-1}$ en zelfs daarboven. Het model profiteert meer van de informatie op kleinere schalen dan de andere twee.
   • In vergelijking met EFT verbetert VDG$\infty$ de beperkingen op kosmologische parameters met een factor tot twee.

3. Kosmologische Parameters:

   • Alle modellen kunnen de Hubble-parameter ($h$) met een precisie van ongeveer 1% herstellen.
   • De dichtheid van koud donkere materie ($\omega_c$) wordt binnen 3% nauwkeurigheid hersteld.
   • De amplitude van de scalarfluctuaties ($A_s$) wordt door VDG$\infty$ het best beperkt (beter dan 5% precisie) zonder significante bias, terwijl EFT hier moeite mee heeft op kleine schalen.

4. Nuisance Parameters en Counterterms:

   • Analyse van de counterterms ($c_0, c_2, c_4$) toont aan dat bij EFT deze parameters bij hoge $k_{max}$ beginnen te "drijven" om residu-effecten op te vangen, wat wijst op een tekortkoming in het fysieke model. Bij VDG$\infty$ convergeren deze parameters naar nul of stabiele waarden, wat suggereert dat de analytische dempingsfunctie de RSD-effecten effectiever modelleert.
   • Het model VDG$\infty$ toont een betere herstelbaarheid van shot-noise parameters ($N_P^0$) dan EFT.

Bijdrage en Significantie

• Optimalisatie voor Euclid: Dit werk biedt een cruciale validatie voor de data-analysestrategie van de Euclid-missie. Het identificeert dat modellen die de real-naar-roodshift mapping beter behandelen (zoals VDG$\infty$ en BACCO), een groter geldigheidsbereik hebben dan traditionele EFT-uitbreidingen.
• Uitbreiding van het Schaalbereik: De studie toont aan dat het gebruik van VDG$\infty$ het bereik van geldige analyse kan uitbreiden van $k_{max} \approx 0.25 \, h \, \text{Mpc}^{-1}$ (EFT) naar $k_{max} \approx 0.40 \, h \, \text{Mpc}^{-1}$. Dit betekent dat Euclid aanzienlijk meer kosmologische informatie kan extraheren uit dezelfde dataset.
• Modelkeuze: Hoewel BACCO zeer nauwkeurig is, is het computatie-intensiever. VDG$\infty$ biedt een uitstekend compromis: het is volledig theoretisch (zoals EFT) maar bereikt een nauwkeurigheid en geldigheidsbereik dat dicht bij de hybride BACCO-aanpak ligt, met een aanzienlijke verbetering ten opzichte van standaard EFT.
• Toekomstige Richtingen: De resultaten benadrukken de noodzaak van fysiek gemotiveerde priors en verdere verfijning van countertermen, vooral voor het modelleren van de sterk niet-lineaire regime. Het werk legt de basis voor de toekomstige analyse van de volledige Euclid-dataset.

Kortom, dit artikel concludeert dat voor de analyse van Euclid H$\alpha$-data de VDG$\infty$-benadering (en in mindere mate BACCO) superieur is aan de standaard EFT-benadering, omdat het een robuustere modellering biedt op niet-lineaire schalen, wat leidt tot nauwkeurigere en minder bevooroordeelde schattingen van de kosmologische parameters.