Numerical simulations of oscillating and differentially rotating neutron stars

De auteurs breiden de pseudospectrale code ROXAS uit om oscillerende, differentieel roterende neutronensterren dynamisch te simuleren, waarmee ze voor het eerst de frequenties van niet-axiale modi in deze configuraties bepalen en zo een stap zetten naar realistischere modellering van post-merger overblijfselen en hun zwaartekrachtsgolven.

De kern

Hoe sterren dansen: Een simpele uitleg van het onderzoek over neutronensterren

Stel je voor dat je een dansvloer hebt met een enorme, zware danser: een neutronenster. Deze ster is zo zwaar dat een theelepel ervan zou wegen als een berg, en hij draait razendsnel om zijn as, soms wel honderden keren per seconde.

Wanneer twee van deze sterren botsen (zoals bij een kosmische danspartij), ontstaat er een nieuw, nog zwaarder monster. Dit nieuwe object draait niet als een strakke ijskoningin die met één snelheid draait. Nee, het gedraagt zich meer als een slappe soep in een blender: de binnenkant draait sneller dan de buitenkant. Dit noemen wetenschappers "differentiële rotatie".

Dit artikel vertelt over een nieuwe computercode genaamd ROXAS die wetenschappers hebben verbeterd om precies te begrijpen hoe deze "slappe soep" trilt en wat voor geluid (of beter gezegd: zwaartekrachtgolven) dat maakt.

Hier is hoe het werkt, vertaald naar alledaagse taal:

1. Het probleem: De dans is te complex

Vroeger konden computers alleen simuleren hoe een ster draait als hij overal even snel draait (zoals een strakke balletdanser). Maar na een botsing is de ster chaotisch: de binnenkant draait sneller dan de buitenkant.

• De analogie: Stel je voor dat je een deegbol rolt. Als je hem strak vasthoudt, draait hij als één stuk. Maar als je de binnenkant harder rolt dan de buitenkant, wordt het deeg rekbaar en onvoorspelbaar.
• Het doel: De wetenschappers wilden hun computerprogramma (ROXAS) leren om deze "rekbare" dans te simuleren, zodat ze kunnen voorspellen wat er gebeurt als zo'n ster trilt.

2. De oplossing: Een slimme rekenmethode

Het programma ROXAS is als een super-snel rekenmachine die de bewegingen van de ster berekent.

• De truc: In plaats van alles vanaf nul te berekenen (wat duizenden supercomputers zou kosten), gebruiken ze een slimme benadering. Ze kijken naar de ster alsof hij op een vlakke vloer staat (een wiskundige truc genaamd "conformal flatness"), maar houden wel rekening met de zwaartekracht.
• Het resultaat: Dankzij deze truc kan de simulatie draaien op een normale kantoorcomputer, in plaats van een enorme supercomputer. Het is alsof je een complexe 3D-film kunt draaien op een oude laptop, terwijl anderen daarvoor een datacenter nodig hebben.

3. Wat hebben ze ontdekt? (De verrassingen)

Toen ze de sterren lieten trillen, kwamen ze twee belangrijke dingen tegen:

• De "Geest" die verdween:
  In eerdere, minder nauwkeurige berekeningen (waar ze de zwaartekracht van de ster zelf negeerden, de zogenaamde "Cowling-benadering"), zagen ze een dubbel geluid: een hoofdton en een "tweede" ton.

  • De analogie: Het was alsof je een gitaar bespeelde en dacht dat er twee snaren tegelijk trilden.
  • De ontdekking: Toen ze de volledige zwaartekracht meerekenden (de dynamische ruimtetijd), bleek die "tweede ton" een fictie te zijn. Het was een artefact van de rekenmethode, net als een spookbeeld in een spiegel. In de echte natuur (zoals de computer het nu simuleert) is er maar één hoofdton.

• Nieuwe danspasjes:
  Voor het eerst hebben ze de trillingen van sterren berekend die niet alleen rondom hun as draaien, maar ook een beetje "waggelen" (niet-axiale trillingen).

  • De analogie: Stel je voor dat de ster niet alleen rond zijn as draait, maar ook een beetje schommelt als een tol die bijna omvalt. Ze hebben nu de exacte frequentie van die schommeling gemeten voor sterren die snel en ongelijkmatig draaien.

4. Waarom is dit belangrijk voor ons?

Je vraagt je misschien af: "Wat heb ik hieraan?"

• De kosmische radio: Wanneer deze sterren trillen, sturen ze zwaartekrachtgolven uit. Dit zijn rimpelingen in de ruimte zelf, zoals golven in een meer.
• De toekomst: Huidige telescopen (zoals LIGO) kunnen deze hoge tonen nog niet goed horen. Maar in de toekomst, met nieuwe telescopen (zoals de "Einstein Telescope"), zullen we deze geluiden kunnen opvangen.
• De boodschap: Als we weten hoe deze sterren klinken (zoals in dit artikel berekend), kunnen we bij het horen van een signaal zeggen: "Ah, dit is een ster die zo en zo snel draait!" Dit helpt ons te begrijpen waaruit sterren zijn opgebouwd en wat er gebeurt als ze botsen.

Samenvatting

De auteurs hebben hun computerprogramma ROXAS opgefrist zodat het nu kan simuleren hoe neutronensterren die ongelijkmatig draaien, trillen. Ze hebben bewezen dat sommige eerdere resultaten "spooktonen" waren en hebben voor het eerst nieuwe trillingen gemeten. Het mooiste is dat ze dit doen met een lichtgewicht programma dat snel werkt, waardoor ze veel meer scenario's kunnen testen dan ooit tevoren.

Het is alsof ze een nieuwe set brillen hebben gemaakt om beter te kijken naar de dans van de zwaarste objecten in het universum.

Technische samenvatting

Probleemstelling

De directe detectie van zwaartekrachtgolven (GW) van samensmeltingen van dubbele neutronensterren (zoals GW170817) heeft de weg vrijgemaakt voor het bestuderen van de overblijfselen van deze gebeurtenissen. Deze overblijfselen worden vaak beschreven als hypermassieve neutronensterren (HMNS), die gedurende enkele tienden van een seconde worden ondersteund door differentieel draaiing voordat ze instorten tot een zwart gat. Tijdens deze fase stralen ze zwaartekrachtgolven uit in het kilohertz-band.

De structuur, stabiliteit en het oscillatiespectrum van deze sterren worden sterk beïnvloed door hun differentieel draaiingsprofiel. Huidige detectoren hebben echter onvoldoende gevoeligheid in het kHz-bereik om deze signalen direct waar te nemen; toekomstige generaties (zoals Einstein Telescope en Cosmic Explorer) zullen dit wel kunnen. Om de interne structuur en de toestandsvergelijking (EoS) van neutronensterren uit toekomstige signalen te kunnen halen, is een nauwkeurig begrip nodig van hoe rotatie en andere parameters de uitgezonden golven vormen. Bestaande numerieke codes zijn vaak computergewijs zeer zwaar, wat parametrische studies bemoeilijkt.

Methodologie

De auteurs hebben de bestaande pseudospectrale code ROXAS (Relativistic Oscillations of non-aXisymmetric neutron stArS) uitgebreid om dynamische evoluties van gestoorde, differentieel draaiende neutronensterren mogelijk te maken.

• Formalisme: De code gebruikt een 3+1-decompositie van de ruimtetijd met de Conformal Flatness Condition (CFC) als benadering voor de volledige Algemene Relativiteitstheorie. Dit vereenvoudigt de vergelijkingen aanzienlijk terwijl de nauwkeurigheid behouden blijft.
• Hydrodynamica: De ster wordt gemodelleerd als een perfect vloeistof met een barotrope toestandsvergelijking (polytroop: $p = \kappa n_B^\gamma$). De evolutie wordt beschreven in termen van primitieve variabelen (dichtheid, druk, snelheid).
• Differentieel draaien: Een cruciale wijziging is de implementatie van differentieel draaien via een rotatiewet (gebaseerd op het Komatsu-Eriguchi-Hachisu (KEH) profiel). Dit vereist een "well-balanced" formulering waarbij de hydrodynamische vergelijkingen worden herschreven rondom een evenwichtstoestand. Hierdoor ontstaan extra termen in de evolutievergelijkingen die specifiek zijn voor differentieel draaien (zoals integralen over de hoeksnelheid $\Omega$).
• Numerieke Implementatie:
  • De code maakt gebruik van de LORENE-infrastructuur om evenwichtsconfiguraties te genereren.
  • Er worden twee roosters gebruikt: een hydrodynamisch rooster (sferische kern en schil) en een metriekrooster (meerdere sferische schillen met een gecomprimeerd extern domein).
  • De simulaties worden gestart met kleine perturbaties (in snelheid of enthalpie) om oscillaties op te wekken.
  • Er wordt zowel in de Cowling-benadering (statische metriek) als in dynamische ruimtetijd (CFC) gesimuleerd.

Belangrijkste Bijdragen

1. Code-uitbreiding: ROXAS is voor het eerst uitgebreid om differentieel draaiende neutronensterren dynamisch te evolueren, inclusief de benodigde wiskundige aanpassingen voor de hydrodynamische vergelijkingen.
2. Validatie: De code is grondig gevalideerd door vergelijking met bestaande literatuur (o.a. Stergioulas et al. 2004, Krüger et al. 2010, Dimmelmeier et al. 2006) voor zowel as-symmetrische als niet-as-symmetrische modi.
3. Ontmaskering van een artefact: De auteurs tonen aan dat de "secundaire fundamentele modus" (FII) die vaak wordt waargenomen in simulaties onder de Cowling-benadering, een artefact is van deze benadering en niet voorkomt in dynamische ruimtetijden.
4. Nieuwe resultaten: Voor het eerst worden frequentiewaarden voor niet-as-symmetrische modi in differentieel draaiende configuraties (binnen CFC) gepubliceerd.

Resultaten

• Overeenkomst met literatuur: De as-symmetrische modi en de resultaten onder de Cowling-benadering tonen uitstekende overeenkomst (fouten vaak < 2%) met eerder gepubliceerde resultaten.
• Het FII-artefact: Simulaties in de Cowling-benadering tonen twee fundamentele pieken (F en FII) rond 2,6 kHz. Wanneer echter de volledige dynamische ruimtetijd (CFC) wordt opgelost, verdwijnt de FII-piek en blijft er slechts één fundamentele modus over (rond 1,3 kHz voor de B4-configuratie). Dit bevestigt dat FII een numeriek artefact is van het verwaarlozen van de metriekdynamica bij differentieel draaien.
• Niet-as-symmetrische modi: De code levert nieuwe frequentiegegevens voor niet-as-symmetrische modi ($2f_2$ en $2f_{-2}$). Bij hoge rotatiesnelheden (vanaf model B7) wordt de $2f_2$-frequentie negatief, wat aangeeft dat deze modus retrograde draait in het meedraaiende referentiestelsel van de ster (CFS-instabiliteit).
• Vergelijking met Dimmelmeier et al. (2006): Voor de dynamische CFC-simulaties is er over het algemeen goede overeenkomst met Dimmelmeier et al., hoewel er enkele discrepanties zijn bij de fundamentele frequentie (F) voor de snelst draaiende sterren. De auteurs wijzen op mogelijke oorzaken zoals verschillende resoluties en simulatietijden, maar benadrukken dat CFC betrouwbaar is voor dit type studie.
• Efficiëntie: ROXAS is een "lightweight" code die aanzienlijk sneller is dan zwaardere codes zoals CoCoNuT (ongeveer 5x sneller voor 1D, en theoretisch 25x en 125x sneller voor respectievelijk as-symmetrische en niet-as-symmetrische simulaties door de gebruikte symmetrieën).

Significantie

Dit werk vormt een belangrijke stap in de modellering van post-samensmeltingsremnants. Door differentieel draaien te integreren in een efficiënte pseudospectrale code, maken de auteurs parametrische studies mogelijk die eerder te rekenintensief waren. Het onthullen van het artefact van de secundaire fundamentele modus is cruciaal voor de interpretatie van toekomstige zwaartekrachtgolfsignalen. De code biedt nu een robuust en snel instrument om de oscillatiespectra van neutronensterren te voorspellen, wat essentieel zal zijn voor de analyse van data van toekomstige detectoren zoals de Einstein Telescope en Cosmic Explorer.