Maximal mass of neutron stars constrained by neutron star observations

Durch die Anwendung eines bayesschen Gewichtsrahmens, der Multimessenger-Beobachtungen (einschließlich GW170817, NICER und Kandidaten für kompakte Objekte) mit hybriden Zustandsgleichungen integriert, bestimmt diese Studie, dass die maximale Neutronensternmasse robust auf etwa 2,2–2,3 Sonnenmassen eingeschränkt ist, während der entsprechende Radius stärker vom zugrunde liegenden hadronischen Modell abhängt und typischerweise in der Nähe von 12 km liegt.

Das Wesentliche

Stellen Sie sich vor, das Universum ist mit einem mysteriösen, superdichten Material gefüllt, das „Neutronensternmaterie" genannt wird. Es ist so schwer, dass ein einziger Teelöffel davon so viel wiegen würde wie ein Berg. Seit langem versuchen Physiker herauszufinden, welche „Spielregeln" für dieses Material gelten – genauer gesagt, wie steif oder quetschbar es ist. Diese Regelmenge wird als Zustandsgleichung (EOS) bezeichnet.

Die große Frage, die dieser Artikel zu beantworten versucht, lautet: Wie schwer kann ein Neutronenstern maximal werden, bevor er zu einem Schwarzen Loch kollabiert?

Hier ist die Geschichte, wie die Autoren dieses Rätsel gelöst haben, einfach erklärt:

1. Die beiden Ausgangspunkte (Die Rezepte)

Um die Regeln herauszufinden, starteten die Wissenschaftler mit zwei verschiedenen „Rezepten" dafür, wie sich diese dichte Materie bei niedrigeren Dichten verhält. Stellen Sie sich diese als zwei verschiedene Theorien darüber vor, wie die Zutaten gemischt werden:

• Rezept A (SFHo): Ein „weicheres" Rezept, was bedeutet, dass die Materie etwas leichter zu quetschen ist.
• Rezept B (DD2): Ein „steiferes" Rezept, was bedeutet, dass die Materie sich stärker gegen das Quetschen wehrt.

Sie wussten, dass diese Rezepte im „Anfangsbereich" der Dichteskala gut funktionierten, aber sie wussten nicht, was bei den extremen, superhohen Dichten im Zentrum eines Neutronensterns passiert. Um diese Lücke zu füllen, nutzten sie eine mathematische „Brücke", um ihre Rezepte mit dem zu verbinden, was wir über die Teilchenphysik bei den höchstmöglichen Energien wissen.

2. Die Detektivarbeit (Echte Hinweise nutzen)

Anstatt nur zu raten, gingen die Autoren wie Detektive vor. Sie nahmen ihre beiden Rezepte und testeten sie gegen echte Hinweise, die von Teleskopen und Gravitationswellendetektoren gesammelt wurden. Sie verwendeten eine spezielle statistische Methode (sogenanntes Bayes'sches Weighting), um zu sehen, welche Versionen ihrer Rezepte den Test bestanden.

Hier sind die Hinweise, die sie verwendeten:

• Der „Große Crash" (GW170817): Als zwei Neutronensterne ineinander krachten, schickten sie Wellen durch den Raum. Das Verhalten dieser Wellen verrät den Wissenschaftlern, wie „quetschbar" die Sterne waren.
• Die „Taschenlampe" (NICER): Ein Weltraumteleskop machte Aufnahmen von heißen Flecken auf rotierenden Neutronensternen. Indem sie maß, wie groß die Sterne aussahen und wie schwer sie waren, erhielten sie ein direktes Verhältnis von Größe zu Gewicht.
• Der „Leichtgewicht"-Kandidat (HESS J1731–347): Ein sehr kleines, leichtes Objekt, das möglicherweise ein Neutronenstern ist.
• Der „Schwergewicht"-Kandidat (GW190814): Ein mysteriöses Objekt, das schwerer ist als die meisten Neutronensterne, aber leichter als die meisten Schwarzen Löcher. Die Wissenschaftler fragten: Könnte dies tatsächlich ein super schwerer Neutronenstern sein?

3. Die Ergebnisse: Was die Hinweise ihnen sagten

Die Wissenschaftler führten ihre beiden Rezepte durch diese Hinweise und betrachteten die Ergebnisse.

Die Gewichtsbeschränkung (Maximale Masse):

• Die Überraschung: Es spielte kaum eine Rolle, welches Startrezept (Weich oder Steif) sie verwendeten. Die realen Hinweise waren so stark, dass sie beide Rezepte zwangen, sich auf dieselbe Antwort zu einigen.
• Das Urteil: Als sie die zuverlässigsten Hinweise verwendeten (den „Großen Crash" und die „Taschenlampe"), beträgt das maximale Gewicht, das ein Neutronenstern tragen kann, etwa 2,2 bis 2,3 Mal die Masse unserer Sonne.
• Die „Schwergewicht"-Wendung: Wenn sie annehmen, dass dieses mysteriöse schwere Objekt (GW190814) ein Neutronenstern ist, springt die Grenze auf etwa 2,6 bis 2,7 Sonnenmassen an. Dies erzeugt jedoch einen Konflikt mit den „Quetschbarkeits"-Hinweisen des Großen Crashes, was die Situation kompliziert macht.

Die Größenbeschränkung (Radius):

• Der Unterschied: Im Gegensatz zum Gewicht hing die Größe des Sterns davon ab, welches Startrezept sie verwendeten.
• Das Urteil: Das „Weiche" Rezept sagte einen Radius von etwa 11,8 km voraus, während das „Steife" Rezept etwa 12,4 km vorhersagte.
• Der Sweet Spot: Wenn alle besten Hinweise kombiniert werden, beträgt die wahrscheinlichste Größe für diese Sterne rund 12 Kilometer (plus oder minus 1 km).

4. Das große Ganze

Der Artikel kommt zu dem Schluss, dass wir durch die Betrachtung der „Endpunkte" (der schwersten und größten möglichen Sterne) und die Verwendung einer Mischung aus echten astronomischen Daten die Regeln für die dichteste Materie des Universums eingrenzen können.

• Das Gewicht: Das Universum scheint eine „Geschwindigkeitsbegrenzung" dafür zu haben, wie schwer ein Neutronenstern sein kann, die bequem bei 2,2 bis 2,3 Sonnenmassen liegt. Dies entspricht dem schwersten Neutronenstern, den wir bisher tatsächlich gesehen haben.
• Die Größe: Sie sind ungefähr so groß wie eine kleine Stadt, etwa 12 km im Durchmesser.
• Die Erkenntnis: Die realen Beobachtungen (die Hinweise) sind viel mächtiger als die theoretischen Startannahmen. Egal, mit welcher Theorie man beginnt, die Daten der Sterne selbst zwingen die Antwort dazu, auf dieselben Zahlen zu konvergieren.

Kurz gesagt, das Universum hat uns eine sehr klare Antwort gegeben: Neutronensterne können unglaublich schwer werden, aber es gibt eine harte Obergrenze, und sie sind für ihr Gewicht überraschend klein.

Technische Zusammenfassung

1. Problemstellung

Die zentrale Herausforderung in der Kern- und Schwerionenphysik besteht darin, die Zustandsgleichung (EOS) stark wechselwirkender Materie bei mittleren Dichten und chemischen Potentialen zu bestimmen. Während Materie bei niedrigen Dichten durch Kernexperimente und Gitter-QCD eingeschränkt ist und Materie bei hohen Dichten durch störungstheoretische QCD (pQCD), bleibt der für Neutronensternkerne relevante Bereich (einige Vielfache der nuklearen Sättigungsdichte, $\rho_0$) schlecht verstanden.

Das spezifische behandelte Problem ist die Bestimmung der maximalen Masse ($M_{TOV}$) und des entsprechenden Radius ($R_{TOV}$) von Neutronensternen. Diese Endpunkte der Masse-Radius-Sequenzen ($M(R)$) sind entscheidend für das Verständnis der Steifigkeit der EOS bei supra-nuklearen Dichten. Die Autoren zielen darauf ab, zu quantifizieren, wie aktuelle Multimessenger-Beobachtungen diese Endpunkte einschränken und in welchem Maße die Ergebnisse von der Wahl des zugrunde liegenden hadronischen Basismodells abhängen.

2. Methodik

A. Konstruktion hybrider EOS

Die Autoren konstruieren Familien kausaler hybrider EOS, indem sie Materie bei niedrigen Dichten (hadronisch) mit Materie bei hohen Dichten (Quarkmaterie) verbinden:

• Bereich niedriger Dichte: Zwei repräsentative hadronische Modelle werden als Basen verwendet:
  • SFHo: Ein relativistisches Mean-Field-Modell (weichere EOS, Inkompressibilität $K=245$ MeV).
  • DD2: Eine steifere hadronische EOS.
• Bereich hoher Dichte: Ein erweitertes Lineares Sigma-Modell (eLSM) basierend auf $U(3) \times U(3)$ chiraler Symmetrie, das Konstituentenquarks und Meson-Nonette integriert. Dieses Modell wird so eingeschränkt, dass es bei niedrigen chemischen Potentialen mit der Thermodynamik der Gitter-QCD und bei hohen chemischen Potentialen mit pQCD-Ergebnissen übereinstimmt.
• Interpolation: Um die beiden Bereiche zu überbrücken, wird eine Polynominterpolation fünften Grades auf die Energiedichte als Funktion der Baryondichte angewendet. Dies gewährleistet thermodynamische Konsistenz und Kontinuität von Druck und Schallgeschwindigkeit ($c_s$) über den Übergangsbereich hinweg und vermeidet große Phasenübergänge erster Ordnung.
• Asymptotische Einschränkung: Die Quark-Meson-EOS muss sich bei asymptotisch hohen Dichten ($\mu \approx 2.6$ GeV) glatt an pQCD-Ergebnisse anschließen.

B. Bayesscher Rahmen

Die Studie verwendet ein Bayessches Gewichtungssystem, um die Wahrscheinlichkeitsverteilungen der Endpunkte der $M(R)$-Sequenz zu analysieren.

1. Prior: Ein breites Ensemble kausaler hybrider EOS wird generiert, das einen Bereich der Steifigkeit bei supra-nuklearen Dichten abdeckt.
2. Likelihood-Gewichtung: Das Ensemble wird gegen Beobachtungsdaten gewichtet, um posteriori-Wahrscheinlichkeitsverteilungen für $M_{TOV}$ und $R_{TOV}$ zu erzeugen.

C. Beobachtungsbeschränkungen

Die folgenden Beschränkungen werden sequenziell und in Kombination angewendet:

1. pQCD-Beschränkung: Stellt sicher, dass sich die EOS dem N3LO pQCD-Referenzpunkt ($\mu=2.6$ GeV, $\rho=6.47$ fm$^{-3}$, $p=3823$ MeV/fm$^3$) nähert, während die Kausalität gewahrt bleibt ($c_s < 1$).
2. GW170817: Gezeitenverformbarkeitsbeschränkung ($\tilde{\Lambda} < 720$) und eine strengere kombinierte Beschränkung („$\Lambda_{190}$": $70 < \Lambda_{1.4} < 580$ und $9.1 < R_{1.4} < 12.8$ km).
3. NICER: Masse-Radius-Messungen für fünf Pulsare (J0030+0451, J0740+6620, J0614–3329, J1231–1411, J0437–4715) mit einer Präzision von $\sim 10\%$.
4. HESS J1731–347: Ein Kandidat für einen sehr massearmen Neutronenstern ($M \approx 0.77 M_\odot$), der aufgrund von Interpretationsunsicherheiten mit Vorsicht behandelt wird.
5. GW190814: Das „Mass Gap"-Objekt ($M \approx 2.59 M_\odot$). Die Autoren testen die Hypothese, dass dieses Objekt ein Neutronenstern ist, indem sie eine Gaußsche Beschränkung anwenden, die bei $2.59 M_\odot$ zentriert ist.

3. Hauptbeiträge

• Vereinheitlichte Diagnose: Zeigt, dass die Wahrscheinlichkeitsverteilungen der Endpunkte der $M(R)$-Sequenz als sensitive, komplementäre Diagnose dienen, um das Verhalten der EOS bei hohen Dichten innerhalb eines Multimessenger-Rahmens einzuschränken.
• Unabhängigkeit vs. Abhängigkeit von der Basis: Trennt systematisch den Einfluss von Beobachtungsdaten von der theoretischen Wahl der hadronischen Basis (SFHo vs. DD2).
• Identifizierung von Spannungen: Identifiziert eine spezifische Spannung zwischen der steifen DD2-EOS, der Hypothese des „Mass Gap"-Neutronensterns und den strengen Gezeitenverformbarkeitsbeschränkungen ($\Lambda_{190}$).

4. Hauptergebnisse

Verteilungen der maximalen Masse ($M_{TOV}$)

• Dominanz der Beobachtungen: Der Peak der $M_{TOV}$-Verteilung wird primär durch Beobachtungsbeschränkungen getrieben und nicht durch die Wahl der hadronischen Basis-EOS (SFHo oder DD2).
• Fortschreiten der Beschränkungen:
  • Nur pQCD: Liefert einen Peak knapp über $1.5 M_\odot$.
  • GW170817/HESS: Keine signifikante Verschiebung der Peak-Position.
  • NICER: Verschiebt den Peak auf 2.1–2.2 $M_\odot$ und schärft die Verteilung.
  • Mass Gap (GW190814): Verschiebt den Peak signifikant auf 2.6–2.7 $M_\odot$.
• Robuste Schlussfolgerung: Bei Anwendung der robustesten Beschränkungen (pQCD + GW170817 + NICER) konvergieren beide EOS-Modelle auf eine bevorzugte maximale Masse von 2.2–2.3 $M_\odot$. Dies ist konsistent mit der beobachteten Masse von PSR J0952–0607 ($2.35 \pm 0.17 M_\odot$).

Verteilungen des Radius ($R_{TOV}$)

• Abhängigkeit von der Basis: Im Gegensatz zur Masse zeigt die Radiusverteilung eine starke Abhängigkeit von der zugrunde liegenden hadronischen EOS.
  • SFHo+eLSM: Peak nahe 11.8 km (mit robusten Beschränkungen).
  • DD2+eLSM: Peak nahe 12.4 km (mit robusten Beschränkungen).
• Allgemeine Präferenz: Kombinierte Beschränkungen bevorzugen einen Radiusbereich von $12 \pm 1$ km.
• Auswirkung der Gezeitenverformbarkeit: Die Anwendung der strengeren $\Lambda_{190}$-Beschränkung verengt die Radiusverteilungen und senkt die bevorzugten Werte, da eine geringere Gezeitenverformbarkeit mit kleineren Radien korreliert.

Die „Mass Gap"-Spannung

• Wenn das GW190814-Objekt als Neutronenstern interpretiert wird, versagt die auf DD2 basierende EOS, eine konsistente Beschreibung zu liefern, wenn sie mit der $\Lambda_{190}$-Gezeitenbeschränkung kombiniert wird. Die Steifigkeit, die erforderlich ist, um einen $2.6 M_\odot$-Stern im DD2-Modell zu tragen, führt zu Gezeitenverformbarkeiten, die zu groß sind, um die GW170817-Grenzen zu erfüllen.
• Die auf SFHo basierende EOS (weicher) kommt mit dieser Kombination etwas besser zurecht, steht jedoch weiterhin vor Herausforderungen.

5. Bedeutung

Diese Arbeit liefert eine rigorose statistische Bewertung der maximalen Masse von Neutronensternen und geht über Einzelwert-Schätzungen hinaus hin zu Wahrscheinlichkeitsverteilungen.

1. Validierung aktueller Modelle: Die Ergebnisse unterstützen die Existenz von Neutronensternen bis zu $\sim 2.2-2.3 M_\odot$ und validieren aktuelle kernphysikalische Modelle gegenüber astrophysikalischen Daten.
2. Einschränkung von Phasenübergängen: Das glatte Anpassungsverfahren und die resultierenden Verteilungen deuten darauf hin, dass große Phasenübergänge erster Ordnung im von aktuellen Beobachtungen abgetasteten Dichtebereich unwahrscheinlich sind oder zumindest stark eingeschränkt werden.
3. Zukünftige Richtungen: Die Studie hebt hervor, dass zukünftige Präzision bei Radiusmessungen (z. B. von NICER oder zukünftigen Röntgenmissionen) und die definitive Natur des GW190814-Objekts (Neutronenstern vs. Schwarzes Loch) entscheidend sind, um die verbleibenden Unsicherheiten in der EOS bei hohen Dichten zu klären. Die Spannung zwischen der Mass-Gap-Hypothese und den Gezeitenverformbarkeitsbeschränkungen legt nahe, dass, falls GW190814 tatsächlich ein Neutronenstern ist, die EOS weicher sein muss als vom DD2-Modell vorhergesagt, oder dass die Gezeitenbeschränkungen einer Neubewertung bedürfen.