Incorporating curved geometry in cosmological simulations

Deze paper introduceert een volledig relativistisch raamwerk voor het uitvoeren van kosmologische simulaties met gekromde ruimtelijke geometrie door een bolvormig stukje gekromde ruimtetijd in te bedden in een gat van een platte buitenwereld, waardoor consistente randvoorwaarden worden geboden en een nauwkeurige voorspelling van waarnemingen mogelijk wordt.

De kern

Kosmologische Simulaties met een Kromme Wereld: Een Simpele Uitleg

Stel je voor dat je een gigantische, driedimensionale puzzel wilt bouwen om te begrijpen hoe het heelal is ontstaan en hoe sterrenstelsels zich vormen. Normaal gesproken bouwen wetenschappers deze puzzels in een perfect platte ruimte, alsof ze op een groot, onuitgerekt laken werken. Maar wat als het laken eigenlijk niet plat is, maar gebogen? Wat als het heelal een beetje op een ballon lijkt (gesloten) of op een zadel (open)?

Dit is precies het probleem dat Julian Adamek en Renan Boschetti in hun nieuwe onderzoek proberen op te lossen. Ze hebben een nieuwe manier bedacht om computersimulaties te draaien voor een heelal dat krom is, in plaats van plat.

Hier is hoe ze dat doen, vertaald naar alledaagse taal:

1. Het Probleem: De "Periode" van de Puzzel

De meeste supercomputers die het heelal simuleren, werken met een trucje: ze nemen een kubus van ruimte en zeggen: "Als je aan de ene kant uitkomt, kom je direct weer aan de andere kant binnen." Dit heet periodieke randvoorwaarden. Het is alsof je in een Pac-Man-spel zit; als je links uit het scherm vliegt, kom je rechts weer binnen.

Dit werkt perfect als het heelal plat is. Maar als het heelal krom is (zoals een bol), werkt deze truc niet meer. Je kunt een bol niet perfect inpakken in een kubus zonder dat het scheurt of vervormt. Als je dit toch doet, krijg je een verkeerd beeld van hoe licht zich verplaatst over enorme afstanden.

2. De Oplossing: De "Kromme Eiland" in een "Platte Zee"

De auteurs hebben een slimme oplossing bedacht die ze de "Einstein-Straus" methode noemen (vernoemd naar de natuurkundigen die dit idee al in de jaren '40 hadden).

Stel je voor:

• Je hebt een groot, plat zwembad (de externe ruimte). Dit is waar de computer zijn rekenkracht op richt en waar de "Pac-Man-regels" gelden.
• In het midden van dit zwembad hakken ze een perfect ronde gat uit.
• In dat gat plaatsen ze een kromme wereldbol (de interne ruimte). Dit is het stukje heelal dat we eigenlijk willen bestuderen.

Deze bol zit precies in het gat en past er perfect in. De randen van de bol sluiten naadloos aan op de randen van het gat in het platte water.

• Voor de computer: Het ziet eruit als een normaal, plat zwembad met een gat. De regels blijven hetzelfde.
• Voor de waarnemer: Als je in het midden van die bol staat, zie je een krom heelal. Licht dat reist, volgt de kromming van die bol, net zoals een vliegtuig dat een route vliegt over de aarde (die bolvormig is) een andere route neemt dan een lijn op een platte kaart.

3. Waarom is dit belangrijk?

Wetenschappers willen weten of het heelal plat is of krom. Ze kijken naar het licht van verre sterrenstelsels.

• Als het heelal krom is, verandert dit hoe we afstanden meten. Het is alsof je een liniaal gebruikt op een ballon: de afstand tussen twee punten is anders dan op een plat vel papier.
• De oude manier om dit te simuleren (de "Scheiding van het Universum"-methode) deed alsof de ruimte plat was, maar deed alsof de uitdijing van het heelal anders was. Dit werkt goed voor kleine stukjes, maar faalt als je kijkt over enorme afstanden (zoals licht dat miljarden jaren reist).

De nieuwe methode van Adamek en Boschetti houdt rekening met de echte geometrie. Ze simuleren niet alleen hoe het heelal groeit, maar ook hoe de ruimte zelf gebogen is.

4. De Resultaten: Alles klopt!

Ze hebben hun methode getest met twee scenario's:

1. Een heelal zonder sterren (alleen materie): Ze keken hoe licht zich verplaatste. De resultaten kwamen exact overeen met de wiskundige theorieën voor een krom heelal.
2. Een heelal met sterrenstelsels: Ze lieten materie klonten vormen (zoals in het echte heelal). Ze ontdekten dat, ongeacht waar je in die kromme bol staat (in het midden of aan de rand), je precies dezelfde patronen ziet. Dit bevestigt een fundamenteel principe van de kosmologie: het heelal ziet er overal hetzelfde uit.

5. De "Relativistische Massadefect" (Een lastig woord, maar simpel idee)

Er is nog een klein, maar belangrijk detail. In een krom heelal is de zwaartekracht anders dan in een plat heelal. De auteurs ontdekten dat ze een beetje extra "massa" moesten toevoegen aan hun simulatie om de zwaartekracht correct te houden.

• Analogie: Stel je voor dat je een koffer vol kleding (deeltjes) in een kofferbak (de ruimte) stopt. Als de kofferbak krom is, past er meer kleding in dan als hij plat is, maar de kleding zelf is niet zwaarder. De "kromming" zorgt voor een extra zwaartekrachts-effect dat je moet compenseren. Ze hebben dit precies berekend zodat hun simulatie eerlijk blijft.

Conclusie: Waarom moeten we hier blij mee zijn?

Deze paper is een grote stap voorwaarts. Het stelt ons in staat om voorspellingen te doen voor toekomstige telescopen (zoals de Euclid-missie of de Rubin-observatorium) die het heelal tot in de kleinste details gaan bekijken.

Als het heelal inderdaad een klein beetje krom is (wat sommige metingen suggereren), dan kunnen we dat alleen maar begrijpen als we simulaties draaien die de kromming van de ruimte zelf in rekening brengen, en niet alleen de snelheid waarmee het heelal uitdijt.

Kortom: Ze hebben een manier gevonden om een kromme wereld te bouwen binnen een platte computer, zodat we eindelijk kunnen zien wat er echt gebeurt als we naar de verre oorden van het heelal kijken.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Incorporating curved geometry in cosmological simulations" van Adamek en Boschetti, geschreven in het Nederlands.

Probleemstelling

Ruimtelijke kromming is een fundamenteel kosmologisch parameter dat strikt wordt beperkt door waarnemingen (zoals CMB-anisotropieën en BAO-metingen). Hoewel de "separate universe"-benadering het effect van kromming op de expansiegeschiedenis in kleine sub-volumes goed kan modelleren, faalt deze methode bij het accuraat weergeven van de niet-Euclidische geometrie op grote schaal.

De kern van het probleem is dat observables (zoals afstanden en hoekverdelingen) worden berekend over grote afstanden waar fotonen reizen van hoge roodverschuivingen naar de waarnemer. In een gebogen ruimte zijn omtrekken en volumes anders dan in een platte ruimte (bijvoorbeeld: de omtrek van een cirkel is niet $\pi$ keer de diameter). Bestaande N-body simulaties gebruiken doorgaans periodieke randvoorwaarden, wat impliceert dat de ruimte op de schaal van de simulatiebox noodzakelijk plat is. Het is technisch uitdagend om een gebogen ruimtetijd te simuleren binnen deze beperkingen zonder de consistentie van de Einstein-vergelijkingen te schenden of de randvoorwaarden te negeren.

Methodologie

De auteurs presenteren een volledig relativistisch raamwerk om kosmologische simulaties uit te voeren met een gebogen ruimtelijke geometrie. De aanpak bestaat uit de volgende stappen:

1. Embedding (Inbedding):

   • Een sferisch stuk van een gebogen Friedmann-Lemaître-Robertson-Walker (FLRW) ruimtetijd (de "curved patch") wordt ingebed in een extern, ruimtelijk vlak FLRW-gebied.
   • Dit wordt gerealiseerd door een exacte oplossing van Einstein's vergelijkingen te gebruiken: een uniforme bol van materie die vrij expandeert in een vacuüm.
   • De overgang tussen de gebogen binnenkant en het vlakke buitenste gebied wordt beheerd via de Israel-junctievoorwaarden.
   • Het vacuümgebied tussen de gebogen patch en de rand van de simulatiebox wordt beschreven door de Schwarzschild-metriek.
   • Door dit systeem in een kubische box te plaatsen, kunnen standaard periodieke randvoorwaarden worden toegepast.

2. Coördinatentransformatie (Gauge Transformatie):

   • De oplossing voor drukloze materie (stof) is analytisch bekend in de comoveerende synchrone coördinaten (Lemaître-Tolman-Bondi of LTB metriek).
   • Simulatiecodes (zoals gevolution) werken echter in de Poisson-gauge.
   • De auteurs leiden een transformatie af van de LTB-metriek naar de Poisson-gauge, tot op tweede orde in perturbatietheorie. Dit is cruciaal voor nauwkeurigheid bij grote krommingen en voor het correct modelleren van de koppeling tussen kromming en materieperturbaties.

3. Initieel Voorwaarde en Massadefect:

   • Er wordt rekening gehouden met het relativistische massadefect: de actieve gravitationele massa is niet simpelweg de som van de deeltjesmassa's vanwege gravitationele bindingsenergie.
   • Om de Hamiltoniaanse constraint te voldoen en periodieke randvoorwaarden te respecteren, wordt een extra massa-term toegevoegd aan de deeltjes (of gelijkmatig verdeeld over de simulatie) om het volume-overschot in de gebogen patch te compenseren.

4. Materieperturbaties:

   • Voor de initiële perturbaties wordt aangenomen dat de modusfuncties op de schaal van de simulatie goed worden benaderd door Fourier-modus van de vlakke geometrie (een benadering van de "separate universe"-aanpak), hoewel de auteurs erkennen dat in een gesloten universum hypersferische harmonischen de correcte basis zouden zijn.
   • De perturbaties worden "geapodiseerd" (afgevlakt) aan de rand van de gebogen patch om overgangen naar het vlakke gebied glad te maken.

Belangrijkste Bijdragen

• Eerste volledig relativistische simulatie: Dit is het eerste werk dat een volledige relativistische behandeling toepast op N-body simulaties met globale ruimtelijke kromming, waarbij zowel de expansiegeschiedenis als de niet-Euclidische geometrie correct worden meegenomen.
• Exacte oplossing voor Schwarzschild in Poisson-gauge: Als bijproduct presenteren de auteurs een exacte uitdrukking voor de Schwarzschild-metriek in Poisson-gauge coördinaten, uitgedrukt via Weierstrass-elliptische functies (zie Appendix B).
• Oplossing voor randvoorwaarden: Het biedt een praktische methode om gebogen ruimtetijden te simuleren binnen de bestaande infrastructuur van periodieke N-body codes zonder de code-architectuur fundamenteel te hoeven veranderen.

Numerieke Resultaten

De auteurs hebben hun methode getest met de relativistische N-body code gevolution in drie experimenten:

1. Einstein-de Sitter simulatie (Hoge kromming):

   • Getest met $\tilde{\Omega}_K = -0.25$ (gesloten universum).
   • Resultaat: De waargenomen afstand-roodverschuivingsrelatie ($d_A(z)$) voor zowel een centrale waarnemer als een randwaarnemer komt overeen met de exacte theoretische voorspelling tot op een nauwkeurigheid van 1 op 1000, zolang het lichtkegel binnen de gebogen patch blijft.
   • De "separate universe"-benadering (die alleen de expansiegeschiedenis aanpast maar de geometrie plat houdt) vertoont fouten van ongeveer 1% bij $z \approx 1$.

2. Gebogen $\Lambda$CDM simulatie (Moderate kromming):

   • Getest met $\tilde{\Omega}_K = -0.1$ en een kosmologische constante.
   • Resultaat: Dezelfde hoge nauwkeurigheid wordt behaald. De Hubble-diagrammen blijven isotroop binnen de patch, wat de consistentie van het model bevestigt.

3. Simulatie met materieperturbaties:

   • Toevoeging van initiële dichtheidsfluctuaties.
   • Meting van de hoekelijke tweepuntsfunctie (angular power spectrum $C_\ell$) in roodverschuivingsbinten.
   • Resultaat: De waargenomen clusteringstatistieken zijn consistent met lineaire theorie (voorspeld door CLASS) en voldoen aan het kosmologische principe (waarnemers op verschillende locaties zien dezelfde statistieken). De "separate universe"-benadering zou hier fouten van ongeveer 10% kunnen introduceren voor deze specifieke krommingswaarde.

Significantie en Conclusie

De studie toont aan dat het negeren van de niet-Euclidische geometrie in grote schaal simulaties leidt tot systematische fouten in de voorspelling van observables, zelfs als de expansiegeschiedenis correct wordt gemodelleerd.

• Toekomstige surveys: Voor toekomstige surveys zoals Euclid, DESI en Rubin Observatory, die precisies onder het procentniveau eisen, is het essentieel om de geometrische effecten van kromming correct te modelleren.
• Newtonse codes: Hoewel de methode hier volledig relativistisch is, suggereert de auteurs dat bestaande Newtonse codes deze resultaten met een goede benadering kunnen reproduceren, mits de initiële voorwaarden en post-processing (lichtkegel reconstructie) correct worden aangepast.
• Generalisatie: De huidige implementatie is beperkt tot gesloten universa ($\Omega_K < 0$), maar de methode kan worden uitgebreid naar open universa ($\Omega_K > 0$) door een compenserende massaschil toe te voegen.

Kortom, dit werk levert een robuust raamwerk voor het forward-modelling van waarnemingen in een gebogen universum, wat cruciaal is voor het testen van inflatietheorieën en het meten van fundamentele kosmologische parameters.