Internal Heating in Magnetars: Role of Electron Captures

Résumé technique : Chauffage interne des magnétars : rôle des captures électroniques

Énoncé du problème
Les magnétars, caractérisés par des champs magnétiques extrêmes ($B \gtrsim 10^{14}$ G), présentent des luminosités de rayons X persistantes ($\sim 10^{33} - 10^{35}$ erg/s) et des éruptions géantes occasionnelles qui dépassent leur production d'énergie de rotation. Bien que les déformations de la croûte induites par les contraintes magnétiques soient un mécanisme de chauffage largement accepté, elles sont plus efficaces dans la croûte interne où les températures de fusion sont élevées, alors qu'aucun source de chaleur n'est attendue dans ces régions profondes. Un mécanisme alternatif postule que l'énergie magnétique est convertie en chaleur via des captures d'électrons par les noyaux dans la croûte externe, pilotées par la décroissance lente du champ magnétique qui comprime la matière. Les estimations précédentes de ce chauffage reposaient sur le modèle de masse atomique HFB-24 ; cette étude examine le rôle des captures électroniques en utilisant des données nucléaires mises à jour et la table de masses atomiques HFB-27 afin d'affiner les prédictions de libération de chaleur et de localisation des sources.

Méthodologie
Les auteurs modélisent la croûte externe d'un magnétar comme un plasma chargé électriquement neutre composé de noyaux $(A, Z)$ et d'électrons libres. L'étude suppose que la croûte est initialement en plein équilibre thermodynamique sous un champ magnétique fort ($B_\star \equiv B/B_{rel} = 2000$, correspondant à $B \approx 8,8 \times 10^{16}$ G). Les étapes méthodologiques clés incluent :

• Effets de quantification : La quantification de Landau-Rabi est appliquée aux électrons, en supposant qu'ils occupent le niveau de Landau le plus bas en raison du champ magnétique élevé et des températures relativement basses ($T \sim 10^8 - 10^9$ K).
• Composition à l'équilibre : La composition initiale est déterminée en minimisant l'énergie libre de Gibbs par nucléon en utilisant une approche itérative.
• Seuils de réaction : Le déclenchement des captures électroniques est calculé de manière analytique. La pression seuil $P_\beta$ pour qu'un noyau $(A, Z)$ capture un électron est dérivée, en tenant compte de l'intensité du champ magnétique et de l'arrangement spatial des noyaux (approximation de Wigner-Seitz).
• Calcul de la chaleur : La chaleur libérée par noyau est estimée en considérant la différence d'énergie entre le noyau parent et le noyau fils, incluant les énergies d'excitation des états fils. L'étude calcule la libération de chaleur pour les transitions état fondamental-vers-état fondamental et état fondamental-vers-état excité.
• Modèles nucléaires : Les calculs utilisent les masses nucléaires expérimentales de l'évaluation des masses atomiques de 2016 et les prédictions théoriques du modèle HFB-27 (basé sur des calculs Hartree-Fock-Bogoliubov déformés et autocohérents). Celles-ci sont comparées aux résultats précédents obtenus avec le modèle HFB-24.
• Fusion pycnonucléaire : Le potentiel de fusion d'éléments légers (ex: Carbone, Oxygène) est considéré comme une limite supérieure de la libération de chaleur dans les régions les plus denses de la croûte externe.

Résultats clés

• Réactions de chauffage dominantes : L'étude identifie des réactions de capture électronique spécifiques qui libèrent les quantités de chaleur les plus significatives. Les transitions vers l'état fondamental incluent $^{12}\text{C} \to {}^{12}\text{Be}$ ($\sim 0,14$ MeV/nucléon), $^{16}\text{O} \to {}^{16}\text{C}$ ($\sim 0,15$ MeV/nucléon) et $^{62}\text{Cr} \to {}^{62}\text{Ti}$ ($\sim 0,09$ MeV/nucléon). Les transitions état fondamental-vers-état excité, telles que $^{82}\text{Ge} \to {}^{82}\text{Zn}$ et diverses transitions du Strontium/Krypton, contribuent à une chaleur supplémentaire ($\sim 0,07 - 0,1$ MeV/nucléon).
• Comparaison de modèles : Bien que les modèles HFB-27 et HFB-24 prédisent des sources de chaleur primaires similaires, des écarts importants apparaissent dans les réactions mineures. Spécifiquement, la chaleur prédite pour les transitions $^{56}\text{Ti} \to {}^{56}\text{Ca}$ et $^{80}\text{Zn} \to {}^{80}\text{Ni}$ est respectivement environ 2 et 3 fois plus élevée dans le modèle HFB-27 en raison des différences dans les valeurs de $Q_{EC}$ calculées.
• Magnitude et localisation de la chaleur : La chaleur maximale libérée par nucléon est estimée à $\sim 0,1$ MeV par les captures électroniques et $\sim 1 - 2$ MeV par les fusions pycnonucléaires. Ces réactions se produisent à des densités $\rho_\beta \sim 10^{10} - 10^{11}$ g cm$^{-3}$ et des pressions $P_\beta \sim 10^{29} - 10^{30}$ dyn cm$^{-2}$.
• Échelles de temps et puissance : L'échelle de temps pour les captures électroniques ($\tau$) à travers les différentes couches de la croûte est de l'ordre de l'âge cinématique des magnétars (quelques milliers d'années). La puissance de chauffage totale résultante est estimée à $W_\infty \sim 10^{35} - 10^{36}$ erg/s, ce qui est cohérent avec les luminosités observées des magnétars et les simulations de refroidissement antérieures.

Signification et affirmations
L'article conclut que les captures électroniques et la fusion pycnonucléaire dans la croûte externe, pilotées par la décroissance du champ magnétique, constituent une source robuste de chauffage interne pour les magnétars. Une conclusion clé est que la libération maximale de chaleur est essentiellement déterminée par les données nucléaires (masses et énergies d'excitation) plutôt que par l'intensité spécifique du champ magnétique ou l'état de phase (solide ou liquide) du plasma électron-ion. Par conséquent, les estimations de chauffage restent des approximations valides même pour les étoiles à neutrons possédant des champs magnétiques plus faibles, bien que la localisation physique des sources de chaleur sera déplacée. Les auteurs notent que ce mécanisme reste efficace même si certaines parties de la croûte sont fondues, offrant un avantage distinct par rapport aux modèles reposant sur la dissipation d'énergie élastique due à la rupture de la croûte. Les résultats sont cohérents avec les données de refroidissement des étoiles à neutrons et les plages de densité déterminées empiriquement pour le dépôt de chaleur.