Strange quark star I: the maximum gravitational mass and deformation of magnetized spinning model

Technische Zusammenfassung: Seltsamer Quark-Stern I: Die maximale Gravitationsmasse und Deformation eines magnetisierten, rotierenden Modells

Problemstellung
Die Arbeit befasst sich mit den strukturellen Parametern von seltsamen Quark-Sternen (Strange Quark Stars, SQS) unter den extremen Bedingungen, die in kompakten astrophysikalischen Objekten herrschen: starken Magnetfeldern (bis zu $5 \times 10^{17}$ G) und schneller Rotation (Frequenzen bis zu 1200 Hz). Während die Existenz von SQS – hypothetische kompakte Objekte, die vollständig aus seltsamer Quark-Materie (SQM) bestehen – bisher unbestätigt bleibt, stellen rezente Beobachtungen massiver kompakter Objekte (z. B. PSR J0952-0607 und das GW190814-Ereignis) Standardmodelle für Neutronensterne in Frage. Die Autoren zielen darauf ab, zu bestimmen, ob ein SQS-Modell, das dichteabhängige Zustandsgleichungen (EOS) einbezieht, starke Magnetfelder und Rotation berücksichtigt, Gravitationsmassen im Einklang mit diesen Beobachtungen unterstützen kann, während gleichzeitig die resultierende stellare Deformation und Energetik charakterisiert wird.

Methodik
Die Studie verwendet einen facettenreichen theoretischen und numerischen Ansatz:

1. Zustandsgleichung (EOS): Die Autoren nutzen das MIT-Bag-Modell mit einer dichteabhängigen Bag-Konstante ($B_{bag}(\rho)$), die durch eine Gauß-Beziehung definiert ist. Um die starken internen Magnetfelder zu berücksichtigen, werden die Fermi-Relationen für Quarks (Up, Down und Strange) modifiziert, um Landau-Quantisierungseffekte einzubeziehen. Das Modell geht von einer Null-Temperatur- und Null-Elektronenanteil-Approximation aus und vernachlässigt Elektronenbeiträge, da diese in der Größenordnung von $10^{-3}$ liegen.
2. Stellare Strukturgleichungen: Die Einsteinschen Feldgleichungen werden in einer axisymmetrischen Raumzeit mittels des 3+1-Formalismus gelöst. Der Energie-Impuls-Tensor enthält Beiträge des perfekten Fluids, der Magnetisierung und des reinen Magnetfeldes. Dies führt zu einem System von vier elliptischen partiellen Differentialgleichungen (PDEs), welche die Metrikfunktionen steuern.
3. Numerische Implementierung: Die Autoren verwenden die LORENE-Bibliothek, einen auf Spektralmethoden basierenden C++ Code, um die PDEs zu lösen. Das Rechengebiet ist in drei Regionen unterteilt: das Innere des Sterns, eine dünne Schale nahe der Oberfläche und ein Äußeres, das bis ins Unendliche reicht.
4. Parameterraum: Die Berechnungen wurden für Rotationsfrequenzen von 0, 400, 800 und 1200 Hz sowie zentralen Magnetfeldstärken von 0 bis $5 \times 10^{17}$ G durchgeführt. Für jede Konfiguration wurden 51 Sternmodelle berechnet, indem die zentrale Entalpie von $0,01c^2$ bis $0,51c^2$ variiert wurde.
5. Validierung: Die Genauigkeit der numerischen Lösungen wurde mithilfe der dreidimensionalen und zweidimensionalen allgemeinen relativistischen Virial-Identitäten (GRV3 und GRV2) verifiziert, wobei Fehler von $\approx 10^{-5}$ für nicht-magnetisierte/nicht-rotierende Modelle und $\approx 10^{-2}$ für magnetisierte/schnell rotierende Modelle erzielt wurden.

Wesentliche Beiträge und Ergebnisse

• Maximale Gravitationsmasse ($M_g^{max}$):

  • Das Modell sagt eine maximale Gravitationsmasse von 2,35 $M_\odot$ für einen nicht-magnetisierten, nicht-rotierenden SQS voraus.
  • Eine Erhöhung der Magnetfeldstärke (bis zu $5 \times 10^{17}$ G) in einem nicht-rotierenden Modell steigert $M_g^{max}$ leicht auf 2,37 $M_\odot$.
  • Rotation verstärkt die maximale Masse signifikant. Im am schnellsten rotierenden Modell (1200 Hz) mit dem stärksten Magnetfeld ($5 \times 10^{17}$ G) erreicht die maximale Gravitationsmasse 2,80 $M_\odot$.
  • Die Autoren stellen eine angepasste Funktion bereit, die $M_g^{max}$ mit dem magnetischen Moment ($\mu$) und der Rotationsfrequenz ($f$) in Beziehung setzt: $M_g^{max}(\mu, f)/M_g^{max}(0,0) \simeq (1 + a\mu^2)(1 + bf^c)$.

• Deformationsparameter ($a$):

  • Der Deformationsparameter ist als das Verhältnis des Äquatorradius zum Polradius definiert ($a = R_{eq}/R_{pol}$).
  • Die Studie stellt fest, dass die Rotation der dominierende Faktor für die Deformation ist, während das Magnetfeld einen geringfügigen Effekt auf die Abplattung hat.
  • Der maximale Deformationsparameter von 1,55 tritt im Modell mit $f = 1200$ Hz und $B_c = 5 \times 10^{17}$ G auf.
  • Eine angepasste Funktion für die Deformation wird ebenfalls bereitgestellt: $a(\mu, f)/a(0,0) \simeq (1 + \tilde{a}\mu^2)(1 + \tilde{b}f^{\tilde{c}})$.

• Energetik und Stabilität:

  • Bindungsenergie: Die gesamte Bindungsenergie pro Baryon ($|E_{BE}|/A$) liegt zwischen 171 MeV und 184 MeV, was die Stabilität der SQS-Modelle bestätigt. Das Verhältnis der Bindungsenergie zur Gravitationsmasse beträgt etwa 0,24 und zeigt eine lineare Beziehung zur Kompaktheit ($\beta = M_g/R_{circ}$).
  • Gesamtenergie: Die Gesamtenergie des SQS wird auf $\approx 5,5 \times 10^{47}$ J ($\approx 10^{54}$ erg) geschätzt. Der Beitrag der externen magnetischen Energie zur Gesamtenergie beträgt weniger als 1 %.
  • Stabilitätskriterium: Die Stabilität gegen axialsymmetrische Störungen wurde dadurch sichergestellt, dass $(dM/dR)_J > 0$ gilt.

• Vergleich mit Beobachtungen:

  • Die vom Modell vorhergesagten Massen sind mit beobachteten kompakten Objekten kompatibel. Insbesondere die nicht-rotierende, nicht-magnetisierte Masse (2,35 $M_\odot$) stimmt mit PSR J0952-0607 überein.
  • Die höheren Massengrenzen (bis zu 2,80 $M_\odot$) in rotierenden/magnetisierten Modellen sind konsistent mit dem Massenbereich der Sekundärkomponente des GW190814-Ereignisses (2,5–2,67 $M_\odot$).
  • Die Modelle decken zudem die Masse- und Radienschätzungen für den Magnetar SGR J1745–2900 und andere massive Pulsare wie PSR J1614-2230 und PSR J0348+0432 ab.

Bedeutung und Ansprüche
Die Arbeit behauptet, dass das vorgeschlagene SQS-Modell, welches ein dichteabhängiges MIT-Bag-Modell mit Landau-Quantisierung nutzt, die hohen beobachteten Gravitationsmassen ohne die Notwendigkeit exotischer Physik jenseits des Standardrahmens für Quark-Materie reproduzieren kann. Die Autoren betonen, dass Magnetfelder unter $10^{19}$ G die Steifigkeit der EOS (wie in früheren Studien von Chatterjee et al. angemerkt) nicht signifikant verändern, aber in Kombination mit der Rotation die strukturellen Parameter (Masse und Form) des Sterns erheblich beeinflussen.

Die Studie etabliert, dass schnell rotierende, stark magnetisierte SQS Massen von bis zu 2,8 $M_\odot$ tragen können, was eine potenzielle Erklärung für die durch LIGO/Virgo detektierten „Mass Gap“-Objekte und schweren Pulsare bietet. Die Autoren geben explizit an, dass die minimalen Gravitationsmasse-Grenzen und Kepler-Konfigurationen für diese Modelle in einer Begleitpublikation untersucht werden, und räumen ein, dass die genaue Berechnung kritischer Konfigurationen für stark magnetisierte, schnell rotierende Sterne eine Limitation für zukünftige Verfeinerungen bleibt.