# Electron Capture Cross Sections for Stellar Nucleosynthesis

Authors: P. G. Giannaka; T. S. Kosmas
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

In the first stage of this work, we perform detailed calculations for the cross sections of the electron capture on nuclei under laboratory conditions. Towards this aim we exploit the advantages of a refined version of the proton-neutron quasiparticle random-phase approximation (pn-QRPA) and carry out state-by-state evaluations of the rates of exclusive processes that lead to any of the accessible transitions within the chosen model space. In the second stage of our present study, we translate the abovementionede
--capture cross sections to the stellar environment ones by inserting the temperature dependence through a Maxwell-Boltzmann distribution describing the stellar electron gas. As a concrete nuclear target we use the 66Zn isotope, which belongs to the iron group nuclei and plays prominent role in stellar nucleosynthesis at core collapse supernovae environment.


---

## Body

## 1. Introduction

Weak interaction processes occurring in the presence of nuclei under stellar conditions play crucial role in the late stages of the evolution of massive stars and in the presupernova stellar collapse [1–6]. As it is known, the core of a massive star, at the end of its hydrostatic burning, is stabilized by electron degeneracy pressure as long as its mass does not exceed an appropriate mass (the Chandrasekhar mass limit, M
Ch) [6–10]. When the core mass exceeds M
Ch, electron degeneracy pressure cannot longer stabilize the center of the star and the collapse starts. In the early stage of collapse electrons are captured by nuclei in the iron group region [6, 10].During the presupernova evolution of core collapse supernova, the Fermi energy (or equivalently the chemical potential) of the degenerate electron gas is sufficiently large to overcome the threshold energyE
thr (E
thr is given by negative Q values of the reactions involved in the interior of the stars) [11] and the nuclear matter in the stellar core is neutronized. This high Fermi energy of the degenerate electron gas leads to enormous e
--capture on nuclei and reduces the electron to baryon ratio Y
e [12, 13]. In this way, the electron pressure is reduced and the energy as well as the entropy drop. One of the important characteristics of the early pre-explosion evolution is the fact that electron capture on nuclei (specifically on nuclei of the pf shell) plays a key role [14, 15].In the early stage of collapse (for densities lower than a few1
0
10g cm−3), the electron chemical potential is of the same order of magnitude as the nuclear Q value, and the e
--capture cross sections are sensitive to the details of GT strength distributions in daughter nuclei. For this reason, some authors restrict the calculations only to the GT strength and evaluate e
--capture rates on the basis of the GT transitions (at these densities, electrons are captured mostly on nuclei with mass number A
≤
60) [9–12, 15, 16]. Various methods, used for calculating e
--capture on nuclei during the collapse phase, have shown that this process produces neutrinos with rather low energies in contrast to the inelastic neutrino-nucleus reactions occurring in supernova [17–20]. These neutrinos escape the star carrying away energy and entropy from the core which is an effective cooling mechanism of the exploding massive star [21]. For higher densities and temperatures, e
- capture occurs on heavier nuclei A
≥
65 [8, 10, 13–15]. As a consequence, the nuclear composition is shifted to more neutron-rich and heavier nuclei (including those with N
>
40) which dominate the matter composition for densities larger than about 1
0
10 g cm−3  [1, 15, 21, 22].The first calculations of stellar electron capture rates for iron group nuclei have been performed by employing the independent particle model (IPM) [2–4]. Recently, similar studies have been addressed by using continuum RPA (CRPA) [25], large scale shell model [26, 27], RPA [11], and so forth [28]. In the present work e
--capture cross sections are obtained within a refined version of the quasiparticle random phase approximation (QRPA) which is reliable for constructing all the accessible final (excited) states of the daughter nuclei in the iron group region of the periodic table [29–38]. For the description of the required correlated nuclear ground states we determine single-particle occupation numbers calculated within the BCS theory [29, 31, 32]. Our nuclear method is tested through the reproducibility of experimental muon capture rates relying on detailed calculations of exclusive, partial, and total muon capture rates [23, 24, 39–42]. The agreement with experimental data provided us with high confidence level of our method and we continued with the calculations of electron capture cross sections in supernova conditions (where the densities and temperatures are high) using the pn-QRPA method. In this paper, we performed calculations for Zn
66 isotope (it belongs to the iron group nuclei) that plays prominent role in core collapse supernovae stellar nucleosynthesis [18, 19, 43].Our strategy in this work is, at first, to perform extensive calculations of the transition rates for all the abovementioned nuclear processes, assuming laboratory conditions, and then to translate these rates to the corresponding quantities within stellar environment through the use of an appropriate convolution procedure [9, 14, 15, 21, 26]. To this purpose, we assume that leptons under such conditions follow Maxwell-Boltzmann energy distribution [9, 26].


## 2. Construction of Nuclear Ground and Excited States

Electrons of energyE
e are captured by nuclei interacting weakly with them via W
- boson exchange as follows:(1)
A
,
Z
+
e
-
⟶
A
,
Z
-
1
∗
+
ν
e
.The outgoing ν
e neutrino carries energy E
ν while the daughter nucleus (
A
,
Z
-
1
) absorbs a part of the incident electron energy given by the difference between the initial E
i and the final E
f nuclear energies as E
ν
=
E
f
-
E
i.The nuclear calculations for the cross sections of reaction (1) start by writing down the weak interaction Hamiltonian H
^
w which is given as a product of the leptonic, j
μ
lept, and the hadronic, J
^
μ, currents (current-current interaction Hamiltonian) as follows:(2)
H
^
w
=
G
2
j
μ
lept
J
^
μ
,where G
=
G
F
cos
⁡
θ
c with G
F and θ
c being the well-known weak interaction coupling constant and the Cabibbo angle, respectively [31, 32, 44].From the nuclear theory point of view, the main task is to calculate the cross sections of reaction (1) which are based on the evaluation of the nuclear transition matrix elements between the initial i and a final f nuclear states of the form(3)
f
H
w
^
i
=
G
2
l
μ
∫
‍
d
3
x
e
-
i
q
·
x
f
J
μ
^
i
.The quantity l
μ
e
-
i
q
·
x stands for the leptonic matrix element written in coordinate space with q being the 3-momentum transfer. For the calculation of these transition matrix elements one may take advantage of the Donnelly-Walecka multipole decomposition which leads to a set of eight independent irreducible tensor multipole operators containing polar-vector and axial-vector components [44] (see Appendix A).In the present work, in (3) the ground state of the parent nucleus i is computed by solving the relevant BCS equations which give us the quasiparticle energies and the amplitudes V and U that determine the probability for each single particle level to be occupied or unoccupied, respectively [31]. Towards this aim, at first, we consider a Coulomb corrected Woods-Saxon potential with a spin orbit part as a mean field for the description of the strong nuclear field [45, 46]. For the latter potential we adopt the parameterization of IOWA group [47]. Then, we use as pairing interaction the monopole part of the Bonn C-D one meson exchange potential. The renormalization of this interaction, to fit in the Zn
66 isotope, is achieved through the two pairing parameters g
pair
p
,
n, p
(
n
) for proton (neutron) pairs, and the values of which are tabulated in Table 1.Table 1
Parameters for the renormalization of the interaction of proton pairs,g
pair
p, and neutron pairs, g
pair
n. They have been fixed in such a way that the corresponding experimental gaps, Δ
p
exp
⁡ and Δ
n
exp
⁡ of 66Zn isotope, are quite accurately reproduced.
Nucleus
g
pair
n
g
pair
p
Δ
n
exp
⁡ (MeV)
Δ
n
theor (MeV)
Δ
p
exp
⁡ (MeV)
Δ
p
theor (MeV)
66Zn
1.0059
0.9271
1.7715
1.7716
1.2815
1.2814As it is well known, the pairing parameters,g
pair
p
,
n, are determined through the reproduction of the energy gaps, Δ
p
,
n
exp
⁡, from neighboring nuclei as follows (3-point formula):(4)
Δ
n
exp
⁡
=
-
1
4
S
n
A
-
1
,
Z
-
2
S
n
A
,
Z
+
S
n
A
+
1
,
Z
,
Δ
p
exp
⁡
=
-
1
4
S
p
A
-
1
,
Z
-
1
Δ
p
exp
⁡
Δ
p
exp
⁡
-
2
S
p
A
,
Z
+
S
p
A
+
1
,
Z
+
1
,where S
p and S
n are the experimental separation energies for protons and neutrons, respectively, of the target nucleus (
A
,
Z
) and the neighboring nuclei (
A
±
1
,
Z
±
1
) and (
A
±
1
,
Z
). For the readers convenience in Table 2 we show the values of experimental separation energies for the target Zn
30
66 and the neighboring nuclei Cu
29
65, Ga
31
67, Zn
30
65, and Zn
30
67.Table 2
The experimental separation energies in MeV for protons and neutrons of the target (A, Z) and neighboring (
A
±
1
,
Z
±
1
) and (
A
±
1
,
Z
) nuclei.
Nucleus
S
n
(
A
-
1
,
Z
)
S
n
(
A
,
Z
)
S
n
(
A
+
1
,
Z
)
S
p
(
A
-
1
,
Z
-
1
)
S
p
(
A
,
Z
)
S
p
(
A
+
1
,
Z
+
1
)
66Zn
7.979
11.059
7.052
7.454
8.924
5.269Subsequently, the excited statesf of the studied daughter nucleus Cu
66 are constructed by solving the pn-QRPA equations [29–38], which are written as follows in matrix form [29]:(5)
A
B
-
B
-
A
X
ν
Y
ν
=
Ω
J
π
ν
X
ν
Y
ν
,
Ω
J
π
ν denotes the excitation energy of the QRPA state J
ν
π with spin J and parity π.The solution of (5) is an eigenvalue problem which provides the amplitudes for forward and backward scattering X and Y, respectively, and the QRPA excitation energies Ω
J
π
ν [31–34]. In our method the solution of the QRPA equations is carried out separately for each multipole set of states J
π.For the renormalization of the residual 2-body interaction (Bonn C-D potential), the strength parameters, for the particle-particle(
g
pp
) and particle-hole (
g
ph
) interaction entering the QRPA matrices A and B, are determined (separately for each multipolarity) from the reproducibility of the low-lying experimental energy spectrum of the final nucleus. The values of these parameters in the case of the spectrum of Cu
66 are listed in Table 3.Table 3
Strength parameters for the particle-particle(
g
pp
) and particle-hole (
g
ph
) interaction for various multipolarities (for the rest of multipolarities, J
π
≤
5
±, the bare 2-body interaction has been used), in the case of the spectrum of 66Cu nucleus.
Positive parity states
Negative parity states
J
π
0
+
1
+
2
+
3
+
4
+
1
-
2
-
3
-
4
-
g
pp
0.827
0.547
0.686
0.854
1.300
0.994
0.200
0.486
0.622
g
ph
0.336
0.200
1.079
0.235
0.200
1.200
0.200
1.200
1.200At this point, it is worth mentioning that, for measuring the excitation energies of the daughter nucleusCu
66 from the ground state of the initial one Zn
66, a shifting of the entire set of QRPA eigenvalues is necessary. Such a shifting is required whenever in the pn-QRPA a BCS ground state is used, a treatment adopted by other groups previously [9, 40, 48, 49]. The shifting, for the spectrum of the daughter nucleus Cu
66, is done in such a way that the first calculated value of each multipole state of Cu
66 (i.e. 1
1
+
,
2
1
+
,
…, etc.) approaches as close as possible the corresponding lowest experimental multipole excitation. Table 4 shows the shifting applied to our QRPA spectrum for each multipolarity of the parent nucleus Zn
66. We note that a similar treatment is required in pn-QRPA calculations performed for double-beta decay studies where the excitations derived for the intermediate odd-odd nucleus (intermediate states) through p-n and n-p reactions from the neighboring nuclei, left or right nuclear isotope, do not match to each other [48, 49]. The resulting low-energy spectrum, after using the parameters of Tables 1 and 3 and the shifting shown in Table 4, agrees well with the experimental one (see Figure 1).Table 4
The shift (in MeV) applied on the spectrum (seperately of each multipole set of states) of66Cu isotope, daughter nucleus of the electron capture on 66Zn.
Positive parity states
Negative parity states
0
+
0.90
0
-
5.00
1
+
2.50
1
-
6.80
2
+
2.55
2
-
3.85
3
+
2.50
3
-
2.60
4
+
1.75
4
-
3.55
5
+
0.55
5
-
3.00Figure 1
Comparison of the theoretical excitation spectrum (resulting from the solution of the QRPA eigenvalue problem) with the low-lying (up to about 3 MeV) experimental one forCu
66 nucleus. As it can be seen, the agreement is very good below 1 MeV but for higher excitation energies it becomes moderate.We must also mention that, usually, in nuclear structure calculations we test a nuclear method in two phases: first through the construction of the excitation spectrum as discussed before and second through the calculations of electron scattering cross sections or muon capture rates. Following the above steps, we test the reproducibility of the relevant experimental data for many nuclear models employed in nuclear applications (nuclear structure and nuclear reactions) and in nuclear astrophysics [1, 9].


## 3. Results and Discussion

In this work we perform detailed cross section calculations for the electron capture onZn
66 isotope on the basis of the pn-QRPA method. The required nuclear matrix elements between the initial J
i and the final J
f states are determined by solving the BCS equations for the ground state [29, 31, 32] and the pn-QRPA equations for the excited states [31–34] (see Section 2). For the calculations of the original cross sections, a quenched value of g
A (see Appendix B) is considered which subsequently modifies all relevant multipole matrix elements [23, 24, 50, 51].At this point of the present work and in order to increase the confidence level of our method, we perform total muon capture rates calculations [23, 24, 39–42]. The comparison with experimental and other theoretical results is shown in Section 3.1. Afterwards, we study in detail the electron capture process as follows. (i) Initially we consider laboratory conditions; that is, the initial (parent) nucleus is considered in the ground state and no temperature dependence is assumed (see Section 3.2.1). (ii) Second, we consider stellar conditions; that is, the parent nucleus is assumed to be in any initial excited state and due to the e
--capture process it goes to any final excited state of the daughter nucleus. At these conditions it is necessary to take into account the temperature dependence of the cross sections (see Section 3.2.2) [10].


### 3.1. Calculations of Muon Capture Rates for66Zn

Despite the fact that the muon capture on nuclei does not play a crucial role in stellar-nucleosynthesis, it is, however, important to start our study from this process since the nuclear matrix elements required for an accurate description of theμ-capture are the same for all semileptonic charge-changing weak interaction processes. In addition, the excitation spectrum of the daughter (
A
,
Z
-
1
) nucleus, as we saw before, is in good agreement with the experimental data.The calculations of the muon capture rates are performed in three steps. In the first step we carry out realistic state-by-state calculations of exclusive ordinary muon capture (OMC) rates inZn
66 isotope for all multipolarities with J
π
≤
5
± (higher multipolarities contribute negligibly). The appropriate expression for the exclusive muon capture rates is written as follows:(6)
Λ
g
s
→
J
f
π
≡
Λ
J
f
π
=
2
G
2
Φ
1
s
2
R
f
q
f
2
×
J
f
π
M
^
J
-
L
^
J
0
g
s
+
2
+
J
f
π
T
^
J
el
-
T
^
J
magn
0
g
s
+
2
,where Φ
1
s represents the muon wave function in the 1
s muonic orbit. The operators in (6) are refered to as Coulomb M
^
J, longitudinal L
^
J, transverse electric T
^
J
el, and transverse magnetic T
^
J
mag multipole operators (see Appendix A). The factor R
f in (6) takes into consideration the nuclear recoil which is written as R
f
=
(
1
+
q
f
/
M
targ
)
-
1, with M
targ being the mass of the target (parent) nucleus.Due to the fact that there are no available data in the literature for exclusive muon capture rates, the test of our method is realized by comparing partial and total muon capture rates with experimental data and other theoretical results [23, 24]. Towards this purpose, our second step includes calculations of the partial μ
--capture rates for various low-spin multipolarities, Λ
J
π (for J
π
≤
4
±), in the studied nucleus. These partial rates are found by summing over the contributions of all the individual multipole states of the studied multipolarity as follows: (7)
Λ
J
π
=
∑
f
‍
Λ
g
s
→
J
f
π
=
2
G
2
Φ
1
s
2
∑
f
‍
q
f
2
R
f
J
f
π
M
^
J
-
L
^
J
0
g
s
+
2
k
k
k
k
k
k
k
k
k
k
k
k
+
∑
f
‍
q
f
2
R
f
J
f
π
T
^
J
el
-
T
^
J
magn
0
g
s
+
2(f runs over all states of the multipolarity J
π). We also estimate the percentage (portion) of their contribution into the total μ-capture rate for the most important multipolarities. In Table 5 we tabulate the individual portions of the low-spin multipole transitions (J
π
=
4
±) into the total muon capture rate. As it can be seen, the contribution of the 1
- multipole transitions is the most important multipolarity exhausting about 44
% of the total muon capture rate. Such an important contribution was found in O
16 and Ca
48 isotopes studied in [41].Table 5
The percentage of each multipolarity into the total muon capture rate evaluated with our pn-QRPA method.
Positive parity transitions
Negative parity transitions
J
π
Portions(
%
)
J
π
Portions(
%
)
0
+
8.22
0
-
7.94
1
+
21.29
1
-
44.21
2
+
2.85
2
-
13.32
3
+
1.58
3
-
0.34
4
+
0.01
4
-
0.23In the last step of testing our method, we evaluate total muon capture rates for theZn
66 isotope. These rates are obtained by summing over all partial multipole transition rates (up to J
π
=
4
±) as follows:(8)
Λ
tot
=
∑
J
π
‍
Λ
J
π
=
∑
J
π
‍
∑
f
‍
Λ
J
f
π
.For the sake of comparison, the abovementionedμ-capture calculations have been carried out using the quenched value g
A
=
1.135 [23, 24]. The results are listed in Table 6, where we also include the experimental total rates and the theoretical ones of [23, 24]. Moreover, in Table 6 we show the individual contribution into the total muon capture rate of the polar-vector (Λ
tot
V), the axial-vector (Λ
tot
A), and the overlap (Λ
tot
V
A) parts. As it can be seen, our results obtained with the quenched g
A coupling constant are in very good agreement with the experimental total muon capture rates (the deviations from the corresponding experimental rates are smaller than 7
%). This agreement provides us with high confidence level for our method.Table 6
Individual contribution of polar-vector, axial-vector, and overlap parts into the total muon-capture rate. The total muon capture rates, obtained by using the pn-QRPA with the quenched value ofg
A
=
1.135 for the medium-weight nucleus 66Zn, are compared with the available experimental data and with the theoretical rates of [23, 24].
Total muon capture ratesΛ
tot(×106 s
-
1)
Present pn-QRPA calculations
Experiment
Other theoretical methods
Nucleus
Λ
tot
V
Λ
tot
A
Λ
tot
V
A
Λ
tot
Λ
tot
exp
⁡
Λ
tot
theor [23]
Λ
tot
theor [24]
66Zn
1.651
4.487
−0.204
5.934
5.809
4.976
5.809


### 3.2. Electron Capture Cross Section

After acquiring a high confidence level for our nuclear method, we proceed with the main goal of the present study which concerns the calculations of the electron capture cross sections. As mentioned before, this includes original (see Section3.2.1) and stellar electron capture investigations (Section 3.2.2).


#### 3.2.1. Original Electron Capture Cross Section on66Zn Isotope

The original cross sections for the electron capture process in theZn
66 isotope are obtained by using the pn-QRPA method considering all the accessible transitions of the final nucleus Cu
66. In the Donnelly-Walecka formalism the expression for the differential cross section in electron capture by nuclei reads [10](9)
d
σ
e
c
d
Ω
=
G
F
2
cos
⁡
2
θ
c
2
π
F
Z
,
E
e
2
J
i
+
1
H
×
∑
J
≥
1
‍
W
E
e
,
E
ν
H
H
i
×
1
-
α
cos
⁡
Φ
+
b
sin
2
Φ
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
×
J
f
T
^
J
mag
J
i
2
+
J
f
T
^
J
el
J
i
2
H
H
H
H
i
-
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
×
2
Re
J
f
T
^
J
mag
J
i
J
f
T
^
J
el
J
i
∗
H
H
i
+
∑
J
≥
0
‍
W
E
e
,
E
ν
H
H
i
i
i
i
i
i
i
×
1
+
α
cos
⁡
Φ
J
f
M
^
J
J
i
2
ω
q
1
+
α
cos
⁡
Φ
+
d
H
H
H
H
i
i
i
i
i
i
+
1
+
α
cos
⁡
Φ
-
2
b
sin
2
Φ
J
f
L
^
J
J
i
2
H
H
H
H
i
i
i
i
i
i
-
ω
q
1
+
α
cos
⁡
Φ
+
d
2
Re
J
f
L
^
J
J
i
∑
J
≥
1
‍
H
H
H
H
i
i
i
i
i
ω
q
1
+
α
cos
⁡
Φ
+
d
×
J
f
M
^
J
J
i
∗
,where F
(
Z
,
E
e
) is the well-known Fermi function [18]. The factor W
(
E
e
,
E
ν
)
=
E
ν
2
/
(
1
+
E
ν
/
M
T
) accounts for the nuclear recoil [8], M
T is the mass of the target nucleus, and parameters α, b, and d are given, for example, in [31]. The nuclear transition matrix elements between the initial state J
i and the final state |
J
f
〉 correspond to the Coulomb M
^
J
M, longitudinal L
^
J
M, transverse electric T
^
J
M
el, and transverse magnetic T
^
J
M
mag multipole operators (discussed in Appendix A).From the energy conservation in the reaction (1), the energy of the outgoing neutrino E
ν is written as follows:(10)
E
ν
=
E
e
-
Q
+
E
i
-
E
f
,which includes the difference between the initial E
i and the final E
f nuclear states. The Q value of the process is determined from the experimental masses of the parent (M
i) and the daughter (M
f) nuclei as Q
=
M
f
-
M
i [9].It is worth mentioning that for low momentum transfer, various authors use the approximationq
→
0 for all multipole operators of (9). Then, the transitions of the Gamow-Teller operator (G
T
+
=
∑
i
‍
τ
i
+
σ
i) provide the dominant contribution to the total cross section [9].While performing detailed calculations for the original electron capture cross sections onZn
66 isotope we assumed that (i) the initial state of the parent nucleus Zn
66 is the ground state 0
+ and (ii) the nuclear system is under laboratory conditions (no temperature dependence of the cross sections is needed). The cross sections as a function of the incident electron energy E
e are calculated with the use of realistic two-body interactions as mentioned before. The obtained total original electron capture cross sections for Zn
66 target nucleus are illustrated in Figure 2 where the individual contributions of various multipole channels (J
π
≤
5
±) are also shown. The electron capture cross sections in Figure 2 exhibit a sharp increase by several orders of magnitude within the first few MeV above energy-threshold, and this reflects the G
T
+ strength distribution. For electron energy E
e
≥
10
MeV the calculated cross sections show a moderate increase. From experimental and astrophysical point of view, the important range of the incident electron energy E
e is up to 30 MeV. At these energies the 1
+ multipolarity has the largest contribution to the total electron capture cross sections [9, 10]. In the present work we have extended the range of E
e up to 50 MeV since at higher energies (around 40 MeV) the contribution of other multipolarities like 1
-, 0
+, and 0
- becomes noticeable and cannot be omitted (see Figure 2).Figure 2
Original total cross sections of electron capture on theZn
66 (parent) nucleus calculated with pn-QRPA method as a function of the incident electron energy E
e. The individual contributions of various multipole channels (for J
π
≤
5
±) are also demonstrated.From the study of the original electron capture cross sections we conclude that the total cross sections can be well approximated with the Gamow-Teller transitions only in the region of low energies [9–12, 15, 16]. For higher incident electron energies the inclusion of the contributions originating from other multipolarities leads to better agreement [10].


#### 3.2.2. Stellar Electron Capture on66Zn Isotope

As it is well known, electron capture process plays a crucial role in late stages of evolution of a massive star, in presupernova and in supernova phases [1–6]. In presupernova collapse, that is, at densities ρ
≤
1
0
10
g
cm
-
3 and temperatures 300
keV
≤
T
≤
800
keV, electrons are captured by nuclei with A
≤
60 [9–12, 15, 16]. During the collapse phase, at higher densities ρ
≥
1
0
10
g
cm
-
3 and temperatures T
≃
1.0
MeV, electron capture process is carried out on heavier and more neutron rich nuclei with Z
<
40 and N
≥
40 [8, 10, 13–15].In an independent particle picture, the Gamow-Teller transitions (which is the most important in the electron capture cross section calculations) are forbidden for these nuclei [2–4]. However, as it has been demonstrated in several studies, GT transitions in these nuclei are unblocked by finite temperature excitations [21, 22]. At high temperatures, T
≃
1.5
MeV, GT transitions are thermally unblocked as a result of the excitation of neutrons from the pf-shell into the g
9
/
2 orbital.For astrophysical environment, where the finite temperature and the matter density effects cannot be ignored (the initial nucleus is at finite temperature), in general, the initial nuclear state needs to be a weighted sum over an appropriate energy distribution. Then, assuming Maxwell-Boltzmann distribution of the initial statei in (9) [9, 26], the total e
--capture cross section is given by the expression [10](11)
σ
E
e
,
T
=
G
F
2
cos
⁡
2
θ
c
2
π
∑
i
‍
F
Z
,
E
e
2
J
i
+
1
e
-
E
i
/
k
T
G
Z
,
A
,
T
×
∑
f
,
J
‍
E
e
-
Q
+
E
i
-
E
f
2
i
O
^
J
f
2
2
J
i
+
1
.The sum over initial states in the latter equation denotes a thermal average of levels, with the corresponding partition function G
(
Z
,
A
,
T
) [10]. The finite temperature induces the thermal population of excited states in the parent nucleus. In the present work we assume that these excited states in the parent nucleus are all the possible states up to about 2.5 MeV. Calculations involving in addition other states lying at higher energies show that they have no sizeable contribution to the total electron capture cross sections. As mentioned before, for the evaluation of the total electron capture cross sections, the use of a quenched value of g
A is necessary [23, 24, 50, 51]. Since the form factor F
A
(
q
2
) multiplies the four components of the axial-vector operator (see (A.2)–(A.5)), a quenched value of g
A must enter the multipole operators generating the pronounced excitations 0
-
,
1
±
,
…, and so forth. For this reason, in our QRPA calculations we multiplied the free nucleon coupling constant g
A
=
1.262 by the factor 0.8 [23, 24, 50, 51].The results coming out of the study of electron capture cross sections under stellar conditions are shown in Figure3, where the same picture as in the original cross section calculations, but now with larger contribution, is observed. As discussed before, the dominant multipolarity is the 1
+, which contributes by more than 40
% to the total cross section. In the region of low energies (up to 30 MeV), the total e
--capture cross section can be described by taking into account only the GT transitions, but at higher incident energies the contributions of other multipolarities become significant and cannot be omitted.Figure 3
Electron capture cross sections for theZn
66 parent nucleus at high temperature (T
=
0.5 MeV) in stellar environment obtained assuming Maxwell-Boltzmann statistics for the incident electrons. The total cross section and the dominant individual multipole channels (J
π
≤
5
±) are demonstrated as functions of the incident electron energy E
e.The percentage contributions of various multipolarities (withJ
π
≤
5
±) into the total e
--capture cross section at T
=
0.5 MeV and for incident electron energy E
e
=
25
MeV are tabulated in Table 7. In addition, in this table we list the values of the individual e
--capture cross sections of each multipolarity with J
π
≤
5
±. More specifically, for E
e
=
25
MeV the 1
+ multipolarity contributes to about 44
%, the 0
+ contributes to about 26
%, and the 1
- contributes to about 11
%. The contributions coming from other multipolarities are less important (smaller than 5
%).Table 7
Totale
--capture cross sections (in 1
0
-
42 MeV−1 cm2) for E
e
=
25 MeV. The percentage of each multipolarity into the total e
--capture cross section evaluated with our pn-QRPA code is also tabulated here.
Positive parity transitions
Negative parity transitions
J
π
σ
e
(×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
J
π
σ
e (×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
0
+
31.164
25.96
0
-
5.288
4.41
1
+
52.779
43.98
1
-
13.409
11.14
2
+
6.921
5.77
2
-
3.262
2.72
3
+
5.499
4.58
3
-
0.905
0.75
4
+
0.244
0.20
4
-
0.299
0.25
5
+
0.208
0.17
5
-
0.042
0.04In performing state-by-state calculations for the electron capture cross sections, our code has the possibility to provide separately the contribution of the polar-vector, the axial-vector, and the overlap parts induced by the corresponding components of the electron capture operators. In Figure4 we illustrate the stellar differential cross sections of each individual transition of the polar-vector and axial-vector components.Figure 4
Individual contributions of the polar-vector, (Λ
V), and axial-vector, (Λ
A), components and total electron capture rate as functions of the excitation energy ω (Zn
66 is the parent nucleus).As mentioned before, our code gives separately the partiale
--capture cross sections of each multipolarity. In order to study the dependence of the differential cross sections on the excitation energy ω throughout the entire pn-QRPA spectrum of the daughter nucleus, a rearrangement of all possible excitations ω in ascending order, with the corresponding cross sections, is required. This was performed by using a special code appropriate for matrices [33]. In the model space chosen for Zn
66 isotope, for all multipolarities up to J
π
=
5
±, we have a number of 447 final states. The differential electron capture cross sections illustrated in Figure 4 present some characteristic clearly pronounced peaks at various excitation energies ω. These peaks correspond mainly to 0
+
,
1
+, and 2
+ transitions. More specifically, in the Cu
66 daughter nucleus the maximum peak corresponds to the 0
1
+ QRPA transition at ω
=
2.538
MeV and other characteristic peaks correspond to 1
7
+, 1
8
+, and 1
10
+ transitions, located at energies ω
=
3.194
MeV, ω
=
3.686
MeV, and ω
=
6.555
MeV, respectively (see Figure 4). The other less important peaks are also shown in Figure 4.Before closing, it should be mentioned that thee-capture cross sections presented in this work may be useful in estimating neutrino-spectra arising from e-capture on nuclei during supernova phase. The knowledge of ν-spectra at every point and time in the core is quite relevant for simulations of the final collapse and explosion phase of a massive star. As it is known [21], in the collapse phase, neutrinos are mainly produced by e-capture on nuclei and on free protons. The energy spectra of the emerging neutrinos from both reactions are important ingredients in stellar modelling and stellar simulations [21, 27].Furthermore, in core collapse simulations one defines the reaction rate of electron capture on nuclei given by(12)
R
h
=
∑
i
‍
Y
i
λ
i
,where the sum runs over all nuclear isotopes present in the astrophysical environment (Y
i denotes the abundance of a given nuclear isotope and λ
i is the calculated electron capture rate for this isotope). The rates of (12) must be known for a wide range of the parameters: T (temperature) and ρ (nuclear density) of the studied star. Thus, for the calculation of the quantity Y
·
λ of a specific nuclear isotope one needs to know in addition to the nuclear composition Y the electron capture rates λ calculated as we have shown in our present work. The rates of electron capture on various nuclear isotopes and the corresponding emitted neutrino spectra in the range of the parameters (T, ρ, Y
e) describing the star until reaching equilibrium during the core collapse are comprehensively studied in [21, 26, 27] for a great number of nuclear isotopes by using the large scale shell model. We are currently performing similar calculations for a set of isotopes by employing the present pn-QRPA method [52].Furthermore, the average neutrino energy,〈
E
ν
〉, of the neutrinos emitted by e-capture on nuclei can be obtained by dividing the neutrino-energy loss rate (defined by an expression similar to (12) by replacing the rate λ
i with the energy loss rate E
j) with the reaction rate for e-capture on nuclei R
h. Assuming, for example, power-law energy distribution for the neutrino spectrum produced by the e-capture in supernova phase, the average neutrino-energy 〈
E
ν
〉 determines a specific supernova-neutrino scenario. In addition, the neutrino emissivity is obtained by multiplying the electron capture rate at nuclear statistical equilibrium with the neutrino-spectra [21, 26, 27]. Finally we note that the rates for the inverse neutrino absorption process are also determined from the electron capture rates obtained as discussed in this section [21].


## 3.1. Calculations of Muon Capture Rates for66Zn

Despite the fact that the muon capture on nuclei does not play a crucial role in stellar-nucleosynthesis, it is, however, important to start our study from this process since the nuclear matrix elements required for an accurate description of theμ-capture are the same for all semileptonic charge-changing weak interaction processes. In addition, the excitation spectrum of the daughter (
A
,
Z
-
1
) nucleus, as we saw before, is in good agreement with the experimental data.The calculations of the muon capture rates are performed in three steps. In the first step we carry out realistic state-by-state calculations of exclusive ordinary muon capture (OMC) rates inZn
66 isotope for all multipolarities with J
π
≤
5
± (higher multipolarities contribute negligibly). The appropriate expression for the exclusive muon capture rates is written as follows:(6)
Λ
g
s
→
J
f
π
≡
Λ
J
f
π
=
2
G
2
Φ
1
s
2
R
f
q
f
2
×
J
f
π
M
^
J
-
L
^
J
0
g
s
+
2
+
J
f
π
T
^
J
el
-
T
^
J
magn
0
g
s
+
2
,where Φ
1
s represents the muon wave function in the 1
s muonic orbit. The operators in (6) are refered to as Coulomb M
^
J, longitudinal L
^
J, transverse electric T
^
J
el, and transverse magnetic T
^
J
mag multipole operators (see Appendix A). The factor R
f in (6) takes into consideration the nuclear recoil which is written as R
f
=
(
1
+
q
f
/
M
targ
)
-
1, with M
targ being the mass of the target (parent) nucleus.Due to the fact that there are no available data in the literature for exclusive muon capture rates, the test of our method is realized by comparing partial and total muon capture rates with experimental data and other theoretical results [23, 24]. Towards this purpose, our second step includes calculations of the partial μ
--capture rates for various low-spin multipolarities, Λ
J
π (for J
π
≤
4
±), in the studied nucleus. These partial rates are found by summing over the contributions of all the individual multipole states of the studied multipolarity as follows: (7)
Λ
J
π
=
∑
f
‍
Λ
g
s
→
J
f
π
=
2
G
2
Φ
1
s
2
∑
f
‍
q
f
2
R
f
J
f
π
M
^
J
-
L
^
J
0
g
s
+
2
k
k
k
k
k
k
k
k
k
k
k
k
+
∑
f
‍
q
f
2
R
f
J
f
π
T
^
J
el
-
T
^
J
magn
0
g
s
+
2(f runs over all states of the multipolarity J
π). We also estimate the percentage (portion) of their contribution into the total μ-capture rate for the most important multipolarities. In Table 5 we tabulate the individual portions of the low-spin multipole transitions (J
π
=
4
±) into the total muon capture rate. As it can be seen, the contribution of the 1
- multipole transitions is the most important multipolarity exhausting about 44
% of the total muon capture rate. Such an important contribution was found in O
16 and Ca
48 isotopes studied in [41].Table 5
The percentage of each multipolarity into the total muon capture rate evaluated with our pn-QRPA method.
Positive parity transitions
Negative parity transitions
J
π
Portions(
%
)
J
π
Portions(
%
)
0
+
8.22
0
-
7.94
1
+
21.29
1
-
44.21
2
+
2.85
2
-
13.32
3
+
1.58
3
-
0.34
4
+
0.01
4
-
0.23In the last step of testing our method, we evaluate total muon capture rates for theZn
66 isotope. These rates are obtained by summing over all partial multipole transition rates (up to J
π
=
4
±) as follows:(8)
Λ
tot
=
∑
J
π
‍
Λ
J
π
=
∑
J
π
‍
∑
f
‍
Λ
J
f
π
.For the sake of comparison, the abovementionedμ-capture calculations have been carried out using the quenched value g
A
=
1.135 [23, 24]. The results are listed in Table 6, where we also include the experimental total rates and the theoretical ones of [23, 24]. Moreover, in Table 6 we show the individual contribution into the total muon capture rate of the polar-vector (Λ
tot
V), the axial-vector (Λ
tot
A), and the overlap (Λ
tot
V
A) parts. As it can be seen, our results obtained with the quenched g
A coupling constant are in very good agreement with the experimental total muon capture rates (the deviations from the corresponding experimental rates are smaller than 7
%). This agreement provides us with high confidence level for our method.Table 6
Individual contribution of polar-vector, axial-vector, and overlap parts into the total muon-capture rate. The total muon capture rates, obtained by using the pn-QRPA with the quenched value ofg
A
=
1.135 for the medium-weight nucleus 66Zn, are compared with the available experimental data and with the theoretical rates of [23, 24].
Total muon capture ratesΛ
tot(×106 s
-
1)
Present pn-QRPA calculations
Experiment
Other theoretical methods
Nucleus
Λ
tot
V
Λ
tot
A
Λ
tot
V
A
Λ
tot
Λ
tot
exp
⁡
Λ
tot
theor [23]
Λ
tot
theor [24]
66Zn
1.651
4.487
−0.204
5.934
5.809
4.976
5.809


## 3.2. Electron Capture Cross Section

After acquiring a high confidence level for our nuclear method, we proceed with the main goal of the present study which concerns the calculations of the electron capture cross sections. As mentioned before, this includes original (see Section3.2.1) and stellar electron capture investigations (Section 3.2.2).


### 3.2.1. Original Electron Capture Cross Section on66Zn Isotope

The original cross sections for the electron capture process in theZn
66 isotope are obtained by using the pn-QRPA method considering all the accessible transitions of the final nucleus Cu
66. In the Donnelly-Walecka formalism the expression for the differential cross section in electron capture by nuclei reads [10](9)
d
σ
e
c
d
Ω
=
G
F
2
cos
⁡
2
θ
c
2
π
F
Z
,
E
e
2
J
i
+
1
H
×
∑
J
≥
1
‍
W
E
e
,
E
ν
H
H
i
×
1
-
α
cos
⁡
Φ
+
b
sin
2
Φ
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
×
J
f
T
^
J
mag
J
i
2
+
J
f
T
^
J
el
J
i
2
H
H
H
H
i
-
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
×
2
Re
J
f
T
^
J
mag
J
i
J
f
T
^
J
el
J
i
∗
H
H
i
+
∑
J
≥
0
‍
W
E
e
,
E
ν
H
H
i
i
i
i
i
i
i
×
1
+
α
cos
⁡
Φ
J
f
M
^
J
J
i
2
ω
q
1
+
α
cos
⁡
Φ
+
d
H
H
H
H
i
i
i
i
i
i
+
1
+
α
cos
⁡
Φ
-
2
b
sin
2
Φ
J
f
L
^
J
J
i
2
H
H
H
H
i
i
i
i
i
i
-
ω
q
1
+
α
cos
⁡
Φ
+
d
2
Re
J
f
L
^
J
J
i
∑
J
≥
1
‍
H
H
H
H
i
i
i
i
i
ω
q
1
+
α
cos
⁡
Φ
+
d
×
J
f
M
^
J
J
i
∗
,where F
(
Z
,
E
e
) is the well-known Fermi function [18]. The factor W
(
E
e
,
E
ν
)
=
E
ν
2
/
(
1
+
E
ν
/
M
T
) accounts for the nuclear recoil [8], M
T is the mass of the target nucleus, and parameters α, b, and d are given, for example, in [31]. The nuclear transition matrix elements between the initial state J
i and the final state |
J
f
〉 correspond to the Coulomb M
^
J
M, longitudinal L
^
J
M, transverse electric T
^
J
M
el, and transverse magnetic T
^
J
M
mag multipole operators (discussed in Appendix A).From the energy conservation in the reaction (1), the energy of the outgoing neutrino E
ν is written as follows:(10)
E
ν
=
E
e
-
Q
+
E
i
-
E
f
,which includes the difference between the initial E
i and the final E
f nuclear states. The Q value of the process is determined from the experimental masses of the parent (M
i) and the daughter (M
f) nuclei as Q
=
M
f
-
M
i [9].It is worth mentioning that for low momentum transfer, various authors use the approximationq
→
0 for all multipole operators of (9). Then, the transitions of the Gamow-Teller operator (G
T
+
=
∑
i
‍
τ
i
+
σ
i) provide the dominant contribution to the total cross section [9].While performing detailed calculations for the original electron capture cross sections onZn
66 isotope we assumed that (i) the initial state of the parent nucleus Zn
66 is the ground state 0
+ and (ii) the nuclear system is under laboratory conditions (no temperature dependence of the cross sections is needed). The cross sections as a function of the incident electron energy E
e are calculated with the use of realistic two-body interactions as mentioned before. The obtained total original electron capture cross sections for Zn
66 target nucleus are illustrated in Figure 2 where the individual contributions of various multipole channels (J
π
≤
5
±) are also shown. The electron capture cross sections in Figure 2 exhibit a sharp increase by several orders of magnitude within the first few MeV above energy-threshold, and this reflects the G
T
+ strength distribution. For electron energy E
e
≥
10
MeV the calculated cross sections show a moderate increase. From experimental and astrophysical point of view, the important range of the incident electron energy E
e is up to 30 MeV. At these energies the 1
+ multipolarity has the largest contribution to the total electron capture cross sections [9, 10]. In the present work we have extended the range of E
e up to 50 MeV since at higher energies (around 40 MeV) the contribution of other multipolarities like 1
-, 0
+, and 0
- becomes noticeable and cannot be omitted (see Figure 2).Figure 2
Original total cross sections of electron capture on theZn
66 (parent) nucleus calculated with pn-QRPA method as a function of the incident electron energy E
e. The individual contributions of various multipole channels (for J
π
≤
5
±) are also demonstrated.From the study of the original electron capture cross sections we conclude that the total cross sections can be well approximated with the Gamow-Teller transitions only in the region of low energies [9–12, 15, 16]. For higher incident electron energies the inclusion of the contributions originating from other multipolarities leads to better agreement [10].


### 3.2.2. Stellar Electron Capture on66Zn Isotope

As it is well known, electron capture process plays a crucial role in late stages of evolution of a massive star, in presupernova and in supernova phases [1–6]. In presupernova collapse, that is, at densities ρ
≤
1
0
10
g
cm
-
3 and temperatures 300
keV
≤
T
≤
800
keV, electrons are captured by nuclei with A
≤
60 [9–12, 15, 16]. During the collapse phase, at higher densities ρ
≥
1
0
10
g
cm
-
3 and temperatures T
≃
1.0
MeV, electron capture process is carried out on heavier and more neutron rich nuclei with Z
<
40 and N
≥
40 [8, 10, 13–15].In an independent particle picture, the Gamow-Teller transitions (which is the most important in the electron capture cross section calculations) are forbidden for these nuclei [2–4]. However, as it has been demonstrated in several studies, GT transitions in these nuclei are unblocked by finite temperature excitations [21, 22]. At high temperatures, T
≃
1.5
MeV, GT transitions are thermally unblocked as a result of the excitation of neutrons from the pf-shell into the g
9
/
2 orbital.For astrophysical environment, where the finite temperature and the matter density effects cannot be ignored (the initial nucleus is at finite temperature), in general, the initial nuclear state needs to be a weighted sum over an appropriate energy distribution. Then, assuming Maxwell-Boltzmann distribution of the initial statei in (9) [9, 26], the total e
--capture cross section is given by the expression [10](11)
σ
E
e
,
T
=
G
F
2
cos
⁡
2
θ
c
2
π
∑
i
‍
F
Z
,
E
e
2
J
i
+
1
e
-
E
i
/
k
T
G
Z
,
A
,
T
×
∑
f
,
J
‍
E
e
-
Q
+
E
i
-
E
f
2
i
O
^
J
f
2
2
J
i
+
1
.The sum over initial states in the latter equation denotes a thermal average of levels, with the corresponding partition function G
(
Z
,
A
,
T
) [10]. The finite temperature induces the thermal population of excited states in the parent nucleus. In the present work we assume that these excited states in the parent nucleus are all the possible states up to about 2.5 MeV. Calculations involving in addition other states lying at higher energies show that they have no sizeable contribution to the total electron capture cross sections. As mentioned before, for the evaluation of the total electron capture cross sections, the use of a quenched value of g
A is necessary [23, 24, 50, 51]. Since the form factor F
A
(
q
2
) multiplies the four components of the axial-vector operator (see (A.2)–(A.5)), a quenched value of g
A must enter the multipole operators generating the pronounced excitations 0
-
,
1
±
,
…, and so forth. For this reason, in our QRPA calculations we multiplied the free nucleon coupling constant g
A
=
1.262 by the factor 0.8 [23, 24, 50, 51].The results coming out of the study of electron capture cross sections under stellar conditions are shown in Figure3, where the same picture as in the original cross section calculations, but now with larger contribution, is observed. As discussed before, the dominant multipolarity is the 1
+, which contributes by more than 40
% to the total cross section. In the region of low energies (up to 30 MeV), the total e
--capture cross section can be described by taking into account only the GT transitions, but at higher incident energies the contributions of other multipolarities become significant and cannot be omitted.Figure 3
Electron capture cross sections for theZn
66 parent nucleus at high temperature (T
=
0.5 MeV) in stellar environment obtained assuming Maxwell-Boltzmann statistics for the incident electrons. The total cross section and the dominant individual multipole channels (J
π
≤
5
±) are demonstrated as functions of the incident electron energy E
e.The percentage contributions of various multipolarities (withJ
π
≤
5
±) into the total e
--capture cross section at T
=
0.5 MeV and for incident electron energy E
e
=
25
MeV are tabulated in Table 7. In addition, in this table we list the values of the individual e
--capture cross sections of each multipolarity with J
π
≤
5
±. More specifically, for E
e
=
25
MeV the 1
+ multipolarity contributes to about 44
%, the 0
+ contributes to about 26
%, and the 1
- contributes to about 11
%. The contributions coming from other multipolarities are less important (smaller than 5
%).Table 7
Totale
--capture cross sections (in 1
0
-
42 MeV−1 cm2) for E
e
=
25 MeV. The percentage of each multipolarity into the total e
--capture cross section evaluated with our pn-QRPA code is also tabulated here.
Positive parity transitions
Negative parity transitions
J
π
σ
e
(×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
J
π
σ
e (×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
0
+
31.164
25.96
0
-
5.288
4.41
1
+
52.779
43.98
1
-
13.409
11.14
2
+
6.921
5.77
2
-
3.262
2.72
3
+
5.499
4.58
3
-
0.905
0.75
4
+
0.244
0.20
4
-
0.299
0.25
5
+
0.208
0.17
5
-
0.042
0.04In performing state-by-state calculations for the electron capture cross sections, our code has the possibility to provide separately the contribution of the polar-vector, the axial-vector, and the overlap parts induced by the corresponding components of the electron capture operators. In Figure4 we illustrate the stellar differential cross sections of each individual transition of the polar-vector and axial-vector components.Figure 4
Individual contributions of the polar-vector, (Λ
V), and axial-vector, (Λ
A), components and total electron capture rate as functions of the excitation energy ω (Zn
66 is the parent nucleus).As mentioned before, our code gives separately the partiale
--capture cross sections of each multipolarity. In order to study the dependence of the differential cross sections on the excitation energy ω throughout the entire pn-QRPA spectrum of the daughter nucleus, a rearrangement of all possible excitations ω in ascending order, with the corresponding cross sections, is required. This was performed by using a special code appropriate for matrices [33]. In the model space chosen for Zn
66 isotope, for all multipolarities up to J
π
=
5
±, we have a number of 447 final states. The differential electron capture cross sections illustrated in Figure 4 present some characteristic clearly pronounced peaks at various excitation energies ω. These peaks correspond mainly to 0
+
,
1
+, and 2
+ transitions. More specifically, in the Cu
66 daughter nucleus the maximum peak corresponds to the 0
1
+ QRPA transition at ω
=
2.538
MeV and other characteristic peaks correspond to 1
7
+, 1
8
+, and 1
10
+ transitions, located at energies ω
=
3.194
MeV, ω
=
3.686
MeV, and ω
=
6.555
MeV, respectively (see Figure 4). The other less important peaks are also shown in Figure 4.Before closing, it should be mentioned that thee-capture cross sections presented in this work may be useful in estimating neutrino-spectra arising from e-capture on nuclei during supernova phase. The knowledge of ν-spectra at every point and time in the core is quite relevant for simulations of the final collapse and explosion phase of a massive star. As it is known [21], in the collapse phase, neutrinos are mainly produced by e-capture on nuclei and on free protons. The energy spectra of the emerging neutrinos from both reactions are important ingredients in stellar modelling and stellar simulations [21, 27].Furthermore, in core collapse simulations one defines the reaction rate of electron capture on nuclei given by(12)
R
h
=
∑
i
‍
Y
i
λ
i
,where the sum runs over all nuclear isotopes present in the astrophysical environment (Y
i denotes the abundance of a given nuclear isotope and λ
i is the calculated electron capture rate for this isotope). The rates of (12) must be known for a wide range of the parameters: T (temperature) and ρ (nuclear density) of the studied star. Thus, for the calculation of the quantity Y
·
λ of a specific nuclear isotope one needs to know in addition to the nuclear composition Y the electron capture rates λ calculated as we have shown in our present work. The rates of electron capture on various nuclear isotopes and the corresponding emitted neutrino spectra in the range of the parameters (T, ρ, Y
e) describing the star until reaching equilibrium during the core collapse are comprehensively studied in [21, 26, 27] for a great number of nuclear isotopes by using the large scale shell model. We are currently performing similar calculations for a set of isotopes by employing the present pn-QRPA method [52].Furthermore, the average neutrino energy,〈
E
ν
〉, of the neutrinos emitted by e-capture on nuclei can be obtained by dividing the neutrino-energy loss rate (defined by an expression similar to (12) by replacing the rate λ
i with the energy loss rate E
j) with the reaction rate for e-capture on nuclei R
h. Assuming, for example, power-law energy distribution for the neutrino spectrum produced by the e-capture in supernova phase, the average neutrino-energy 〈
E
ν
〉 determines a specific supernova-neutrino scenario. In addition, the neutrino emissivity is obtained by multiplying the electron capture rate at nuclear statistical equilibrium with the neutrino-spectra [21, 26, 27]. Finally we note that the rates for the inverse neutrino absorption process are also determined from the electron capture rates obtained as discussed in this section [21].


## 3.2.1. Original Electron Capture Cross Section on66Zn Isotope

The original cross sections for the electron capture process in theZn
66 isotope are obtained by using the pn-QRPA method considering all the accessible transitions of the final nucleus Cu
66. In the Donnelly-Walecka formalism the expression for the differential cross section in electron capture by nuclei reads [10](9)
d
σ
e
c
d
Ω
=
G
F
2
cos
⁡
2
θ
c
2
π
F
Z
,
E
e
2
J
i
+
1
H
×
∑
J
≥
1
‍
W
E
e
,
E
ν
H
H
i
×
1
-
α
cos
⁡
Φ
+
b
sin
2
Φ
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
×
J
f
T
^
J
mag
J
i
2
+
J
f
T
^
J
el
J
i
2
H
H
H
H
i
-
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
H
H
H
H
i
ɛ
i
+
ɛ
f
q
1
-
α
cos
⁡
Φ
-
d
×
2
Re
J
f
T
^
J
mag
J
i
J
f
T
^
J
el
J
i
∗
H
H
i
+
∑
J
≥
0
‍
W
E
e
,
E
ν
H
H
i
i
i
i
i
i
i
×
1
+
α
cos
⁡
Φ
J
f
M
^
J
J
i
2
ω
q
1
+
α
cos
⁡
Φ
+
d
H
H
H
H
i
i
i
i
i
i
+
1
+
α
cos
⁡
Φ
-
2
b
sin
2
Φ
J
f
L
^
J
J
i
2
H
H
H
H
i
i
i
i
i
i
-
ω
q
1
+
α
cos
⁡
Φ
+
d
2
Re
J
f
L
^
J
J
i
∑
J
≥
1
‍
H
H
H
H
i
i
i
i
i
ω
q
1
+
α
cos
⁡
Φ
+
d
×
J
f
M
^
J
J
i
∗
,where F
(
Z
,
E
e
) is the well-known Fermi function [18]. The factor W
(
E
e
,
E
ν
)
=
E
ν
2
/
(
1
+
E
ν
/
M
T
) accounts for the nuclear recoil [8], M
T is the mass of the target nucleus, and parameters α, b, and d are given, for example, in [31]. The nuclear transition matrix elements between the initial state J
i and the final state |
J
f
〉 correspond to the Coulomb M
^
J
M, longitudinal L
^
J
M, transverse electric T
^
J
M
el, and transverse magnetic T
^
J
M
mag multipole operators (discussed in Appendix A).From the energy conservation in the reaction (1), the energy of the outgoing neutrino E
ν is written as follows:(10)
E
ν
=
E
e
-
Q
+
E
i
-
E
f
,which includes the difference between the initial E
i and the final E
f nuclear states. The Q value of the process is determined from the experimental masses of the parent (M
i) and the daughter (M
f) nuclei as Q
=
M
f
-
M
i [9].It is worth mentioning that for low momentum transfer, various authors use the approximationq
→
0 for all multipole operators of (9). Then, the transitions of the Gamow-Teller operator (G
T
+
=
∑
i
‍
τ
i
+
σ
i) provide the dominant contribution to the total cross section [9].While performing detailed calculations for the original electron capture cross sections onZn
66 isotope we assumed that (i) the initial state of the parent nucleus Zn
66 is the ground state 0
+ and (ii) the nuclear system is under laboratory conditions (no temperature dependence of the cross sections is needed). The cross sections as a function of the incident electron energy E
e are calculated with the use of realistic two-body interactions as mentioned before. The obtained total original electron capture cross sections for Zn
66 target nucleus are illustrated in Figure 2 where the individual contributions of various multipole channels (J
π
≤
5
±) are also shown. The electron capture cross sections in Figure 2 exhibit a sharp increase by several orders of magnitude within the first few MeV above energy-threshold, and this reflects the G
T
+ strength distribution. For electron energy E
e
≥
10
MeV the calculated cross sections show a moderate increase. From experimental and astrophysical point of view, the important range of the incident electron energy E
e is up to 30 MeV. At these energies the 1
+ multipolarity has the largest contribution to the total electron capture cross sections [9, 10]. In the present work we have extended the range of E
e up to 50 MeV since at higher energies (around 40 MeV) the contribution of other multipolarities like 1
-, 0
+, and 0
- becomes noticeable and cannot be omitted (see Figure 2).Figure 2
Original total cross sections of electron capture on theZn
66 (parent) nucleus calculated with pn-QRPA method as a function of the incident electron energy E
e. The individual contributions of various multipole channels (for J
π
≤
5
±) are also demonstrated.From the study of the original electron capture cross sections we conclude that the total cross sections can be well approximated with the Gamow-Teller transitions only in the region of low energies [9–12, 15, 16]. For higher incident electron energies the inclusion of the contributions originating from other multipolarities leads to better agreement [10].


## 3.2.2. Stellar Electron Capture on66Zn Isotope

As it is well known, electron capture process plays a crucial role in late stages of evolution of a massive star, in presupernova and in supernova phases [1–6]. In presupernova collapse, that is, at densities ρ
≤
1
0
10
g
cm
-
3 and temperatures 300
keV
≤
T
≤
800
keV, electrons are captured by nuclei with A
≤
60 [9–12, 15, 16]. During the collapse phase, at higher densities ρ
≥
1
0
10
g
cm
-
3 and temperatures T
≃
1.0
MeV, electron capture process is carried out on heavier and more neutron rich nuclei with Z
<
40 and N
≥
40 [8, 10, 13–15].In an independent particle picture, the Gamow-Teller transitions (which is the most important in the electron capture cross section calculations) are forbidden for these nuclei [2–4]. However, as it has been demonstrated in several studies, GT transitions in these nuclei are unblocked by finite temperature excitations [21, 22]. At high temperatures, T
≃
1.5
MeV, GT transitions are thermally unblocked as a result of the excitation of neutrons from the pf-shell into the g
9
/
2 orbital.For astrophysical environment, where the finite temperature and the matter density effects cannot be ignored (the initial nucleus is at finite temperature), in general, the initial nuclear state needs to be a weighted sum over an appropriate energy distribution. Then, assuming Maxwell-Boltzmann distribution of the initial statei in (9) [9, 26], the total e
--capture cross section is given by the expression [10](11)
σ
E
e
,
T
=
G
F
2
cos
⁡
2
θ
c
2
π
∑
i
‍
F
Z
,
E
e
2
J
i
+
1
e
-
E
i
/
k
T
G
Z
,
A
,
T
×
∑
f
,
J
‍
E
e
-
Q
+
E
i
-
E
f
2
i
O
^
J
f
2
2
J
i
+
1
.The sum over initial states in the latter equation denotes a thermal average of levels, with the corresponding partition function G
(
Z
,
A
,
T
) [10]. The finite temperature induces the thermal population of excited states in the parent nucleus. In the present work we assume that these excited states in the parent nucleus are all the possible states up to about 2.5 MeV. Calculations involving in addition other states lying at higher energies show that they have no sizeable contribution to the total electron capture cross sections. As mentioned before, for the evaluation of the total electron capture cross sections, the use of a quenched value of g
A is necessary [23, 24, 50, 51]. Since the form factor F
A
(
q
2
) multiplies the four components of the axial-vector operator (see (A.2)–(A.5)), a quenched value of g
A must enter the multipole operators generating the pronounced excitations 0
-
,
1
±
,
…, and so forth. For this reason, in our QRPA calculations we multiplied the free nucleon coupling constant g
A
=
1.262 by the factor 0.8 [23, 24, 50, 51].The results coming out of the study of electron capture cross sections under stellar conditions are shown in Figure3, where the same picture as in the original cross section calculations, but now with larger contribution, is observed. As discussed before, the dominant multipolarity is the 1
+, which contributes by more than 40
% to the total cross section. In the region of low energies (up to 30 MeV), the total e
--capture cross section can be described by taking into account only the GT transitions, but at higher incident energies the contributions of other multipolarities become significant and cannot be omitted.Figure 3
Electron capture cross sections for theZn
66 parent nucleus at high temperature (T
=
0.5 MeV) in stellar environment obtained assuming Maxwell-Boltzmann statistics for the incident electrons. The total cross section and the dominant individual multipole channels (J
π
≤
5
±) are demonstrated as functions of the incident electron energy E
e.The percentage contributions of various multipolarities (withJ
π
≤
5
±) into the total e
--capture cross section at T
=
0.5 MeV and for incident electron energy E
e
=
25
MeV are tabulated in Table 7. In addition, in this table we list the values of the individual e
--capture cross sections of each multipolarity with J
π
≤
5
±. More specifically, for E
e
=
25
MeV the 1
+ multipolarity contributes to about 44
%, the 0
+ contributes to about 26
%, and the 1
- contributes to about 11
%. The contributions coming from other multipolarities are less important (smaller than 5
%).Table 7
Totale
--capture cross sections (in 1
0
-
42 MeV−1 cm2) for E
e
=
25 MeV. The percentage of each multipolarity into the total e
--capture cross section evaluated with our pn-QRPA code is also tabulated here.
Positive parity transitions
Negative parity transitions
J
π
σ
e
(×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
J
π
σ
e (×1
0
-
42 cm2 MeV−1 )
Portions(
%
)
0
+
31.164
25.96
0
-
5.288
4.41
1
+
52.779
43.98
1
-
13.409
11.14
2
+
6.921
5.77
2
-
3.262
2.72
3
+
5.499
4.58
3
-
0.905
0.75
4
+
0.244
0.20
4
-
0.299
0.25
5
+
0.208
0.17
5
-
0.042
0.04In performing state-by-state calculations for the electron capture cross sections, our code has the possibility to provide separately the contribution of the polar-vector, the axial-vector, and the overlap parts induced by the corresponding components of the electron capture operators. In Figure4 we illustrate the stellar differential cross sections of each individual transition of the polar-vector and axial-vector components.Figure 4
Individual contributions of the polar-vector, (Λ
V), and axial-vector, (Λ
A), components and total electron capture rate as functions of the excitation energy ω (Zn
66 is the parent nucleus).As mentioned before, our code gives separately the partiale
--capture cross sections of each multipolarity. In order to study the dependence of the differential cross sections on the excitation energy ω throughout the entire pn-QRPA spectrum of the daughter nucleus, a rearrangement of all possible excitations ω in ascending order, with the corresponding cross sections, is required. This was performed by using a special code appropriate for matrices [33]. In the model space chosen for Zn
66 isotope, for all multipolarities up to J
π
=
5
±, we have a number of 447 final states. The differential electron capture cross sections illustrated in Figure 4 present some characteristic clearly pronounced peaks at various excitation energies ω. These peaks correspond mainly to 0
+
,
1
+, and 2
+ transitions. More specifically, in the Cu
66 daughter nucleus the maximum peak corresponds to the 0
1
+ QRPA transition at ω
=
2.538
MeV and other characteristic peaks correspond to 1
7
+, 1
8
+, and 1
10
+ transitions, located at energies ω
=
3.194
MeV, ω
=
3.686
MeV, and ω
=
6.555
MeV, respectively (see Figure 4). The other less important peaks are also shown in Figure 4.Before closing, it should be mentioned that thee-capture cross sections presented in this work may be useful in estimating neutrino-spectra arising from e-capture on nuclei during supernova phase. The knowledge of ν-spectra at every point and time in the core is quite relevant for simulations of the final collapse and explosion phase of a massive star. As it is known [21], in the collapse phase, neutrinos are mainly produced by e-capture on nuclei and on free protons. The energy spectra of the emerging neutrinos from both reactions are important ingredients in stellar modelling and stellar simulations [21, 27].Furthermore, in core collapse simulations one defines the reaction rate of electron capture on nuclei given by(12)
R
h
=
∑
i
‍
Y
i
λ
i
,where the sum runs over all nuclear isotopes present in the astrophysical environment (Y
i denotes the abundance of a given nuclear isotope and λ
i is the calculated electron capture rate for this isotope). The rates of (12) must be known for a wide range of the parameters: T (temperature) and ρ (nuclear density) of the studied star. Thus, for the calculation of the quantity Y
·
λ of a specific nuclear isotope one needs to know in addition to the nuclear composition Y the electron capture rates λ calculated as we have shown in our present work. The rates of electron capture on various nuclear isotopes and the corresponding emitted neutrino spectra in the range of the parameters (T, ρ, Y
e) describing the star until reaching equilibrium during the core collapse are comprehensively studied in [21, 26, 27] for a great number of nuclear isotopes by using the large scale shell model. We are currently performing similar calculations for a set of isotopes by employing the present pn-QRPA method [52].Furthermore, the average neutrino energy,〈
E
ν
〉, of the neutrinos emitted by e-capture on nuclei can be obtained by dividing the neutrino-energy loss rate (defined by an expression similar to (12) by replacing the rate λ
i with the energy loss rate E
j) with the reaction rate for e-capture on nuclei R
h. Assuming, for example, power-law energy distribution for the neutrino spectrum produced by the e-capture in supernova phase, the average neutrino-energy 〈
E
ν
〉 determines a specific supernova-neutrino scenario. In addition, the neutrino emissivity is obtained by multiplying the electron capture rate at nuclear statistical equilibrium with the neutrino-spectra [21, 26, 27]. Finally we note that the rates for the inverse neutrino absorption process are also determined from the electron capture rates obtained as discussed in this section [21].


## 4. Summary and Conclusions

The electron capture on nuclei plays crucial role during the presupernova and collapse phase (in the late stagee
--capture on free protons is also significant). It becomes increasingly possible as the density in the star’s center is enhanced and it is accompanied by an increase of the chemical potential (Fermi energy) of the degenerate electron gas. This process reduces the electron-to-baryon ratio Y
e of the matter composition.In this work, by using our numerical approach based on a refinement of the pn-QRPA that describes reliably all the semileptonic weak interaction processes in nuclei, we studied in detail the electron capture process onZn
66 isotope and calculated original and stellar e
--capture cross sections. We tested our nuclear model (the pn-QRPA) through the reproducibility of orbital muon capture rates for this isotope. The agreement with experimental data and other reliable theoretical results of partial and total μ-capture rates and of the percentage contributions of various low-lying excitations is quite good which provides us with high confidence level for the obtained cross sections.Our future plans are to extend the application of this method and make similar calculations for other interesting nuclei [52]. Also this method could be applied to other semileptonic nuclear processes like beta-decay and charged-current neutrino-nucleus processes important in nuclear astrophysics and neutrino nucleosynthesis.



---

*Source: 398796-2015-05-13.xml*