Rotating Fermion-Boson Stars in $R$-squared Gravity

Technische Zusammenfassung: Rotierende Fermionen-Bosonen-Sterne in der $R^2$-Gravitation

Problemstellung
Kompakte Objekte, insbesondere Neutronensterne (NSs), dienen als entscheidende Laboratorien zur Untersuchung der Gravitation im starken Feldregime. Während die Allgemeine Relativitätstheorie (GR) der Standard bleibt, führen Erweiterungen wie die $f(R)$-Gravitation zusätzliche Skalarfreiheitsgrade (Skalaronen) ein, die die Sternstruktur verändern können. Gleichzeitig bietet die potenzielle Akkumulation von Dunkler Materie innerhalb von NSs – modelliert als bosonische Komponente, die nur gravitativ mit fermionischer Materie interagiert – einen Mechanismus, um ungewöhnlich massive kompakte Objekte oder solche zu erklären, die in die „Low-Mass-Gap“ zwischen Neutronensternen und Schwarzen Löchern fallen (z. B. die Sekundärkomponente von GW190814). Vorherige Studien haben entweder statische Mischsterne in modifizierter Gravitation oder rotierende Sterne in der GR getrennt behandelt. Dieser Arbeit fehlte eine konsistente Behandlung, die modifizierte Gravitation (speziell $R^2$-Gravitation), gleichmäßige Rotation und einen gemischten Fermionen-Bosonen-Sektor kombiniert. Diese Arbeit adressiert diese Lücke durch die Konstruktion statischer und gleichmäßig rotierender Gleichgewichtmodelle von Fermionen-Bosonen-Sternen im Rahmen des $f(R) = R + aR^2$-Modells, um zu untersuchen, wie der Skalarfreiheitsgrad und die Rotation gemeinsam die Masse-Radius-Beziehungen und die interne Struktur beeinflussen.

Methodik
Die Autoren konstruieren Gleichgewichtslösungen unter Verwendung eines numerischen Rahmens, der den RNS-Code (basierend auf der Komatsu-Eriguchi-Hachisu Self-Consistent-Field-Methode) erweitert, um $f(R)$-Gravitation und einen gemischten Materiesektor einzubeziehen.

1. ** Theoretischer Rahmen:**

   • Gravitation: Die Studie verwendet das $R^2$-Modell, $f(R) = R + aR^2$, wobei $a$ ein positiver Parameter ist, der die Abweichung von der GR kontrolliert. Diese Theorie ist dynamisch äquivalent zu einer Skalar-Tensor-Theorie mit einem massiven Skalarfeld. Die Gleichungen werden im Einstein-Frame über eine konforme Transformation formuliert, was ein Skalarfeld $\phi$ mit einem spezifischen Potenzial $V(\phi)$ einführt.
   • Materiesektor: Das System besteht aus zwei Komponenten:
     • Fermionischer Sektor: Modelliert als perfekte Flüssigkeit mit einer Nulltemperatur-Zustandsgleichung (EOS) unter Verwendung der tabellierten AkmalPR-EOS.
     • Bosonischer Sektor: Repräsentiert durch ein selbstinteragierendes komplexes Skalarfeld mit einem Potenzial $V_s(|\Phi|^2) = \mu_b^2|\Phi|^2 + \frac{1}{2}\lambda|\Phi|^4$.
   • Approximation: Der bosonische Sektor wird unter Verwendung der RPV-Approximation (Rigidly Rotating Potential/High-Coupling) behandelt. Im Regime starker Selbstwechselwirkung (großes $\lambda$) wird das bosonische Feld als effektive perfekte Flüssigkeit mit kompakter Unterstützung approximiert. Der exponentiell abfallende „Schwanz“ des Feldes wird vernachlässigt, und das Feld wird außerhalb eines durch die Feldparameter bestimmten Radius auf Null gesetzt. Dies ermöglicht es, den bosonischen Beitrag als Quellterm ähnlich der fermionischen Flüssigkeit zu behandeln, was die numerische Implementierung vereinfacht.

2. Numerische Implementierung:

   • Die Feldgleichungen für die Metrikpotentiale und das Skalarfeld werden iterativ unter Verwendung von Green'schen Funktionen gelöst.
   • Der bosonische Feldansatz für rotierende Konfigurationen beinhaltet eine harmonische Zeit- und Azimutabhängigkeit: $\Phi(t, r, \theta, \psi) = \phi(r, \theta)e^{i(\eta\psi - wt)}$.
   • Die fermionische Flüssigkeit wird als starr rotierend angenommen, während der bosonische Sektor eine effektive differenzielle Rotation aufweist, die durch die Feldparameter bestimmt wird, obwohl die RPV-Approximation dieses Profil näher an eine starre Rotation rückt.
   • Globale Größen (Masse $M$, Drehimpuls $J$) werden über Komar-Integrale im Einstein-Frame berechnet, welche mit den Jordan-Frame-Größen im Unendlichen übereinstimmen, da das Skalarfeld dort verschwindet.

Zentrale Beiträge und Ergebnisse
Die Arbeit präsentiert Gleichgewichtsequenzen für statische und gleichmäßig rotierende Fermionen-Bosonen-Sterne für drei Gravitationsmodelle: GR ($a=0$) sowie zwei $R^2$-Fälle ($a=10$ und $a=10^4$).

1. Strukturelle Modifikationen:

   • Einfluss des Skalarfeldes: Die Anwesenheit des Skalarfreiheitsgrades modifiziert systematisch die räumliche Verteilung beider Materiesektoren. Im Vergleich zur GR führt die $R^2$-Gravitation zu einer Unterdrückung der Amplitude des bosonischen Feldes (insbesondere in den inneren Regionen) und einer Umverteilung der fermionischen Energiedichte hin zu räumlich ausgedehnteren Profilen.
   • Massenumverteilung: Für eine feste Gesamtmasse erfordern die modifizierten Gravitationsmodelle geringere Zentraldichten. Die Skalarwechselwirkung reduziert effektiv die lokalisierte Materiekonzentration, die zur Unterstützung einer gegebenen Masse benötigt wird, was zu größeren Sternradien und einer „weicheren“ effektiven Gravitationskopplung im Kern führt.
   • Rotationseffekte: Rotierende Konfigurationen weisen eine toroidale bosonische Struktur auf (verschwindet auf der Rotationsachse). Die dimensionslose interne Frequenz $\bar{w}$ steuert, ob die bosonische Komponente eine kernartige Struktur oder einen Halo bildet. Eine Erhöhung von $a$ vergrößert den Parameterraum, der kompakte, kernartige Mischkonfigurationen zugänglich macht, selbst bei niedrigeren internen Frequenzen, in denen solche Lösungen in der GR nicht existieren würden.

2. Masse-Radius-Beziehungen und maximale Massen:

   • Erhöhte maximale Masse: Der $R^2$-Term erhöht die maximal unterstützte Masse sowohl für statische als auch für Kepler-Limits (Massenverlustgrenze) relativ zur GR.
     • Für $a=10$ steigt die statische maximale Masse um $\sim 2,9\,\%$ und das Kepler-Limit um $\sim 6,2\,\%$.
     • Für $a=10^4$ steigt die statische maximale Masse um $\sim 7,5\,\%$ und das Kepler-Limit um $\sim 13,9\,\%$.
   • Kepler vs. Statisch: Der relative Zuwachs der maximalen Masse ist am Kepler-Limit konsistent größer als im statischen Limit.

3. Beobachtungsbedingte Einschränkungen:

   • Die Modelle werden gegen aktuelle astrophysikalische Einschränkungen getestet, einschließlich der NICER-Messungen von PSR J0030+0451 und PSR J0740+6620, der GW170817-Gezeitenverformbarkeit und der Masse der Sekundärkomponente von GW190814 ($\sim 2,5\text{--}2,67 M_\odot$).
   • Interpretation von GW190814: In der GR erfordert die Reproduktion der GW190814-Sekundärmasse mit der AkmalPR-EOS eine nahezu Kepler-Rotation, was den Radius über die durch PSR J0740+6620 gesetzten Grenzen hinaus aufbläht. Im Gegensatz dazu können die $R^2$-Modelle (insbesondere $a=10^4$) diese Masse selbst in nahezu statischen Sequenzen unter Einhaltung der Radienbeschränkungen aufnehmen. Dies deutet darauf hin, dass gemischte Fermionen-Bosonen-Sterne in der $R^2$-Gravitation eine praktikable Erklärung für Objekte in der Low-Mass-Gap zwischen Neutronensternen und Schwarzen Löchern bieten, ohne extreme nukleare Steifigkeit oder schnelle Rotation vorauszusetzen.

Bedeutung und Ansprüche
Das Papier behauptet, dass Fermionen-Bosonen-Sterne in der $R^2$-Gravitation einen vielversprechenden Rahmen zur Untersuchung der kombinierten Effekte von dunkler bosnischer Materie, Rotation und Modifikationen der starken Gravitation bilden. Die primäre Bedeutung liegt in der Demonstration, dass der Skalarfreiheitsgrad in der $f(R)$-Gravitation das Domänenintervall zulässiger Gleichgewichtslösungen vergrößert und die maximal unterstützten Massen im Vergleich zur GR erhöht.

Die Autoren schließen bescheiden, dass ihre Modelle zwar mit aktuellen Beobachtungseinschränkungen kompatibel bleiben, aber die spezifische Fähigkeit der $R^2$-Gravitation, die GW190814-Sekundärmasse in einem nicht-rotierenden oder langsam rotierenden Zustand zu unterstützen, sie von Standard-GR-Modellen mit konventioneller nuklearer Materie unterscheidet. Sie betonen, dass dieses Ergebnis spezifisch für die gewählte EOS und die Approximation (RPV) ist und dass der Rahmen eine testbare Alternative für kompakte Objekte im Low-Mass-Gap darstellt. Die Arbeit behauptet nicht, die Existenz solcher Sterne endgültig bewiesen zu haben, sondern etabliert sie als eine physikalisch konsistente und beobachtungstechnisch lebensfähige Kandidatenklasse innerhalb modifizierter Gravitationstheorien.