On the treatment of thermal effects in the equation of state on neutron star merger remnants

Probleemstelling

De detectie van gravitationele golven (GW) van binaire neutronenster-samensmeltingen (BNS), zoals GW170817, heeft een nieuw venster geopend voor het bestuderen van de toestand van materie bij supranucleaire dichtheden. Een cruciale uitdaging bij het interpreteren van het signaal na de samensmelting (post-merger) is de behandeling van thermische effecten in de Toestandvergelijking (EOS).
In numerieke relativiteitssimulaties worden twee hoofdmethoden gebruikt:

1. Hybride EOS: Een veelgebruikte benadering die een koude, stuksgewijs polytropische component combineert met een thermische drukterm die lijkt op een ideaal gas ($P_{th} \propto \rho^{\Gamma_{th}}$). Dit is computatie-efficiënt maar vereenvoudigt de complexe thermodynamica.
2. Getabuleerde EOS: Volledige microfysische tabellen die afhankelijkheden van temperatuur, samenstelling en dichtheid bevatten. Deze zijn nauwkeuriger maar computatie-intensiever.

Het probleem is dat het niet duidelijk is in welke mate de vereenvoudigde hybride benadering de dynamica van het overblijfsel (de hypermassieve neutronenster of HMNS) en de daaruit voortvloeiende gravitationele golven correct weergeeft, vooral op lange tijdschalen (>100 ms) waar convectie en instabiliteiten een rol spelen.

Methodologie

De auteurs hebben een reeks lange-termijnsimulaties (tot >150 ms) uitgevoerd van binaire neutronenster-samensmeltingen met behulp van de Einstein Toolkit en de ILLINOISGRMHD code.

• EOS-modellen: Er zijn vier verschillende microfysische EOS-tabels gebruikt (SLy4, DD2, HShen, LS220). Voor elke tabel werden twee representaties vergeleken:
  • De volledige getabuleerde (tab) versie met eindige temperatuur.
  • De hybride (hyb) versie, waarbij de koude component is benaderd met een stuksgewijs polytroop (7 en 10 stukken) en de thermische component als een ideaal gas met een vaste thermische index ($\Gamma_{th} = 1.8$).
• Initiële Data: Er zijn gelijkmassige binairsystemen ($M_0 = 1.4 M_\odot$) geconstrueerd. De auteurs analyseerden de verschillen in evenwichtsconfiguraties tussen de twee EOS-benaderingen en hoe deze de initiële data beïnvloeden.
• Analyse van Convectie en Stabiliteit: Om de stabiliteit van de HMNS tegen convectie te beoordelen, werden twee criteria toegepast:
  • Het Ledoux-criterium (voor thermische en samenstellingsgradiënten).
  • Het Solberg-Høiland-criterium (een generalisatie voor axiale perturbaties die de Brunt-Väisälä-frequentie ($N^2$) en de epicyclische frequentie ($K^2$) combineert). Dit is essentieel omdat rotatie de stabiliteit beïnvloedt.
• Gravitationele Golven: De GW-signaal spectra ($h_{22}$) werden geanalyseerd met behulp van Prony's methode (ESPRIT) om de frequenties van oscillatiemodi (zoals de fundamentele $f_2$-modus en inertiale modi) te extraheren.

Belangrijkste Bijdragen

1. Vergelijking van EOS-benaderingen: De studie biedt een directe, systematische vergelijking tussen volledig getabuleerde en hybride EOS-modellen in lange-termijnsimulaties, waarbij wordt aangetoond dat de keuze van de EOS de evolutie van het overblijfsel fundamenteel beïnvloedt.
2. Validatie van Inertiale Modi: De auteurs bevestigen dat inertiale modi (die worden aangedreven door convectie) ook worden opgewekt in simulaties met volledig getabuleerde EOS, niet alleen in hybride modellen. Dit ondersteunt eerdere bevindingen maar met een meer consistente thermische behandeling.
3. Kritische Analyse van Temperatuurdefinitie: In Appendix A wordt aangetoond dat de gebruikelijke definitie van temperatuur in hybride modellen ($T_{th} \propto \epsilon_{th}$) de werkelijke temperatuur in de degeneratieve regime (hoge dichtheid) significant onderschat.
4. Nieuwe Quasi-universele Relaties: Er worden nieuwe empirische relaties afgeleid voor de post-merger GW-frequenties, waarbij wordt aangetoond dat afwijkingen van bestaande relaties significant worden in de late fase van de samensmelting.

Resultaten

• Dynamica en Instabiliteit:

  • Hybride modellen produceren over het algemeen compactere overblijfsels met een hogere kans op vroege ineenstorting tot een zwart gat (BH), vooral bij "zachte" EOS (zoals LS220).
  • Getabuleerde modellen tonen een betere thermische drukondersteuning, waardoor sommige overblijfsels (zoals bij LS220) stabiel blijven langer dan 100 ms, terwijl hun hybride tegenhangers instorten.
  • Er zijn duidelijke verschillen in de maximale dichtheid en de evolutie van de lapse-functie tussen de twee benaderingen.

• Gravitationele Golven en Spectra:

  • In de vroege fase (<10 ms) lijken de GW-signaalspectra (vooral de $f_{2,i}$-modus) tussen de modellen redelijk overeen te komen.
  • In de late fase (>30-40 ms) ontstaan significante verschillen. De frequenties van de $f_2$-modus en inertiale modi verschuiven tussen de hybride en getabuleerde modellen (verschillen tot wel 390 Hz voor HShen).
  • Inertiale modi worden opgewekt door convectieve patronen die worden aangedreven door thermische stratificatie en differentieel draaiing. Deze modi hebben lagere frequenties dan de fundamentele kwadrupolaire modus en zijn potentieel detecteerbaar met derde-generatie detectors (zoals Einstein Telescope en Cosmic Explorer).

• Convectieve Stabiliteit:

  • Analyse met het Solberg-Høiland-criterium toont aan dat er lokale, maar aanhoudende convectieve instabiliteiten ontstaan in de HMNS.
  • Hoewel hybride modellen soms meer convectieve regio's lijken te vertonen, zijn de patronen in getabuleerde modellen complexer en worden ze beter gekarakteriseerd door de nauwkeurigere thermodynamica.
  • De rotatieprofielen verschillen: hybride modellen vertonen vaak steilere rotatiegradiënten, wat kan leiden tot schuifinstabiliteiten.

• Quasi-universele Relaties:

  • De relatie tussen de piekfrequentie bij samensmelting en de getijdenpolariseerbaarheid ($\kappa^T_2$) blijft grotendeels geldig.
  • Echter, voor de frequenties in de late post-merger fase wijken de resultaten van hybride modellen significant af van de getabuleerde resultaten en van eerdere quasi-universele relaties. Dit suggereert dat het gebruik van hybride EOS voor het afleiden van EOS-eigenschappen uit late GW-signalen tot fouten kan leiden.

Betekenis en Conclusie

De studie concludeert dat de behandeling van thermische effecten in de EOS van cruciaal belang is voor het begrijpen van de langetermijnevolutie van neutronenster-samensmeltingen.

• Nauwkeurigheid: De hybride benadering, hoewel efficiënt, faalt in het volledig vastleggen van de complexe thermodynamica, wat leidt tot onnauwkeurigheden in de voorspelling van de levensduur van het overblijfsel, de frequentie van GW-modi en de stabiliteit tegen instorting.
• Toekomstige Detecties: Voor de interpretatie van signalen van toekomstige, gevoelige detectors (3e generatie) is het essentieel om volledig getabuleerde EOS te gebruiken of hybride modellen te kalibreren, omdat de late-fase inertiale modi een waardevolle bron van informatie kunnen zijn over de thermische en rotatie-eigenschappen van neutronensterren (GW-asteroseismologie).
• Temperatuurdefinitie: De auteurs waarschuwen dat de standaarddefinitie van temperatuur in hybride modellen de werkelijke temperatuur in degeneratief materie onderschat, wat de fysica van convectie en koeling beïnvloedt.

Kortom, voor een betrouwbare interpretatie van toekomstige GW-observaties van BNS-mergers is een consistente behandeling van thermische effecten via getabuleerde EOS noodzakelijk.