Maximal mass of neutron stars constrained by neutron star observations

Door een Bayesiaans wegingskader te hanteren dat multimessenger-observaties (waaronder GW170817, NICER en kandidaat-compacte objecten) integreert met hybride toestandsvergelijkingen, concludeert deze studie dat de maximale neutronenstermassa robuust beperkt is tot ongeveer 2,2–2,3 zonnemassa's, terwijl de bijbehorende straal sterker afhankelijk is van het onderliggende hadronische model en doorgaans rond de 12 km ligt.

De kern

Stel je voor dat het heelal gevuld is met een mysterieus, superdicht materiaal dat "neutronenstermateriaal" wordt genoemd. Het is zo zwaar dat een enkele theelepel evenveel zou wegen als een berg. Al lang proberen natuurkundigen de "spelregels" voor dit materiaal te achterhalen—specifiek, hoe stijf of kneedbaar het is. Deze set regels wordt de Toestandsvergelijking (EOS) genoemd.

De grote vraag die dit artikel probeert te beantwoorden is: Wat is het absolute zwaarste gewicht dat een neutronenster kan bereiken voordat het instort tot een zwart gat?

Hier is het verhaal van hoe de auteurs deze puzzel oplossen, uitgelegd in eenvoudige bewoordingen:

1. De Twee Startpunten (De Recepten)

Om de regels te achterhalen, begonnen de wetenschappers met twee verschillende "recepten" voor hoe dit dichte materiaal zich gedraagt bij lagere dichtheden. Denk hierbij aan twee verschillende theorieën over hoe de ingrediënten mengen:

• Recept A (SFHo): Een "zachter" recept, wat betekent dat het materiaal wat makkelijker te comprimeren is.
• Recept B (DD2): Een "stijver" recept, wat betekent dat het materiaal meer weerstand biedt tegen comprimeren.

Ze wisten dat deze recepten goed werkten aan het "begin" van de dichtheidsschaal, maar ze wisten niet wat er gebeurde bij de extreme, superhoge dichtheden die in het centrum van een neutronenster voorkomen. Om het gat op te vullen, gebruikten ze een wiskundige "brug" om hun recepten te verbinden met wat we weten over deeltjesfysica bij de hoogst mogelijke energieën.

2. Het Detectivewerk (Het Gebruik van Echte Aanwijzingen)

In plaats van alleen maar te gissen, traden de auteurs op als detectives. Ze namen hun twee recepten en testten ze tegen echte aanwijzingen verzameld door telescopen en zwaartekrachtsgolf-detectoren. Ze gebruikten een speciale statistische methode (genaamd Bayesiaanse weging) om te zien welke versies van hun recepten de test doorstonden.

Hier zijn de aanwijzingen die ze gebruikten:

• De "Grote Crash" (GW170817): Toen twee neutronensterren tegen elkaar botsten, stuurden ze rimpelingen door de ruimte. De manier waarop deze rimpelingen zich gedroegen, vertelde de wetenschappers hoe "kneedbaar" de sterren waren.
• De "Zaklamp" (NICER): Een ruimtetelescoop maakte foto's van hete vlekken op draaiende neutronensterren. Door te meten hoe groot de sterren eruit zagen en hoe zwaar ze waren, kregen ze een directe verhouding tussen grootte en gewicht.
• De "Lichte" Kandidaat (HESS J1731–347): Een zeer klein, licht object dat een neutronenster zou kunnen zijn.
• De "Zware" Kandidaat (GW190814): Een mysterieus object dat zwaarder is dan de meeste neutronensterren, maar lichter dan de meeste zwarte gaten. De wetenschappers vroegen zich af: Zou dit eigenlijk een superzware neutronenster kunnen zijn?

3. De Resultaten: Wat de Aanwijzingen Ze Vertelden

De wetenschappers voerden hun twee recepten door deze aanwijzingen en keken naar de resultaten.

Het Gewichtsbeperking (Maximale Massa):

• De Verrassing: Het maakte niet veel uit welk startrecept (Zacht of Stijf) ze gebruikten. De echte aanwijzingen waren zo sterk dat ze beide recepten dwongen om overeen te komen met hetzelfde antwoord.
• Het Vonnis: Toen ze de meest betrouwbare aanwijzingen gebruikten (de "Grote Crash" en de "Zaklamp"), is het maximale gewicht dat een neutronenster kan dragen ongeveer 2,2 tot 2,3 keer de massa van onze Zon.
• De "Zware" Twist: Als ze aannemen dat dat mysterieuze zware object (GW190814) een neutronenster is, springt de limiet omhoog naar ongeveer 2,6 tot 2,7 keer de massa van de Zon. Dit creëert echter een conflict met de "kneedbaarheids"-aanwijzingen van de Grote Crash, wat het een lastige situatie maakt.

De Groottebeperking (Straal):

• Het Verschil: In tegenstelling tot het gewicht, hing de grootte van de ster wel af van welk startrecept ze gebruikten.
• Het Vonnis: Het "Zachte" recept voorspelde een straal van ongeveer 11,8 km, terwijl het "Stijve" recept ongeveer 12,4 km voorspelde.
• Het Sweet Spot: Wanneer alle beste aanwijzingen worden gecombineerd, is de meest waarschijnlijke grootte voor deze sterren ongeveer 12 kilometer (plus of minus 1 km).

4. Het Grote Plaatje

Het artikel concludeert dat door te kijken naar de "eindpunten" (de zwaarste en grootste mogelijke sterren) en gebruik te maken van een mix van echte astronomische data, we de regels voor het dichtste materiaal van het heelal kunnen verengen.

• Het Gewicht: Het heelal lijkt een "snelheidslimiet" te hebben voor hoe zwaar een neutronenster kan zijn, comfortabel rond de 2,2 tot 2,3 zonsmassa's. Dit komt overeen met de zwaarste neutronenster die we tot nu toe daadwerkelijk hebben gezien.
• De Grootte: Ze zijn ongeveer zo groot als een kleine stad, ongeveer 12 km in doorsnee.
• De Kernboodschap: De waarnemingen uit de echte wereld (de aanwijzingen) zijn veel krachtiger dan de theoretische startgissingen. Welke theorie je ook als startpunt neemt, de data van de sterren zelf dwingt het antwoord om te convergeren naar dezelfde cijfers.

Kortom, het heelal heeft ons een zeer duidelijk antwoord gegeven: Neutronensterren kunnen ongelooflijk zwaar worden, maar er is een harde bovengrens, en ze zijn verrassend klein voor hoeveel ze wegen.

Technische samenvatting

1. Probleemstelling

De centrale uitdaging in de kern- en zware-ionenfysica is het bepalen van de toestandsvergelijking (EOS) van sterk wisselwerkende materie bij intermediaire dichtheden en chemische potentialen. Waar materie bij lage dichtheid wordt beperkt door kernexperimenten en rooster-QCD, en materie bij hoge dichtheid door perturbatieve QCD (pQCD), blijft het gebied dat relevant is voor de kernen van neutronensterren (enkele malen de nucleaire verzadigingsdichtheid, $\rho_0$) slecht begrepen.

Het specifieke probleem dat wordt aangepakt, is het bepalen van de maximale massa ($M_{TOV}$) en de bijbehorende straal ($R_{TOV}$) van neutronensterren. Deze eindpunten van de massa-straalreeks ($M(R)$) zijn cruciaal voor het begrijpen van de stijfheid van de EOS bij supra-nucleaire dichtheden. De auteurs beogen te kwantificeren hoe huidige multimessenger-waarnemingen deze eindpunten beperken en in welke mate de resultaten afhankelijk zijn van de keuze van het onderliggende hadronische basismodel.

2. Methodologie

A. Constructie van Hybride EOS's

De auteurs construeren families van causale hybride EOS's door materie met lage dichtheid (hadronisch) te verbinden met materie met hoge dichtheid (quark):

• Regime met lage dichtheid: Twee representatieve hadronische modellen worden gebruikt als basissen:
  • SFHo: Een relativistisch gemiddeld-veldmodel (zachtere EOS, oncompressibiliteit $K=245$ MeV).
  • DD2: Een stijvere hadronische EOS.
• Regime met hoge dichtheid: Een uitgebreid Lineair Sigma-model (eLSM) gebaseerd op $U(3) \times U(3)$ chiraal symmetrie, met inbegrip van constituent-quarks en meson-nonetten. Dit model is beperkt om te overeenkomen met rooster-QCD-thermodynamica bij lage chemische potentialen en pQCD-resultaten bij hoge chemische potentialen.
• Interpolatie: Om de twee regimes te overbruggen, wordt een polynoominterpolatie van de vijfde graad toegepast op de energiedichtheid als functie van de baryondichtheid. Dit zorgt voor thermodynamische consistentie en continuïteit van de druk en de geluidssnelheid ($c_s$) over het overgangsgebied heen, waardoor grote fase-overgangen van de eerste orde worden vermeden.
• Asymptotische Beperking: De quark-meson EOS moet soepel aansluiten op pQCD-resultaten bij asymptotisch hoge dichtheden ($\mu \approx 2.6$ GeV).

B. Bayesiaans Kader

De studie maakt gebruik van een Bayesiaans wegingskader om de waarschijnlijkheidsverdelingen van de eindpunten van de $M(R)$-reeks te analyseren.

1. Prior: Een breed ensemble van causale hybride EOS's wordt gegenereerd dat een reeks stijfheden bij supra-nucleaire dichtheden beslaat.
2. Likelihood-weging: Het ensemble wordt gewogen tegen waarnemingsdata om posterior-waarschijnlijkheidsverdelingen voor $M_{TOV}$ en $R_{TOV}$ te genereren.

C. Waarnemingsbeperkingen

De volgende beperkingen worden sequentieel en in combinatie toegepast:

1. pQCD-beperking: Zorgt ervoor dat de EOS nadert naar het N3LO pQCD-referentiepunt ($\mu=2.6$ GeV, $\rho=6.47$ fm$^{-3}$, $p=3823$ MeV/fm$^3$) terwijl causaliteit wordt gehandhaafd ($c_s < 1$).
2. GW170817: Beperking op getijdenvervorming ($\tilde{\Lambda} < 720$) en een strengere gecombineerde beperking ("$\Lambda_{190}$": $70 < \Lambda_{1.4} < 580$ en $9.1 < R_{1.4} < 12.8$ km).
3. NICER: Massa-straalmetingen voor vijf pulsars (J0030+0451, J0740+6620, J0614–3329, J1231–1411, J0437–4715) met een precisie van $\sim 10\%$.
4. HESS J1731–347: Een kandidaat voor een neutronenster met zeer lage massa ($M \approx 0.77 M_\odot$), behandeld met voorzichtigheid vanwege interpretatieonzekerheden.
5. GW190814: Het object in de "massagap" ($M \approx 2.59 M_\odot$). De auteurs testen de hypothese dat dit object een neutronenster is door een Gaussische beperking toe te passen die gecentreerd is op $2.59 M_\odot$.

3. Belangrijkste Bijdragen

• Gecombineerde Diagnose: Demonstreert dat de waarschijnlijkheidsverdelingen van de eindpunten van de $M(R)$-reeks fungeren als een gevoelige, complementaire diagnose voor het beperken van het gedrag van de EOS bij hoge dichtheid binnen een multimessenger-kader.
• Onafhankelijkheid versus Afhankelijkheid van de Basis: Splits systematisch de invloed van waarnemingsdata los van de theoretische keuze van het hadronische basismodel (SFHo versus DD2).
• Identificatie van Spanning: Identificeert een specifieke spanning tussen de stijve DD2 EOS, de hypothese van de neutronenster in de "massagap", en de strikte beperkingen op getijdenvervorming ($\Lambda_{190}$).

4. Belangrijkste Resultaten

Verdelingen van de Maximale Massa ($M_{TOV}$)

• Dominantie van Waarnemingen: De piek van de $M_{TOV}$-verdeling wordt voornamelijk gedreven door waarnemingsbeperkingen en niet door de keuze van de hadronische basiseos (SFHo of DD2).
• Progressie van Beperkingen:
  • Alleen pQCD: Levert een piek iets boven $1.5 M_\odot$ op.
  • GW170817/HESS: Geen significante verschuiving in de piekpositie.
  • NICER: Verplaatst de piek naar 2.1–2.2 $M_\odot$ en verscherpt de verdeling.
  • Massagap (GW190814): Verplaatst de piek significant naar 2.6–2.7 $M_\odot$.
• Robuuste Conclusie: Bij toepassing van de meest robuuste beperkingen (pQCD + GW170817 + NICER) convergeren beide EOS-modellen naar een voorkeursmaximale massa van 2.2–2.3 $M_\odot$. Dit is consistent met de waargenomen massa van PSR J0952–0607 ($2.35 \pm 0.17 M_\odot$).

Verdelingen van de Straal ($R_{TOV}$)

• Afhankelijkheid van de Basis: In tegenstelling tot de massa, toont de straalverdeling een sterke afhankelijkheid van de onderliggende hadronische EOS.
  • SFHo+eLSM: Piekt rond 11.8 km (met robuuste beperkingen).
  • DD2+eLSM: Piekt rond 12.4 km (met robuuste beperkingen).
• Algemene Voorkeur: Gecombineerde beperkingen geven de voorkeur aan een straalbereik van $12 \pm 1$ km.
• Impact van Getijdenvervorming: Het toepassen van de strengere $\Lambda_{190}$-beperking versmalt de straalverdelingen en verlaagt de voorkeurswaarden, aangezien lagere getijdenvervorming correleert met kleinere stralen.

De Spanning rond de "Massagap"

• Als het GW190814-object wordt geïnterpreteerd als een neutronenster, faalt de op DD2 gebaseerde EOS om een consistente beschrijving te bieden wanneer deze wordt gecombineerd met de $\Lambda_{190}$-getijdenbeperking. De stijfheid die nodig is om een ster van $2.6 M_\odot$ te ondersteunen in het DD2-model, resulteert in getijdenvervormingen die te groot zijn om te voldoen aan de GW170817-grenzen.
• De op SFHo gebaseerde EOS (zachter) gaat iets beter om met deze combinatie, maar staat nog steeds voor uitdagingen.

5. Betekenis

Dit werk biedt een rigoureuze statistische beoordeling van de maximale massa van neutronensterren, die verder gaat dan enkel-punts schattingen naar waarschijnlijkheidsverdelingen.

1. Validatie van Huidige Modellen: De resultaten ondersteunen het bestaan van neutronensterren tot $\sim 2.2-2.3 M_\odot$, en valideren huidige kernfysica-modellen tegen astrofysische data.
2. Beperking op Fase-overgangen: De gladde matchprocedure en de resulterende verdelingen suggereren dat grote, fase-overgangen van de eerste orde onwaarschijnlijk zijn in het dichtheidsbereik dat wordt onderzocht door huidige waarnemingen, of ten minste sterk beperkt zijn.
3. Toekomstige Richtingen: De studie benadrukt dat toekomstige precisie in straalmetingen (bijvoorbeeld van NICER of toekomstige röntgenmissies) en de definitieve aard van het GW190814-object (neutronenster versus zwart gat) cruciaal zijn voor het oplossen van de resterende onzekerheden in de EOS bij hoge dichtheid. De spanning tussen de massagap-hypothese en beperkingen op getijdenvervorming suggereert dat, als GW190814 inderdaad een neutronenster is, de EOS zachter moet zijn dan het DD2-model voorspelt, of dat de getijdenbeperkingen opnieuw moeten worden geëvalueerd.