Evolution of proto-neutron stars to pulsars, magnetars and central compact objects

Technische Zusammenfassung: Evolution von Proto-Neutronensternen zu Pulsaren, Magnetaren und Central Compact Objects

Problemstellung
Der Ursprung der ultra-starken Magnetfelder, die bei Magnetaren beobachtet werden ($B \sim 10^{14}–10^{15}$ G), bleibt Gegenstand der Debatte. Während die „Fossil Field“-Hypothese (Flusserhaltung aus den Vorläufersternen) die Felder konventioneller Neutronensterne und Central Compact Objects (CCOs) erklären kann, hat sie Schwierigkeiten, die Magnetarfelder zu erklären, ohne extrem seltene Progenitor-Bedingungen vorauszusetzen. Umgekehrt gilt die Dynamo-Hypothese, die eine schnelle Rotation und kräftige Konvektion während der Phase des Proto-Neutronensterns (PNS) erfordert, als weitgehend akzeptiert für Magnetare, steht jedoch vor einer Herausforderung: Die abgeleitete Geburtsrate von Magnetaren ist vergleichbar mit der Rate der Kernkollaps-Supernovae, was darauf hindeutet, dass die Entstehung von Magnetaren nicht so selten ist, wie es extreme Parameteranforderungen implizieren würden. Zudem legen Standardberechnungen zur Flusserhaltung nahe, dass selbst bei einem kollabierenden Eisenkern die resultierenden Felder ($B \sim 10^{10}$ G) nicht ausreichen, um die typischen Pulsarfelder ($B \sim 10^{12}$ G) zu erklären, was impliziert, dass selbst für gewöhnliche Pulsare ein Dynamo-Mechanismus erforderlich ist. Das zentrale Problem besteht darin, die hohen Geburtsraten verschiedener Neutronenstern-Typen mit ihren vielfältigen Magnetfeldstärken unter einem einheitlichen Evolutionsrahmen in Einklang zu bringen.

Methodik
Die Autoren verwenden ein reduziertes, eindimensionales schergetriebenes $\alpha-\Omega$-Dynamomodell, das von Wickramasinghe et al. (2014) adaptiert wurde, welche es ursprünglich auf Verschmelzungen von Weißen Zwergen angewendet haben. Das Modell verfolt die Entwicklung der poloidalen ($B_p$) und toroidalen ($B_\phi$) Magnetfelder während der PNS-Phase. Zu den wesentlichen Merkmalen der Methodik gehören:

• Gleichungsordnung: Das Modell löst gekoppelte Differentialgleichungen für $B_p$ und $B_\phi$. Das toroidale Feld wird durch den $\Omega$-Effekt (Differenzialrotation, die auf das poloidale Feld wirkt) erzeugt, während das poloidale Feld über den $\alpha$-Effekt (turbulente Konvektion, die auf das toroidale Feld wirkt) regeneriert wird.
• PNS-Evolution: Im Gegensatz zu statischen Modellen beinhaltet dieses Framework die dynamische Kontraktion des PNS von einem Anfangsradius $R_0 \sim 40$ km zu einem Endradius $R_\infty \approx 12$ km. Die Radius-Evolution wird basierend auf Camelio et al. (2017) für verschiedene Massen und Zustandsgleichungen (EoS) parametrisiert.
• Masse und Drehimpuls: Das Modell berücksichtigt die Abnahme der Gravitationsmasse durch Neutrino-Kühlung sowie die Erhaltung des Drehimpulses, was den Spin-up des PNS antreibt.
• Viskosität und Scherungsdämpfung: Die Dämpfung der Differenzialrotation (Scherrate $q$) wird über eine Kombination aus magnetorotationsinstabilitäts-getriebener (MRI) Viskosität und konvektiver Viskosität modelliert. Die Scherrate nimmt über die Zeit ab, was den Dynamo schließlich beendet.
• Flussverlust und Stabilität: Das Modell enthält einen Rekonnektions-Term, um den magnetischen Flussverlust zu berücksichtigen (initial auf 20 % gesetzt), und integriert Stabilitätskriterien (Braithwaite 2009), wobei Feldinstabilitäten die Zerfallszeiten von $B_p$ und $B_\phi$ bestimmen.
• Parameterraum-Scan: Die Autoren integrieren die Gleichungen numerisch mittels der Runge-Kutta-Methode und variieren systematisch die Anfangsparameter, einschließlich der Anfangsscherrate ($q_0$), der Anfangsmasse ($M_0$), des Anfangsradius ($R_0$), der Anfangs-Seed-Felder und der finalen Rotationsperiode ($P_\infty$).

Zentrale Ergebnisse
Die Simulationen zeigen, dass die finale Magnetfeldkonfiguration hochsensibel auf die Anfangsrotationsrate und die Scherung reagiert, was es einem einzigen Dynamo-Mechanismus ermöglicht, das volle Spektrum beobachteter Neutronenstern-Typen zu erzeugen:

1. Scherraten-Schwellenwert: Der Sättigungswert des toroidalen Feldes ($B_{\phi, \infty}$) hängt stark von der Anfangsscherrate $q_0$ ab. Es existiert ein kritischer Schwellenwert bei $q_0 \simeq 0.23$.

   • Unterhalb des Schwellenwerts ($q_0 < 0.23$): Der $\alpha$-Effekt wird nicht ausreichend ausgelöst. Das poloidale Feld bleibt nahe dem Wert, der durch die Flusserhaltung diktiert wird ($B_p \approx 2,5 \times 10^{10}$ G), und das toroidale Feld wächst nicht signifikant an. Dieses Regime entspricht den Central Compact Objects (CCOs).
   • Oberhalb des Schwellenwerts: Der Dynamo wird aktiv. Höhere Scherraten führen zu stärkeren toroidalen Feldern.

2. Bildung verschiedener Populationen:

   • Magnetare: Für die kürzesten Anfangsperioden (schnelle Rotation, z. B. $P_\infty \sim 2,5$ ms) produziert das Modell Magnetar-ähnliche Stärken mit $B_p \simeq 10^{15}$ G und $B_\phi \simeq 10^{16}$ G.
   • Gewöhnliche Pulsare: Für langsamere Rotatoren liefert das Modell typische Pulsarfelder ($B_p \simeq 10^{12}$ G, $B_\phi \simeq 10^{14}$ G).
   • Low-B Magnetare / High-B Pulsare: Mittlere Rotationsraten erzeugen Felder im Bereich von Low-B Magnetaren oder High-B Pulsaren.

3. Rolle der Flusserhaltung: Die Anfangsstärke des poloidalen Feldes wird weitgehend durch die Flusserhaltung während des Kollapses bestimmt. Der Dynamo-Prozess (speziell der $\alpha$-Effekt) ist jedoch erforderlich, um diese Felder auf das Niveau von Pulsaren und Magnetaren zu verstärken. Das toroidale Feld wird generisch vom $\Omega$-Effekt dominiert.

4. Central Compact Objects (CCOs): Das Modell legt nahe, dass CCOs keinen kräftigen Dynamo-Prozess durchlaufen, da sie eine langsame Anfangsrotation ($P_0 \sim 100$ ms) aufweisen. Ihre poloidalen Felder werden durch die Flusserhaltung bestimmt, aber sie könnten dennoch toroidale Felder erwerben, die allein durch den $\Omega$-Effekt verstärkt werden, was ihre kleinen Dipolfelder, aber potenziell signifikanten internen toroidalen Komponenten erklärt.

Bedeutung und Behauptungen
Das Paper behauptet, die scheinbare Diskrepanz zwischen den hohen Geburtsraten von Magnetaren und den extremen Bedingungen, die oft für Formationsszenarien gefordert werden, aufzulösen. Indem sie zeigen, dass ein während der PNS-Phase operierender Dynamo-Prozess natürlich ein Kontinuum von Magnetfeldstärken je nach Anfangsrotationsrate erzeugen kann, schlagen die Autoren vor, dass:

• Sowohl gewöhnliche Pulsare als auch Magnetare wahrscheinlich eine Dynamo-Phase durchlaufen; der Unterschied liegt in der Effizienz dieses Prozesses, der durch die Anfangsrotation angetrieben wird.
• Die „Fossil Field“-Hypothese unzureichend ist, um typische Pulsarfelder zu erklären, was einen Dynamo selbst für Nicht-Magnetar-Neutronensterne notwendig macht.
• CCOs das Niedrig-Rotations-Ende dieses Spektrums darstellen, in dem der Dynamo unterdrückt wird, wodurch die Felder von der Flusserhaltung dominiert werden.

Die Autoren betonen, dass ihr eindimensionales Modell, obwohl vereinfacht, die wesentlichen Abhängigkeiten von Rotation, Scherung und Auftrieb erfasst und qualitativ mit komplexeren dreidimensionalen Simulationen (z. B. Raynaud et al. 2020) übereinstimmt. Sie kommen zu dem Schluss, dass der $\alpha-\Omega$-Dynamo einen einheitlichen Rahmen bietet, um die Vielfalt der jungen Neutronenstern-Magnetfelder ohne die Anforderung seltener, extremer Progenitor-Bedingungen für Magnetare zu verstehen.