Convective stability analysis of massive neutron stars formed in binary mergers

Technische Samenvatting: Convectieve Stabiliteitsanalyse van Massieve Neutronensterren gevormd in Binaire Mergers

Probleemstelling
De post-merger evolutie van binaire neutronenster (BNS) mergers wordt beheerst door de dynamica van de resulterende massieve neutronenster (NS). Een cruciaal aspect van deze evolutie is de potentiële ontwikkeling van hydrodynamische instabiliteiten, met name convectieve instabiliteit, die grootschalige fluïdumstromingen kan aandrijven, de thermische structuur kan veranderen en gravitatiegolfmodi (GW-modi) kan exciteren. Eerdere studies hebben tegenstrijdige conclusies opgeleverd: sommige simulaties die het Schwarzschild-criterium gebruikten, suggereerden dat massieve NS's grootschalige convectieve instabiliteit ontwikkelen, wat vervolgens observeerbare inertiële modi exciteert [74, 75]. Daarentegen vonden andere studies die neutrino-straling incorporeerden een stabiele stratificatie [55]. Een belangrijke beperking in eerdere analyses is de frequente verwaarlozing van rotationele herstellende krachten. Gezien het feit dat merger-remnants snelle en differentiële rotatie vertonen, is de rotationele bijdrage aan de stabiliteit (vaak beschreven door het Solberg-Høiland criterium) essentieel. Dit artikel beoogt deze discrepanties op te lossen door volledig algemeen-relativistische hydrodynamische simulaties uit te voeren en een uitgebreid stabiliteitscriterium toe te passen dat zowel rekening houdt met thermische drijfvermogen als rotationele effecten.

Methodologie
De auteurs maken gebruik van de SACRA-MPI numerieke relativiteitscode om de merger van gelijk massa BNS-systemen te simuleren over een duur van ongeveer 100 ms post-merger.

• Toestandsvergelijking (EOS): Drie stuksgewijze polytopische EOS's (APR4, SLy4, MPA1) worden gebruikt, aangevuld met een $\Gamma$-wet thermisch component ($\Gamma_{\text{th}} = 1.8$) om schokverhitting te modelleren.
• Numerieke Opstelling: Simulaties maken gebruik van een Cartesiaans genest rooster met een 2:1 verfijning en spiegelymmetrie over het orbitale vlak. Drie resoluties worden getest voor het APR4-model om convergentie te waarborgen.
• Stabiliteitsanalyse: De auteurs leiden en passen de relativistische Solberg-Høiland stabiliteitscriteria toe (Criterium I en Criterium II) voor hete, differentiële roterende sterren. Deze criteria incorporeren de interactie tussen thermisch drijfvermogen (Schwarzschild-discriminant) en rotationele herstellende krachten (specifieke impulsmomentgradiënten). De analyse wordt uitgevoerd op de hoekgemiddelde profielen van specifieke entropie en specifiek impulsmoment.
• Modianalyse: De studie analyseert de excitatie van fluïdummodi en GW-emissie, waarbij specifiek wordt gefocust op de $m=1$ één-armige modus en de potentiële excitatie van $m=2$ inertiële modi volgend op de demping van de kwadrupolaire $f$-modus.

Belangrijkste Bijdragen

1. Afleiding van Relativistische Stabiliteitscriteria: Het artikel biedt een rigoureuze toepassing van de relativistische Solberg-Høiland criteria (Vergelijkingen 46 en 47) op post-merger massieve NS's, waarbij expliciet rekening wordt gehouden met de rotationele herstellende kracht die werd verwaarloosd in studies die uitsluitend vertrouwden op het Schwarzschild-criterium.
2. Onafhankelijk Simulatiekader: De auteurs voeren een onafhankelijke analyse uit met numerieke opstellingen (EOS-keuzes, resoluties, hydrodynamische schema's) die verschillen van eerdere tegenstrijdige studies, om directe één-op-één vergelijkingen te vermijden en een zelfconsistente onderzoek van de fysica te waarborgen.
3. Symmetrie- en Impulsmomentbehoudstests: De studie introduceert een symmetrie-behoudende test om te verifiëren dat de numerieke code de $\pi$-symmetrie niet kunstmatig verbreekt. Bovendien onderzoeken zij de correlatie tussen de groei van de $m=1$ modusamplitude en de schending van het lineaire impulsmomentbehoud, een kritische controle om onderscheid te maken tussen fysieke instabiliteiten en numerieke artefacten.

Resultaten

• Convectieve Stabiliteit: De simulaties tonen geen bewijs van grootschalige convectieve instabiliteit. Zowel de hoekgemiddelde specifieke entropie als het specifieke impulsmoment nemen naar buiten toe toe binnen de massieve NS's. Hoewel lokale regio's onder het Schwarzschild-criterium alleen instabiel kunnen lijken, stabiliseren deze regio's door de inclusie van rotationele effecten (met name de toename van het specifieke impulsmoment naar buiten toe). De algemene stabiliteitscriteria (Criterium I en II) worden voor het grootste deel van de ster voldaan.
• Inertiële Modi: In tegen tegenstelling tot de bevindingen in Refs. [74, 75], wordt er geen excitatie van observeerbare inertiële modi gedetecteerd na de demping van de kwadrupolaire $f$-modi. Het GW-signaal blijft gedurende de simulatie gedomineerd door de $f$-modus, zonder de opkomst van discrete inertiële modi op kilohertz-frequenties.
• $m=1$ Instabiliteit: Een $m=1$ één-armige modus wordt waargenomen terwijl deze groeit en verzadigt, wat consistent is met eerdere literatuur. De auteurs vinden echter een sterke correlatie tussen de groei van deze modusamplitude en de schending van het behoud van lineair impulsmoment. Dit suggereert dat de excitatie van de één-armige modus in deze specifieke simulaties een numerieke oorsprong kan hebben in plaats van een puur fysieke instabiliteit te zijn.
• Thermische en Rotationele Profielen: De massieve NS bereikt een quasi-evenwichtstoestand binnen 20–30 ms. De kern is koud en roteert traag, terwijl de envelop wordt verhit door spiraalarm-interacties en differentiële rotatie vertoont. Het specifieke impulsmoment neemt naar buiten toe toe, wat een sterke stabiliserende centrifugale kracht biedt.

Betekenis en Claims
Het artikel beweert dat de inclusie van rotationele herstellende krachten de stabiliteitsbeoordeling van post-merger massieve NS's fundamenteel verandert. De auteurs stellen dat hun resultaten, die een stabiele configuratie laten zien, incompatibel zijn met eerdere studies die grote convectie en daaropvolgende inertiële modus-excitatie voorspelden. Zij argumenteren dat het gebrek aan convectieve instabiliteit in hun simulaties de afwezigheid van geëxciteerde inertiële modi verklaart.

Wat betreft de $m=1$ modus nemen de auteurs een gematigd standpunt in door te stellen dat ze "niet de mogelijkheid kunnen uitsluiten dat de excitatie van de één-armige modus een numerieke oorsprong heeft" vanwege de correlatie met het niet-behoud van lineair impulsmoment. Zij benadrukken dat hoewel de code symmetrie-behoudende tests doorstaat, het behoud van lineair impulsmoment in simulaties met mesh-verfijning een uitdaging blijft die verder onderzoek vereist om de fysieke aard van de één-armige instabiliteit definitief vast te stellen.

De studie concludeert dat toekomstig onderzoek dat microfysica (samenstellingsgradiënten) en magnetische velden incorporeert noodzakelijk is om een definitief oordeel te vellen, maar binnen de reikwijdte van hun huidige hydrodynamische simulaties lijken de massieve NS's gevormd in BNS-mergers convectief stabiel te zijn.