Magnetar crust electron capture for 55\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{55}$$\end{document}Co and 56\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{56}$$\end{document}Ni

Based on the relativistic mean-field effective interaction principle and random phase approximation theory in superstrong magnetic fields (SMFs), we present an analysis of the influence of SMFs on the electron Fermi energy, nuclear blinding energy, single-particle level structure and electron capture for 55\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{55}$$\end{document}Co, and 56\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{56}$$\end{document}Ni by the shell-model Monte Carlo method in the magnetar’s crust. The electron capture rates increase by two orders of magnitude due to an increase in the electron Fermi energy and a change in single-particle level structure by SMFs. Then the rates decrease by more than two orders of magnitude due to an increase in the nuclear binding energy and a reduction in the electron Fermi energy by SMFs.

Very recently, Mereghetti et al. [22] reviewed explicitly the magnetars' properties, origin and evolution. Some interesting behaviors such as a wide array of X-ray activity including short bursts, large outbursts, giant flares, quasiperiodic oscillations, enhanced spin-down, glitches and antiglitches are displayed in magnetars and discussed by Kaspi and Beloborodov [23]. The electron Fermi energy, and the a e-mail: liujingjing68@126.com electron fraction in SMFs and the dipole magnetic field and spin-down evolutions were explored in detail by Gao et al. [24][25][26][27][28] and Zhu et al. [29].
A new improved relativistic mean-field effective interaction theory (RMFEIT) with explicit density dependence of the meson-nucleon couplings is used to calculate the energy level for nuclear ground state and properties of excited states [30,31]. RMFEIT has recently gained considerable success in describing various facets of nuclear structure properties. RMFEIT also can give a quantitative description of ground state properties of spherical and deformed nuclei and has been used to discuss the properties of nuclear matter by some authors (e.g., [30][31][32][33]).
Our work differs from previous work [1][2][3]5,6,9]. FFN, and AUFD based their work on the theory of the Brink hypothesis. The Brink hypothesis is a very crude approximation, which assumes that the Gamow-Teller (GT) strength distribution on excited states is the same as that on the ground state. The GT strength distribution is only shifted by the excitation energy of the states. NKK [9] investigated the EC issue by using the quasi-particle random phase approximation and addressed the problem only in the case without SMFs. We analyze the EC in SMFs according to the shell-model Monte Carlo (SMMC) method and random phase approximation (RPA) theory [7,8]. Although Refs. [7,8] investigated in detail the EC process, they have lost sight on the influence of SMFs on the EC. Our discussions also differ from [34], which analyzed the EC process by using the method of the Brink hypothesis only in non-zero temperature crusts of neutron stars.
In this paper, based on the Lai dong model [19][20][21] and the RMFEIT [30][31][32][33], we improve our previous work (i.e., [12,13]) and study the electron Fermi energy, nuclear blinding energy, single-particle level structure and the EC problem in SMFs. In Sect. 2, we describe the theoretical framework and analyze the effect of SMFs on the electron property and nuclear blinding energy in magnetars. The SMMC method, Gamow-Teller response functions, and the EC in SMFs are also discussed. In Sect. 3 we present numerical results and discussion and summarize our conclusions in Sect. 4.

The electron properties in magnetar's crust
The properties of matter are significantly modified by strong magnetic fields, such as the equation of state, the electron energy, the outer crust structure and the composition in neutron stars. Some work presented the motions of charged particle in a SMF with the quantum mechanism in detail (e.g., [35][36][37]). At first, we consider the non-relativistic motion of a charged particle (charge e i and mass m i ) in a uniform magnetic field, which is assumed to be along the z-axis, the circular orbit radius and the cyclotron frequency in the process of particle gyrates are given by r = m i cv ⊥ /|e i |B, ω c = |e i |B/m i c, respectively, here v ⊥ is the velocity perpendicular to the magnetic field. The kinetic energy for the electron (m i → m e , e i → −e) transverse motion is quantized and can be written where n l = 0, 1, 2..., is the Landau level number. The cyclotron energy for an electron is given by where B 12 = B/10 12 G is the magnetic field strength in units of 10 12 G, The total electron energy, which includes the kinetic energy associated with the z-momentum ( p z ) and the spin energy in non-relativistic quantum environment can be written as [35][36][37] where ν = n l + (1 + σ z )/2. σ z = −1, σ z = ±1 are the spin degeneracy for the ground Landau level (n l = 0), and excited levels, respectively. We define a critical magnetic field strength B cr from the relation ofhω c = m e c 2 (i.e. B cr = m 2 e c 3 /eh = 4.414×10 13 G). The transverse motion of the electron becomes relativistic whenhω c ≥ m e c 2 (i.e. B ≥ B cr ) for extremely strong magnetic fields. The energy eigenstates of electrons must obey the relativistic Dirac equation and are given by [35][36][37] E n = c 2 p 2 z + m 2 e c 4 1 + 2ν where the shape of the Landau wavefunction in the relativistic theory is the same as that in the non-relativistic theory due to the fact that the cyclotron radius is independent of the particle mass.
In SMFs the number density of electrons n e is related to the chemical potential U F by [19][20][21] where ρ = (hc/eB) 1/2 = 2.5656 × 10 −10 B 1/2 12 cm is the cyclotron radius (the characteristic size of the wave packet), and g n0 is the spin degeneracy of the Landau level, g 00 = 1 and g n0 = 2 for n ≥ 1, and According to the relation of the usual relativistic energy and momentum from Eq. (4), the interaction energy term, which is proportional to the quantum number ν, cannot exceed the electron chemical potential due to the electron interaction with the magnetic field. Thus the maximum number of Landau levels ν max , related to the highest value of the allowed interaction energy, should be satisfied with E n = U F when ν max , and p z = 0. So we have . However, in the general case (i.e., 0 ≤ ν ≤ ν max ), when the maximum electron momentum is equaled to the Fermi momentum P F for different Landau level value ν, the electron chemical potential from Eq. (4) can be computed as follows: If we define a non-dimensional Fermi momentum x e (ν) = p F /m e c, and a non-dimensional Fermi energy γ e = U F / m e c 2 , the electron density, the electron energy, and the electron pressure can be written as [19][20][21] and , λ e =h/m e c is the electron Compton wavelength.

The nuclear energy in magnetar's crust
The matter in the outer crust of a cold (T = 0 K) magnetar consists of a Coulomb lattice of completely ionized atoms and a uniform Fermi gas of relativistic electrons. The Gibbs free energy per baryon g(A, z, P) at a constant pressure and zero temperature is given by [38] g where ε is the corresponding energy per nucleon, n = A/V is the baryon density in a cell, and V is the volume occupied by a unit cell of the Coulomb lattice. The energy per nucleon ε, which consists of three different contributions from nucleus, electrons, and lattice is given by where the nuclear contribution to the total energy per nucleus is simple and independent of the density and is written as where M(A, z) is the nucleus mass, ε b (A, z) is the corresponding binding energy, and m n and m p are neutron and proton masses, respectively. Based on RMFEIT [32,33], Pena Arteaga et al. [39] discussed the influence of SMFs on the nuclear binding energies. An effective Lagrangian with nucleons and mesons is given by the simple and independent density function [39][40][41] where L N , L m , and L int are the Lagrangian of the free nucleus, the free meson fields and the electromagnetic field generated by protons, and the Lagrangian describing the interactions, respectively. These Lagrangians are represented as and respectively, where ψ is the Dirac spinor. m nu , and m σ , m ω , m ρ , are the nucleon and meson masses, respectively. U (σ ) = (g 2 /3)σ 3 + g 3 /4)σ 4 is the standard form for the nonlinear coupling of the σ meson field. g σ , g ω , g ρ are the coupling constants for the σ, ω, ρ, respectively. e is photon fields which vanishes for neutrons. The coupling of the proton orbital motion with the external magnetic field, and the coupling of proton and neutron intrinsic dipole magnetic moments with the external magnetic field can be expressed as [42] and is the nuclear magneton, κ n = g n /2, κ p = g p /2 − 1 (here g n = −3.8263, g p = 5.5856) are the intrinsic magnetic moments of protons and neutrons, respectively.
The contribution of electrons, which are treated as a degenerate free Fermi gas, is given by The lattice energy per baryon ε l (A, Z , n) can be written as [42] where C bcc = 3.40665 × 10 −3 , and a is the lattice constant. Similar calculations such as faced-centered cubic or simple cubic ones can be carried out by evaluating different lattice configurations.

The SMMC method and Gamow-Teller response functions
Here we neglect the effect of SMFs on the GT properties here because the GT transition matrix elements for electron capture do not depend on the magnetic fields [35,43]. A detailed discussion of the SMMC method can be found in [44]. Based on a statistical formulation of the nuclear many-body problem in the finite-temperature version of this approach, an observable is calculated as the canonical expectation value of a corresponding operatorÂ by the SMMC method at a given temperature T , and is written by [45][46][47][48] The problem of the shell-model HamiltonianĤ has been investigated in detail by Ref. [48]. When a certain manybody HamiltonianĤ is given, a tractable expression for an imaginary-time evolution operator is written byÛ = exp −βĤ , where β = 1/T N , T N is the nuclear temperature in units of MeV. Tr AÛ is the canonical partition function for A nucleons. In terms of a spectral expansion, the total strength of a transition operatorÂ is then given by the following expectation value: here |i , | f are the many-body states of the initial and final nucleus with energy E i , E f , respectively. The SMMC method is used to calculate the response function R A (τ ) of an operatorÂ at an imaginary time τ . By using a spectral distribution of initial and final states |i and | f with energies E i and E f . R A (τ ) is given by [4,7,8] Note that the total strength for the operator is given by R(τ = 0). S GT + is the total amount of the GT strength available for an initial state given by summing over a complete set a final states in GT transition matrix elements |M GT | 2 i f . The strength distribution is given by [4] which is related to R A (τ ) by a Laplace transform, Note that here E is the energy transfer within the parent nucleus, and the strength distribution S GT + (E) has units of MeV −1 .

The study of EC
The weak interactions were investigated by some authors, e.g., Refs. [8,10], by using the RPA theory and SMMC method with a global parameterizations of the single-particle numbers, The EC rate in the absence of a SMFs is related to the electron capture cross-section and can be given by Eq. (4) in Ref. [4]. We can compute the EC rates from one of the initial states to all possible final states in a SMFs and in a weak magnetic field approximation (B << B cr ). The EC rates in both cases are expressed by [12,13] and respectively, where θ n = g no the ε n is the total rest mass and kinetic energies, and F(Z , ε n ) is the Coulomb wave correction.
In the case without and with SMFs, the p 0 and p B 0 are defined as and respectively, where = m 2 e c 4 (1 + 2ν B/B cr ) = m 2 e c 4 (1 + 2νb).

Results and discussions
The magnetar surface thermal temperature is high as (4−6)× 10 6 K, and the crust temperature may be higher, about several 10 8 K. For example, the maximum of inner crust temperature for some light elements may be as high as 2 × 10 8 K (e.g., [23,49]), but the maximum of inner crust temperature for some iron elements may be as high as 10 9 K (e.g., [50][51][52]). In this paper, we discuss the electron capture for iron group nuclei, so we select several typical inner crust temperature value at the order of 10 9 K.
In considering relativistic electrons in magnetar crust, the electron chemical potential should range from 10 m e c 2 to 100 m e c 2 . Figure 1 shows the EC rates of 55 Co, and 56 Ni as a function of U F (note ρ 7 and T 9 are in units of 10 7 g/cm 3 and 10 9 K, respectively). The electron chemical potential effects greatly on EC process in SMFs due to the different Q-value and transition orbits. When density and temperature are given (e.g., ρ 7 = 5.86, Y e = 0.47, T 9 = 3.40), the rates decrease when 1 m e c 2 < U F < 5 m e c 2 , then increase by about more than one order of magnitude as electron chemical potential increases when 5 m e c 2 < U F ≤ 100 m e c 2 . According to the discussions in Refs. [19][20][21], the electron chemical potential is strongly depended on the factor B −2 12 (Y e ρ) 2 . When magnetic field strength is given, the higher the density, the larger the electron chemical potential. When the density and temperature are given, the smaller the magnetic fields, the larger the electron chemical potential becomes. Thus, a large number of electrons involve an EC reaction because their chemical potential exceeds the Q-value.
A SMFs has a great influence on EC rates. From Fig. 2, one sees that the higher the density, the larger the rates become. Because the electron kinetic energy and electron chemical potential are so high that the influence of SMFs on rates is dramatic at relativity higher density. As the SMFs increases, the rates increase by about two orders of magnitude. Then the rates decrease by more than two orders magnitude when 10 14 G < B < 10 15 G. But then another increase of the EC rates appears when B > 10 15 G.
Form Fig. 2, we find that the EC rates decrease when B 12 > 100. According to Eqs. (4)-(8), when ν max 1, for electron Fermi energy, the Landau energy level spacing becomes a very small fraction. When the electron gas is in a mildly relativistic state, we have ν max → U e /hω c . But for relativistic electron gas, the ν max ≥ 100/B 12 when the U e ≥ 1 MeV. Thus, the ν max tends to the order of unity, and the electron Landau levels will be termed strongly quantized as the magnetic field strength increases (i.e. when B 12 ≥ 10 ∼ 100).
Based on the RMFEFT models of [32,33], we study the effect of SMFs on the binding energy per particle, which has a parabolic increase as the SMFs increases. According to Eqs. (12)-(18), we calculate the binding energy for 56 Ni, and 55 Co, and find that the binding energies increase by 0.601, and 0.402 MeV, respectively when the SMFs increases from 10 16 to 10 18 G. This is equivalent to significantly raise the threshold energy of EC due to an increase in nuclear binding energy. Thus, the EC is crippled greatly by a SMFs. Meanwhile, because of the interaction between the electrons and SMFs, the electron Fermi energy decreases as a SMFs increases. This actually discourages the EC reaction.
Based on Eqs. (17) and (18), and Refs. [32,33], we find that an abrupt increase is shown for EC rates when 10 14 G ≤ B ≤ 10 16 G shown as in Fig. 2. Such a jump shows that the underlying shell structure may be changed in a fundamental way. This jump in nuclear properties can be traced to the single-particle behavior. A particle will move from a level going upwards and go to a level downward with increasing spin by SMFs. Because of the two levels have opposite angular momentum along the symmetry axis, the nucleus becomes spin-polarized. The single-particle structure for protons and neutrons is strongly modified by a SMFs. Firstly, the strong interaction between the magnetic field and the neutron (proton) magnetic dipole moment will lead to the nucleon paramagnetism. Secondly, the coupling of the orbital motion of protons with SMFs also make the proton orbital magnetism. The interaction between the nucleus and SMFs maybe removes all degeneracies in the single-particle spectrum, and it breaks significantly the formerly degenerate levels with opposing signs of the angular momentum projec-   . 3 The GT + distributions S GT + for nuclei 55 Co and 56 Ni as a function of the energy transfer from the electron, E = E f − E i at the temperature of T = 0.8 MeV tion. A reduction of the neutron and proton pairing gaps will show because of splitting from single-particle energy as the SMFs increases. Finally, it will lead to their disappearance.
Here we neglect the effect of SMFs on the GT because the GT transition matrix elements for electron capture do not depend on the SMFs [35,43]. Figure 3 shows the strength distributions S GT + for 55 Co and 56 Ni as a function of the daughter state excitation energy. According to Eqs. (21)-(23), we calculate the S GT + and find the peak of it reaches to 2.211, and 4.89 MeV −1 at 3.5, 1.45 MeV of daughter nuclei 55 Fe, and 56 Co, respectively. The total GT strength for 56 Ni in a full p f -shell may be calculated by B(GT) = 10.1g 2 A according to the data of experiment. For 55 Co, we discuss the total GT strength and the distribution in a truncated calculation, in which five particles in the final nucleus are maximally allowed to be excited out of the state of f 7/2 . We obtain a Table 1 Comparisons of the maximum value of our calculations λ B max (LJ) with those of minimum value λ B min (LJ) for 55 Co, and 56 Ni when 10 13 G ≤ B ≤ 10 18 G under the conditions CD i , i = 1-6, corresponding to ρ 7 = 5.86, Y e = 0.47, T 9 = 3.40; ρ 7 = 6.63, Y e = 0.46, T 9 = 3.51; ρ 7 = 11.7, Y e = 0.45, T 9 = 3.73; ρ 7 = 14.5, Y e = 0.44, T 9 = 3.80; ρ 7 = 106, Y e = 0.43, T 9 = 4.93, and ρ 7 = 4010, Y e = 0.41, T 9 = 7.33. The ratio is computed as  total GT strength of 8.7g 2 A from the ground states of 55 Co, and 8.9g 2 A from both of the excited states of J = 3/2. According to Eqs. (6)- (8) and (25), we calculate the maximum and the minimum value of the EC rates in SMFs for the two typical nuclei by considering the effect of all evaluating factors, and list them in Table 1. The maximum rates reach 5.875×10 5 s −1 , and 3.062×10 5 s −1 when B = 3.432×10 16 G for 55 Co, 56 Ni, at ρ 7 = 4010, Y e = 0.41, T 9 = 7.33, respectively. However, the minimum EC rates are 7.207 s −1 , and 2.873 s −1 when B = 3.678 × 10 14 G for 55 Co, 56 Ni, at ρ 7 = 5.86, Y e = 0.47, T 9 = 3.40, respectively. The rates increase by about three orders magnitude as the SMFs increases.
FFN [1][2][3], AUFD [5,6], and NKK [9] studied EC rates in the case that there are no SMFs. In Tables 2 and 3, we present the comparisons of our results in SMFs with those of FFN, AUFD, and NKK. Our results in the case that there are no SMFs are about 14.7, and 26.07% lower than FFN, and AUFD, respectively, but they are about one order of magnitude larger than that of NKK (e.g. at ρ 7 = 5.86, T 9 = 3.40, Y e = 0.47). Our SMFs rates are by about six orders of magnitude higher than those of FFN (AUFD), and NKK when 10 < B 12 < 10 6 for 55 Co, and 56 Ni. FFN used the semiempirical atomic mass formula from Ref. [53] to estimate the Q-value and the EC rates, thus the Q-value used in the effective rates are quite different. Based on the nuclear shell model and the Brink hypothesis method, AUFD expanded the FFN work and analyzed the nuclear excited level by a simple calculation on the nuclear excitation levels transitions. The Brink hypothesis is a very crude approximation. Thus the calculation method is a little rough. Using the quasi-particle random phase approximation theory, NKK expanded the nuclear excitation energy distribution by considering the particle emission processes, which constrained the parent excitation energies. However, only low angular momentum states are considered. The SMMC method actually draws an average of GT intensity distribution of the EC process, and the calculated results are in good agreement with experiments. Thus the method is relatively accurate.

Summary
In this work, we investigate the EC for 55 Co, and 56 Ni and the influences of SMFs on electron Fermi energy, binding energy per nuclei, and single-particle level structure in magnetar crust based on the RMFEIT, and Lai dong model; we compare our results with those of FFN, AUFD, and NKK. Our results increase by about two orders of magnitude as SMFs increase, and then decrease by more than two orders magnitude. There is an abrupt increase in EC rates around B = 3 × 10 14 G for ρ 7 < 100 (but around B = 3 × 10 15 G for ρ 7 > 100). Such a jump may be an indication that the underlying shell structure has changed in a fundamental way because of single-particle behavior by SMFs. We find that the EC rates in the case that there are no SMFs are about 14.70, and 26.07 % lower than those of FFN, and AUFD (e.g. for 55 Co at ρ 7 = 4.32, T 9 = 3.26, Y e = 0.47), but they are about one order of magnitude larger than that of NKK (e.g., at ρ 7 = 5.86, T 9 = 3.40, Y e = 0.47). But our EC rates in SMFs are by about six orders of magnitude higher than those of FFN, AUFD, and NKK when 10 < B 12 < 10 6 for 55 Co, and 56 Ni.
The observations show that the persistent X-ray emission of a magnetar could originate from magnetic field decay or heating from magnetospheric current and/or from the EC reaction in the crust. The research on EC reaction in magnetar's crust is still a long-term and arduous task. Our results may be helpful to the future study of the magnetar thermal evolution, and soft X-ray emission mechanism.