Heating in Magnetar Crusts from Electron Captures

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}-10^{11}$ g cm$^{-3}$ with heat power $W^\infty\sim 10^{35}-10^{36}$ 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.


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 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]. 1 http://www.physics.mcgill.ca/~pulsar/magnetar/main.html -accessed on 7 May 2021 arXiv:2106.05040v1 [astro-ph.HE] 9 Jun 2021 In particular, these mechanisms are only effective deep enough beneath the surface of the star where the temperature lies below the melting temperature 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][12][13][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.

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

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]): and provided the mass density ρ does not exceed: 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 λ e =h/(m e c) is the electron Compton wavelength (h being the Dirac-Planck constant), and we introduced the dimensionless magnetic field strength B ≡ B/B rel with: where α = e 2 /(hc) is the fine-structure constant (e being the elementary electric charge). The condition T < T B is fulfilled in magnetars since B 1, and the temperature is typically of order 10 8 − 10 9 K (see, e.g., [8]). The threshold density ρ B coincides with the density associated with the neutron drip transition (delimiting the outer and inner parts of the crust) for B ≈ 5.72 × 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]: where γ e is the electron Fermi energy in units of the electron mass energy m e c 2 , and ρ 6 = ρ/10 6 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 ≈ 8.83 × 10 16 G) and B = 3000 (B ≈ 1.32 × 10 17 G). As we show, electron captures occur at pressures well above 10 28 dyn cm −2 (densities above 10 10 g cm −3 ). As can be seen in Figures 1 and 2, T Fe exceeds 10 10 K for the crustal regions of interest; therefore, electrons are highly degenerate.    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]: where a N = (3Z/(4πn e )) 1/3 is the ion sphere's radius, n e the electron number density, and Γ m the Coulomb coupling parameter at melting. As shown in Figures 1 and 2, T m is of order 10 9 K (assuming Γ 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 Γ 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 = Γ m T m , i.e., the Coulomb potential energy Z 2 e 2 /a N is much larger than the thermal energy k B T.
In the following, we thus neglect thermal effects.

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 µ (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µ ≈ 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 β (depending on the nucleus under consideration, as well as the magnetic field strength): Recalling that the magnetic field in the crust of a neutron star roughly evolves on timescales τ B of order ∼ Myr (see, e.g., [5,[25][26][27][28][29][30][31][32][33]), the above reaction essentially proceeds in quasiequilibrium 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: 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]:

Baryon Chemical Potential and Matter Pressure
The model we adopted here was described in [40]. The baryon chemical potential is given by: 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]: For nuclei in excited states, the baryon chemical potential, which we denote by µ * , takes a similar form except that M (A, Z)c 2 must be replaced by 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., Chap. 4 of [19]): where we introduced the dimensionless parameter: 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.

Onset of Electron Captures
The onset of electron captures by nuclei (A, Z) is determined by the condition µ(A, Z, P β , B ) = µ(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 α as [24]: where we introduced the Q-value (in vacuum) associated with electron capture by nuclei (A, Z): These Q-values can be obtained from the tabulated Q-values of β decay by the following relation: 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 γ e 1 in Equation (12): Equation (13) thus reduces to: Introducing the dimensionless parameter: the solutions (real roots) of Equation (20) for γ e are given by the following formulas [16]: Using Equation (11), the threshold pressure and baryon number density are thus given by: respectively. The transition is accompanied by a discontinuous change of density given by: The threshold condition for transitions to daughter nuclei in excited states is . The corresponding pressure P * β and density n * β (A, Z) can be obtained from the previous formulas by merely substituting . As discussed in [39,42], the first electron capture by the nucleus (A, Z) may be accompanied by the emission of ∆N > 0 neutrons. The corresponding pressure P drip and baryon density n drip are given by similar expressions as Equations (23) 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 . 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 γ e at the neutron drip point must obey the following inequality (see, e.g., [40]):

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 µ(A, 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: in the first case and: where we used Equation (11) and the relation 2 dP e = n e m e c 2 dγ e . Expanding the electron Fermi energy γ e (n e + δn e ) to first order in δn e /n e , substituting in the expression of µ(A, Z − 2, P β , B ), and using Equations (13)- (15) and (30), we finally obtain: where γ e is calculated from Equation (22). Equation (31) , which is generally satisfied for even A nuclei, but not necessarily for odd A nuclei. In the latter case, we typically have . In other words, as the pressure reaches P β (A, Z, B ), 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 ) = 0. The daughter nucleus sinks deeper in the crust and only captures a second electron in quasi- 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 ≈ −1.4621 for the liquid (this value was calculated using the constant A 1 on page 75 of [19] as C = A 1 (4π/3) 1/3 ) to C ≈ −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 ) released from the ground-state-to-excited-state transitions can be obtained from Equation (31) 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: 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.

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 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 β (A, Z) or P pyc = P * β (A, Z), respectively: or: 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: with γ 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: The expressions for q * pyc (A, Z) and q * fus (A, Z) follow after substituting M (A, Z − 1)c 2 with M (A, Z − 1)c 2 + E ex (A, Z − 1). Numerical results are presented and discussed in the next section.

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 3 [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 4 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 ±h (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 ρ = 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 γ e 1, 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].    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 = 2000 and B = 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 1. Maximum possible heat released q per nucleon (in MeV) from electron captures in the outer crust of a magnetar with B = 2000. The pressure P β (in dyn cm −2 ) and density ρ β (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.  30 6.34 × 10 11 122 Zr → 122 Sr − 2e − + 2ν e ( ) 0.007 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 = 2000. The pressure P β (in dyn cm −2 ) and density ρ β (in g cm −3 ) at which electron capture occurs are given considering ground-state-to-ground-state transitions.
8.46 × 10 29 2.47 × 10 11 12 C → 12 Be − 2e − + 2ν e 0.143 5.38 × 10 29 1.98 × 10 11 The heating resulting from pycnonuclear fusions of carbon and oxygen is summarized in Tables 5 and 6 for B = 2000 and B = 3000, respectively. Contrary to our previous study of accreting neutron stars [16], we did not consider here fusions of oxygen occurring at pressure P * β because we realized that µ β e > µ drip e : 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 pyc could be much smaller, especially if fusions occur at densities substantially lower than n β . 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 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 = 2000. The numbers in parentheses indicate the contributions from the sole fusion. The pressure P pyc (in dyn cm −2 ) and density ρ pyc (in g cm −3 ) at which fusion occurs were fixed by the onset of electron captures considering ground-state-toground-state transitions. 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 β (A, Z, B ) and density n β (A, Z, B ). 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 = 3000, the results obtained from Equations (23), (24), and (31) differed from the exact values by −7.6 × 10 −2 %, −3.6 × 10 −2 %, and 5.0 × 10 −3 %, respectively.

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): where Φ is the gravitational potential obeying Poisson's equation (G being the gravitational constant): and the free electric charge current j j j free is given by Maxwell's equation: The magnetization of a relativistic electron gas in a strongly quantizing magnetic field is negligibly small [23]; therefore, we can safely replace H H H ≈ B B B. Expanding the Lorentz force term in the right-hand side of Equation (40) using Equation (42) yields: 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 from the spin down, as measured by timing analyses). These theoretical expectations were corroborated by astrophysical observations [32,[54][55][56][57][58]. Assuming for simplicity that the magnetic field is axially symmetric and purely toroidal, i.e., B B B = B( , z)1 φ 1 φ 1 φ in cylindrical coordinates [59] ( is the radial distance, z the height, φ the azimuthal angle, and 1 φ 1 φ 1 φ the associated unit vector), the last term in Equation (43) 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 −2 R and δM ∼ 10 −5 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 2 /(4πR) is much smaller than the second term in the left-hand side of order B 2 /(8π∆R). We thus dropped the last term in Equation (43) as in [9]. Introducing the local gravitational field g g g = −∇ ∇ ∇Φ, which is essentially uniform in the outer crust and given by g ≈ GM/R 2 , Equation (43) finally reduces to: 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 2 ) < 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.

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: must remain constant. Initially, the nuclei (A, Z) are present in the crust in a layer delimited by the matter pressures P min (A, Z, B + δB ) and P max (A, Z, B + δB ), 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 such that the shallowest layer containing this element can be compressed up to the threshold pressure P β (A, Z, B ), thus maintaining the same total pressure. Using Equation (45) and expanding to first order in δB , we find: where the pressures are evaluated here for the same magnetic field strength B . Results are summarized in Tables 7 and 8 for B = 2000 and B = 3000, respectively, using the results of [21] for the initial pressures P min (note that the initial pressure of 56 Fe is simply P min = 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 min = 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 β . The changes of magnetic field required to trigger the various chains of reactions depend on the layer, but are typically of order δB /B ≈ 10 −3 − 10 −4 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, δ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 δB /B of the magnetic field thus drops from 2.09 × 10 −2 to 2.21 × 10 −3 when the magnetic field is varied from B = 100 to B = 1000. In the ultrarelativistic regime (γ e 1), both P min (A, Z, B ) and P β (A, Z, B ) increase linearly with B so that δB /B ∝ 1/B . The time τ ∼ τ B δB /B after which the magnetic field has changed by an amount δB is thus approximately inversely proportional to the magnetic field strength (τ B ∼ 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 δB /B ∼ τ/τ B ∼ 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 Σ A q A,tot (∆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 ∆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 ∼0.1 MeV for all layers, the heat power can be estimated as W ∞ ∼ q A,tot (∆M/m u )/τ, where ∆M = ∑ A ∆M A is the mass of the outer crust. With ∆M ∼ 10 −5 − 10 −4 M [60] with M the mass of the Sun, and τ ∼ 1 kyr, we found W ∞ ∼ 10 35 − 10 36 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 δRδB/τ ∼ 10 42 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 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.    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 = 2000, from zero pressure (at the stellar surface) to the pressure at which the last capture occurs.

Heating Induced by Spin-Down
Observed isolated magnetars rotate at frequencies f ∼ 0.1 − 0.5 Hz and exhibit a very high spin-down rate. They are believed to be born with initial frequencies f i ∼ 1 kHz. The time to spin-down to a frequency f f i can be estimated as: 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 ∼ 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.

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 δB/B ∼ 10 −3 − 10 −4 for B ∼ 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 ∞ ∼ 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 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.