COSMOS: A numerical relativity code specialized for PBH formation

COSMOS is een standalone, met OpenMP geparallelliseerde C++ numerieke relativiteitscode die is ontworpen om de vorming van oerzwarte gaten te simuleren door de Einstein-vergelijkingen in 3+1 dimensies op te lossen met behulp van gespecialiseerde niet-Cartesiaanse opschalingscoördinaten en vaste mesh-verfijning om niet-lineaire gravitationele dynamica te verwerken.

De kern

Het Grote Plaatje: Zwart Gat Bouwen vanaf de Basis

Stel je het vroege universum niet voor als een gladde, kalme oceaan, maar als een onstuimige zee met gigantische golven. Soms worden deze golven zo hoog en zwaar dat ze in elkaar klappen, waardoor Primaire Zwarte Gaten (PBH's) ontstaan. In tegen tegenstelling tot de zwarte gaten waar we meestal over horen — die de overblijfselen zijn van dode sterren — werden PBH's direct geboren uit de "rimpelingen" in het weefsel van ruimte en tijd vlak na de oerknal.

Wetenschappers willen begrijpen hoe deze zwarte gaten ontstaan, maar de wiskunde is ongelooflijk moeilijk. Het is alsof je precies probeert te voorspellen hoe een specifieke druppel water zal spatten wanneer deze een plas raakt, maar de plas is aan het uitdijen en het water is gemaakt van pure zwaartekracht.

De Oplossing: COSMOS (Het Digitale Laboratorium)

Dit artikel introduceert COSMOS, een computerprogramma (geschreven in C++) dat specifiek is ontworpen om deze geboorten van zwarte gaten te simuleren. Denk aan COSMOS als een hoogtechnologische, digitale windtunnel voor zwaartekracht. Net zoals ingenieurs schaalmodellen van auto's bouwen om te testen hoe lucht eromheen stroomt, gebruiken natuurkundigen COSMOS om "schaalmodellen" van het vroege universum te bouwen om te zien hoe zwaartekracht zich gedraagt wanneer de boel chaotisch wordt.

Hoe het Werkt: De "Zoomlens"-truc

Een van de grootste uitdagingen in deze simulatie is dat er twee heel verschillende groottes tegelijkertijd spelen:

1. De Kleine Plek: De specifieke regio waar het zwarte gat instort (zeer klein).
2. Het Grote Plaatje: Het gehele universum dat eromheen uitdijt (zeer groot).

Als je probeert naar het hele universum te kijken met een microscoop, verlies je het grote plaatje. Als je naar het hele universum kijkt met een groothoeklens, kun je de kleine details van de instorting niet zien.

COSMOS lost dit op met een "Smart Zoom"-functie.
Stel je voor dat je een film kijkt van een instortende ster. Het grootste deel van het scherm toont het uitgestrekte, lege universum. Maar zodewegs de ster begint te krimpen, zoomt de camera automatisch super dicht in op die ene plek, waarbij het meer "pixels" (resolutie) toevoegt precies waar ze nodig zijn. Hierdoor kan de computer de kleine, gewelddadige instorting aan, zonder dat er een supercomputer ter grootte van een stad nodig is om alleen al de lege ruimte eromheen te berekenen.

De Ingrediënten: Wat zit er in de Simulatie?

Om de simulatie te draaien, mengt COSMOS twee hoofdingrediënten, zoals een bakker die een specifiek soort deeg maakt:

1. Het Fluïdum: Een perfect, glad fluïdum dat de materie in het vroege universum vertegenwoordigt. Het volgt eenvoudige regels (zoals een ballon die opblaast of leegloopt).
2. Het Scalaire Veld: Een spookachtig energieveld dat door de ruimte rimpelt.

Het programma lost de Einstein-vergelijkingen op (de regelset voor hoe zwaartekracht werkt) om te zien hoe deze ingrediënten met elkaar interageren. Het vraat: "Als ik begin met een specifieke bult in het universum, zal deze dan gladstrijken, of zal deze in zichzelf imploderen tot een zwart gat?"

Het "Geen-Voorbereiding"-Voordeel

Normaal gesproken is het opzetten van een natuurkundesimulatie als het bakken van een cake waarbij je eerst het beslag moet mengen (complexe wiskundige vergelijkingen oplossen) voordat je de oven überhaupt aanzet.

COSMOS is anders. Omdat het een universum simuleert dat al vol zit met fluïdum (in plaats van lege ruimte), kan het starten met het beslag dat al gemengd is. De begincondities zijn zo perfect ingesteld dat de computer geen tijd hoeft te verspillen aan het oplossen van moeilijke "elliptische" vergelijkingen om te kunnen starten. Het drukt gewoon op "run" en kijkt hoe het verhaal zich ontvouwt. Dit maakt de code lichter, sneller en gemakkelijker voor andere wetenschappers om te installeren en te gebruiken.

Wat het Papier Laat Zien (De Voorbeelden)

Het artikel bevat drie "proefritten" om te laten zien hoe de code werkt:

1. De Simpele Golf: Het simuleert een kleine, zachte rimpeling in de ruimte om te bewijzen dat de code overeenkomt met bekende, eenvoudige wiskunde.
2. De Perfecte Bol (Adiabatisch): Het simuleert een perfect ronde, instortende bubbel van de ruimte. Het artikel toont een afbeelding van de "lapse function" (een maatstaf voor hoe de tijd stroomt) en hoe de computer inzoomt op het centrum terwijl het zwarte gat vormt.
3. De Ghost Wave (Iso-curvature): Het simuleert een instorting veroorzaakt door dat "spookachtige" energieveld waar eerder over werd gesproken.

In al deze gevallen vindt de code succesvol het moment waarop een Zwart Gat wordt geboren (technisch genoemd een "apparent horizon") en brengt het de vorm ervan in kaart.

Waarom Dit Er Toe Doet

De auteurs maken niet alleen een speeltje; ze bouwen een gespecialiseerd hulpmiddel voor een specifieke taak. Hoewel er andere instrumenten bestaan voor algemene botsingen van zwarte gaten (zoals die van samensmeltende sterren), is COSMOS de gespecialiseerde moersleutel die specifiek is ontworpen voor de unieke, chaotische en uitdijende omgeving van het vroege universum.

Door deze code open en gemakkelijk toegankelijk te maken, hopen de auteurs dat andere wetenschappers het kunnen gebruiken om nieuwe vragen te stellen over de verborgen geschiedenis van het universum, wat potentieel kan verklaren waar donkere materie uit bestaat of waarom we vandaag de dag zwaartekrachtgolven waarnemen.

---

In een notendop: Dit artikel presenteert COSMOS, een gespecialiseerd computerprogramma dat fungeert als een "smart zoom"-camera voor het vroege universum. Het stelt wetenschappers in staat om te simuleren hoe minuscule rimpelingen in de ruimtetijd instorten om de allereerste zwarte gaten te vormen, waarbij het efficiënt is door de rekenkracht precies daar te richten waar de actie plaatsvindt.

Technische samenvatting

Technische Samenvatting van COSMOS: Een Numerieke Relativiteitty-code voor PBH-vorming

Probleemstelling
Primordiale zwarte gaten (PBH's) worden verondersteld te ontstaan in het vroege universum uit super-horizonale primordiale fluctuaties met niet-lineair grote initiële amplitudes, waarbij zij de stellaire evolutie omzeilen. Het simuleren van hun vorming vereist het oplossen van de Einstein-vergelijkingen onder omstandigheden waarin meerdere lengteschalen gelijktijdig voorkomen: de microscopische grootte van de instortende regio en de macroscopische schaal van de kosmologische expansie. Standaard numerieke benaderingen hebben vaak moeite met de computationele kosten om deze uiteenlopende schalen efficiënt op te lossen. Bovendien brengt het verkrijgen van initiële data voor ruimtetijden gevuld met fluïda en scalaire velden, in plaats van asymptotisch vlakke vacuümregio's, specifieke uitdagingen met betrekking tot de beperkingsvergelijkingen (constraint equations) met zich mee.

Methodologie
Het artikel introduceert COSMOS, een C++ pakket ontworpen om de Einstein-vergelijkingen in 3+1 dimensies specifiek voor PBH-vorming op te lossen. De code is een C++ vertaling en evolutie van de SACRA-code, die oorspronkelijk voor andere doeleinden is ontwikkeld.

• Fysisch Model: COSMOS modelleert materie als een perfect fluïdum met een lineaire toestandsvergelijking ($P = w\rho$) en een massaloos scalair veld. De energie-impuls-tensor is de som van de bijdragen van het fluïdum ($T^{FL}_{\mu\nu}$) en het scalaire veld ($T^{SC}_{\mu\nu}$). Het systeem lost de Einstein-vergelijkingen op ($G_{\mu\nu} = 8\pi G/c^4 T_{\mu\nu}$), de scalaire veldvergelijking ($\nabla_\mu\nabla^\mu\phi = 0$), en de fluïdum-conserveringsvergelijkingen ($\nabla_\mu T^{FL}_{\mu\nu} = 0$). De evolutie van het fluïdum volgt het schema dat is vastgesteld door Kurganov & Tadmor (2000) en Shibata & Font (2005).
• Numerieke Technieken: Om het multi-schaalprobleem aan te pakken, implementeert COSMOS:
  • Niet-Cartesiaanse schaal-op coördinaten: Om de instortende regio efficiënt op te lossen.
  • Fixed Mesh Refinement (FMR): Om resolutie waar nodig dynamisch aan te passen.
  • Parallellisatie: OpenMP wordt gebruikt om praktisch aanvaardbare computationele snelheden te bereiken.
  • Afhankelijkheden: De code is ontworpen zonder externe afhankelijkheden om installatie en uitbreiding te vergemakkelijken.
• Constructie van Initiële Data: In tegen tegenover de standaardmethoden voor asymptotisch vlakke vacuümruimtetijden die het oplossen van elliptische differentiële vergelijkingen vereisen, genereert COSMOS initiële data door langgolvige groeiende-modus oplossingen te adopteren (tot de volgende leidende orde in de expansieparameter $\epsilon = k/(aH) \ll 1$).
  • Voor adiabatische fluctuaties wordt de initiële data gekenmerkt door de kromming-perturbatie $\zeta(\vec{x})$.
  • Voor iso-kromatische perturbaties wordt deze gekenmerkt door $\Upsilon(\vec{x})$ geassocieerd met het massaloze scalaire veld.
  • De fluïdumdistributie wordt gegenereerd om de Einstein-beperkingsvergelijkingen tot op machineprecisie te voldoen, waardoor de noodzaak voor elliptische solvers wordt geëlimineerd.

Belangrijkste Bijdragen

1. Gespecialiseerde Softwarerelease: Het artikel presenteert de publieke release van COSMOS, een hulpmiddel dat specifiek is afgestemd op het simuleren van PBH-vorming, wat het onderscheidt van algemene numerieke relativiteit-codes.
2. Efficiënte Generatie van Initiële Data: De code implementeert een methode om de beperkingsvergelijkingen analytisch te voldoen via langgolvige benaderingen, waardoor de computationele overhead van elliptische solvers voor initiële data in kosmische fluïdumachtergronden wordt vermeden.
3. Documentatie en Voorbeelden: Het pakket bevat drie demonstratievoorbeelden bedoeld voor instructiedoeleinden:
   • Evolutie van een enkelvoudige sinusvormige perturbatie (validatie tegen lineaire perturbatietheorie).
   • Adiabatische sferisch symmetrische initiële fluctuaties (demonstratie van de vorming van een schijnbare horizon).
   • Sferisch symmetrische iso-kromatische perturbaties gedreven door een massaloos scalair veld.
   • Deze voorbeelden maken gebruik van een lage resolutie om functionaliteit te demonstreren in plaats van hoge-precisie resultaten.

Resultaten
Het artikel biedt visuele en kwalitatieve resultaten van de opgenomen voorbeelden:

• Lineaire Validatie: De tijdsevolutie van de spoor van de extrinsieke kromming ($trK$) voor een kleine sinusvormige fluctuatie komt overeen met de oplossing van de lineaire perturbatievergelijking.
• Horizonvorming: Simulaties van adiabatische en iso-kromatische perturbaties demonstreren succesvol de vorming van schijnbare horizonten. De code produceert databestanden (bijv. ini_all.dat) waarmee gebruikers de geometrie en de materiedistributie op het moment dat een schijnbare horizon wordt gedetecteerd, kunnen reconstrueren.
• Mesh Refinement: Visualisaties bevestigen de werking van de Fixed Mesh Refinement-lagen, waarbij duidelijke regio's van resolutie (laagste, 1e en 2e laag) rond de instortende regio worden getoond.

Betekenis en Claims
De auteurs positioneren COSMOS als een gespecialiseerd instrument om de groeiende belangstelling voor PBH's te ondersteunen, met name in de context van gravitatiegolf-observaties. Het artikel beweert dat de code een gestroomlijnde aanpak biedt voor het simuleren van niet-lineaire gravitationele dynamica in het vroege universum door:

• Een zelfvoorzienende C++ omgeving te bieden die gemakkelijk te installeren en uit te breiden is voor diverse wetenschappelijke settings zodra de broncode begrepen is.
• Een robuuste methode te bieden voor het afhandelen van het specifieke initiële beginselprobleem van kosmische fluïda zonder elliptische solvers.
• Te dienen als fundament voor eerder en toekomstig onderzoek, zoals blijkt uit het gebruik in talrijke publicaties door de auteurs en hun collaborateurs (bijv. Yoo et al., 2013; Okawa et al., 2014; Escrivà & Yoo, 2025).

Het artikel merkt bescheiden op dat de opgenomen voorbeelden voor demonstratie en instructie zijn, waarbij de resolutie opzettelijk laag is gehouden. Het verduidelijkt ook dat sommige geavanceerde functies die in eerdere publicaties van de auteurs werden gebruikt, mogelijk niet aanwezig zijn in de publieke release, wat impliceert dat het reproduceren van alle eerdere resultaten aanvullende, niet-publieke implementaties kan vereisen.