# Heating in magnetar crusts from electron captures

Nicolas Chamel, Anthea Francesca Fantina, Lami Suleiman, Julian-Leszek Zdunik, Pawel Haensel

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Nicolas Chamel, Anthea Francesca Fantina, Lami Suleiman, Julian-Leszek Zdunik, Pawel Haensel. Heating in magnetar crusts from electron captures. Universe, 2021, 7 (6), pp.193. 10.3390/universe7060193. [hal-03260055]

HAL Id: hal-03260055

https://hal.science/hal-03260055

Submitted on 21 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Article

# Heating in Magnetar Crusts from Electron Captures

Nicolas Chamel <fnref n="1" />*, Anthea Francesca Fantina <fnref n="1" />,<fnref n="2" /> ID, Lami Suleiman <fnref n="3" />,<fnref n="4" /> ID, Julian-Leszek Zdunik <fnref n="3" /> ID
and Pawel Haensel <fnref n="3" /> ID

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Chamel, N.; Fantina, A. F.; Suleiman, L.; Zdunik, J.-L.; Haensel, P. Heating in Magnetar Crusts from Electron Captures. *Universe* **2021**, *7*, 193. https://doi.org/10.3390/universe7060193

Academic Editors: Isaac Vidaña and
Marcello Baldo

Received: 7 May 2021

Accepted: 4 June 2021

Published: 8 June 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The persistent thermal luminosity of magnetars and their outbursts suggest the existence of some internal heat sources located in their outer crust. The compression of matter accompanying the decay of the magnetic field may trigger exothermic electron captures and, possibly, pycnonuclear fusions of light elements that may have been accreted onto the surface from the fallback of supernova debris, from a disk or from the interstellar medium. This scenario bears some resemblance to deep crustal heating in accreting neutron stars, although the matter composition and the thermodynamic conditions are very different. The maximum possible amount of heat that can be released by each reaction and their locations are determined analytically taking into account the Landau–Rabi quantization of electron motion. Numerical results are also presented using experimental, as well as theoretical nuclear data. Whereas the heat deposited is mainly determined by atomic masses, the locations of the sources are found to be very sensitive to the magnetic field strength, thus providing a new way of probing the internal magnetic field of magnetars. Most sources are found to be concentrated at densities 10¹⁰–10¹<fnref n="1" /> g cm⁻<fnref n="3" /> with heat power W∞ ~ 10³⁵–10³⁶ erg/s, as found empirically by comparing cooling simulations with observed thermal luminosity. The change of magnetic field required to trigger the reactions is shown to be consistent with the age of known magnetars. This suggests that electron captures and pycnonuclear fusion reactions may be a viable heating mechanism in magnetars. The present results provide consistent microscopic inputs for neutron star cooling simulations, based on the same model as that underlying the Brussels-Montreal unified equations of state.

**Keywords:** neutron star; magnetar; outburst; magnetic field; heating; cooling; electron capture; pycnonuclear fusion

## 1. Introduction

Magnetars form a subclass of known neutron stars exhibiting various astrophysical phenomena powered by their extreme magnetic field [1] (see, e.g., [2] for a recent review) and currently consisting of 12 confirmed soft gamma-ray repeaters and 12 anomalous X-ray pulsars [3] according to the McGill Online Magnetar Catalog<fnref n="1" />. In particular, enhancements of the X-ray flux by several orders of magnitude lasting for weeks or even years have been observed in these stars [4]. It is widely thought that these “outbursts” are caused by some sort of internal heat deposition; however, a detailed understanding is still lacking. Internal heating may also explain why magnetars are hotter than their weakly magnetized relatives. Heat may come from the dissipation of mechanical energy during crustquakes (see, e.g., [5]) or from Ohmic dissipation. However, both scenarios have shortcomings, as discussed in detail in [6]. In particular, these mechanisms are only effective deep enough beneath the surface of the star where the temperature lies below the melting temperature

<footnote n="1">¹ Institute of Astronomy and Astrophysics, Université Libre de Bruxelles, CP 226, Boulevard du Triomphe, 1050 Brussels, Belgium</footnote>

<footnote n="2">² Grand Accélérateur National d'Ions Lourds (GANIL), CEA/DRF-CNRS/IN2P3, Boulevard Henri Becquerel, 14076 Caen, France; anthea.fantina@ganil.fr</footnote>

<footnote n="3">³ N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, 00-716 Warszawa, Poland; lsuleiman@camk.edu.pl (L.S.); jlz@camk.edu.pl (J.-L.Z.); haensel@camk.edu.pl (P.H.)</footnote>

<footnote n="4">⁴ Laboratoire Univers et Théories, Observatoire de Paris, CNRS, Université de Paris, 92195 Meudon, France</footnote>

<footnote>* Correspondence: nicolas.chamel@ulb.be</footnote>
or where the electric conductivity is sufficiently low. However, it has been shown that the sources should be located in the shallow region of the crust to avoid excessive neutrino losses [7,8].

In this paper, we focused on the scenario proposed in [9], who pointed out that the decay of the magnetic field may trigger exothermic nuclear reactions (see, e.g., [6] for a critical review of other scenarios). This mechanism is similar to deep crustal heating in accreting neutron stars (see, e.g., [10] and the references therein), the compression of matter being induced here by the loss of magnetic pressure (and more generally by the local rearrangement of the magnetic field lines), rather than from the accumulation of accreted material onto the stellar surface (the existence of accreting magnetars—conceivably in the form of ultraluminous X-ray sources—remains a matter of debate; see, e.g., [11–14] and the references therein). The authors of [9] implicitly assumed that the same amount of heat is released in both cases and at similar locations, although the conditions prevailing in these two classes of neutron stars are very different. Finally, nuclear reactions may also be triggered by the spin-down of the star due to electromagnetic radiation [15].

The crustal heating scenario of [9,15] was further examined here from the nuclear physics perspective. Our microscopic model of magnetar crusts and our analysis of electron captures by nuclei and pycnonuclear fusion reactions caused by matter compression are described in Section 2. Our results for the calculated heat and its location are presented in Section 3. The astrophysical implications are discussed in Section 4.

## 2. Microphysics of Magnetar Crusts

We closely followed our recent analysis for (unmagnetized) accreting neutron stars [16], which we adapted to the present context.

### 2.1. Thermodynamic Conditions

Let us consider an electrically charged neutral matter element composed of fully ionized atomic nuclei (A, Z) with proton number Z and mass number A embedded in a relativistic electron gas.

In the presence of a magnetic field, the electron motion perpendicular to the field is quantized into Landau–Rabi levels [17,18]. All electrons are confined to the lowest level whenever the temperature T lies below (see, e.g., Chapter 4 of [19]):

$$T_B = \frac{m_e c^2}{k_B} B_* \approx 5.93 \times 10^9 B_* \text{ K}, \quad (1)$$

and provided the mass density ρ does not exceed:

$$\rho_B = \frac{A}{Z} m_u \frac{B_*^{3/2}}{\sqrt{2\pi^2\lambda_e^3}} \approx 2.07 \times 10^6 \frac{A}{Z} B_*^{3/2} \text{ g cm}^{-3}, \quad (2)$$

where $k_B$ denotes Boltzmann's constant, c is the speed of light, $m_e$ is the electron mass, $m_u$ is the unified atomic mass unit, and $\lambda_e = \hbar/(m_e c)$ is the electron Compton wavelength ($\hbar$ being the Dirac–Planck constant), and we introduced the dimensionless magnetic field strength $B_* \equiv B/B_{rel}$ with:

$$B_{rel} = \left( \frac{m_e c^2}{\alpha \lambda_e^3} \right)^{1/2} \approx 4.41 \times 10^{13} \text{ G}, \quad (3)$$

where $\alpha = e^2/(\hbar c)$ is the fine-structure constant ($e$ being the elementary electric charge). The condition $T < T_B$ is fulfilled in magnetars since $B_* \gg 1$, and the temperature is typically of order $10^8–10^9$ K (see, e.g., [8]). The threshold density $\rho_B$ coincides with the density associated with the neutron drip transition (delimiting the outer and inner parts of the crust) for $B \approx 5.72 \times 10^{16}$ G [20]. For the densities and the magnetic field strengths that we considered, electrons remain in the strongly quantizing regime.
Whether the electron gas is degenerate or not is determined by the Fermi temperature given by [19]:

$$T_{\text{Fe}} = \frac{m_e c^2}{k_B} (\gamma_e - 1) \approx 5.93 \times 10^9 \left( \sqrt{1 + 1.02 \left[ \frac{4 Z}{3 A} \rho_6 \left( \frac{\rho}{\rho_B} \right)^2 \right]^{2/3}} - 1 \right) \text{K}, \quad (4)$$

where $\gamma_e$ is the electron Fermi energy in units of the electron mass energy $m_e c^2$, and $\rho_6 = \rho/10^6 \text{ g cm}^{-3}$. We estimated this temperature in the different layers of the outer crust of a magnetar using the composition calculated in [21] for $B_* = 2000$ ($B \approx 8.83 \times 10^{16} \text{ G}$) and $B_* = 3000$ ($B \approx 1.32 \times 10^{17} \text{ G}$). As we show, electron captures occur at pressures well above $10^{28} \text{ dyn cm}^{-2}$ (densities above $10^{10} \text{ g cm}^{-3}$). As can be seen in Figures 1 and 2, $T_{\text{Fe}}$ exceeds $10^{10} \text{ K}$ for the crustal regions of interest; therefore, electrons are highly degenerate.

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Characteristic temperatures in the outer crust of a magnetar with a magnetic field strength $B_* = 2000$ ($B \approx 8.83 \times 10^{16} \text{ G}$) as a function of the pressure $P$ (in dyn cm$^{-2}$): electron Fermi temperature $T_{\text{Fe}}$ (solid line) separating the quantum (degenerate gas) and classical regimes and melting temperature $T_m$ (dotted line) separating the liquid and solid phases (assuming $\Gamma_m = 175$). The Coulomb plasma temperature $T_\ell$ separating the weakly and strongly coupled plasma regimes is much higher and is not shown. The magnetic field is strongly quantizing for all pressures and temperatures shown (in the present case, $T_B = 1.19 \times 10^{13} \text{ K}$): electrons are confined to the lowest Landau-Rabi level. Composition taken from [21].</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Same as Figure 1 for $B_{\star} = 3000$ ($B \approx 1.32 \times 10^{17}$ G). The magnetic field is strongly quantizing for all pressures and temperatures shown ($T_B = 1.78 \times 10^{13}$ K).</figcaption></figure>

At these temperatures, the electron-ion plasma may not form a crust, but may be in a liquid state, especially in the shallow layers. The crystallization temperature is given by [19]:

$$T_m = \frac{Z^2 e^2}{a_N k_B \Gamma_m} \approx 1.3 \times 10^5 Z^2 \left(\frac{175}{\Gamma_m}\right) \left(\frac{\rho_6}{A}\right)^{1/3} \text{ K}, \quad (5)$$

where $a_N = (3Z/(4\pi n_e))^{1/3}$ is the ion sphere’s radius, $n_e$ the electron number density, and $\Gamma_m$ the Coulomb coupling parameter at melting. As shown in Figures 1 and 2, $T_m$ is of order $10^9$ K (assuming $\Gamma_m = 175$ as for unmagnetized matter), which is comparable to the temperature prevailing in magnetar crusts. However, the situation remains uncertain since the presence of a strong magnetic field tends to stabilize the solid phase by lowering $\Gamma_m$, therefore increasing $T_m$; see, e.g., [22]. In either case, the Coulomb plasma remains strongly coupled since the temperature $T$ is much smaller than $T_{\ell} = \Gamma_m T_m$, i.e., the Coulomb potential energy $Z^2e^2/a_N$ is much larger than the thermal energy $k_B T$.

In the following, we thus neglect thermal effects.

## 2.2. Compression-Induced Nuclear Processes

Ignoring the negligible magnetic susceptibility of the outer crust [23] and recalling that the electron gas is highly degenerate, any increase of the matter pressure $P$ (caused by the decay of the magnetic field or the spin-down of the star) must necessarily be accompanied by an increase of the baryon chemical potential $\mu$ (as shown in Appendix A of Ref. [24], the baryon chemical potential coincides with the Gibbs free energy per nucleon of the electron-ion plasma) according to the Gibbs-Duhem relation $nd\mu \approx dP$, where $n$ is the baryon number density, and we neglected the thermal contribution $sdT$ (s being the entropy density), see, e.g., Appendix A of [24]. The compression of a matter element may thus trigger electron captures by nuclei ($A, Z$) with the emission of neutrinos at some pressure $P_{\beta}$ (depending on the nucleus under consideration, as well as the magnetic field strength):

$$(A, Z) + e^{-} \rightarrow (A, Z-1) + v_{e}. \quad (6)$$

Recalling that the magnetic field in the crust of a neutron star roughly evolves on timescales $\tau_B$ of order $\sim$ Myr (see, e.g., [5,25–33]), the above reaction essentially proceeds in quasi-equilibrium with no energy release. The threshold pressure of the reaction (6) can thus
be determined by comparing the baryon chemical potentials of the parent and daughter nuclei. The latter are generally highly unstable, and the capture of a second electron (at the same pressure) may thus be accompanied by an energy release Q per nucleus:

$$ (A, Z - 1) + e^{-} \rightarrow (A, Z - 2) + v_{e} + Q. \qquad (7) $$

We assumed that the temperature is sufficiently low that the parent nucleus is in its ground state. According to the Fermi and Gamow-Teller selection rules, the reaction (6) will generally involve a transition to the lowest excited state of the daughter nucleus. Alternatively, a transition to the ground state may actually be a more realistic scenario in this astrophysical context given the very long time scales (~Myr) for the compression of matter. We studied both types of reactions. Since we were interested in the maximum possible amount of heat that can be potentially released, we calculated Q considering that the daughter nuclei (A, Z - 2) are in their ground state. We also considered pycnonuclear fusion reactions of light elements [34,35]. Such elements may have been accreted onto the neutron star surface from the fallback of supernova debris after the explosion, from a remnant disk around the star [36,37], as observed in AXP 4U 0142+61 [38], or from the interstellar medium.

As nuclei sink deeper into the crust, further compression may give rise to delayed neutron emission, thus marking the transition to the inner crust [20,39]:

$$ (A, Z) + e^{-} \rightarrow (A - \Delta N, Z - 1) + \Delta N n + v_{e}. \qquad (8) $$

## 2.3. Baryon Chemical Potential and Matter Pressure

The model we adopted here was described in [40]. The baryon chemical potential is given by:

$$ \mu = \frac{M'(A, Z)c^2}{A} + \frac{Z}{A}m_e c^2 \left[ \gamma_e - 1 + \frac{4}{3}C\alpha\lambda_e n_e^{1/3}Z^{2/3} \right], \qquad (9) $$

where $M'(A, Z)$ denotes the nuclear mass (including the rest mass of Z protons, A - Z neutrons, and Z electrons—the reason for including the electron rest mass in $M(A, Z)$ is that atomic masses are generally tabulated rather than nuclear masses). We approximated the constant C of the electrostatic correction by the Wigner-Seitz estimate [41]:

$$ C = -\frac{9}{10} \left( \frac{4\pi}{3} \right)^{1/3} \approx -1.4508. \qquad (10) $$

For nuclei in excited states, the baryon chemical potential, which we denote by $\mu^*$, takes a similar form except that $M'(A, Z)c^2$ must be replaced by $M'(A, Z)c^2 + E_{ex}(A, Z)$, where $E_{ex}(A, Z)$ is the excitation energy.

Ignoring the small anomalous magnetic moment of electrons, the thermodynamic matter pressure, which consists of the pressure $P_e$ of the electron Fermi gas and the lattice contribution $P_L$, is expressible as (see, e.g., Chp. 4 of [19]):

$$ P = \frac{B_{\star}m_{e}c^{2}}{(2\pi)^{2}\lambda_{e}^{3}} \left[ x_{e}\sqrt{1+x_{e}^{2}} - \ln\left(x_{e}+\sqrt{1+x_{e}^{2}}\right) \right] + \frac{1}{3}C\alpha\hbar c n_{e}^{4/3}Z^{2/3}, \qquad (11) $$

where we introduced the dimensionless parameter:

$$ x_e = \sqrt{\gamma_e^2 - 1} = \frac{2\pi^2\lambda_e^3 n_e}{B_*}. \qquad (12) $$

Since the thermal pressure of nuclei is neglected, the pressure is the same whether nuclei are in their ground state or in an excited state.
## 2.4. Onset of Electron Captures

The onset of electron captures by nuclei (A, Z) is determined by the condition $\mu(A, Z, P_β, B_★) = \mu(A, Z - 1, P_β, B_★)$, assuming that the daughter nucleus is in its ground state. This condition takes the same form as in the absence of magnetic fields and can be expressed to first order in $\alpha$ as [24]:

$$ \gamma_e + C\alpha\lambda_e n_e^{1/3}F(Z) = \gamma_e^β(A, Z), \quad (13) $$

$$ F(Z) = Z^{5/3} - (Z - 1)^{5/3} + \frac{1}{3}Z^{2/3}, \quad (14) $$

$$ \gamma_e^β(A, Z) = -\frac{Q_{EC}(A, Z)}{m_e c^2} + 1, \quad (15) $$

where we introduced the Q-value (in vacuum) associated with electron capture by nuclei (A, Z):

$$ Q_{EC}(A, Z) = M'(A, Z)c^2 - M'(A, Z - 1)c^2. \quad (16) $$

These Q-values can be obtained from the tabulated Q-values of $β$ decay by the following relation:

$$ Q_{EC}(A, Z) = -Q_β(A, Z - 1). \quad (17) $$

Note that if $Q_{EC}(A, Z) > 0$, the nucleus (A, Z) is unstable against electron captures at any density.

The threshold condition (13) is amenable to analytical solutions if the electron density in the second term of the left-hand side (electrostatic correction) is expressed in terms of the electron Fermi energy using the ultrarelativistic approximation $\gamma_e \gg 1$ in Equation (12):

$$ n_e \approx \frac{B_★}{2\pi^2\lambda_e^3}\gamma_e. \quad (18) $$

Introducing

$$ \bar{F}(Z, B_★) = \frac{1}{3}C\alpha F(Z) \left( \frac{B_★}{2\pi^2} \right)^{1/3} < 0, \quad (19) $$

Equation (13) thus reduces to:

$$ \gamma_e + 3\bar{F}(Z, B_★)\gamma_e^{1/3} = \gamma_e^β. \quad (20) $$

Introducing the dimensionless parameter:

$$ v \equiv \frac{\gamma_e^β}{2|\bar{F}(Z, B_★)|^{3/2}}, \quad (21) $$

the solutions (real roots) of Equation (20) for $\gamma_e$ are given by the following formulas [16]:

$$ \gamma_e = \begin{cases} 8|\bar{F}(Z, B_★)|^{3/2} \cosh^3\left(\frac{1}{3}\operatorname{arccosh} v\right) & \text{if } v \ge 1, \\ 8|\bar{F}(Z, B_★)|^{3/2} \cos^3\left(\frac{1}{3}\operatorname{arccos} v\right) & \text{if } 0 \le v < 1. \end{cases} \quad (22) $$

Using Equation (11), the threshold pressure and baryon number density are thus given by:

$$ P_β(A, Z, B_★) = \frac{B_★ m_e c^2}{4\pi^2 \lambda_e^3} \left[ \gamma_e \sqrt{\gamma_e^2 - 1} - \ln\left(\sqrt{\gamma_e^2 - 1} + \gamma_e\right) + \frac{C\alpha}{3} \left(\frac{4B_★ Z^2}{\pi^2}\right)^{1/3} \left(\gamma_e^2 - 1\right)^{2/3} \right], \quad (23) $$

$$ n_β(A, Z, B_★) = \frac{B_★ A}{2\pi^2 \lambda_e^3 Z} \sqrt{\gamma_e^2 - 1}, \quad (24) $$
respectively. The transition is accompanied by a discontinuous change of density given by:

$$ \frac{\Delta n}{n_{\beta}} = \frac{Z}{Z-1} \left\{ 1 + \frac{1}{3} C \alpha \left( \frac{B_{\star}}{2\pi^2} \right)^{1/3} \left[ Z^{2/3} - (Z-1)^{2/3} \right] \frac{\gamma_e}{(\gamma_e^2 - 1)^{5/6}} \right\} - 1. \quad (25) $$

The threshold condition for transitions to daughter nuclei in excited states is $g(A, Z, P_{\beta}^*, B_*) = g^*(A, Z-1, P_{\beta}^*, B_*)$. The corresponding pressure $P_{\beta}^*$ and density $n_{\beta}^*(A, Z)$ can be obtained from the previous formulas by merely substituting $M'(A, Z-1)c^2$ with $M'(A, Z-1)c^2 + E_{ex}(A, Z-1)$. Because $E_{ex}(A, Z-1) \ge 0$, such transitions occur at higher pressure $P_{\beta}^* \ge P_{\beta}$ and density $n_{\beta}^*(A, Z, B_*) \ge n_{\beta}(A, Z, B_*)$.

As discussed in [39,42], the first electron capture by the nucleus (A, Z) may be accompanied by the emission of $\Delta N > 0$ neutrons. The corresponding pressure $P_{drip}$ and baryon density $n_{drip}$ are given by similar expressions as Equations (23) and (24), respectively, except that the threshold electron Fermi energy $\gamma_e^\beta$ is now replaced by $\gamma_e^\text{drip} \equiv \mu_e^\text{drip} / (m_e c^2)$ with:

$$ \mu_e^{\text{drip}}(A, Z) = M'(A - \Delta N, Z - 1)c^2 - M'(A, Z)c^2 + \Delta N m_n c^2 + m_e c^2, \quad (26) $$

assuming that the daughter nucleus is in the ground state. Transitions to excited states can be taken into account by adding the suitable excitation energy $E_{ex}(A - \Delta N, Z - 1)$ to $M'(A - \Delta N, Z - 1)c^2$. Neutron emission will thus occur whenever $\mu_e^{\text{drip}}(A, Z) < \mu_e^\beta(A, Z)$. All the expressions provided here assume that electrons are confined to the lowest Landau-Rabi level in all layers of the outer crust, which translates into the condition that the solution for $\gamma_e$ at the neutron drip point must obey the following inequality (see, e.g., [40]):

$$ \gamma_e \le \sqrt{1 + 2B_*}. \quad (27) $$

### 3. Crustal Heating

#### 3.1. Heat Released by Electron Captures

The heat released by electron captures can be determined analytically as follows. The reaction (6) will be generally almost immediately followed by a second electron capture (7) on the daughter nucleus provided $\mu(A, Z-2, P_{\beta}, B_*) < \mu(A, Z, P_{\beta}, B_*)$ or $\mu(A, Z-2, P_{\beta}^*, B_*) < \mu(A, Z, P_{\beta}^*, B_*)$ depending on whether the daughter nucleus after the first capture is in its ground state or in an excited state, respectively. The maximum possible amount of heat deposited in matter per one nucleus is given by:

$$ Q(A, Z, B_*) = A \left[ \mu(A, Z, P_{\beta}, B_*) - \mu(A, Z-2, P_{\beta}, B_*) \right] \quad (28) $$

in the first case and:

$$ Q^*(A, Z, B_*) = A \left[ \mu(A, Z, P_{\beta}^*, B_*) - \mu(A, Z-2, P_{\beta}^*, B_*) \right], \quad (29) $$

in the second case. The corresponding amounts of heat per one nucleon are given by $q(A, Z, B_*) \equiv Q(A, Z, B_*)/A$ and $q^*(A, Z, B_*) \equiv Q^*(A, Z, B_*)/A$, respectively. These estimates represent upper limits since part of the energy is radiated away by neutrinos.

Let us first consider ground-state-to-ground-state transitions. The two successive electron captures are accompanied by a small discontinuous change $\delta n_e$ of the electron density $n_e$. Requiring the pressure to remain fixed $P_{\beta}(A, Z, B_*) = P(A, Z-2, B_*)$ leads to:

$$ \delta n_e \approx \frac{C \alpha \lambda_e}{3} \left[ Z^{2/3} - (Z-2)^{2/3} \right] \frac{dn_e}{d\gamma_e} n_e^{1/3}, \quad (30) $$
where we used Equation (11) and the relation² $dP_e = n_e m_e c^2 d\gamma_e$. Expanding the electron Fermi energy $\gamma_e(n_e + \delta n_e)$ to first order in $\delta n_e / n_e$, substituting in the expression of $\mu(A, Z - 2, P_\beta, B_\star)$, and using Equations (13)–(15) and (30), we finally obtain:

$$Q(A, Z, B_*) = Q_{EC}(A, Z-1) - Q_{EC}(A, Z) - m_e c^2 C \alpha \left(\frac{B_*}{2\pi^2}\right)^{1/3} (\gamma_e^2 - 1)^{1/6} \left[Z^{5/3} + (Z-2)^{5/3} - 2(Z-1)^{5/3}\right], \quad (31)$$

where $\gamma_e$ is calculated from Equation (22). Equation (31) is only valid if $\mu(A, Z - 2, P_\beta, B_\star) < \mu(A, Z, P_\beta, B_\star)$, which is generally satisfied for even A nuclei, but not necessarily for odd A nuclei. In the latter case, we typically have $Q_\beta(A, Z - 1, B_\star) < Q_\beta(A, Z - 2, B_\star)$. Using Equation (15), this implies that $\gamma_e^2(A, Z) < \gamma_e^2(A, Z - 1)$. In other words, as the pressure reaches $P_\beta(A, Z, B_\star)$, the nucleus (A, Z) decays, but the daughter nucleus (A, Z - 1) is actually stable against electron capture, and therefore, no heat is released $Q(A, Z, B_\star) = 0$. The daughter nucleus sinks deeper in the crust and only captures a second electron in quasi-equilibrium at pressure $P_\beta(A, Z - 1, B_\star) > P_\beta(A, Z, B_\star)$.

Equation (31) shows that the heating associated with electron captures is very weakly dependent on the spatial arrangement of ions since the electron-ion and ion-ion interactions only enter through the very small electrostatic correction proportional to the structure constant C, as confirmed by numerical calculations, which will be presented in Section 3.3. Moreover, it is worth recalling that the variations of the constant C between the liquid and solid phases are very small, from $C \approx -1.4621$ for the liquid (this value was calculated using the constant $A_1$ on page 75 of [19] as $C = A_1(4\pi/3)^{1/3}$) to $C \approx -1.4442$ for a perfect body-centered cubic crystal [19]. Whether the plasma is liquid or solid is therefore not expected to have any significant impact on the heating, contrary to the more popular mechanism involving crustal failures.

The heat $Q^*(A, Z, B_\star)$ released from the ground-state-to-excited-state transitions can be obtained from Equation (31) by substituting $M'(A, Z - 1)c^2$ with $M'(A, Z - 1)c^2 + E_{ex}(A, Z - 1)$, leading to:

$$Q^*(A, Z, B_\star) \approx Q(A, Z, B_\star) + 2E_{ex}(A, Z - 1) > Q(A, Z, B_\star). \quad (32)$$

For odd A nuclei, heat can thus only possibly be released if the first electron capture proceeds via excited states of the daughter nucleus and provided:

$$2E_{ex}(A, Z - 1) - Q_{\beta}(A, Z - 2, B_{\star}) + Q_{\beta}(A, Z - 1, B_{\star}) > 0. \quad (33)$$

It should be remarked that once the excited nuclei have decayed to their ground state and the nuclear equilibrium has been attained, the threshold densities and pressures between the different crustal layers will be given by Equations (23) and (24), respectively.

## 3.2. Heat Released by Pycnonuclear Fusions

In the densest regions of the outer crust, light nuclei may undergo pycnonuclear fusion reactions above some pressure $P_{pyc}$. The daughter nuclei are usually highly unstable against electron captures at that pressure, and thus decay by releasing some additional heat provided $P_\beta(2A, 2Z) < P_\beta(A, Z)$ or $P_\beta^*(2A, 2Z) < P_\beta^*(A, Z)$ depending on whether transitions to the ground state or excited state are considered. The fusion rates remain highly uncertain [43], and for this reason, the pressure $P_{pyc}$ is very difficult to estimate. Still, an upper limit on the heat released per nucleon can be obtained by setting $P_{pyc} = P_\beta(A, Z)$ or $P_{pyc} = P_\beta^*(A, Z)$, respectively:

$$q_{pyc}(A, Z) = \mu(A, Z, P_\beta) - \mu(2A, 2Z - 2, P_\beta), \quad (34)$$

or:
$$q_{\text{pyc}}^{*}(A, Z) = \mu(A, Z, P_{\beta}^{*}) - \mu(2A, 2Z - 2, P_{\beta}^{*}). \qquad (35)$$

Expanding to first order the electron density $n_e$ associated with the layer containing nuclei $(2A, 2Z - 2)$, using Equations (20) and (30), we find:

$$q_{\text{pyc}}(A, Z) = \frac{M'(A, Z - 1)c^2}{A} - \frac{M'(2A, 2Z - 2)c^2}{2A} + \frac{m_e c^2}{A} C\alpha(1 - 2^{2/3}) \left(\frac{B_* \gamma_e}{2\pi^2}\right)^{1/3} (Z - 1)^{5/3}, \qquad (36)$$

with $\gamma_e$ from Equation (22). The heat actually consists of two contributions: the first from the fusion itself and the second from the subsequent electron captures by the highly unstable nuclei $(2A, 2Z)$. The heat released by the sole fusion is given by:

$$q_{\text{fus}}(A, Z) = \frac{M'(A, Z)c^2}{A} - \frac{M'(2A, 2Z)c^2}{2A} + \frac{m_e c^2}{A} C\alpha(1 - 2^{2/3}) \left(\frac{B_* \gamma_e}{2\pi^2}\right)^{1/3} Z^{5/3}. \qquad (37)$$

The expressions for $q_{\text{pyc}}^*(A, Z)$ and $q_{\text{fus}}^*(A, Z)$ follow after substituting $M'(A, Z - 1)c^2$ with $M'(A, Z - 1)c^2 + E_{\text{ex}}(A, Z - 1)$. Numerical results are presented and discussed in the next section.

### 3.3. Results and Discussions

We estimated the amount of heat deposited in the outer crust of a magnetar and the locations of the sources from the analytical formulas presented in the previous sections using the experimental atomic masses and the $Q_β$ values from the 2016 Atomic Mass Evaluation [44] supplemented with the microscopic atomic mass table HFB-24 [45] from the BRUSLIB database³ [46]. These same nuclear data have already been used to determine the equilibrium composition of the outer crust of a magnetar [21]. In particular, we took as the input for our calculations the results given in Tables II and III of [21]. We also considered light elements that may have been accreted onto the neutron star surface. We focused in particular on carbon and oxygen. As a matter of fact, the presence of carbon has been inferred in the atmosphere of the isolated neutron star in Cassiopeia A [47]. The nuclear masses $M'(A, Z)$ were obtained from measured atomic masses after subtracting out the binding energy of atomic electrons according to Equation (A4) of [48]. Excitation energies were taken from the Nuclear Data section of the International Atomic Energy Agency website⁴ following the Gamow–Teller selection rules, namely that the parity of the final state is the same as that of the initial state, whereas the total angular momentum $J$ can either remain unchanged or vary by $\pm\hbar$ (excluding transitions from $J = 0$ to $J = 0$).

The main heat sources associated with electron captures are plotted in Figures 3 and 4 as a function of the pressure $P$ and the mass density $\rho = nm_u$ for two different values of the magnetic field strength, namely $B_* = 2000$ and $B_* = 3000$, respectively. We checked that no neutron emission occurred for the transitions that are shown. Most of the heat is released at densities and pressures that are substantially higher than in accreting neutron stars [16]. Quite remarkably, heat sources are not as uniformly distributed as in accreting neutron stars (compare with Figure 1 of [16]), but are concentrated at densities $10^{10}–10^{11}$ g cm$^{-3}$ (corresponding to pressures $10^{29}–10^{30}$ dyn cm$^{-2}$). Comparing Figures 3 and 4 shows that the locations of heat sources are systematically shifted to higher densities with increasing magnetic field strength, whereas the sources remain essentially unchanged. These results can be easily understood from Equations (24) and (31). In particular, the amount of heat deposited is mainly determined by nuclear masses (more precisely by the relevant Q-values); the magnetic field only enters through the small electrostatic correction proportional to the fine structure constant. On the contrary, the
threshold pressure and the density both increase linearly with $B_*$ in the ultrarelativistic regime $\gamma_e \gg 1$,

$$P_{\beta}(A, Z, B_{*}) \approx \frac{B_{*}m_{e}c^{2}}{4\pi^{2}\lambda_{e}^{3}}\gamma_{e}^{2}, \qquad (38)$$

$$n_{\beta}(A, Z, B_{*}) \approx \frac{B_{*} A}{2\pi^{2}\lambda_{e}^{3} Z} \gamma_{e}^{\beta}, \qquad (39)$$

where we neglected the small electrostatic correction. These expressions also show that knowing the depth where the main heat sources should be located could thus be used to directly probe the internal magnetic field in the outer crust. For the magnetic fields expected in magnetars, the range of densities where most of the heat from nuclear reactions is deposited corresponds to that determined empirically by comparing cooling simulations with observations [7,8].

<figure><img src="image_8.png" /><figcaption>Figure 3. Heat $q$ released per nucleon (in MeV) from selected electron captures by nuclei in the outer crust of a magnetar with a magnetic field strength $B_* = 2000$ in a pressure $P$ (in dyn cm$^{-2}$)-mass density $\rho$ (in g cm$^{-3}$) diagram, considering transitions from the ground state of the parent nucleus to either the ground state (squares) or the first excited state (triangles) of the daughter nucleus. The size of each symbol is proportional to the amount of heat deposited as indicated.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>Figure 4. Same as Figure 3 for $B_* = 3000$.</figcaption></figure>

As expected, transitions from the ground state of the parent nucleus to the first excited state of the daughter nucleus are more exothermic than ground-state-to-ground-state transitions. Moreover, the heat released from electron captures by odd nuclei is negligible.
All reactions are listed in Tables 1 and 2 for $B_{\star} = 2000$ and $B_{\star} = 3000$, respectively. We also indicated in Tables 3 and 4 the reactions involving light elements that may have been accreted. In the shallow regions of the outer crust that are most relevant for the present discussion, results are completely determined by experimental measurements. It is only in the densest parts that theoretical atomic masses are needed. For the model adopted here, being microscopically grounded, its predictions for nuclei far from the stability valley are expected to remain fairly accurate (see Table 2 of [49], where the predictions of HFB-24 were compared to recent atomic mass measurements not included in the original fit). The sensitivity with respect to different mass models was discussed in [10] in the context of accreted neutron star crusts.

<table><caption>**Table 1.** Maximum possible heat released *q* per nucleon (in MeV) from electron captures in the outer crust of a magnetar with $B_{\star} = 2000$. The pressure $P_{\beta}$ (in dyn cm$^{-2}$) and density $\rho_{\beta}$ (in g cm$^{-3}$) at which electron capture occurs are given considering ground-state-to-ground-state transitions (first line of reaction) and ground-state-to-excited-state transitions (second line). The initial composition was taken from [21]. The symbol (⋆) is used to distinguish reactions for which theoretical atomic masses were needed.</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>6.62 × 10<sup>28</sup></td><td>6.31 × 10<sup>10</sup></td><td><sup>56</sup>Fe → <sup>56</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.037</td></tr><tr><td>6.94 × 10<sup>28</sup></td><td>6.46 × 10<sup>10</sup></td><td><sup>56</sup>Fe → <sup>56</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.041</td></tr><tr><td>3.06 × 10<sup>29</sup></td><td>1.44 × 10<sup>11</sup></td><td><sup>56</sup>Cr → <sup>56</sup>Ti − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.041</td></tr><tr><td>7.03 × 10<sup>29</sup></td><td>2.36 × 10<sup>11</sup></td><td><sup>56</sup>Ti → <sup>56</sup>Ca − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.023</td></tr><tr><td>1.23 × 10<sup>29</sup></td><td>8.75 × 10<sup>10</sup></td><td><sup>62</sup>Ni → <sup>62</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.045</td></tr><tr><td>1.43 × 10<sup>29</sup></td><td>9.43 × 10<sup>10</sup></td><td><sup>62</sup>Ni → <sup>62</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.061</td></tr><tr><td>3.88 × 10<sup>29</sup></td><td>1.65 × 10<sup>11</sup></td><td><sup>62</sup>Fe → <sup>62</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.044</td></tr><tr><td>8.94 × 10<sup>29</sup></td><td>2.70 × 10<sup>11</sup></td><td><sup>62</sup>Cr → <sup>62</sup>Ti − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.086</td></tr><tr><td>1.37 × 10<sup>30</sup></td><td>3.64 × 10<sup>11</sup></td><td><sup>62</sup>Ti → <sup>62</sup>Ca − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.041</td></tr><tr><td>1.31 × 10<sup>29</sup></td><td>9.50 × 10<sup>10</sup></td><td><sup>88</sup>Sr → <sup>88</sup>Kr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.027</td></tr><tr><td>2.36 × 10<sup>29</sup></td><td>1.26 × 10<sup>11</sup></td><td><sup>88</sup>Sr → <sup>88</sup>Kr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.078</td></tr><tr><td>3.16 × 10<sup>29</sup></td><td>1.54 × 10<sup>11</sup></td><td><sup>88</sup>Kr → <sup>88</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.025</td></tr><tr><td>4.43 × 10<sup>29</sup></td><td>1.81 × 10<sup>11</sup></td><td><sup>88</sup>Kr → <sup>88</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.068</td></tr><tr><td>6.31 × 10<sup>29</sup></td><td>2.28 × 10<sup>11</sup></td><td><sup>88</sup>Se → <sup>88</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.031</td></tr><tr><td>1.00 × 10<sup>30</sup></td><td>3.05 × 10<sup>11</sup></td><td><sup>88</sup>Ge → <sup>88</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.031</td></tr><tr><td>2.38 × 10<sup>29</sup></td><td>1.31 × 10<sup>11</sup></td><td><sup>86</sup>Kr → <sup>86</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.029</td></tr><tr><td>3.87 × 10<sup>29</sup></td><td>1.66 × 10<sup>11</sup></td><td><sup>86</sup>Kr → <sup>86</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.086</td></tr><tr><td>4.89 × 10<sup>29</sup></td><td>1.97 × 10<sup>11</sup></td><td><sup>86</sup>Se → <sup>86</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.023</td></tr><tr><td>8.93 × 10<sup>29</sup></td><td>2.81 × 10<sup>11</sup></td><td><sup>86</sup>Ge → <sup>86</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.037</td></tr><tr><td>1.25 × 10<sup>30</sup></td><td>3.53 × 10<sup>11</sup></td><td><sup>86</sup>Zn → <sup>86</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.027</td></tr><tr><td>3.85 × 10<sup>29</sup></td><td>1.71 × 10<sup>11</sup></td><td><sup>84</sup>Se → <sup>84</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.029</td></tr><tr><td>7.13 × 10<sup>29</sup></td><td>2.46 × 10<sup>11</sup></td><td><sup>84</sup>Ge → <sup>84</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.029</td></tr><tr><td>1.09 × 10<sup>30</sup></td><td>3.24 × 10<sup>11</sup></td><td><sup>84</sup>Zn → <sup>84</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (⋆)</td><td>0.029</td></tr><tr><td>5.59 × 10<sup>29</sup></td><td>2.13 × 10<sup>11</sup></td><td><sup>82</sup>Ge → <sup>82</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.023</td></tr><tr><td>8.26 × 10<sup>29</sup></td><td>2.58 × 10<sup>11</sup></td><td><sup>82</sup>Ge → <sup>82</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.096</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>9.09 × 10<sup>29</sup></td><td>2.88 × 10<sup>11</sup></td><td><sup>82</sup>Zn → <sup>82</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.019</td></tr><tr><td>1.61 × 10<sup>30</sup></td><td>4.10 × 10<sup>11</sup></td><td><sup>82</sup>Ni → <sup>82</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.026</td></tr><tr><td>7.44 × 10<sup>29</sup></td><td>2.54 × 10<sup>11</sup></td><td><sup>132</sup>Sn → <sup>132</sup>Cd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.022</td></tr><tr><td>8.61 × 10<sup>29</sup></td><td>2.73 × 10<sup>11</sup></td><td><sup>132</sup>Sn → <sup>132</sup>Cd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.040</td></tr><tr><td>9.01 × 10<sup>29</sup></td><td>2.90 × 10<sup>11</sup></td><td><sup>132</sup>Cd → <sup>132</sup>Pd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.017</td></tr><tr><td>1.25 × 10<sup>30</sup></td><td>3.56 × 10<sup>11</sup></td><td><sup>132</sup>Pd → <sup>132</sup>Ru − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.022</td></tr><tr><td>7.00 × 10<sup>29</sup></td><td>2.47 × 10<sup>11</sup></td><td><sup>80</sup>Zn → <sup>80</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.005</td></tr><tr><td>1.44 × 10<sup>30</sup></td><td>3.79 × 10<sup>11</sup></td><td><sup>80</sup>Ni → <sup>80</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.026</td></tr><tr><td>1.04 × 10<sup>30</sup></td><td>3.15 × 10<sup>11</sup></td><td><sup>128</sup>Pd → <sup>128</sup>Ru − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.018</td></tr><tr><td>1.16 × 10<sup>30</sup></td><td>3.42 × 10<sup>11</sup></td><td><sup>126</sup>Ru → <sup>126</sup>Mo − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.015</td></tr><tr><td>1.48 × 10<sup>30</sup></td><td>3.97 × 10<sup>11</sup></td><td><sup>124</sup>Mo → <sup>124</sup>Zr − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.024</td></tr><tr><td>1.53 × 10<sup>30</sup></td><td>4.17 × 10<sup>11</sup></td><td><sup>122</sup>Zr → <sup>122</sup>Sr − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.007</td></tr></tbody></table>

<table><caption>**Table 2.** Same as Table 1, but for $B_{\star} = 3000$.</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>1.04 × 10<sup>29</sup></td><td>9.74 × 10<sup>10</sup></td><td><sup>56</sup>Fe → <sup>56</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.037</td></tr><tr><td>1.09 × 10<sup>29</sup></td><td>9.96 × 10<sup>10</sup></td><td><sup>56</sup>Fe → <sup>56</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.041</td></tr><tr><td>4.72 × 10<sup>29</sup></td><td>2.19 × 10<sup>11</sup></td><td><sup>56</sup>Cr → <sup>56</sup>Ti − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.042</td></tr><tr><td>1.08 × 10<sup>30</sup></td><td>3.58 × 10<sup>11</sup></td><td><sup>56</sup>Ti → <sup>56</sup>Ca − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.023</td></tr><tr><td>1.92 × 10<sup>29</sup></td><td>1.34 × 10<sup>11</sup></td><td><sup>62</sup>Ni → <sup>62</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.045</td></tr><tr><td>2.23 × 10<sup>29</sup></td><td>1.45 × 10<sup>11</sup></td><td><sup>62</sup>Ni → <sup>62</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.062</td></tr><tr><td>5.98 × 10<sup>29</sup></td><td>2.52 × 10<sup>11</sup></td><td><sup>62</sup>Fe → <sup>62</sup>Cr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.044</td></tr><tr><td>1.37 × 10<sup>30</sup></td><td>4.09 × 10<sup>11</sup></td><td><sup>62</sup>Cr → <sup>62</sup>Ti − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.086</td></tr><tr><td>2.09 × 10<sup>30</sup></td><td>5.50 × 10<sup>11</sup></td><td><sup>62</sup>Ti → <sup>62</sup>Ca − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.041</td></tr><tr><td>2.06 × 10<sup>29</sup></td><td>1.47 × 10<sup>11</sup></td><td><sup>88</sup>Sr → <sup>88</sup>Kr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.027</td></tr><tr><td>3.68 × 10<sup>29</sup></td><td>1.94 × 10<sup>11</sup></td><td><sup>88</sup>Sr → <sup>88</sup>Kr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.078</td></tr><tr><td>4.90 × 10<sup>29</sup></td><td>2.35 × 10<sup>11</sup></td><td><sup>88</sup>Kr → <sup>88</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.025</td></tr><tr><td>6.86 × 10<sup>29</sup></td><td>2.77 × 10<sup>11</sup></td><td><sup>88</sup>Kr → <sup>88</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.068</td></tr><tr><td>9.72 × 10<sup>29</sup></td><td>3.48 × 10<sup>11</sup></td><td><sup>88</sup>Se → <sup>88</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.031</td></tr><tr><td>1.54 × 10<sup>30</sup></td><td>4.63 × 10<sup>11</sup></td><td><sup>88</sup>Ge → <sup>88</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.032</td></tr><tr><td>3.72 × 10<sup>29</sup></td><td>2.01 × 10<sup>11</sup></td><td><sup>86</sup>Kr → <sup>86</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.029</td></tr><tr><td>6.00 × 10<sup>29</sup></td><td>2.54 × 10<sup>11</sup></td><td><sup>86</sup>Kr → <sup>86</sup>Se − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.086</td></tr><tr><td>7.55 × 10<sup>29</sup></td><td>3.00 × 10<sup>11</sup></td><td><sup>86</sup>Se → <sup>86</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.023</td></tr><tr><td>1.37 × 10<sup>30</sup></td><td>4.28 × 10<sup>11</sup></td><td><sup>86</sup>Ge → <sup>86</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.037</td></tr><tr><td>1.91 × 10<sup>30</sup></td><td>5.36 × 10<sup>11</sup></td><td><sup>86</sup>Zn → <sup>86</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.027</td></tr><tr><td>5.96 × 10<sup>29</sup></td><td>2.61 × 10<sup>11</sup></td><td><sup>84</sup>Se → <sup>84</sup>Ge − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.029</td></tr><tr><td>1.10 × 10<sup>30</sup></td><td>3.74 × 10<sup>11</sup></td><td><sup>84</sup>Ge → <sup>84</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.029</td></tr><tr><td>1.67 × 10<sup>30</sup></td><td>4.91 × 10<sup>11</sup></td><td><sup>84</sup>Zn → <sup>84</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.029</td></tr><tr><td>8.61 × 10<sup>29</sup></td><td>3.24 × 10<sup>11</sup></td><td><sup>82</sup>Ge → <sup>82</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.023</td></tr><tr><td>1.27 × 10<sup>30</sup></td><td>3.92 × 10<sup>11</sup></td><td><sup>82</sup>Ge → <sup>82</sup>Zn − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.096</td></tr></tbody></table>
<table><caption>**Table 2.** *Cont.*</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>1.39 × 10<sup>30</sup></td><td>4.38 × 10<sup>11</sup></td><td><sup>82</sup>Zn → <sup>82</sup>Ni − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.019</td></tr><tr><td>2.46 × 10<sup>30</sup></td><td>6.21 × 10<sup>11</sup></td><td><sup>82</sup>Ni → <sup>82</sup>Fe − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.026</td></tr><tr><td>1.15 × 10<sup>30</sup></td><td>3.88 × 10<sup>11</sup></td><td><sup>132</sup>Sn → <sup>132</sup>Cd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.022</td></tr><tr><td>1.33 × 10<sup>30</sup></td><td>4.17 × 10<sup>11</sup></td><td><sup>132</sup>Sn → <sup>132</sup>Cd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.040</td></tr><tr><td>1.39 × 10<sup>30</sup></td><td>4.44 × 10<sup>11</sup></td><td><sup>132</sup>Cd → <sup>132</sup>Pd − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.017</td></tr><tr><td>1.93 × 10<sup>30</sup></td><td>5.42 × 10<sup>11</sup></td><td><sup>132</sup>Pd → <sup>132</sup>Ru − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.023</td></tr><tr><td>1.60 × 10<sup>30</sup></td><td>4.80 × 10<sup>11</sup></td><td><sup>128</sup>Pd → <sup>128</sup>Ru − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.018</td></tr><tr><td>1.79 × 10<sup>30</sup></td><td>5.21 × 10<sup>11</sup></td><td><sup>126</sup>Ru → <sup>126</sup>Mo − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.015</td></tr><tr><td>2.27 × 10<sup>30</sup></td><td>6.04 × 10<sup>11</sup></td><td><sup>124</sup>Mo → <sup>124</sup>Zr − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.024</td></tr><tr><td>2.35 × 10<sup>30</sup></td><td>6.34 × 10<sup>11</sup></td><td><sup>122</sup>Zr → <sup>122</sup>Sr − 2e<sup>−</sup> + 2v<sub>e</sub> (★)</td><td>0.007</td></tr></tbody></table>

<table><caption>**Table 3.** Maximum possible heat released *q* per nucleon (in MeV) from electron captures by carbon and oxygen in the outer crust of a magnetar with $B_{\star} = 2000$. The pressure $P_{\beta}$ (in dyn cm$^{-2}$) and density $\rho_{\beta}$ (in g cm$^{-3}$) at which electron capture occurs are given considering ground-state-to-ground-state transitions.</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>5.59 × 10<sup>29</sup></td><td>1.64 × 10<sup>11</sup></td><td><sup>12</sup>C → <sup>12</sup>Be − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.143</td></tr><tr><td>3.55 × 10<sup>29</sup></td><td>1.31 × 10<sup>11</sup></td><td><sup>16</sup>O → <sup>16</sup>C − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.153</td></tr></tbody></table>

<table><caption>**Table 4.** Same as in Table 3, but for $B_{\star} = 3000$.</caption><thead><tr><th>P<sub>β</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>β</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q (MeV)</th></tr></thead><tbody><tr><td>8.46 × 10<sup>29</sup></td><td>2.47 × 10<sup>11</sup></td><td><sup>12</sup>C → <sup>12</sup>Be − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.143</td></tr><tr><td>5.38 × 10<sup>29</sup></td><td>1.98 × 10<sup>11</sup></td><td><sup>16</sup>O → <sup>16</sup>C − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>0.154</td></tr></tbody></table>

The heating resulting from pycnonuclear fusions of carbon and oxygen is summarized in Tables 5 and 6 for $B_{\star} = 2000$ and $B_{\star} = 3000$, respectively. Contrary to our previous study of accreting neutron stars [16], we did not consider here fusions of oxygen occurring at pressure $P_{\beta}^*$ because we realized that $\mu_e^\beta > \mu_e^{\text{drip}}$: oxygen is unstable against electron-capture-induced neutron emission. The heat from pycnonuclear fusion reactions actually consists of two contributions: the first from the fusion itself and the second from the subsequent electron captures. The two processes were found to be roughly equally exothermic. If they occurred, the reactions would contribute as much heat per nucleon as all electron captures by heavier nuclei combined. However, the total amount of heat actually deposited depends on the unknown abundance of light elements. Let us stress that the amount of heat estimated here is a conservative upper limit. In reality, $q_{\text{pyc}}$ could be much smaller, especially if fusions occur at densities substantially lower than $n_{\beta}$. This may actually be the case since the fusion rates are expected to be thermally enhanced in magnetars due to their higher temperatures than in ordinary neutron stars [43].
<table><caption>**Table 5.** Maximum possible heat released $q_{pyc}$ per nucleon (in MeV) from pycnonuclear fusion followed by electron captures in the outer crust of a magnetar with $B_{\star} = 2000$. The numbers in parentheses indicate the contributions from the sole fusion. The pressure $P_{pyc}$ (in dyn cm$^{-2}$) and density $\rho_{pyc}$ (in g cm$^{-3}$) at which fusion occurs were fixed by the onset of electron captures considering ground-state-to-ground-state transitions.</caption><thead><tr><th>P<sub>pyc</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>pyc</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q<sub>pyc</sub> (MeV)</th></tr></thead><tbody><tr><td>5.59 × 10<sup>29</sup></td><td>1.64 × 10<sup>11</sup></td><td><sup>12</sup>C + <sup>12</sup>C → <sup>24</sup>Ne - 2e<sup>-</sup> + 2v<sub>e</sub></td><td>1.42 (0.65)</td></tr><tr><td>3.55 × 10<sup>29</sup></td><td>1.31 × 10<sup>11</sup></td><td><sup>16</sup>O + <sup>16</sup>O → <sup>32</sup>Si - 2e<sup>-</sup> + 2v<sub>e</sub></td><td>1.17 (0.60)</td></tr></tbody></table>

<table><caption>**Table 6.** Same as Table 5 for $B_{\star} = 3000$.</caption><thead><tr><th>P<sub>pyc</sub> (dyn cm<sup>-2</sup>)</th><th>ρ<sub>pyc</sub> (g cm<sup>-3</sup>)</th><th>Reactions</th><th>q<sub>pyc</sub> (MeV)</th></tr></thead><tbody><tr><td>8.46 × 10<sup>29</sup></td><td>2.47 × 10<sup>11</sup></td><td><sup>12</sup>C + <sup>12</sup>C → <sup>24</sup>Ne - 2e<sup>-</sup> + 2v<sub>e</sub></td><td>1.42 (0.67)</td></tr><tr><td>5.38 × 10<sup>29</sup></td><td>1.98 × 10<sup>11</sup></td><td><sup>16</sup>O + <sup>16</sup>O → <sup>32</sup>Si - 2e<sup>-</sup> + 2v<sub>e</sub></td><td>1.18 (0.61)</td></tr></tbody></table>

To assess the validity of the analytical treatment, we numerically solved the threshold conditions $g(A, Z, n_e) = g(A, Z - 1, n_{e1})$ and $P(Z, n_e) = P(Z - 1, n_{e1}) = P(Z - 2, n_{e2})$ without any approximation to determine the exact values for the transition pressure $P_{\beta}(A, Z, B_{\star})$ and density $n_{\beta}(A, Z, B_{\star})$. The heat released was then calculated from Equation (28). The relative deviations for the transition pressure, densities, and amount of heat are typically of order $10^{-3}\%$. In view of the ultrarelativistic approximation (18) we made for the electrostatic correction in Equation (13), the largest errors were found for the shallowest transitions and for the strongest magnetic field. In particular, considering the ground-state-to-ground-state transition from $^{56}$Fe to $^{56}$Cr for $B_{\star} = 3000$, the results obtained from Equations (23), (24), and (31) differed from the exact values by $-7.6 \times 10^{-2}\%$, $-3.6 \times 10^{-2}\%$, and $5.0 \times 10^{-3}\%$, respectively.

# 4. Astrophysical Implications

## 4.1. Equilibrium of Self-Gravitating Magnetized Stars

For the sake of the argument, we restricted ourselves to Newtonian gravity, as in [9]. The magneto-hydrodynamic equilibrium equation for an ideally conducting self-gravitating fluid star reads [50] (see [51] for the general relativistic equations for stationary axisymmetric magnetized stars and [52] for the discussion of the Newtonian limit):

$$n \nabla \mu = -\rho \nabla \Phi + \frac{1}{c} j_{\text{free}} \times \mathbf{B}, \qquad (40)$$

where $\Phi$ is the gravitational potential obeying Poisson's equation ($G$ being the gravitational constant):

$$\nabla^2 \Phi = 4\pi G \rho, \qquad (41)$$

and the free electric charge current $j_{\text{free}}$ is given by Maxwell's equation:

$$\nabla \times H = \frac{4\pi}{c} j_{\text{free}}. \qquad (42)$$

The magnetization of a relativistic electron gas in a strongly quantizing magnetic field is negligibly small [23]; therefore, we can safely replace $H \approx B$. Expanding the Lorentz force term in the right-hand side of Equation (40) using Equation (42) yields:

$$n \nabla \mu + \nabla \left( \frac{B^2}{8\pi} \right) = -\rho \nabla \Phi + \frac{1}{4\pi} (\mathbf{B} \cdot \nabla) \mathbf{B}. \qquad (43)$$

Besides, the toroidal component of the magnetic field inside a neutron star is expected to be much stronger than the poloidal component [53] (the strength of which can be estimated
mated from the spin down, as measured by timing analyses). These theoretical expectations were corroborated by astrophysical observations [32,54–58]. Assuming for simplicity that the magnetic field is axially symmetric and purely toroidal, i.e., **B** = B(ω, z)1<sub>φ</sub> in cylindrical coordinates [59] (ω is the radial distance, z the height, φ the azimuthal angle, and 1<sub>φ</sub> the associated unit vector), the last term in Equation (43) reduces to −B<sup>2</sup>/(4πω)1<sub>φ</sub>. Recalling that the thickness δR and the (gravitational or baryonic) mass δM of the outer crust are very small compared to the neutron star radius R and mass M, respectively, typically, δR ~ 10<sup>−2</sup>R and δM ~ 10<sup>−5</sup>M (see, e.g., [60]), we adopted the plane parallel approximation. Since the toroidal component of the magnetic field is expected to be confined in the crust of thickness ΔR ~ 0.1R (see, e.g., [61]), the last term in the right-hand side of Equation (43) of order B<sup>2</sup>/(4πR) is much smaller than the second term in the left-hand side of order B<sup>2</sup>/(8πΔR). We thus dropped the last term in Equation (43) as in [9]. Introducing the local gravitational field **g** = −∇Φ, which is essentially uniform in the outer crust and given by **g** ~ GM/R<sup>2</sup>, Equation (43) finally reduces to:

$$nd\mu + \frac{1}{8\pi}d(B^2) = gd\Sigma. \qquad (44)$$

where Σ is the column density as measured from the surface and defined by dΣ = −ρdr (r being the radial distance in spherical coordinates).

Let us focus on a crustal layer located at a given column density (dΣ = 0). A decrease of the magnetic pressure d(B<sup>2</sup>) < 0 must thus be compensated by an increase of the baryon chemical potential dμ > 0 for the layer to remain in mechanical equilibrium. Due to this excess energy, nuclei may become unstable against electron captures and pycnonuclear fusion reactions, as discussed in Section 2.

## 4.2. Heating Induced by Magnetic Field Decay

The change of magnetic field strength required to trigger electron captures is very small and can be estimated as follows. Ignoring magnetization effects and recalling that the electron gas is highly degenerate, the Gibbs–Duhem relation reduces to ndμ ~ dP (see, e.g., Appendix A of [24]). The stellar equilibrium condition (44) thus entails that the total pressure (matter plus magnetic), given by:

$$P_{\text{tot}} = P + \frac{B^2}{8\pi}, \qquad (45)$$

must remain constant. Initially, the nuclei (A, Z) are present in the crust in a layer delimited by the matter pressures P<sub>min</sub>(A, Z, B<sub>+</sub> + δB<sub>+</sub>) and P<sub>max</sub>(A, Z, B<sub>+</sub> + δB<sub>+</sub>), as determined by the minimization of the baryon chemical potential. For all the nuclei (A, Z) to capture electrons, the magnetic field strength must therefore decay by an amount δB<sub>+</sub> such that the shallowest layer containing this element can be compressed up to the threshold pressure P<sub>β</sub>(A, Z, B<sub>+</sub>), thus maintaining the same total pressure. Using Equation (45) and expanding to first order in δB<sub>+</sub>, we find:

$$\frac{\delta B_{\star}}{B_{\star}} \approx \frac{4\pi}{B^2} \left[ P_{\beta}(A, Z, B_{\star}) - P_{\min}(A, Z, B_{\star}) \right], \qquad (46)$$

where the pressures are evaluated here for the same magnetic field strength B<sub>+</sub>. Results are summarized in Tables 7 and 8 for B<sub>+</sub> = 2000 and B<sub>+</sub> = 3000, respectively, using the results of [21] for the initial pressures P<sub>min</sub> (note that the initial pressure of <sup>56</sup>Fe is simply P<sub>min</sub> = 0 since this element is predicted to be present at the stellar surface). Results for carbon and oxygen supposedly initially accumulated onto the surface are presented in Tables 9 and 10 (in this case, we set P<sub>min</sub> = 0). Note that the required magnetic field changes for the pycnonuclear fusions of these elements are the same since we assumed that these reactions occur at the pressure P<sub>β</sub>. The changes of magnetic field required to trigger the various chains of reactions depend on the layer, but are typically of order δB<sub>+</sub>/B<sub>+</sub> ~ 10<sup>−3</sup>–10<sup>−4</sup>
for the magnetic field considered and are largest in the shallow layers where the magnetic field is the most strongly quantizing. Inspecting the numerical results more closely, $\delta B_*/B_*$ is found to decrease systematically with increasing $B_*$. The trend becomes more significant if one considers much weaker magnetic fields. Since we are limited by our assumption of a strongly quantizing magnetic field, let us focus on the electron captures by $^{56}$Fe. The relative change $\delta B_*/B_*$ of the magnetic field thus drops from $2.09 \times 10^{-2}$ to $2.21 \times 10^{-3}$ when the magnetic field is varied from $B_* = 100$ to $B_* = 1000$. In the ultrarelativistic regime ($\gamma_e \gg 1$), both $P_{\min}(A, Z, B_*)$ and $P_\beta(A, Z, B_*)$ increase linearly with $B_*$ so that $\delta B_*/B_* \propto 1/B_*$.

The time $\tau \sim \tau_B \delta B_*/B_*$ after which the magnetic field has changed by an amount $\delta B_*$ is thus approximately inversely proportional to the magnetic field strength ($\tau_B \sim \text{Myr}$ being the characteristic time scale of magnetic field dissipation). Crustal heating due to electron captures and pycnonuclear fusion reactions is thus likely to be relevant for young and middle-age magnetars. This conclusion is consistent with cooling simulations [7,8]. The relative change of the magnetic field since the birth of magnetars can thus potentially be inferred from their age, as estimated from the kinematics of the expanding supernova remnant. With ages typically of the order of a few kyr according to the Magnetar Outburst Online Catalog$^5$, we found that the magnetic fields have decayed by $\delta B_*/B_* \sim \tau/\tau_B \sim 10^{-3}$, which is comparable to the change of magnetic field required by nuclear reactions.

The total energy deposited in the outer crust due to the magnetic field decay is given by $\Sigma_A q_{A,tot}(\Delta M_A/m_u)$ where $q_{A,tot}$ is the total energy per baryon released for a given initial nucleus defined by A in Table 7 and $\Delta M_A$ the mass of this layer (we ignored here the difference between the gravitational and baryonic masses). Because $q_{A,tot}$ is of the same order $\sim 0.1 \text{ MeV}$ for all layers, the heat power can be estimated as $W^\infty \sim q_{A,tot}(\Delta M/m_u)/\tau$, where $\Delta M = \sum_A \Delta M_A$ is the mass of the outer crust. With $\Delta M \sim 10^{-5}-10^{-4} M_\odot$ [60] with $M_\odot$ the mass of the Sun, and $\tau \sim 1 \text{ kyr}$, we found $W^\infty \sim 10^{35}-10^{36} \text{ erg/s}$, which is comparable to the heat power obtained empirically by fitting cooling curves to observational data [7,8]. For comparison, the power resulting from the decay of the magnetic field amounts roughly to $BR^2\delta R\delta B/\tau \sim 10^{42} \text{ erg/s}$. This energy is expected to be the main source of the magnetar activity such as X-ray bursts and giant flares. However, some fraction could also contribute to the thermal emission.

These simple estimates suggest that nuclear reactions may be a viable source of heating in mature magnetars. However, numerical simulations of the full magneto-thermal evolution are needed to make more reliable predictions.

<table><caption>**Table 7.** Relative decay of the magnetic field required to entirely replace the given element through the indicated chain of electron captures by the nuclei listed in Table 1 in the outer crust of a magnetar for $B_* = 2000$, from the initial pressure given in [21] to the pressure at which the last capture occurs.</caption><thead><tr><th>Reactions</th><th>δB<sub>*</sub>/B<sub>*</sub></th></tr></thead><tbody><tr><td><sup>56</sup>Fe → <sup>56</sup>Ca - 6e<sup>-</sup> + 6v<sub>e</sub></td><td>1.13 × 10<sup>-3</sup></td></tr><tr><td><sup>62</sup>Ni → <sup>62</sup>Ca - 8e<sup>-</sup> + 8v<sub>e</sub></td><td>2.21 × 10<sup>-3</sup></td></tr><tr><td><sup>88</sup>Sr → <sup>88</sup>Zn - 8e<sup>-</sup> + 8v<sub>e</sub></td><td>1.58 × 10<sup>-3</sup></td></tr><tr><td><sup>86</sup>Kr → <sup>86</sup>Ni - 8e<sup>-</sup> + 8v<sub>e</sub></td><td>1.94 × 10<sup>-3</sup></td></tr><tr><td><sup>84</sup>Se → <sup>84</sup>Ni - 6e<sup>-</sup> + 6v<sub>e</sub></td><td>1.58 × 10<sup>-3</sup></td></tr><tr><td><sup>82</sup>Ge → <sup>82</sup>Fe - 6e<sup>-</sup> + 6v<sub>e</sub></td><td>2.23 × 10<sup>-3</sup></td></tr><tr><td><sup>132</sup>Sn → <sup>132</sup>Ru - 6e<sup>-</sup> + 6v<sub>e</sub></td><td>1.39 × 10<sup>-3</sup></td></tr><tr><td><sup>80</sup>Zn → <sup>80</sup>Fe - 4e<sup>-</sup> + 4v<sub>e</sub></td><td>1.66 × 10<sup>-3</sup></td></tr><tr><td><sup>128</sup>Pd → <sup>128</sup>Ru - 2e<sup>-</sup> + 2v<sub>e</sub></td><td>7.91 × 10<sup>-4</sup></td></tr></tbody></table>
<table><caption>**Table 7.** *Cont.*</caption><thead><tr><th>Reactions</th><th>δB<sub>*</sub>/B<sub>*</sub></th></tr></thead><tbody><tr><td><sup>126</sup>Ru → <sup>126</sup>Mo − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>5.22 × 10<sup>−4</sup></td></tr><tr><td><sup>124</sup>Mo → <sup>124</sup>Zr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>8.81 × 10<sup>−4</sup></td></tr><tr><td><sup>122</sup>Zr → <sup>122</sup>Sr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>4.81 × 10<sup>−4</sup></td></tr></tbody></table>

<table><caption>**Table 8.** Same as Table 7, but for $B_{\star} = 3000$.</caption><thead><tr><th>Reactions</th><th>δB<sub>*</sub>/B<sub>*</sub></th></tr></thead><tbody><tr><td><sup>56</sup>Fe → <sup>56</sup>Ca − 6e<sup>−</sup> + 6v<sub>e</sub></td><td>7.70 × 10<sup>−4</sup></td></tr><tr><td><sup>62</sup>Ni → <sup>62</sup>Ca − 8e<sup>−</sup> + 8v<sub>e</sub></td><td>1.49 × 10<sup>−3</sup></td></tr><tr><td><sup>88</sup>Sr → <sup>88</sup>Zn − 8e<sup>−</sup> + 8v<sub>e</sub></td><td>1.08 × 10<sup>−3</sup></td></tr><tr><td><sup>86</sup>Kr → <sup>86</sup>Ni − 8e<sup>−</sup> + 8v<sub>e</sub></td><td>1.32 × 10<sup>−3</sup></td></tr><tr><td><sup>84</sup>Se → <sup>84</sup>Ni − 6e<sup>−</sup> + 6v<sub>e</sub></td><td>1.06 × 10<sup>−3</sup></td></tr><tr><td><sup>82</sup>Ge → <sup>82</sup>Fe − 6e<sup>−</sup> + 6v<sub>e</sub></td><td>1.49 × 10<sup>−3</sup></td></tr><tr><td><sup>132</sup>Sn → <sup>132</sup>Ru − 6e<sup>−</sup> + 6v<sub>e</sub></td><td>1.02 × 10<sup>−3</sup></td></tr><tr><td><sup>128</sup>Pd → <sup>128</sup>Ru − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>5.40 × 10<sup>−4</sup></td></tr><tr><td><sup>126</sup>Ru → <sup>126</sup>Mo − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>3.41 × 10<sup>−4</sup></td></tr><tr><td><sup>124</sup>Mo → <sup>124</sup>Zr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>5.94 × 10<sup>−4</sup></td></tr><tr><td><sup>122</sup>Zr → <sup>122</sup>Sr − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>3.10 × 10<sup>−4</sup></td></tr></tbody></table>

<table><caption>**Table 9.** Relative decay of the magnetic field required to trigger the chains of electron captures by light nuclei listed in Table 3 in the outer crust of a magnetar for $B_{\star} = 2000$, from zero pressure (at the stellar surface) to the pressure at which the last capture occurs.</caption><thead><tr><th>Reactions</th><th>δB<sub>*</sub>/B<sub>*</sub></th></tr></thead><tbody><tr><td><sup>12</sup>C → <sup>12</sup>Be − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>9.02 × 10<sup>−4</sup></td></tr><tr><td><sup>16</sup>O → <sup>16</sup>C − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>5.71 × 10<sup>−4</sup></td></tr></tbody></table>

<table><caption>**Table 10.** Same as Table 9 for $B_{\star} = 3000$.</caption><thead><tr><th>Reactions</th><th>δB<sub>*</sub>/B<sub>*</sub></th></tr></thead><tbody><tr><td><sup>12</sup>C → <sup>12</sup>Be − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>6.05 × 10<sup>−4</sup></td></tr><tr><td><sup>16</sup>O → <sup>16</sup>C − 2e<sup>−</sup> + 2v<sub>e</sub></td><td>3.86 × 10<sup>−4</sup></td></tr></tbody></table>

### 4.3. Heating Induced by Spin-Down

Observed isolated magnetars rotate at frequencies $f \sim 0.1–0.5$ Hz and exhibit a very high spin-down rate. They are believed to be born with initial frequencies $f_i \sim 1$ kHz. The time to spin-down to a frequency $f \ll f_i$ can be estimated as:

$$t \sim 6 \times 10^{-3} \left( \frac{B_p}{10^{14} G} \right)^{-2} \left( \frac{f}{1 \text{ kHz}} \right)^{-2} \text{yr}, \quad (47)$$

where $B_p$ is the value of the poloidal magnetic field at the pole and the time dependence of $B_p$ is neglected. Initially, the crust is rotationally flattened, with the equatorial radius larger than the polar one. Using Equation (47), we found that a newborn magnetar spins down to $f = 10$ Hz in less than $t \sim 100 (B_p/10^{14} G)^{-2}$ yr. According to [15], the decrease of centrifugal force could induce non-equilibrium processes and heating. Some results of [15] were questioned by the authors of [62], who studied the spin-down heating in millisecond pulsars with a fully accreted crust. Notice however that the authors of [15] assumed that initially the crust is composed of catalyzed matter, which is appropriate for a hot newly born magnetar. Anyway, spin-down heating could be efficient only during
the first few days of magnetars' lives, because at later times, spin-down compression is no longer significant. Therefore, the spin-down heating does not contribute significantly to the observed thermal luminosity of known magnetars.

**5. Conclusions**

We investigated electron captures and pycnonuclear fusion reactions in the outer crust of a magnetar induced by the compression of matter accompanying the decay of the magnetic field and the spin-down of the star. Taking into account Landau-Rabi quantization of electron motion and focusing on the strongly quantizing regime, we derived very accurate analytical formulas for the maximum amount of heat that can possibly be released by each reaction and their location. Making use of essentially all available experimental data supplemented with the predictions from the HFB-24 atomic mass model, we found that the energy release was ~0.02-0.1 MeV per nucleon in two subsequent electron captures, similarly to accreting systems, although the composition and the thermodynamic conditions are very different. For the initial composition of the magnetar crust, we considered the sequence of equilibrium nuclei previously calculated in [21] using the same atomic mass model HFB-24. We also studied the possibility that light elements such as carbon and oxygen might have been accreted onto the surface of the star from the fallback of supernovas debris, from a disk, or from the interstellar medium. The pycnonuclear fusions of these elements could potentially release as much heat per nucleon as the electron captures by all the other nuclei.

The maximum amount of heat released by each individual reaction is found to be essentially independent of the magnetic field and is mainly determined by the relevant Q-values. On the contrary, the pressure and the density at which heat is deposited both increase almost linearly with the magnetic field strength. For internal magnetic fields of order $10^{16}-10^{17}$ G, heat sources are found in deeper layers than in accreting neutron stars. Moreover, they are not uniformly distributed, but are concentrated at densities of order $10^{10}-10^{11}$ g cm$^{-3}$ (pressures $10^{29}-10^{30}$ dyn cm$^{-2}$). Quite remarkably, a similar range of densities is supported by the adjustment of cooling simulations to the observed thermal luminosity of magnetars [7,8].

We also showed that the relative change of magnetic field required to trigger the various reactions is approximately inversely proportional to the magnetic field strength and is typically of order $\delta B/B \sim 10^{-3}-10^{-4}$ for $B \sim 10^{16}-10^{17}$ G. Such variations are comparable to those expected from the decay of the magnetic field since the birth of currently known magnetars. Moreover, the heat power $W \sim 10^{35}-10^{36}$ erg/s is found to be consistent with values inferred empirically from cooling simulations varying the composition of the envelope [7,8]. Electron captures and pycnonuclear fusion reactions induced by the decay of the magnetic field may thus potentially explain the origin of internal heating in magnetars. Although nuclear processes could also potentially be triggered by the spin-down of the star, we showed that this mechanism only operates during the early life of newborn magnetars.

Unlike other mechanisms involving crust quakes, the heating induced by nuclear reactions is essentially independent of the detailed crustal structure and will still remain viable if some regions are actually liquid. Indeed, the amount of heat released by each individual reaction mainly depends on nuclear masses, while the density and the pressure at which they occur is governed by the electron gas. Together with the unified equations of state calculated for various magnetic field strengths in [21], the present results calculated using the same nuclear model provide consistent and realistic microscopic inputs for cooling simulations of magnetars.

Our analysis can be easily extended to lower magnetic fields (involving summations over several Landau-Rabi levels) and finite temperatures. However, in this more general situation, the calculations are not amenable to analytical solutions. In our calculations, we ignored the effects of the magnetic field on nuclei. Although the change of nuclear masses is expected to be small for $B \lesssim 10^{17}$ G [63,64], its impact on the Q-values, hence also on the heat sources and their location, may be more significant and deserves further investigation.
We also implicitly assumed that the compression of matter occurs continuously. However, if the layers of interest are in a solid phase, the loss of magnetic pressure may be compensated by the build-up of elastic stresses until the crust fails; in this case, the compression will occur suddenly during crust quakes so that the electron capture (6) will proceed off-equilibrium with the release of heat. The maximum possible total amount of heat deposited in the outer crust could thus be potentially even higher than our present estimate. This alternative scenario, which is more likely to occur in the deepest layers of the outer crust where the melting temperature is highest, requires the detailed knowledge of the cooling and magneto-elastohydrodynamic evolution of the liquid ocean and solid layers beneath.

**Author Contributions:** Conceptualization, N.C.; methodology, N.C.; software, A.F.F. and N.C.; validation, N.C., A.F.F., L.S., J.-L.Z., and P.H.; formal analysis, N.C. and A.F.F.; investigation, N.C.; writing—original draft preparation, N.C.; writing—review and editing, N.C., A.F.F., L.S., J.-L.Z., and P.H.; visualization, A.F.F. and N.C.; supervision, N.C.; project administration, N.C. All authors read and agreed to the published version of the manuscript.

**Funding:** The work of N.C. was funded by Fonds de la Recherche Scientifique-FNRS (Belgium) under Grant Number IISN 4.4502.19. L.S., P.H., and J.-L.Z. acknowledge the financial support from the National Science Centre (Poland) Grant Number 2018/29/B/ST9/02013. This work was also partially supported by the European Cooperation in Science and Technology Action CA16214 and the CNRS International Research Project (IRP) "Origine des éléments lourds dans l'univers: Astres Compacts et Nucléosynthèse (ACNu)".

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data analyzed in this paper can be found in the McGill Online Magnetar Catalog (http://www.physics.mcgill.ca/~pulsar/magnetar/main.html, accessed on 7 May 2021, see [3]), the Magnetar Outburst Online Catalog (http://magnetars.ice.csic.es, accessed on 7 May 2021), the 2016 Atomic Mass Evaluation (see [44]), the BRUSLIB database (http://www.astro.ulb.ac.be/bruslib/, accessed on 7 May 2021, see [46]), and the Nuclear Data section of the International Atomic Energy Agency website (https://www-nds.iaea.org/relnsd/NdsEnнд/QueryForm.html, accessed on 7 May 2021).

**Acknowledgments:** The authors thank Michal Bejger for proofreading the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Notes**

1. http://www.physics.mcgill.ca/~pulsar/magnetar/main.html, accessed on 7 May 2021.

2. The demonstration follows directly from the definitions $P_e = n_e^2 d(\epsilon_e/n_e)/dn_e$ and $\gamma_e = (d\epsilon_e/dn_e)/(m_e c^2)$, where $\epsilon_e$ denotes the energy density of the electron Fermi gas.

3. http://www.astro.ulb.ac.be/bruslib/, accessed on 7 May 2021.

4. https://www-nds.iaea.org/relnsd/NdsEnнд/QueryForm.html, accessed on 7 May 2021.

5. http://magnetars.ice.csic.es, accessed on 7 May 2021.

**References**

1. Duncan, R.C.; Thompson, C. Formation of Very Strongly Magnetized Neutron Stars: Implications for Gamma-ray Bursts. Astrophys. J. Lett. **1992**, *392*, L9. [CrossRef]

2. Esposito, P.; Rea, N.; Israel, G.L. Magnetars: A Short Review and Some Sparse Considerations. In Astrophysics and Space Science Library; Belloni, T.M., Méndez, M., Zhang, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2021; Volume 461, pp. 97–142. [CrossRef]

3. Olausen, S.A.; Kaspi, V.M. The McGill Magnetar Catalog. Astrophys. J. Suppl. Ser. **2014**, *212*, 6. [CrossRef]

4. Coti Zelati, F.; Rea, N.; Pons, J.A.; Campana, S.; Esposito, P. Systematic study of magnetar outbursts. Mon. Not. R. Astron. Soc. **2018**, *474*, 961–1017. [CrossRef]

5. De Grandis, D.; Turolla, R.; Wood, T.S.; Zane, S.; Taverna, R.; Gourgouliatos, K.N. Three-dimensional Modeling of the Magnetothermal Evolution of Neutron Stars: Method and Test Cases. Astrophys. J. **2020**, *903*, 40. [CrossRef]

6. Beloborodov, A.M.; Li, X. Magnetar Heating. Astrophys. J. **2016**, *833*, 261. [CrossRef]
7. Kaminker, A.D.; Yakovlev, D.G.; Potekhin, A.Y.; Shibazaki, N.; Shternin, P.S.; Gnedin, O.Y. Magnetars as cooling neutron stars with internal heating. *Mon. Not. R. Astron. Soc.* **2006**, *371*, 477–483. [CrossRef]

8. Kaminker, A.D.; Potekhin, A.Y.; Yakovlev, D.G.; Chabrier, G. Heating and cooling of magnetars with accreted envelopes. *Mon. Not. R. Astron. Soc.* **2009**, *395*, 2257–2267. [CrossRef]

9. Cooper, R.L.; Kaplan, D.L. Magnetic Field-Decay-Induced Electron Captures: A Strong Heat Source in Magnetar Crusts. *Astrophys. J. Lett.* **2010**, *708*, L80–L83. [CrossRef]

10. Fantina, A.F.; Zdunik, J.L.; Chamel, N.; Pearson, J.M.; Haensel, P.; Goriely, S. Crustal heating in accreting neutron stars from the nuclear energy-density functional theory. I. Proton shell effects and neutron-matter constraint. *Astron. Astrophys.* **2018**, *620*, A105. [CrossRef]

11. Koliopanos, F.; Vasilopoulos, G.; Godet, O.; Bachetti, M.; Webb, N.A.; Barret, D. ULX spectra revisited: Accreting, highly magnetized neutron stars as the engines of ultraluminous X-ray sources. *Astron. Astrophys.* **2017**, *608*, A47. [CrossRef]

12. Tong, H.; Wang, W. Accreting magnetars: Linking ultraluminous X-ray pulsars and the slow pulsation X-ray pulsars. *Mon. Not. R. Astron. Soc.* **2019**, *482*, 4956–4964. [CrossRef]

13. Doroshenko, V.; Santangelo, A.; Suleimanov, V.F.; Tsygankov, S.S. An observational argument against accretion in magnetars. *Astron. Astrophys.* **2020**, *643*, A173. [CrossRef]

14. Brice, N.; Zane, S.; Turolla, R.; Wu, K. Super-Eddington Emission from Accreting, Highly Magnetised Neutron Stars with a Multipolar Magnetic Field. *Mon. Not. R. Astron. Soc.* **2021**, *504*, 701–715. [CrossRef]

15. Iida, K.; Sato, K. Spin Down of Neutron Stars and Compositional Transitions in the Cold Crustal Matter. *Astrophys. J.* **1997**, *477*, 294–312. [CrossRef]

16. Chamel, N.; Fantina, A.F.; Zdunik, J.L.; Haensel, P. Experimental constraints on shallow heating in accreting neutron star crusts. *Phys. Rev. C* **2020**, *102*, 015804. [CrossRef]

17. Rabi, I.I. Das freie Elektron im homogenen Magnetfeld nach der Diracschen Theorie. *Z. Phys.* **1928**, *49*, 507–511. [CrossRef]

18. Landau, L. Diamagnetismus der Metalle. *Z. Phys.* **1930**, *64*, 629–637. [CrossRef]

19. Haensel, P.; Potekhin, A.Y.; Yakovlev, D.G. *Neutron Stars. 1. Equation of State and Structure*; Springer: New York, NY, USA, 2007.

20. Chamel, N.; Stoyanov, Z.K.; Mihailov, L.M.; Mutafchieva, Y.D.; Pavlov, R.L.; Velchev, C.J. Role of Landau quantization on the neutron-drip transition in magnetar crusts. *Phys. Rev. C* **2015**, *91*, 065801. [CrossRef]

21. Mutafchieva, Y.D.; Chamel, N.; Stoyanov, Z.K.; Pearson, J.M.; Mihailov, L.M. Role of Landau–Rabi quantization of electron motion on the crust of magnetars within the nuclear energy density functional theory. *Phys. Rev. C* **2019**, *99*, 055805. [CrossRef]

22. Potekhin, A.Y.; Chabrier, G. Equation of state for magnetized Coulomb plasmas. *Astron. Astrophys.* **2013**, *550*, A43. [CrossRef]

23. Blandford, R.D.; Hernquist, L. Magnetic susceptibility of a neutron star crust. *J. Phys. C Solid State Phys.* **1982**, *15*, 6233–6243. [CrossRef]

24. Chamel, N.; Fantina, A.F. Electron capture instability in magnetic and nonmagnetic white dwarfs. *Phys. Rev. D* **2015**, *92*, 023008. [CrossRef]

25. Hollerbach, R.; Rüdiger, G. Hall drift in the stratified crusts of neutron stars. *Mon. Not. R. Astron. Soc.* **2004**, *347*, 1273–1278. [CrossRef]

26. Pons, J.A.; Geppert, U. Magnetic field dissipation in neutron star crusts: From magnetars to isolated neutron stars. *Astron. Astrophys.* **2007**, *470*, 303–315. [CrossRef]

27. Viganò, D.; Rea, N.; Pons, J.A.; Perna, R.; Aguilera, D.N.; Miralles, J.A. Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models. *Mon. Not. R. Astron. Soc.* **2013**, *434*, 123–141. [CrossRef]

28. Geppert, U.; Viganò, D. Creation of magnetic spots at the neutron star surface. *Mon. Not. R. Astron. Soc.* **2014**, *444*, 3198–3208. [CrossRef]

29. Gourgouliatos, K.N.; Cumming, A. Hall Attractor in Axially Symmetric Magnetic Fields in Neutron Star Crusts. *Phys. Rev. Lett.* **2014**, *112*, 171101. [CrossRef] [PubMed]

30. Wood, T.S.; Hollerbach, R. Three Dimensional Simulation of the Magnetic Stress in a Neutron Star Crust. *Phys. Rev. Lett.* **2015**, *114*, 191101. [CrossRef]

31. Bransgrove, A.; Levin, Y.; Beloborodov, A. Magnetic field evolution of neutron stars-I. Basic formalism, numerical techniques and first results. *Mon. Not. R. Astron. Soc.* **2018**, *473*, 2771–2790. [CrossRef]

32. Igoshev, A.P.; Hollerbach, R.; Wood, T.; Gourgouliatos, K.N. Strong toroidal magnetic fields required by quiescent X-ray emission of magnetars. *Nat. Astron.* **2021**, *5*, 145–149. [CrossRef]

33. Kojima, Y.; Kisaka, S.; Fujisawa, K. Evolution of magnetic deformation in neutron star crust. *Mon. Not. R. Astron. Soc.* **2021**, *502*, 2097–2104. [CrossRef]

34. Horowitz, C.J.; Dussan, H.; Berry, D.K. Fusion of neutron-rich oxygen isotopes in the crust of accreting neutron stars. *Phys. Rev. C* **2008**, *77*, 045807. [CrossRef]

35. Rashdan, M. Sub-barrier fusion calculations for the neutron star crust using the microscopic Brueckner G-matrix and Skyrme energy density functionals. *Phys. Rev. C* **2015**, *91*, 054613. [CrossRef]

36. Chatterjee, P.; Hernquist, L.; Narayan, R. An Accretion Model for Anomalous X-ray Pulsars. *Astrophys. J.* **2000**, *534*, 373–379. [CrossRef]

37. Alpar, M.A. On Young Neutron Stars as Propellers and Accretors with Conventional Magnetic Fields. *Astrophys. J.* **2001**, *554*, 1245–1254. [CrossRef]
38. Ertan, Ü.; Erkut, M.H.; Ekşi, K.Y.; Alpar, M.A. The Anomalous X-Ray Pulsar 4U 0142+61: A Neutron Star with a Gaseous Fallback Disk. Astrophys. J. **2007**, *657*, 441–447. [CrossRef]

39. Fantina, A.F.; Chamel, N.; Mutafchieva, Y.D.; Stoyanov, Z.K.; Mihailov, L.M.; Pavlov, R.L. Role of the symmetry energy on the neutron-drip transition in accreting and nonaccreting neutron stars. Phys. Rev. C **2016**, *93*, 015801. [CrossRef]

40. Chamel, N.; Pavlov, R.L.; Mihailov, L.M.; Velchev, C.J.; Stoyanov, Z.K.; Mutafchieva, Y.D.; Ivanovich, M.D.; Pearson, J.M.; Goriely, S. Properties of the outer crust of strongly magnetized neutron stars from Hartree-Fock-Bogoliubov atomic mass models. Phys. Rev. C **2012**, *86*, 055804. [CrossRef]

41. Salpeter, E.E. Electrons Screening and Thermonuclear Reactions. Aust. J. Phys. **1954**, *7*, 373. [CrossRef]

42. Chamel, N.; Fantina, A.F.; Zdunik, J.L.; Haensel, P. Neutron drip transition in accreting and nonaccreting neutron star crusts. Phys. Rev. C **2015**, *91*, 055803. [CrossRef]

43. Yakovlev, D.G.; Gasques, L.R.; Afanasjev, A.V.; Beard, M.; Wiescher, M. Fusion reactions in multicomponent dense matter. Phys. Rev. C **2006**, *74*, 035803. [CrossRef]

44. Wang, M.; Audi, G.; Kondev, F.G.; Huang, W.J.; Naimi, S.; Xu, X. The AME2016 atomic mass evaluation (II). Tables, graphs and references. Chin. Phys. C **2017**, *41*, 030003. [CrossRef]

45. Goriely, S.; Chamel, N.; Pearson, J.M. Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. XIII. The 2012 atomic mass evaluation and the symmetry coefficient. Phys. Rev. C **2013**, *88*, 024308. [CrossRef]

46. Xu, Y.; Goriely, S.; Jorissen, A.; Chen, G.L.; Arnould, M. Databases and tools for nuclear astrophysics applications. BRUSsels Nuclear LIBrary (BRUSLIB), Nuclear Astrophysics Compilation of REactions II (NACRE II) and Nuclear NETwork GENerator (NETGEN). Astron. Astrophys. **2013**, *549*, A106. [CrossRef]

47. Ho, W.C.G.; Heinke, C.O. A neutron star with a carbon atmosphere in the Cassiopeia A supernova remnant. Nature **2009**, *462*, 71–73. [CrossRef]

48. Lunney, D.; Pearson, J.M.; Thibault, C. Recent trends in the determination of nuclear masses. Rev. Mod. Phys. **2003**, *75*, 1021–1082. [CrossRef]

49. Pearson, J.M.; Chamel, N.; Potekhin, A.Y.; Fantina, A.F.; Ducoin, C.; Dutta, A.K.; Goriely, S. Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals-I. Role of symmetry energy. Mon. Not. R. Astron. Soc. **2018**, *481*, 2994–3026. [CrossRef]

50. Ingraham, R.L.; Wilkes, J.M. Inhibition of Degeneracy by Intense Magnetic Fields-Derivation and Astrophysical Application. Astrophys. Space Sci. **1987**, *135*, 87–104. [CrossRef]

51. Chatterjee, D.; Elghozi, T.; Novak, J.; Oertel, M. Consistent neutron star models with magnetic field-dependent equations of state. Mon. Not. R. Astron. Soc. **2015**, *447*, 3785–3796. [CrossRef]

52. Chatterjee, D.; Fantina, A.F.; Chamel, N.; Novak, J.; Oertel, M. On the maximum mass of magnetized white dwarfs. Mon. Not. R. Astron. Soc. **2017**, *469*, 95–109. [CrossRef]

53. Braithwaite, J. Axisymmetric magnetic fields in stars: Relative strengths of poloidal and toroidal components. Mon. Not. R. Astron. Soc. **2009**, *397*, 763–774. [CrossRef]

54. Tiengo, A.; Esposito, P.; Mereghetti, S.; Turolla, R.; Nobili, L.; Gastaldello, F.; Götz, D.; Israel, G.L.; Rea, N.; Stella, L.; et al. A variable absorption feature in the X-ray spectrum of a magnetar. Nature **2013**, *500*, 312–314. [CrossRef]

55. Makishima, K.; Enoto, T.; Hiraga, J.S.; Nakano, T.; Nakazawa, K.; Sakurai, S.; Sasano, M.; Murakami, H. Possible Evidence for Free Precession of a Strongly Magnetized Neutron Star in the Magnetar 4U 0142+61. Phys. Rev. Lett. **2014**, *112*, 171102. [CrossRef]

56. Borghese, A.; Rea, N.; Coti Zelati, F.; Tiengo, A.; Turolla, R. Discovery of a Strongly Phase-variable Spectral Feature in the Isolated Neutron Star RX J0720.4-3125. Astrophys. J. Lett. **2015**, *807*, L20. [CrossRef]

57. Rodríguez Castillo, G.A.; Israel, G.L.; Tiengo, A.; Salvetti, D.; Turolla, R.; Zane, S.; Rea, N.; Esposito, P.; Mereghetti, S.; Perna, R.; et al. The outburst decay of the low magnetic field magnetar SWIFT J1822.3-1606: Phase-resolved analysis and evidence for a variable cyclotron feature. Mon. Not. R. Astron. Soc. **2016**, *456*, 4145–4155. [CrossRef]

58. Makishima, K.; Enoto, T.; Yoneda, H.; Odaka, H. A NuSTAR confirmation of the 36 ks hard X-ray pulse-phase modulation in the magnetar 1E 1547.0-5408. Mon. Not. R. Astron. Soc. **2021**, *502*, 2266–2284. [CrossRef]

59. Chandrasekhar, S.; Prendergast, K.H. The Equilibrium of Magnetic Stars. Proc. Natl. Acad. Sci. USA **1956**, *42*, 5–9. [CrossRef]

60. Pearson, J.M.; Goriely, S.; Chamel, N. Properties of the outer crust of neutron stars from Hartree-Fock-Bogoliubov mass models. Phys. Rev. C **2011**, *83*, 065810. [CrossRef]

61. Sur, A.; Haskell, B. The impact of superconductivity and the Hall effect in models of magnetized neutron stars. arXiv **2021**, arXiv:2104.14908.

62. Gusakov, M.E.; Kantor, E.M.; Reisenegger, A. Rotation-induced deep crustal heating of millisecond pulsars. Mon. Not. R. Astron. Soc. Lett. **2015**, *453*, L36–L40. [CrossRef]

63. Peña Arteaga, D.; Grasso, M.; Khan, E.; Ring, P. Nuclear structure in strong magnetic fields: Nuclei in the crust of a magnetar. Phys. Rev. C **2011**, *84*, 045806. [CrossRef]

64. Stein, M.; Maruhn, J.; Sedrakian, A.; Reinhard, P.G. Carbon-oxygen-neon mass nuclei in superstrong magnetic fields. Phys. Rev. C **2016**, *94*, 035802. [CrossRef]