#### Abstract

Neutrinos from supernovae constitute important probes of both the currently unknown supernova mechanisms and of neutrino properties. Reliable information about the nuclear responses to supernova neutrinos is therefore crucial. In this work, we compute the cross sections for the charged-current neutrino-nucleus scattering off the even-even molybdenum isotopes. The nuclear responses to supernova neutrinos are subsequently calculated by folding the cross sections with a Fermi-Dirac distribution.

#### 1. Introduction

Studies of neutrino-nucleus interactions with neutrinos of low and intermediate incoming energies (i.e., neutrino energies in the range up to a few hundred MeV) are of great importance for the disentanglement of various unresolved questions in astrophysics, particle physics, and nuclear physics [1, 2]. Supernovae (type II) constitute the inevitable deaths of very massive stars initiated by the collapse of their iron cores. For the details of supernova physics, see for example, [3–5]. In a supernova event, about 99% of the explosion energy is emitted as neutrinos of all flavors. A future detection of neutrinos from a coming supernova would therefore provide a wealth of valuable information both on the currently unknown supernova mechanisms and the associated nucleosynthesis of heavy elements. A high-statistics observation of the neutrino signal from the next nearby supernova could in addition provide important information on unknown neutrino properties (see e.g., [6]). It has for example, been shown by simulations [7] that the signal produced by supernova neutrinos in a large-scale detector such as the IceCube [8] probably could be used to disentangle the important question of normal-or-inverted neutrino-mass hierarchy. As has been proposed by Volpe [9, 10], charged-current neutrino-nucleus scattering with neutrinos from low-energy neutrino beams could in the future be exploited for spectroscopy of the virtual states involved in neutrinoless double-beta decay (see e.g., [11]) and consequently constrain theoretical predictions for the associated nuclear matrix elements.

Neutrinos from astrophysical sources can be detected by Earth-bound detectors by using charged-current and/or neutral-current interactions in nuclei [12]. For a recent review on experimental methods for detection of supernova neutrinos, see [13]. One planned possibility for measurements of astrophysical neutrinos (solar and supernova neutrinos) is the MOON (Mo Observatory Of Neutrinos) experiment [14].

In this paper, we perform theoretical calculations of the cross sections for the charged-current neutrino-nucleus scattering off the stable even-even molybdenum isotopes for neutrino energies, which are relevant for supernova neutrinos. We also present results for the averaged cross sections obtained by folding the cross sections with a Fermi-Dirac distribution. The calculations are based on the general theory for the treatment of semileptonic processes in nuclei, which was first introduced by O'Connell et al. [15]. This theory is discussed comprehensively in [16]. In the present calculations, the initial and final nuclear states are constructed by using the proton-neutron quasiparticle random-phase approximation (pnQRPA, see e.g., [17]).

The nuclear-structure dependence of the -nucleus cross sections is contained in the reduced nuclear matrix elements of various one-body operators, which depend on the momentum transfer. With an increasing number of nuclear final states, the calculations of these matrix elements are obviously increasingly time consuming. Therefore, in [18] we introduced a fast method for the calculation of the involved matrix elements, which is based on the barycentric Lagrange interpolation [19]. This method is adopted in this work and it is shown that the obtained results are in very good agreement with exact calculations.

This paper is organized as follows. In Section 2 we outline the theoretical framework for calculations of charged-current neutrino-nucleus cross sections. First the pnQRPA is introduced. We then briefly summarize the formalism for computations of charged-current neutrino-nucleus scattering off nuclei. Then, in Section 3 we present our results. Finally, in Section 4 we draw the conclusions.

#### 2. Theory

##### 2.1. pnQRPA

In this work, the proton-neutron quasiparticle random-phase approximation (pnQRPA) is adopted to construct the nuclear final and initial states. In this section, we therefore briefly summarize the formalism of the pnQRPA. For a more detailed treatment, the reader is referred to [17, 20].

The nuclear Hamiltonian for a general two-body residual interaction can be expressed in the form where the index contains the single-particle quantum numbers , , and and holds in addition the magnetic quantum number . Here is the particle creation operator and denotes the corresponding particle annihilation operator. In (2.1) are single-particle energies and the antisymmetric two-body matrix elements are defined as .

The pairing correlations in this work are taken into account by adopting the BCS (Bardeen-Cooper-Schrieffer) theory. The starting point in the BCS approach is the variational ansatz for the BCS vacuum where is the Hartree-Fock vacuum and denotes the time-reversed particle annihilation operator, which is defined as where ). The variational parameters and in (2.2) are obtained by solving the BCS equations [17].

The quasiparticles subsequently are defined via the Bogoliubov-Valatin transformation By using this transformation the nuclear Hamiltonian (2.1) can be cast into the form where the term is proportional to a normal ordered product of quasiparticle creation operators and quasiparticle annihilation operators.

In practice, the parameters of the BCS calculation are adjusted to reproduce the experimental pairing gaps for the reference even-even nucleus under consideration. In this work, the empirical pairing gaps have been deduced from the three-point formulae (see e.g., [17]) where () is the proton (neutron) separation energy for the nucleus having mass number and proton number .

In the pnQRPA, the states of the odd-odd nucleus subsequently are constructed by coupling two-quasiparticle operators to good angular momentum and parity . The state vector corresponding to the excitation then is defined by where denotes the pnQRPA vacuum and the additional quantum number enumerates the states with the same angular momentum and parity. Here the pnQRPA creation operator is given by where the sum runs over all possible proton-neutron configurations in the adopted valence space. The pnQRPA equations can then be written in the matrix form [17] where denotes the pnQRPA energy of the phonon . In (2.8) is the well-known pnQTDA matrix and the matrix contains the induced ground-state correlations. Matrices and contain both a particle-particle and a particle-hole channel of the proton-neutron two-body interaction. As is customary [11] also in this work the particle-particle and particle-hole channels are scaled by phenomenological constants and , respectively. This is done for each multipole separately. For more details on the scaling of the pnQRPA Hamiltonian, see [20, 21]. In the present computations, these parameters are tuned for the multipole only by using experimental observables such as beta-decay rates, the energy of the giant Gamow-Teller resonance, and experimental energies of the low-lying states in the odd-odd nucleus under consideration.

##### 2.2. Charged-Current Neutrino-Nucleus Scattering

In this work, we consider charged-current neutrino and antineutrino scatterings off a nucleus () with mass number and proton number , that is, neutrino-induced reactions of either the form or for the lepton flavors . The Feynman diagram for the -nucleus scattering that proceed, via an exchange of a boson is shown in Figure 1. In the figure , () represents the four-momentum of the incoming (outgoing), lepton and and are the four-momenta of the initial and final nuclei. The antineutrino reaction (2.10), which is mediated by a boson, can be obtained via charge-conjugation operation on Figure 1 and therefore is not displayed here.

In this work, we employ conventions that are similar to the ones used in [22]. We consequently define the covariant and contravariant four-vectors as and , respectively, with the metric . We thus write and , where and are the energies and and denote the three-momenta of the incoming and outgoing leptons. The four-momenta of the initial and final nuclear states take similarly the forms and , respectively.

For low-energy neutrinos, such as the supernova neutrinos that are considered in this work, the transferred four momentum is small compared to the mass of the gauge boson ( and resp.), that is, . The process in Figure 1 can consequently be treated in lowest order in the coupling constant . The effective Hamiltonian density can then be written for the neutrino-induced reaction (2.9) in the current-current form where represents the charged-current lepton current and is the hadron current. For charge-changing interactions, the coupling constant is given by where , GeV denotes the Fermi constant, and is the Cabibbo angle. For the antineutrino scattering, the effective Hamiltonian density similarly is given by where and . The nuclear matrix elements for both the considered processes then take the form where and for neutrino (antineutrino) scattering. We have here defined , and . This should be understood such that the appropriate hadron and lepton currents for the process under consideration are combined according to (2.11) and (2.12), respectively. In (2.13) we have moreover introduced the lepton matrix element Here the hadron currents can be decomposed into vector (V) and axial-vector (A) pieces according to At the origin , the single-nucleon matrix elements of the currents and are given by where denotes the state vector of a proton (neutron) having the three momentum and spin projection and is the nucleon mass. The corresponding one-nucleon matrix elements for the nuclear current are related to the ones of through the charge symmetries [23]

By the conservation of the vector current (CVC), the vector form factors in (2.16) can be written in terms of the proton and neutron electromagnetic form factors and as (see e.g., [23]) In the present work, we use the electromagnetic form factors of [22]. The axial-vector form factor in (2.16) is assumed to be of dipole form with the quenched static value . Moreover, the pseudoscalar form factor in (2.16) is obtained from the Goldberger-Treiman relation (see e.g., [16]) where [5] represents the mass of the charged pion.

We assume in this work that the final and initial nuclear states have well-defined angular momenta and and parities and , respectively. The double-differential cross section for the charged-current neutrino-nucleus scattering with transition from an initial nuclear state to a final nuclear state can then be written in the form (see e.g., [24]) where Here the excitation energy of the final nuclear state is defined with respect to the initial nuclear state, that is, . In the above expressions, we have introduced where the magnitude of the three-momentum transfer is given by and denotes the rest mass of the outgoing lepton. In (2.20), the function , where is the proton number of the final nucleus, takes into account the distortion of the outgoing lepton wave function due to the interaction with the (final) nucleus and the + (−) sign is used in the case of neutrinos (antineutrinos). Furthermore, in (2.22) the minus sign is used for neutrinos and the plus sign for antineutrinos.

For the treatment of final-state interactions, we use the method introduced by Engel in [25]. Consequently, we define the effective momentum of the outgoing lepton within the nucleus as where the effective energy is given by Here is the value at the origin of the Coulomb potential produced by the final nucleus. In this work, we approximate the Coulomb field with that of a uniformly charged sphere and therefore where for leptons and for antileptons and denotes the nuclear radius. In the region where is small, that is, , one can, just as in beta-decay calculations, adopt a Fermi function of the form [26] where For larger values of the energy of the final lepton, we employ in the present work the so-called modified effective momentum approximation (EMA) introduced in [25]. In this approximation, the absolute value of the three momentum and the energy of the outgoing lepton are replaced by the effective values given by (2.25) and (2.26), respectively. The effective energy and momentum are also used in the calculation of the three-momentum transfer defined by (2.24).

Guided by [10, 25], we adopt in the present calculations for the neutrino interaction (2.9) the Coulomb treatment which, gives the smallest cross section. For the antineutrino reaction (2.10), the final-state interactions decrease the cross section. Therefore, we choose in this case the Coulomb treatment (Fermi function or modified EMA), which gives the largest cross section.

#### 3. Results and Discussion

In this section, we adopt the formalism presented in Section 2 to compute the cross sections for the charged-current neutrino and antineutrino scatterings off the even-even molybdenum isotopes.

The initial and final nuclear states are first constructed by using the pnQRPA (proton-neutron quasiparticle random-phase approximation) discussed in Section 2.1. In the calculations, the single-particle energies were generated by using the Coulomb-corrected Woods-Saxon potential with the parameters of [27]. For the nucleus, some of the single-particle energies close to the respective Fermi surfaces were then adjusted according to [28].

We employ the formalism discussed in Section 2.2 to calculate the double-differential cross sections (2.20) for all the final nuclear states separately. The total cross section as a function of the energy of the impinging neutrino then is computed by integrating over the scattering angle and adding up all the individual contributions coming from the final states. For the calculations of the involved nuclear matrix elements, we use the efficient method introduced in [18]. For more details, we refer to [29, 30].

In Table 1, we show our calculated cross sections for the charged-current neutrino-nucleus scattering off the even molybdenum isotopes as functions of the energy of the incoming neutrino. In the table, the cross sections are tabulated for a large set of neutrino energies which, are relevant for supernova neutrinos, that is, . The results for the antineutrino-induced reactions are similarly shown in Table 2. As is seen in the tables, both the neutrino and antineutrino cross sections increase strongly with increasing neutrino energy. We can also conclude that the neutrino cross sections increase significantly with increasing neutron number . For low energies of the incoming neutrino, the cross sections are dominated by Gamow-Teller-like transitions mediated by the operator and Fermi-like ones, which proceed via the operator . At zero momentum transfer (), these operators reduce to the usual Gamow-Teller and Fermi () operators. It is well known that the Gamow-Teller operator obeys the Ikeda sum rule [31] where () is the total Gamow-Teller strength for the () channel. The nuclei in this work have large neutron excess and hence the branch is small and therefore . Similarly, for the transitions to final states one has that . Furthermore, the energy threshold, , is decreasing from MeV () to MeV (). This explains the increase of the neutrino cross sections with increasing neutron number at low neutrino energies. Contrary to this, the antineutrino cross sections decrease with increasing . This is explained by the fact that the values are increasing significantly with increasing neutron number.

In the last column of Table 1, we compare our results for with those of [10]. It is seen in the table that our results agree well with the ones published in the aforementioned reference. Our computed neutrino cross sections for are also in qualitative agreement with the results displayed in Figure 1 of [14].

As already mentioned above, we use in the present calculations the method of barycentric Lagrange interpolation for the computations of the required nuclear matrix elements. Therefore, in Figure 2 we compare the cross sections computed with this method with exact calculations in the case of the neutrino scattering off . The results for the antineutrino reaction are similarly shown in Figure 3. It is seen in the figures that the results computed by the approximate method (based on numerical interpolation) are accurate to 4-5 significant digits. Consequently, the numerical errors are very small compared to the errors coming from, for example, uncertainties in the nuclear wave functions. The results for the other isotopes are similar and are thus not shown here.

From the experimental point-of-view, an interesting quantity is the flux-averaged cross section , which is obtained by folding the cross sections with an appropriate energy profile of the incoming neutrinos. In this work, we adopt a two-parameter Fermi-Dirac distribution where denotes the effective neutrino temperature and is the so-called pinching parameter. Here the constant normalizes the total flux to unity.

Because of interactions with the matter, the neutrinos undergo flavor conversions when they propagate through the star (see e.g., [32]). Consequently, the energy profiles of the neutrinos (or antineutrinos) which arrive at the Earth-bound detector are different from the initial ones. For simplicity, we assume in this work that the oscillations occur in the outer layers of the star. Furthermore, we neglect collective neutrino-neutrino interactions. The neutrino profiles for the neutrinos, which reach the detector can then be written in the forms [32] for neutrinos and antineutrinos, respectively. Here () is the energy profile of electron neutrinos (non-electron neutrinos) and correspondingly for the antineutrinos. We also point out that due to the large rest masses of the muon and tau (and their antiparticles), only supernova and can be detected by charged-current neutrino-nucleus scattering.

In Table 3 are shown the computed averaged cross sections for the supernova and induced scatterings off . In the table, () represents nonoscillating electron neutrinos (electron antineutrinos). We moreover denote with () and () the oscillating electron neutrinos (electron antineutrinos) for the normal and inverted hierarchies, respectively. According to (3.3) (see (3.4)), these neutrinos have energy profiles that are altered with respect to the initial ones by the flavor conversions . Here results are shown for the two sets of neutrino parameters shown in Table 4. Averaged cross sections for other energy profiles of the incoming neutrinos (antineutrinos) can be easily computed by using the original cross sections tabulated in Table 1 (Table 2). We conclude from Table 3 that the flavor transformations significantly increase the averaged cross sections both for neutrinos and antineutrinos. The antineutrino cross sections are also much smaller than the corresponding ones for neutrinos. This is explained by the large neutron excess for () since then the antineutrino-induced reactions are suppressed because of Pauli blocking. In the table, we also compare our results with the ones of [33] for the parameter set II of Table 4. The numbers computed in the aforementioned reference have here been transformed to take into account the oscillation effect (3.3) with the updated value of [34] adopted in the present calculations. The discrepancy between the two calculations for is about 16% and for the flavor converted neutrinos in the normal hierarchy the results are almost exactly the same. Here it should be noted that in the aforementioned reference the measured , , and strength distributions of [35] were adopted instead of those calculated by the use of the pnQRPA. We therefore conclude that the results are in very good agreement with each other despite the rather different methods used in the calculations.

The proposed MOON detector can be realized by either using enriched or natural molybdenum [33, 36]. Consequently, accurate estimates of the nuclear responses for all the stable molybdenum isotopes are important. Therefore, we show in Table 5 our calculated averaged cross sections for all the stable even-even () Mo nuclei. It is seen in the table that the values of the neutrino cross sections are increasing significantly with increasing neutron number, while for the antineutrino reaction (2.10) the trend is opposite. We also conclude again that the cross sections for the neutrino scattering are notably larger than the ones for the antineutrino channel. Computations of the cross sections of the charged-current neutrino-nucleus scattering off the odd isotopes, and , are still in progress and the results will be published elsewhere.

#### 4. Conclusions

In this work, we have computed the cross sections for the charged-current neutrino and antineutrino scatterings off the stable molybdenum isotopes for neutrino energies, which are relevant for supernova neutrinos. In the calculations, the initial and final nuclear states have been constructed by using the pnQRPA. We have tabulated the cross sections for an extensive set of neutrino energies. The nuclear responses for the aforementioned nuclei have subsequently been calculated by using a two-parameter Fermi-Dirac distribution. However, averaged cross sections for other energy profiles of the incoming neutrino (or antineutrino) can be easily computed by using the cross sections given in this paper.

Our results for the neutrino-nucleus scattering off agree well with previously published results. We have found for the studied chain of nuclei that the neutrino-scattering cross sections increase significantly with the neutron number . Contrary to this, the cross sections for the antineutrino scattering decrease notably with increasing neutron number.

#### Acknowledgments

This work was supported by the Academy of Finland under the Finnish Center Of Excellence Program 2012–2017 (Nuclear and Accelerator Based Program at JYFL). Dr. E. Ydrefors would like to thank Professor H. Ejiri for fruitful discussions.