Next Article in Journal
The Role of Radio Observations in Studies of Infrared-Bright Galaxies: Prospects for a Next-Generation Very Large Array
Next Article in Special Issue
The Elasticity of the Neutron Star Mantle: The Improved Compressible Liquid Drop Model for Cylindrical Phases
Previous Article in Journal
Estimating the Transit Speed and Time of Arrival of Interplanetary Coronal Mass Ejections Using CME and Solar Flare Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Onset of Electron Captures and Shallow Heating in Magnetars

by
Nicolas Chamel
1,* and
Anthea Francesca Fantina
1,2,†
1
Institute of Astronomy and Astrophysics, Université Libre de Bruxelles, CP 226, Boulevard du Triomphe, B-1050 Brussels, Belgium
2
Grand Accélérateur National d’Ions Lourds (GANIL), CEA/DRF-CNRS/IN2P3, Boulevard Henri Becquerel, 14076 Caen, France
*
Author to whom correspondence should be addressed.
Anthea Francesca Fantina's main affiliation is the affiliation 2.
Submission received: 13 May 2022 / Revised: 3 June 2022 / Accepted: 8 June 2022 / Published: 11 June 2022
(This article belongs to the Special Issue Advances in Magnetars)

Abstract

:
The loss of magnetic pressure accompanying the decay of the magnetic field in a magnetar may trigger exothermic electron captures by nuclei in the shallow layers of the stellar crust. Very accurate analytical formulas are obtained for the threshold density and pressure, as well as for the maximum amount of heat that can be possibly released, taking into account the Landau–Rabi quantization of electron motion. These formulas are valid for arbitrary magnetic field strengths, from the weakly quantizing regime to the most extreme situation in which electrons are all confined to the lowest level. Numerical results are also presented based on experimental nuclear data supplemented with predictions from the Brussels-Montreal model HFB-24. This same nuclear model has been already employed to calculate the equation of state in all regions of magnetars.

1. Introduction

Soft gamma-ray repeaters and anomalous X-ray pulsars are two facets of a very active subclass of neutron stars, called magnetars, exhibiting outbursts and less frequently giant flares that release huge amounts of energy up to ∼ 10 46 erg within a second (see e.g., [1] for a recent review). These phenomena are thought to be powered by internal magnetic fields exceeding 10 14 10 15 G [2]. At the date of this writing, 24 such objects have been discovered and six more candidates remain to be confirmed according to the McGill Online Magnetar Catalog [3]. Their persistent X-ray luminosity ∼ 10 33 10 35 erg/s, which is well in excess of their rotational energy and which implies a higher surface temperature than in weakly magnetized neutron stars of the same age [4], provides further evidence for extreme magnetic fields [5]. It is widely thought that heat is generated by the deformations of the crust beyond the elastic limit due to magnetic stresses (see, e.g., [6]). This mechanism is most effective in the inner region of the crust, where crystallization first occurs [7,8]. However, it has been demonstrated that heat sources should be located in the shallow region of the crust to avoid excessive neutrino losses [9,10]. Alternatively, the magnetic energy in the outer crust may be dissipated into heat through electron captures by nuclei triggered by the magnetic field evolution [11]. This mechanism is analogous to crustal heating in accreting neutron stars [12], the matter compression being induced here by the loss of magnetic support rather than accretion from a stellar companion.
We have recently estimated the maximum amount of heat that could be possibly released by electron captures and the location of the heat sources taking into account Landau–Rabi quantization of electron motion induced by the magnetic field [13]. For simplicity, we focused on the strongly quantizing regime in which only the lowest Landau–Rabi level is occupied, thus allowing for a simple analytical treatment. Results are extended here to arbitrary magnetic fields. We demonstrate that the weakly quantizing regime is also amenable to accurate analytical approximations. Results are presented based on experimental nuclear data supplemented with the Brussels-Montreal atomic mass table HFB-24 [14]. The underlying nuclear energy-density functional BSk24 has been already applied to construct unified equations of state for both unmagnetized neutron stars [15,16,17] and magnetars [18].
The paper is organized as follows. In Section 2, we present the equation of state of the outer crust of a magnetar and the approximations we made. In Section 3, we give the equations to determine the initial composition of the outer crust and the boundaries delimiting different layers. Our analytical treatment of the heating from electron captures is described in Section 4. Numerical results including detailed error estimates are presented and discussed in Section 5.

2. Equation of State of Magnetar Crusts

In the following, we shall consider the crustal region at densities above the ionization threshold and below the neutron-drip point. We assume that each crustal layer is made of fully ionized atomic nuclei ( A , Z ) with proton number Z and mass number A embedded in a relativistic electron gas.

2.1. Main Equations

Whereas nuclei with number density n N exert a negligible pressure P N 0 , they contribute to the mass-energy density
E N = n N M ( A , Z , B ) c 2 ,
where M ( A , Z , B ) denotes the ion mass including the rest mass of Z electrons. In principle, M ( A , Z , B ) may also depend on the magnetic field, which will be conveniently measured in terms of the dimensionless ratio B B / B rel with
B rel = m e 2 c 3 e 4.41 × 10 13 G ,
where m e is the electron mass, c is the speed of light, is the Planck–Dirac constant and e is the elementary electric charge.
To a very good approximation, electrons can be treated as an ideal Fermi gas. In the presence of a magnetic field, the electron motion perpendicular to the field is quantized into Landau–Rabi levels [19,20]. The observed surface magnetic field on a magnetar is typically B s 10 14 10 15 G [3,21,22]. The internal magnetic field B is expected to be even stronger and could potentially reach 10 17 10 18 G (see, e.g., [23]). In our previous study [13], we assumed for simplicity that the magnetic field is strongly quantizing, meaning that electrons remain all confined to the lowest level throughout the outer crust, thus requiring B 5.72 × 10 16 G [24]. Even if weaker fields are considered, quantization effects are not expected to be completely washed out by thermal effects. Indeed, the temperatures T 10 8 10 9 K prevailing in a magnetar for which B 1 (see e.g., [10]) are much lower than the characteristic temperature
T B = m e c 2 k B B 5.93 × 10 9 B K ,
where k B denotes Boltzmann’s constant. Strictly speaking, Equation (3) is only relevant in the strongly quantizing regime. If several Landau–Rabi levels are populated, the characteristic temperature is reduced but only by a factor of a few at most at the bottom of the outer crust (see Chap. 4 in [25]). Neglecting the small electron anomalous magnetic moment and ignoring thermal effects, the electron energy density (with the rest-mass excluded) and electron pressure are given by
E e = B m e c 2 ( 2 π ) 2 λ e 3 ν = 0 ν max g ν ( 1 + 2 ν B ) ψ + x e ( ν ) 1 + 2 ν B n e m e c 2 ,
P e = B m e c 2 ( 2 π ) 2 λ e 3 ν = 0 ν max g ν ( 1 + 2 ν B ) ψ x e ( ν ) 1 + 2 ν B ,
respectively, where we have introduced the electron Compton wavelength λ e = / ( m e c ) , g ν = 1 for ν = 0 and g ν = 2 for ν 1 ,
ψ ± ( x ) = x 1 + x 2 ± ln ( x + 1 + x 2 ) ,
x e ( ν ) = γ e 2 1 2 ν B ,
and ν max is fixed by the electron number density n e given by
n e = 2 B ( 2 π ) 2 λ e 3 ν = 0 ν max g ν x e ( ν ) .
Here γ e denotes the electron Fermi energy in units of m e c 2 . The index ν max is the highest integer for which γ e 2 1 2 ν max B 0 , i.e.
ν max = γ e 2 1 2 B ,
where [ . ] denotes the integer part. The mean baryon number density follows from the requirement of electric charge neutrality
n ¯ = A Z n e = A n N .
The main correction to the ideal electron Fermi gas arises from the electron-ion interactions. According to the Bohr-van Leeuwen theorem [26], the electrostatic corrections to the energy density and to the pressure are independent of the magnetic field apart from a negligibly small contribution due to quantum zero-point motion of ions about their equilibrium position [27]. For pointlike ions embedded in a uniform electron gas, the corresponding energy density is given by (see e.g., Chap. 2 of [25])
E L = C M 4 π 3 1 / 3 e 2 n e 4 / 3 Z 2 / 3 ,
where C M is the Madelung constant. The contribution to the pressure is thus given by
P L = n e 2 d ( E L / n e ) d n e = E L 3 .
The pressure of the Coulomb plasma finally reads P = P e + P L , whereas the energy density is given by E = E N + E e + E L .
For ions arranged in a body-centered cubic lattice, the Madelung constant is given by C M = 0.895929255682 [28]. However, the electron-ion plasma may not necessarily be in a solid state, especially in the shallow layers, which are the main focus of this work. The crystallization temperature can be estimated as [7]:
T m 1.3 × 10 5 Z 2 175 Γ m ρ 6 A 1 / 3 K ,
where ρ 6 is the density in units of 10 6 g cm 3 , and Γ m is the Coulomb coupling parameter at melting. In the absence of magnetic field, Γ m 175 and T m is typically of order 10 9 K. The presence of a magnetic field tends to lower Γ m , thus increasing T m [29]. In any case, the Madelung constant in the liquid phase remains very close to that of the solid phase. In the following, we will adopt the Wigner–Seitz estimate C M = 9 / 10 for the Madelung constant [30]. Thermal effects on thermodynamic quantities are small and will be neglected.

2.2. Weakly Quantizing Magnetic Field

The magnetic field is weakly quantizing if many Landau–Rabi levels are filled: ν max 1 . Using the expansions (41) obtained in [31] for the electron density n e leads to the following estimate for the mean baryon number density:
n ¯ A 2 π 2 Z λ e 3 2 3 γ e 2 1 3 / 2 + ( 2 B ) 3 / 2 ζ 1 2 , γ e 2 1 2 B + B 2 6 γ e 2 1 ,
where ζ ( z , q ) is the Hurwitz zeta function defined by
ζ ( z , q ) = ν = 0 + 1 ( ν + q ) z
for ( z ) > 1 and by analytic continuation to other z 1 (excluding poles ν + q = 0 ). The first term in Equation (14) represents the mean baryon number density in the absence of magnetic field. The second term accounts for quantum oscillations due to the filling of Landau–Rabi levels, while the last term is a higher-order magnetic correction.
The expression for the associated expansion of the pressure is more involved. In the notations of [31], the electron contribution P e to the pressure can be directly obtained from the grand potential density by P e = ω 0 ( mon ) ω 0 ( osc ) . Using Equations (41), (43) and (44) of [31] yields ( α = e 2 / ( c ) is the fine-structure constant):
P e m e c 2 4 π 2 λ e 3 { 1 2 1 2 B + 2 B 2 3 log γ e + 2 B + γ e 2 1 1 + 2 B 1 2 γ e 2 B + γ e 2 1 2 B + 1 3 γ e 2 B + γ e 2 1 3 2 B 3 + B arccosh γ e γ e γ e 2 1 ( 2 B ) 5 / 2 0 + ζ ˜ 3 ( 1 / 2 , q + 1 ) 1 + 2 B q d q + 2 3 ( 2 B ) 5 / 2 γ e ζ 3 2 , γ e 2 1 2 B + 2 15 ( 2 B ) 7 / 2 γ e 3 ζ 5 2 , γ e 2 1 2 B + 1 240 B γ e 4 + 4 B 2 0 1 ζ 1 2 , q ζ 1 2 , q + 1 2 B d q } ,
with
ζ ˜ 3 ( z , q ) = ζ ( z , q ) 1 z 1 q z + 1 1 2 q z z 12 q z 1 .
The total pressure is found by adding the electrostatic correction (12) using the expansion for the electron density:
P m e c 2 4 π 2 λ e 3 { 1 2 1 2 B + 2 B 2 3 log γ e + 2 B + γ e 2 1 1 + 2 B 1 2 γ e 2 B + γ e 2 1 2 B + 1 3 γ e 2 B + γ e 2 1 3 2 B 3 + B arccosh γ e γ e γ e 2 1 ( 2 B ) 5 / 2 0 + ζ ˜ 3 ( 1 / 2 , q + 1 ) 1 + 2 B q d q + 2 3 ( 2 B ) 5 / 2 γ e ζ 3 2 , γ e 2 1 2 B + 2 15 ( 2 B ) 7 / 2 γ e 3 ζ 5 2 , γ e 2 1 2 B + 1 240 B γ e 4 + 4 B 2 0 1 ζ 1 2 , q ζ 1 2 , q + 1 2 B d q + 2 3 2 3 π 1 / 3 C M α Z 2 / 3 [ 2 3 γ e 2 1 3 / 2 + ( 2 B ) 3 / 2 ζ 1 2 , γ e 2 1 2 B + B 2 6 γ e 2 1 ] 4 / 3 } .
In the absence of magnetic field B = 0 (corresponding to the limit ν max + ), the mean baryon number density and the pressure reduce, respectively, to
n ¯ = A Z x r 3 3 π 2 λ e 3 ,
P = m e c 2 8 π 2 λ e 3 x r 2 3 x r 2 1 1 + x r 2 + ln ( x r + 1 + x r 2 ) + C M α 3 4 243 π 7 1 / 3 x r 4 m e c 2 λ e 3 Z 2 / 3 .
Here, x r denotes the relativity parameter defined by x r = λ e k e and k e = ( 3 π 2 n e ) 1 / 3 is the electron Fermi wave number. The Fermi energy is then given by
γ e = 1 + x r 2 .
As shown in Figure 1 and Figure 2 for two representative values γ e = 10 (shallow region of the outer crust) and γ e = 50 (bottom of the outer crust), respectively, the expansions (14) and (16) are surprisingly precise throughout the outer crust. In the limit of vanishingly small magnetic field ( ν max + ), Equations (14) and (16) converge toward the exact results, (19) and (20), respectively. Although the errors increase with the magnetic field as expected, they remain very small in the intermediate regime ν max 1 for which the field is no longer weakly quantizing. When electrons start to be all confined to the lowest Landau–Rabi level, i.e., when ν max = 0 , the error on n ¯ amounts to 0.1% only. The approximate formula for the pressure is found to be more reliable, with errors not exceeding 0.02% and fluctuating.
The expansions (14) and (16) can thus be confidently applied for arbitrary magnetic field strengths, from B = 0 ( ν max + ) up to the threshold magnetic field B = ( γ e 2 1 ) / 2 at the onset of the strongly quantizing regime ( ν max = 0 ) discussed in the next subsection.

2.3. Strongly Quantizing Magnetic Field

The magnetic field is strongly quantizing if all electrons lie in the lowest Landau-Rabi level. Setting ν max = 0 in Equations (5) and (8) leads to the following expression for the pressure and mean baryon number density:
P = B m e c 2 4 π 2 λ e 3 [ γ e γ e 2 1 ln γ e 2 1 + γ e + C M α 3 16 B Z 2 3 π 1 / 3 γ e 2 1 2 / 3 ] ,
n ¯ = B 2 π 2 λ e 3 A Z γ e 2 1 .
These expressions hold if
γ e < 1 + 2 B ,
or equivalently if
n e < B 3 / 2 2 π 2 λ e 3 ,
using Equations (7) and (8).

3. Initial Composition of Magnetar Crusts

Assuming the crust is initially in a full thermodynamic equilibrium in the presence of some magnetic field, the composition is found by minimizing the Gibbs free energy per nucleon, which coincides with the baryon chemical potential (see, e.g., Appendix A in [32]):
μ = E + P n ¯ = M ( A , Z , B ) c 2 A + Z A m e c 2 γ e 1 + 4 3 C M 4 π 3 1 / 3 α λ e n e 1 / 3 Z 2 / 3 .
This minimization can be performed very efficiently following the iterative approach proposed in [33,34]. Freely available computer codes in the limiting cases ν max + and ν max = 0 can be found in [35,36].

3.1. Interface between Adjacent Crustal Layers

The pressure P 1 2 associated with the transition from a crustal layer made of nuclei ( A 1 , Z 1 ) to a denser layer made of nuclei ( A 2 , Z 2 ) is determined by the equilibrium condition
μ ( A 1 , Z 1 , P 1 2 ) = μ ( A 2 , Z 2 , P 1 2 ) ,
which can be approximately written as [33]
γ e + C M α λ e F ( Z 1 , A 1 ; Z 2 , A 2 ) 4 π n e 3 1 / 3 = γ e 1 2 ,
F ( Z 1 , A 1 ; Z 2 , A 2 ) 4 3 Z 1 2 / 3 Z 1 A 1 1 3 Z 1 2 / 3 Z 2 A 2 Z 2 2 / 3 Z 2 A 2 Z 1 A 1 Z 2 A 2 1 ,
γ e 1 2 M ( A 2 , Z 2 , B ) A 2 m e M ( A 1 , Z 1 , B ) A 1 m e Z 1 A 1 Z 2 A 2 1 + 1 .
The bottom of the outer crust is defined by the depth at which neutrons start to drip out of nuclei. The corresponding electron Fermi energy obeys an equation similar to Equation (28), the function F ( Z 1 , A 1 ; Z 2 , A 2 ) being replaced by ( 4 / 3 ) Z 2 / 3 and γ e 1 2 by [24]
γ e drip M ( A , Z , B ) c 2 + A m n c 2 Z m e c 2 + 1 ,
where m n is the neutron mass.
The threshold condition (28) takes formally the same form with and without magnetic fields. However, the solutions do depend on B through the relation between γ e and n e , and potentially also through the ion masses.

3.2. No Magnetic Field

In the absence of magnetic field B = 0 , the solution of Equation (28) reads [33]
x r = γ e 1 2 1 1 F ˜ ( Z 1 , A 1 ; Z 2 , A 2 ) 2 / ( γ e 1 2 ) 2 F ˜ ( Z 1 , A 1 ; Z 2 , A 2 ) × 1 F ˜ ( Z 1 , A 1 ; Z 2 , A 2 ) 2 1 ,
with
F ˜ ( Z 1 , A 1 ; Z 2 , A 2 ) 4 9 π 1 / 3 C M α F ( Z 1 , A 1 ; Z 2 , A 2 ) .
This solution exists only if F ˜ ( Z 1 , A 1 ; Z 2 , A 2 ) 1 .
The mean baryon number density n ¯ 1 max up to which nuclei ( A 1 , Z 1 ) are present and the transition pressure P 1 2 are then given by Equations (19) and (20) respectively.

3.3. Strongly Quantizing Magnetic Field

The solution of Equation (28) was also found in the limit of a strongly quantizing magnetic field [34]. Introducing
F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) 1 3 C M α F ( Z 1 , A 1 ; Z 2 , A 2 ) 2 B 3 π 1 / 3 ,
υ γ e 1 2 2 | F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) | 3 / 2 ,
the electron Fermi energy at the crustal interface is given by the following formulas:
  • γ e 1 2 > 0 and F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) > 0
    γ e = 8 F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) 3 / 2 sinh 3 1 3 arcsinh υ ,
  • γ e 1 2 > 0 and F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) < 0
    γ e = 8 | F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) | 3 / 2 cosh 3 1 3 arccosh υ if υ 1 , 8 | F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) | 3 / 2 cos 3 1 3 arccos υ if 0 υ < 1 ,
  • γ e 1 2 < 0 and F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) < 0
    γ e = 8 | F ¯ ( Z 1 , A 1 ; Z 2 , A 2 ; B ) | 3 / 2 cos 3 θ k if 1 < υ 0 ,
    θ k 1 3 arccos υ + 2 π k 3 and k = 0 , 2 .
The mean baryon number density n ¯ 1 max up to which nuclei ( A 1 , Z 1 ) are present and the transition pressure P 1 2 are then given by Equations (23) and (22), respectively. In the case γ e 1 2 < 0 , the physically admissible solution among k = 0 and k = 2 is the one yielding the lowest transition pressure P 1 2 satisfying the conditions γ e > 1 and n ¯ 2 min n ¯ 1 max .
Let us recall that these solutions are only valid under the assumption ν max = 0 , which translates into a lower bound for the magnetic field B B 1 2 . To find B 1 2 , we substitute Equations (24) and (25) in Equation (28). This leads to
B 1 2 = ( γ e 1 2 ) 2 2 1 + C M α ( 3 π ) 1 / 3 F ( Z 1 , A 1 ; Z 2 , A 2 ) 2 .
We have assumed B 1 2 1 so that γ e 2 B 1 2 .

3.4. Intermediate Magnetic Fields

Approximate analytical solutions can also be found in the intermediate regime. Remarking that the magnetic field enters explicitly in Equation (28) only through the small electrostatic correction, the threshold electron Fermi energy γ e is still approximately given by the solution in the absence of magnetic fields, Equations (21) and (32). However, the density and the pressure are now given by Equations (14) and (18), respectively. As shown in Section 2.2, these expansions in the weakly quantizing limit ν max 1 (including the absence of magnetic field as a limiting case) actually remain very precise for ν max 1 and even at the onset of the strongly quantizing regime ν max = 0 . Combining the solutions thus obtained with those presented in Section 3.3, the full range of possible initial magnetic field strengths can be treated analytically.

4. Magnetic Field Decay and Electron Captures

4.1. Onset of Electron Captures

The initial magnetic field decays on a very long time scale, say of the typical order of millions of years [37]. The compression of the crust thus occurs very slowly. When the pressure of a matter element reaches some value P β ( A , Z , B ) , the capture of an electron by nuclei ( A , Z ) (in their ground state) opens. The daughter nuclei may be in an excited state.
The onset of electron captures by nuclei ( A , Z ) is formally determined by the same condition irrespective of the magnetic field strength by requiring the constancy of the Gibbs free energy per nucleon at fixed temperature and pressure [32]. The threshold electron Fermi energy is found to the first order in the fine-structure constant α from the condition:
γ e + C M 4 π n e 3 1 / 3 α λ e F ( Z ) = γ e β ,
F ( Z ) Z 5 / 3 ( Z 1 ) 5 / 3 + 1 3 Z 2 / 3 ,
γ e β Q EC ( A , Z , B ) m e c 2 + 1 ,
where we have introduced the Q-value (in vacuum) associated with electron capture by nuclei ( A , Z ):
Q EC ( A , Z , B ) = M ( A , Z , B ) c 2 M ( A , Z 1 , B ) c 2 E ex ( A , Z 1 , B ) .
These Q-values can be obtained from the tabulated Q-values of β decay by the following relation:
Q EC ( A , Z , B ) = Q β ( A , Z 1 , B ) E ex ( A , Z 1 , B ) .
Here, E ex ( A , Z 1 , B ) denotes the excitation energy of the daughter nucleus. Transitions to the ground state can be considered by setting E ex ( A , Z 1 , B ) = 0 .

4.2. No Magnetic Field

In the absence of magnetic fields, the threshold condition (41) can be solved exactly [38]:
γ e = 1 + ( x r β ) 2 ,
x r β = γ e β 1 1 F ˜ ( Z ) 2 / ( γ e β ) 2 F ˜ ( Z ) 1 F ˜ ( Z ) 2 1 ,
F ˜ ( Z ) C M 4 9 π 1 / 3 α F ( Z ) .
The pressure P β ( A , Z ) and the density n ¯ β ( A , Z ) at the onset of electron captures are then given by Equations (20) and (19), respectively.

4.3. Intermediate Magnetic Field

For ν max > 0 , the density and the pressure are obtained by substituting the solution (46) for the threshold electron Fermi energy in the absence of magnetic field in Equations (14) and (18), respectively.

4.4. Strongly Quantizing Magnetic Field

In the strongly quantizing regime ( ν max = 0 ), Equation (41) can be solved exactly from the general analytical solutions given in Section 3.3. Introducing
F ¯ ( Z , B ) 1 3 C M α F ( Z ) 2 B 3 π 1 / 3 < 0 ,
υ γ e β 2 | F ¯ ( Z , B ) | 3 / 2 ,
remarking that γ e β > 1 , the solutions are given by the following formulas:
γ e = 8 | F ¯ ( Z , B ) | 3 / 2 cosh 3 1 3 arccosh υ if υ 1 , 8 | F ¯ ( Z , B ) | 3 / 2 cos 3 1 3 arccos υ if 0 υ < 1 .
The threshold pressure and density are respectively given by Equations (22) and (23).
Let us recall that this solution is only valid under the assumption ν max = 0 , which translates into a lower bound for the magnetic field B B β . To find B β , we substitute Equations (24) and (25) in Equation (41). This leads to
B β = ( γ e β ) 2 2 1 + C M α ( 3 π ) 1 / 3 F ( Z ) 2 .
We have assumed B β 1 so that γ e 2 B β .

4.5. Heat Released

The first electron capture does not release any significant heat since it essentially proceeds in quasiequilibrium. However, the daughter nuclei (possibly in some excited state) are generally unstable and capture a second electron off-equilibrium thus depositing some heat at the same pressure P β ( A , Z , B ) . Ignoring the fraction of energy carried away by neutrinos, the maximum amount of heat per nucleus ( A , Z ) is given by
Q ( A , Z , B ) = μ ( A , Z , B ) μ ( A , Z 2 , B ) .
It is to be understood that the baryon chemical potentials must be evaluated at the same pressure. Expressing the electron Fermi energy associated with nuclei ( A , Z 2 ) as γ e + δ γ e with γ e given by the solution of Equation (41), and expanding the pressure to the first order in δ γ e leads to
δ γ e 1 3 C M α λ e 4 π n e 3 1 / 3 Z 2 / 3 ( Z 2 ) 2 / 3 .
We have made use of the Gibbs–Duhem relation d P e = n e m e c 2 d γ e and we have neglected terms of order α δ γ e . Substituting Equation (54) in Equation (53) and eliminating γ e using Equation (41) lead to the following expression for the heat released per nucleus (keeping as before first-order terms):
Q ( A , Z , B ) Q ( 0 ) ( A , Z , B ) C M α m e c 2 λ e 4 π A 3 Z n ¯ 1 / 3 Z 5 / 3 + ( Z 2 ) 5 / 3 2 ( Z 1 ) 5 / 3 ,
where the zeroth-order term is determined by nuclear data alone
Q ( 0 ) ( A , Z , B ) 2 M ( A , Z 1 , B ) c 2 + 2 E ex ( A , Z 1 , B ) M ( A , Z 2 , B ) c 2 M ( A , Z , B ) c 2 .
Apart from the small electrostatic correction (the term proportional to the structure constant C M ), the maximum heat released by electron captures is thus independent of whether the crust is solid or liquid. Unless B 10 17 G [39], the structure of nuclei remains essentially unchanged in the presence of a magnetic field so that Q ( A , Z , B ) Q ( A , Z ) .
To estimate the heat in Equation (53), we implicitly assumed μ ( A , Z 2 , P β , B ) < μ ( A , Z , P β , B ) , which generally holds for even A nuclei, but not necessarily for odd A nuclei. In the latter case, we typically have Q β ( A , Z 1 , B ) < Q β ( A , Z 2 , B ) . Using Equation (43), this implies that γ e β ( A , Z ) < γ e β ( A , Z 1 ) . 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-equilibrium at pressure P β ( A , Z 1 , B ) > P β ( A , Z , B ) .

4.6. Neutron Delayed Emission

As discussed in [24,40], the first electron capture by the nucleus ( A , Z ) may be accompanied by the emission of Δ N > 0 neutrons. The corresponding pressure P β n and baryon density n β n are obtained from similar expressions as for electron captures except that the threshold electron Fermi energy γ e β is now replaced by
γ e β n = M ( A Δ N , Z 1 ) c 2 + E ex ( A Δ N , Z 1 ) M ( A , Z ) c 2 + Δ N m n c 2 m e c 2 + 1 .
Neutron emission will thus occur whenever γ e β n ( A , Z ) < γ e β ( A , Z ) .

5. Results and Discussions

5.1. Initial Composition of the Outer Crust

The initial composition of the outer crust of a magnetar was determined in [18] but only for a few selected magnetic field strengths, namely B = 1000 , 2000, and 3000. We have extended the calculations to the whole range of magnetic field strengths ranging from B = 0 to B = 10 17 G. To this end, we have used the experimental atomic masses from the 2016 Atomic Mass Evaluation [41] supplemented with the same microscopic atomic mass table HFB-24 [14] from the BRUSLIB database1 [42]. The functional BSk24 underlying the model HFB-24 was also adopted to calculate the equation of state of the inner crust of a magnetar [18]. This same functional was also applied to construct a unified equation of state for unmagnetized neutron stars [15,16,17], and to calculate superfluid properties [43]. Results are publicly available on CompOSE2. This equation of state is consistent with the constraints inferred from analyses of the gravitational-wave signal from the binary neutron-star merger GW170817 and of its electromagnetic counterpart [44]. As shown in [18], the magnetic field has a negligible impact on the equation of the state of the inner crust and core of magnetars unless it exceeds about 10 17 G.
Depending on the strength of the magnetic field when the magnetar was born, different nuclides are expected to be produced in the outer crust. Changes in the composition compared to that obtained in [15] in the absence of the magnetic field are summarized in Table 1 and Table 2. The mean baryon number densities and the pressure at the interface between adjacent crustal layers can be calculated for any magnetic field using the analytical formulas given in Section 3. The only nuclear inputs are embedded in the parameter γ e 1 2 defined by Equation (30). Values for this parameter are indicated in Table 3 for all possible transitions.

5.2. Heating

We have estimated the heat released by electron captures and their location using the experimental atomic masses and the Q β values (including the recommended ones) from the 2016 Atomic Mass Evaluation [41] supplemented with the atomic mass model HFB-24 [14]. We have taken excitation energies from the Nuclear Data section of the International Atomic Energy Agency website3 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 ± (excluding transitions from J = 0 to J = 0 ).
The threshold density and pressure for the onset of each electron capture, as well as the amount of heat deposited, can be calculated for any magnetic field strength from the analytical formulas presented in Section 4 using the parameters indicated in Table 4 for ground-state to ground-state transitions, in Table 5 for ground-state to excited state transitions, and in Table 6 for transitions involving light elements that could have been accreted from the interstellar medium. Full numerical results are freely available in [45].
To assess the reliability of our analytical treatment, we have numerically solved the exact threshold conditions μ ( A , Z , n e ) = μ ( A , Z 1 , n e 1 ) and P ( Z , n e ) = P ( Z 1 , n e 1 ) = P ( Z 2 , n e 2 ) without any approximation, i.e., using Equations (5) and (8), to determine the exact values for the threshold pressure P β and baryon number density n ¯ β . The heat deposited is then calculated as Q = μ ( A , Z , n e ) μ ( A , Z 2 , n e 2 ) . We have compared these exact results with the approximate analytical formulas. As an example, we focus on the electron capture by 56 Fe, considering the ground-state to ground-state transition. As shown in Figure 3, the quantum oscillations of the threshold density are correctly reproduced. The errors are found to be the largest for specific values of the magnetic field strength corresponding to exact fillings of Landau–Rabi levels and amount to a few percents, but drop by about an order of magnitude in the strongly quantizing regime, as shown in Figure 4. As previously discussed in Section 2.2, the expansion in the weakly quantizing regime is more reliable for the pressure than for the density. Indeed, the overall errors on the threshold pressure are significantly smaller, as can be observed in Figure 5. In the strongly quantizing regime, the errors on P β are of the same order as those on n ¯ β and are displayed in Figure 6. The heat released, plotted in Figure 7 and Figure 8, also exhibits quantum oscillations, though the amplitude is very small. To a first approximation, the heat is therefore essentially given by that in the absence of magnetic fields, as anticipated. The analytical formula for Q is found to be even more precise than formulas for the density and pressure. To check that these error estimates are not specific to the reaction considered, we have also analyzed the electron capture by 122 Zr, which is present in much deeper layers of the outer crust for all the magnetic field strengths. As can be observed in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, the analytical formulas remain very precise in this other case. We have examined other reactions and reached similar conclusions.

6. Conclusions

We have derived accurate analytical formulas (with typical errors below 1%) for calculating the threshold density n ¯ β and pressure P β for the onset of electron captures by nuclei in the shallow layers of magnetar crusts, as well as the maximum amount of heat released Q taking into account the Landau–Rabi quantization of electron motion induced by the magnetic field. We have also obtained formulas for determining the initial constitution of the outer crust. These formulas are applicable over the whole range of magnetic fields encountered in neutron stars, from the weakly quantizing regime to the most extreme situation in which the electrons all lie in the lowest Landau–Rabi level.
Using experimental nuclear data supplemented with predictions from the atomic mass model HFB-24, we have calculated all the necessary nuclear parameters to calculate the shallow heating for any given magnetic field considering both ground-state to ground-state and ground-state to excited-state transitions. Full numerical results can be found in [45]. Together with the results for the equation of state and superfluid properties published in [15,16,17,18,43], they provide consistent microscopic inputs for modelling the magneto-thermal evolution of neutron stars.

Author Contributions

Conceptualization, N.C.; methodology, N.C.; software, N.C. and A.F.F.; validation, N.C. and A.F.F.; formal analysis, N.C.; investigation, N.C.; writing—original draft preparation, N.C.; writing—review and editing, N.C. and A.F.F.; visualization, N.C.; supervision, N.C.; project administration, N.C. All authors have 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. 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 2016 Atomic Mass Evaluation [41], the BRUSLIB database (http://www.astro.ulb.ac.be/bruslib/, accessed on 9 June 2022, see [42]), and the Nuclear Data section of the International Atomic Energy Agency website (https://www-nds.iaea.org/relnsd/NdsEnsdf/QueryForm.html, accessed on 9 June 2022). The results presented in this study are openly available on the Zenodo repository [45].

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
http://www.astro.ulb.ac.be/bruslib/, accessed on 9 June 2022.
2
https://compose.obspm.fr, accessed on 9 June 2022.
3

References

  1. 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. [Google Scholar] [CrossRef]
  2. Duncan, R.C.; Thompson, C. Formation of Very Strongly Magnetized Neutron Stars: Implications for Gamma-Ray Bursts. Astrophys. J. Lett. 1992, 392, L9. [Google Scholar] [CrossRef]
  3. Olausen, S.A.; Kaspi, V.M. The McGill Magnetar Catalog. Astrophys. J. Suppl. 2014, 212, 6. [Google Scholar] [CrossRef]
  4. 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. [Google Scholar] [CrossRef]
  5. Beloborodov, A.M.; Li, X. Magnetar Heating. Astrophys. J. 2016, 833, 261. [Google Scholar] [CrossRef]
  6. 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. [Google Scholar] [CrossRef]
  7. Fantina, A.F.; De Ridder, S.; Chamel, N.; Gulminelli, F. Crystallization of the outer crust of a non-accreting neutron star. Astron. Astrophys. 2020, 633, A149. [Google Scholar] [CrossRef]
  8. Carreau, T.; Gulminelli, F.; Chamel, N.; Fantina, A.F.; Pearson, J.M. Crystallization of the inner crust of a neutron star and the influence of shell effects. Astron. Astrophys. 2020, 635, A84. [Google Scholar] [CrossRef]
  9. 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. [Google Scholar] [CrossRef]
  10. 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. [Google Scholar] [CrossRef]
  11. 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. [Google Scholar] [CrossRef]
  12. Haensel, P.; Zdunik, J.L. Non-equilibrium processes in the crust of an accreting neutron star. Astron. Astrophys. 1990, 227, 431–436. [Google Scholar]
  13. Chamel, N.; Fantina, A.F.; Suleiman, L.; Zdunik, J.L.; Haensel, P. Heating in Magnetar Crusts from Electron Captures. Universe 2021, 7, 193. [Google Scholar] [CrossRef]
  14. 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. [Google Scholar] [CrossRef]
  15. 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. [Google Scholar] [CrossRef]
  16. Pearson, J.M.; Chamel, N.; Potekhin, A.Y. Unified equations of state for cold nonaccreting neutron stars with Brussels-Montreal functionals. II. Pasta phases in semiclassical approximation. Phys. Rev. C 2020, 101, 015802. [Google Scholar] [CrossRef]
  17. Pearson, J.M.; Chamel, N. Unified equations of state for cold nonaccreting neutron stars with Brussels-Montreal functionals. III. Inclusion of microscopic corrections to pasta phases. Phys. Rev. C 2022, 105, 015803. [Google Scholar] [CrossRef]
  18. 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. [Google Scholar] [CrossRef]
  19. Rabi, I.I. Das freie Elektron im homogenen Magnetfeld nach der Diracschen Theorie. Z. Fur. Phys. 1928, 49, 507–511. [Google Scholar] [CrossRef]
  20. Landau, L. Diamagnetismus der Metalle. Z. Fur Phys. 1930, 64, 629–637. [Google Scholar] [CrossRef]
  21. 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. [Google Scholar] [CrossRef]
  22. An, H.; Kaspi, V.M.; Beloborodov, A.M.; Kouveliotou, C.; Archibald, R.F.; Boggs, S.E.; Christensen, F.E.; Craig, W.W.; Gotthelf, E.V.; Grefenstette, B.W.; et al. NuSTAR Observations of X-Ray Bursts from the Magnetar 1E 1048.1-5937. Astrophys. J. 2014, 790, 60. [Google Scholar] [CrossRef]
  23. Uryū, K.; Yoshida, S.; Gourgoulhon, E.; Markakis, C.; Fujisawa, K.; Tsokaros, A.; Taniguchi, K.; Eriguchi, Y. New code for equilibriums and quasiequilibrium initial data of compact objects. IV. Rotating relativistic stars with mixed poloidal and toroidal magnetic fields. Phys. Rev. D 2019, 100, 123019. [Google Scholar] [CrossRef]
  24. 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. [Google Scholar] [CrossRef]
  25. Haensel, P.; Potekhin, A.Y.; Yakovlev, D.G. Neutron Stars. 1. Equation of State and Structure; Springer: New York, NY, USA, 2007. [Google Scholar]
  26. Van Vleck, J.H. The Theory of Electric and Magnetic Susceptibilities; Oxford University Press: London, UK, 1932. [Google Scholar]
  27. Baiko, D.A. Coulomb crystals in the magnetic field. Phys. Rev. E 2009, 80, 046405. [Google Scholar] [CrossRef]
  28. Baiko, D.A.; Potekhin, A.Y.; Yakovlev, D.G. Thermodynamic functions of harmonic Coulomb crystals. Phys. Rev. E 2001, 64, 057402. [Google Scholar] [CrossRef]
  29. Potekhin, A.Y.; Chabrier, G. Equation of state for magnetized Coulomb plasmas. Astron. Astrophys. 2013, 550, A43. [Google Scholar] [CrossRef]
  30. Salpeter, E.E. Electrons Screening and Thermonuclear Reactions. Aust. J. Phys. 1954, 7, 373. [Google Scholar] [CrossRef]
  31. Dib, C.O.; Espinosa, O. The magnetized electron gas in terms of Hurwitz zeta functions. Nucl. Phys. B 2001, 612, 492–518. [Google Scholar] [CrossRef]
  32. Chamel, N.; Fantina, A.F. Electron capture instability in magnetic and nonmagnetic white dwarfs. Phys. Rev. D 2015, 92, 023008. [Google Scholar] [CrossRef]
  33. Chamel, N. Analytical determination of the structure of the outer crust of a cold nonaccreted neutron star. Phys. Rev. C 2020, 101, 032801. [Google Scholar] [CrossRef]
  34. Chamel, N.; Stoyanov, Z.K. Analytical determination of the structure of the outer crust of a cold nonaccreted neutron star: Extension to strongly quantizing magnetic fields. Phys. Rev. C 2020, 101, 065802. [Google Scholar] [CrossRef]
  35. Chamel, N. Equilibrium Structure of the Outer Crust of a Cold Nonaccreted Neutron Star; Zenodo: Geneva, Switzerland, 2020. [Google Scholar] [CrossRef]
  36. Chamel, N.; Stoyanov, Z. Equilibrium Structure of the Outer Crust of a Magnetar; Zenodo: Geneva, Switzerland, 2020. [Google Scholar] [CrossRef]
  37. Pons, J.A.; Viganò, D. Magnetic, thermal and rotational evolution of isolated neutron stars. Living Rev. Comput. Astrophys. 2019, 5, 3. [Google Scholar] [CrossRef]
  38. Chamel, N.; Fantina, A.F. Binary and ternary ionic compounds in the outer crust of a cold nonaccreting neutron star. Phys. Rev. C 2016, 94, 065802. [Google Scholar] [CrossRef]
  39. 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. [Google Scholar] [CrossRef]
  40. 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. [Google Scholar] [CrossRef]
  41. 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. [Google Scholar] [CrossRef]
  42. 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. [Google Scholar] [CrossRef]
  43. Allard, V.; Chamel, N. 1S0 Pairing Gaps, Chemical Potentials and Entrainment Matrix in Superfluid Neutron-Star Cores for the Brussels-Montreal Functionals. Universe 2021, 7, 470. [Google Scholar] [CrossRef]
  44. Perot, L.; Chamel, N.; Sourie, A. Role of the symmetry energy and the neutron-matter stiffness on the tidal deformability of a neutron star with unified equations of state. Phys. Rev. C 2019, 100, 035801. [Google Scholar] [CrossRef]
  45. Chamel, N.; Fantina, A.F. Onset of Electron Captures and Shallow Heating in Magnetars [Data Set]; Zenodo: Geneva, Switzerland, 2022. [Google Scholar] [CrossRef]
Figure 1. Relative errors (in %) of the approximate analytical expansions (14) and (16) for the baryon number density n ¯ (dashed line) and electron pressure P e (solid line), respectively, as a function of the magnetic field strength B = B / B rel with γ e = 10 . The errors are obtained by taking the difference between the approximate and exact results and dividing by the exact result. All electrons are confined to the lowest Landau–Rabi level at B = ( γ e 2 1 ) / 2 = 49.5 .
Figure 1. Relative errors (in %) of the approximate analytical expansions (14) and (16) for the baryon number density n ¯ (dashed line) and electron pressure P e (solid line), respectively, as a function of the magnetic field strength B = B / B rel with γ e = 10 . The errors are obtained by taking the difference between the approximate and exact results and dividing by the exact result. All electrons are confined to the lowest Landau–Rabi level at B = ( γ e 2 1 ) / 2 = 49.5 .
Universe 08 00328 g001
Figure 2. Same as Figure 1 for γ e = 50 . All electrons are confined to the lowest Landau–Rabi level at B = ( γ e 2 1 ) / 2 = 1249.5 .
Figure 2. Same as Figure 1 for γ e = 50 . All electrons are confined to the lowest Landau–Rabi level at B = ( γ e 2 1 ) / 2 = 1249.5 .
Universe 08 00328 g002
Figure 3. Top panel: Exact threshold density (in cgs units) for the onset of electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 3. Top panel: Exact threshold density (in cgs units) for the onset of electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g003
Figure 4. Same as Figure 3 in the strongly quantizing regime.
Figure 4. Same as Figure 3 in the strongly quantizing regime.
Universe 08 00328 g004
Figure 5. Top panel: Exact threshold pressure (in cgs units) for the onset of electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 5. Top panel: Exact threshold pressure (in cgs units) for the onset of electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g005
Figure 6. Same as Figure 5 in the strongly quantizing regime.
Figure 6. Same as Figure 5 in the strongly quantizing regime.
Universe 08 00328 g006
Figure 7. Top panel: Exact heat released per nucleon from electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 7. Top panel: Exact heat released per nucleon from electron captures by 56 Fe as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g007
Figure 8. Same as Figure 7 in the strongly quantizing regime.
Figure 8. Same as Figure 7 in the strongly quantizing regime.
Universe 08 00328 g008
Figure 9. Top panel: Exact threshold density (in cgs units) for the onset of electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 9. Top panel: Exact threshold density (in cgs units) for the onset of electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g009
Figure 10. Same as Figure 9 in the strongly quantizing regime.
Figure 10. Same as Figure 9 in the strongly quantizing regime.
Universe 08 00328 g010
Figure 11. Top panel: Exact threshold pressure (in cgs units) for the onset of electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 11. Top panel: Exact threshold pressure (in cgs units) for the onset of electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g011
Figure 12. Same as Figure 11 in the strongly quantizing regime.
Figure 12. Same as Figure 11 in the strongly quantizing regime.
Universe 08 00328 g012
Figure 13. Top panel: Exact heat released per nucleon from electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Figure 13. Top panel: Exact heat released per nucleon from electron captures by 122 Zr as a function of the magnetic field strength B = B / B rel up to the onset of the strongly quantizing regime. Bottom panel: Relative error (in %) of the approximate analytical expression.
Universe 08 00328 g013
Figure 14. Same as Figure 13 in the strongly quantizing regime.
Figure 14. Same as Figure 13 in the strongly quantizing regime.
Universe 08 00328 g014
Table 1. Magnetic field strength B = B / B rel for the appearance (+) or the disappearance (−) of a nuclide in the outer crust of a magnetar.
Table 1. Magnetic field strength B = B / B rel for the appearance (+) or the disappearance (−) of a nuclide in the outer crust of a magnetar.
Nuclide B
66 Ni(−)67
88 Sr(+)858
126 Ru(+)1023
80 Ni(−)1072
128 Pd(+)1249
78 Ni(−)1416
64 Ni(−)1669
124 Zr(+)1872
121 Y(−)1907
132 Sn(+)1986
80 Ni(−)2087
Table 2. Sequence of equilibrium nuclides with increasing depth (from top to bottom) in the outer crust of a magnetar for different magnetic field strengths. The first row indicates values of B = B / B rel associated with a change of composition. Results are valid up to B = 2087 .
Table 2. Sequence of equilibrium nuclides with increasing depth (from top to bottom) in the outer crust of a magnetar for different magnetic field strengths. The first row indicates values of B = B / B rel associated with a change of composition. Results are valid up to B = 2087 .
6785810231072124914161669187219071986
56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe 56 Fe
62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni 62 Ni
64 Ni 64 Ni 64 Ni 64 Ni 64 Ni 64 Ni 64 Ni
66 Ni
88 Sr 88 Sr 88 Sr 88 Sr 88 Sr 88 Sr 88 Sr 88 Sr 88 Sr
86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr 86 Kr
84 Se 84 Se 84 Se 84 Se 84 Se 84 Se 84 Se 84 Se 84 Se 84 Se 84 Se
82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge 82 Ge
132 Sn
80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn 80 Zn
78 Ni 78 Ni 78 Ni 78 Ni 78 Ni 78 Ni
80 Ni 80 Ni 80 Ni 80 Ni
128 Pd 128 Pd 128 Pd 128 Pd 128 Pd 128 Pd
126 Ru 126 Ru 126 Ru 126 Ru 126 Ru 126 Ru 126 Ru 126 Ru
124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo 124 Mo
122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr 122 Zr
124 Zr 124 Zr 124 Zr
121 Y 121 Y 121 Y 121 Y 121 Y 121 Y 121 Y 121 Y 121 Y
120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr 120 Sr
122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr 122 Sr
124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr 124 Sr
Table 3. Values of the nuclear parameter γ e 1 2 from which the pressure and the densities at the boundary between adjacent layers in the outer crust of a magnetar can be calculated.
Table 3. Values of the nuclear parameter γ e 1 2 from which the pressure and the densities at the boundary between adjacent layers in the outer crust of a magnetar can be calculated.
γ e 1 2 Interface
1.8908 56 Fe 62 Ni
4.8972 62 Ni 64 Ni
8.6863 62 Ni 88 Sr
8.1312 64 Ni 66 Ni
9.3317 64 Ni 86 Kr
18.098 64 Ni 88 Sr
12.155 66 Ni 86 Kr
10.044 86 Kr 84 Se
5.5622 88 Sr 86 Kr
15.330 84 Se 82 Ge
20.519 82 Ge 80 Zn
38.310 82 Ge 132 Sn
25.926 80 Zn 78 Ni
37.083 80 Zn 128 Pd
33.289 132 Sn 80 Zn
32.978 78 Ni 80 Ni
409.21 78 Ni 128 Pd
48.055 78 Ni 126 Ru
45.122 80 Ni 124 Mo
218.54 80 Ni 126 Ru
30.039 128 Pd 126 Ru
32.010 126 Ru 124 Mo
37.441 124 Mo 122 Zr
40.352 122 Zr 121 Y
42.293 122 Zr 124 Zr
1.8541 124 Zr 121 Y
39.031 124 Zr 120 Sr
40.786 121 Y 120 Sr
44.857 120 Sr 122 Sr
47.747 122 Sr 124 Sr
Table 4. Values of the nuclear parameters γ e β and Q ( 0 ) (considering ground-state to ground-state transitions) from which the threshold density and pressure for the onset of electron captures, as well as the heat released can be calculated. The symbol (★) is used to distinguish reactions for which theoretical atomic masses were needed.
Table 4. Values of the nuclear parameters γ e β and Q ( 0 ) (considering ground-state to ground-state transitions) from which the threshold density and pressure for the onset of electron captures, as well as the heat released can be calculated. The symbol (★) is used to distinguish reactions for which theoretical atomic masses were needed.
γ e β Reaction Q ( 0 ) (MeV)
8.232 56 Fe 56 Cr 2 e + 2 ν e 2.069
18.867 56 Cr 56 Ti 2 e + 2 ν e 2.295
29.313 56 Ti 56 Ca 2 e + 2 ν e 3.514
43.710 56 Ca 56 Ar 2 e + 2 ν e ( ) 2.045
11.415 62 Ni 62 Fe 2 e + 2 ν e 2.776
21.262 62 Fe 62 Cr 2 e + 2 ν e 2.725
31.174 62 Cr 62 Ti 2 e + 2 ν e 2.442
41.470 62 Ti 62 Ca 2 e + 2 ν e ( ) 2.490
15.299 64 Ni 64 Fe 2 e + 2 ν e 2.484
24.445 64 Fe 64 Cr 2 e + 2 ν e 2.471
34.581 64 Cr 64 Ti 2 e + 2 ν e 1.865
19.782 66 Ni 66 Fe 2 e + 2 ν e 3.257
27.062 66 Fe 66 Cr 2 e + 2 ν e 1.287
38.397 66 Cr 66 Ti 2 e + 2 ν e ( ) 3.540
15.938 86 Kr 86 Se 2 e + 2 ν e 2.504
23.585 86 Se 86 Ge 2 e + 2 ν e 1.979
30.980 86 Ge 86 Zn 2 e + 2 ν e ( ) 2.310
38.926 86 Zn 86 Ni 2 e + 2 ν e ( ) 2.320
20.754 84 Se 84 Ge 2 e + 2 ν e 2.389
28.517 84 Ge 84 Zn 2 e + 2 ν e 1.903
36.362 84 Zn 84 Ni 2 e + 2 ν e ( ) 2.390
25.431 82 Ge 82 Zn 2 e + 2 ν e 1.868
34.256 82 Zn 82 Ni 2 e + 2 ν e ( ) 2.154
44.640 82 Ni 82 Fe 2 e + 2 ν e ( ) 2.080
31.233 80 Zn 80 Ni 2 e + 2 ν e 1.879
39.435 78 Ni 78 Fe 2 e + 2 ν e ( ) 2.070
41.842 124 Mo 124 Zr 2 e + 2 ν e ( ) 2.920
42.722 122 Zr 122 Sr 2 e + 2 ν e ( ) 0.790
11.397 88 Sr 88 Kr 2 e + 2 ν e 2.395
18.564 88 Kr 88 Se 2 e + 2 ν e 2.144
26.761 88 Se 88 Ge 2 e + 2 ν e 2.582
34.601 88 Ge 88 Zn 2 e + 2 ν e ( ) 2.740
36.695 126 Ru 126 Mo 2 e + 2 ν e ( ) 1.860
34.444 128 Pd 128 Ru 2 e + 2 ν e ( ) 2.290
28.662 132 Sn 132 Cd 2 e + 2 ν e 1.987
33.237 132 Cd 132 Pd 2 e + 2 ν e ( ) 2.863
38.084 132 Pd 132 Ru 2 e + 2 ν e ( ) 2.940
Table 5. Same as in Table 4 but considering ground-state to excited state transitions (for which the excitation energy is experimentally known and E exc > 0 ).
Table 5. Same as in Table 4 but considering ground-state to excited state transitions (for which the excitation energy is experimentally known and E exc > 0 ).
γ e β Reaction Q ( 0 ) (MeV)
8.448 56 Fe 56 Cr 2 e + 2 ν e 2.290
12.405 62 Ni 62 Fe 2 e + 2 ν e 3.788
20.726 86 Kr 86 Se 2 e + 2 ν e 7.397
31.260 82 Ge 82 Zn 2 e + 2 ν e 7.825
15.764 88 Sr 88 Kr 2 e + 2 ν e 6.858
22.290 88 Kr 88 Se 2 e + 2 ν e 5.951
31.010 132 Sn 132 Cd 2 e + 2 ν e 4.387
Table 6. Same as in Table 4 but considering reactions involving carbon and oxygen.
Table 6. Same as in Table 4 but considering reactions involving carbon and oxygen.
γ e β Reaction Q ( 0 ) (MeV)
21.393 16 O 16 C 2 e + 2 ν e 2.411
27.163 12 C 12 Be 2 e + 2 ν e 1.661
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chamel, N.; Fantina, A.F. Onset of Electron Captures and Shallow Heating in Magnetars. Universe 2022, 8, 328. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060328

AMA Style

Chamel N, Fantina AF. Onset of Electron Captures and Shallow Heating in Magnetars. Universe. 2022; 8(6):328. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060328

Chicago/Turabian Style

Chamel, Nicolas, and Anthea Francesca Fantina. 2022. "Onset of Electron Captures and Shallow Heating in Magnetars" Universe 8, no. 6: 328. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8060328

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop