Internal Heating in Magnetars: Role of Electron Captures

Technische Zusammenfassung: Interne Erwärmung in Magnetaren: Die Rolle von Elektroneninfängen

Problemstellung
Magnetare, die durch extreme Magnetfelder ($B \gtrsim 10^{14}$ G) charakterisiert sind, weisen persistente Röntgenleuchtkräfte ($\sim 10^{33} - 10^{35}$ erg/s) und gelegentliche Riesenflares auf, die deren Rotationsenergie übersteigen. Während krustale Deformationen, induziert durch magnetische Spannungen, ein weit akzeptierter Erwärmungsmechanismus sind, sind sie im inneren Krustenbereich, in dem die Schmelztemperaturen hoch sind, am effektivsten, obwohl dort keine Wärmequellen zu erwarten sind. Ein alternativer Mechanismus postuliert, dass magnetische Energie durch Elektroneninfang an Kernen in der äußeren Kruste in Wärme umgewandelt wird, angetrieben durch den langsamen Zerfall des Magnetfeldes, der die Materie komprimiert. Frühere Schätzungen dieser Erwärmung stützten sich auf das HFB-24-Atommassenmodell; diese Studie untersucht die Rolle des Elektroneninfangs unter Verwendung aktualisierter Nuklear-Daten und der HFB-27-Atommasse-Tabelle, um Vorhersagen zur Wärmefreisetzung und zu den Quellenlokalisationen zu präzisieren.

Methodik
Die Autoren modellieren die äußere Kruste eines Magnetars als elektrisch neutrale Plasma aus Kernen $(A, Z)$ und freien Elektronen. Die Studie geht davon aus, dass die Kruste initial in vollem thermodynamischem Gleichgewicht unter einem starken Magnetfeld ($B_\star \equiv B/B_{rel} = 2000$, entsprechend $B \approx 8,8 \times 10^{16}$ G) ist. Wesentliche methodische Schritte umfassen:

• Quantisierungseffekte: Die Landau-Rabi-Quantisierung wird auf Elektronen angewendet, wobei angenommen wird, dass sie aufgrund des hohen Magnetfeldes und der relativ niedrigen Temperaturen ($T \sim 10^8 - 10^9$ K) das unterste Landau-Niveau besetzen.
• Gleichgewichtszusammensetzung: Die initiale Zusammensetzung wird durch Minimierung der Gibbs-Energie pro Nukleon mittels eines iterativen Ansatzes bestimmt.
• Reaktionsschwellen: Der Beginn der Elektroneninfänge wird analytisch berechnet. Der Schwellendruck $P_\beta$ für einen Kern $(A, Z)$, der ein Elektron einfängt, wird abgeleitet, wobei die Stärke des Magnetfeldes und die räumliche Anordnung der Kerne (Wigner-Seitz-Näherung) berücksichtigt werden.
• Wärmeberechnung: Die freigesetzte Wärme pro Kern wird geschätzt, indem die Energiedifferenz zwischen dem Mutter- und Tochterkern einschließlich der Anregungsenergien der Tochterzustände betrachtet wird. Die Studie berechnet die Wärmefreisetzung sowohl für Grundzustand-zu-Grundzustand- als auch für Grundzustand-zu-Angeregter-Zustand-Übergänge.
• Nuklearmodelle: Die Berechnungen nutzen experimentelle Kernmassen aus der Atomic Mass Evaluation 2016 sowie theoretische Vorhersagen aus dem HFB-27-Modell (basierend auf selbstkonsistenten deformierten Hartree-Fock-Bogoliubov-Berechnungen). Diese werden mit früheren Ergebnissen verglichen, die mit dem HFB-24-Modell erzielt wurden.
• Pyronukleare Fusion: Das Potenzial für die Fusion leichter Elemente (z. B. Kohlenstoff, Sauerstoff) wird als obere Grenze der Wärmefreisetzung in den dichtesten Regionen der äußeren Kruste betrachtet.

Wesentliche Ergebnisse

• Dominante Erwärmungsreaktionen: Die Studie identifiziert spezifische Elektroneninfang-Reaktionen, die die signifikantesten Mengen an Wärme freisetzen. Grundzustandsübergänge umfassen $^{12}\text{C} \to {}^{12}\text{Be}$ ($\sim 0,14$ MeV/Nukleon), $^{16}\text{O} \to {}^{16}\text{C}$ ($\sim 0,15$ MeV/Nukleon) und $^{62}\text{Cr} \to {}^{62}\text{Ti}$ ($\sim 0,09$ MeV/Nukleon). Grundzustand-zu-Angeregter-Zustand-Übergänge, wie $^{82}\text{Ge} \to {}^{82}\text{Zn}$ und verschiedene Strontium/Krypton-Übergänge, tragen zusätzliche Wärme ($\sim 0,07 - 0,1$ MeV/Nukleon) bei.
• Modellvergleich: Während die HFB-27- und HFB-24-Modelle ähnliche primäre Wärmequellen vorhersagen, treten signifikante Diskrepanzen bei Nebenreaktionen auf. Insbesondere ist die vorhergesagte Wärme für $^{56}\text{Ti} \to {}^{56}\text{Ca}$ und $^{80}\text{Zn} \to {}^{80}\text{Ni}$ Übergänge im HFB-27-Modell etwa um das 2- bzw. 3-fache höher, bedingt durch Unterschiede in den berechneten $Q_{EC}$-Werten.
• Wärmemagnitudo und Lokalisation: Die maximale freigesetzte Wärme pro Nukleon wird auf $\sim 0,1$ MeV durch Elektroneninfänge und $\sim 1 - 2$ MeV durch pyronukleare Fusionen geschätzt. Diese Reaktionen treten bei Dichten $\rho_\beta \sim 10^{10} - 10^{11}$ g cm$^{-3}$ und Drücken $P_\beta \sim 10^{29} - 10^{30}$ dyn cm$^{-2}$ auf.
• Zeitskalen und Leistung: Die Zeitskala für Elektroneninfänge ($\tau$) über verschiedene Krustenschichten hinweg liegt in der Größenordnung des kinematischen Alters von Magnetaren (einige tausend Jahre). Die resultierende gesamte Heizleistung wird auf $W_\infty \sim 10^{35} - 10^{36}$ erg/s geschätzt, was konsistent mit beobachteten Magnetar-Leuchtkräften und früheren Abkühlungssimulationen ist.

Bedeutung und Behauptungen
Das Paper kommt zu dem Schluss, dass Elektroneninfänge und pyronukleare Fusion in der äußeren Kruste, angetrieben durch den Zerfall des Magnetfeldes, eine robuste Quelle der internen Erwärmung für Magnetare darstellen. Eine zentrale Erkenntnis ist, dass die maximale Wärmefreisetzung im Wesentlichen durch Nukleardaten (Massen und Anregungsenergien) bestimmt wird und nicht durch die spezifische Stärke des Magnetfeldes oder den Phasenzustand (fest oder flüssig) des Elektron-Ion-Plasmas. Folglich bleiben die Schätzungen der Erwärmung auch für Neutronensterne mit geringeren Magnetfeldern eine gültige Annäherung, wennglich sich der physikalische Ort der Wärmequellen verschieben wird. Die Autoren merken an, dass dieser Mechanismus auch dann effektiv bleibt, wenn Teile der Kruste geschmolzen sind, was einen deutlichen Vorteil gegenüber Modellen bietet, die auf der Dissipation elastischer Energie durch Krustenversagen beruhen. Die Ergebnisse stehen im Einklang mit Neutronenstern-Abkühlungsdaten und den empirisch ermittelten Dichtebereichen der Wärmedeposition.