Large-scale calculations of supernova neutrino-induced reactions in Z=8-82 target nuclei
Abstract
Background: In the environment of high neutrino-fluxes provided in core-collapse supernovae or neutron star mergers, neutrino-induced reactions with nuclei contribute to the nucleosynthesis processes. A number of terrestrial neutrino detectors are based on inelastic neutrino-nucleus scattering and modeling of the respective cross sections allow predictions of the expected detector reaction rates. Purpose: To provide a self-consistent microscopic description of neutrino-nucleus cross sections involving a large pool of nuclei for the implementation in models of nucleosynthesis and neutrino detector simulations. Methods: Self-consistent theory framework based on relativistic nuclear energy density functional is employed to determine the nuclear structure of the initial state and relevant transitions to excited states induced by neutrinos. The weak neutrino-nucleus interaction is employed in the current-current form and a complete set of transition operators is taken into account. Results: We perform large-scale calculations of charged-current neutrino-nucleus cross sections, including those averaged over supernova neutrino fluxes, for the set of even-even target nuclei from oxygen toward lead (), spanning (OPb pool). The model calculations include allowed and forbidden transitions up to multipoles. Conclusions: The present analysis shows that the self-consistent calculations result in considerable differences in comparison to previously reported cross sections, and for a large number of target nuclei the cross sections are enhanced. Revision in modeling r-process nucleosynthesis based on a self-consistent description of neutrino-induced reactions would allow an updated insight into the origin of elements in the Universe and it would provide the estimate of uncertainties in the calculated element abundance patterns.
pacs:
21.30.Fe, 21.60.Jz, 23.40.Bw, 25.30.-cI Introduction
Neutrino-induced reactions on nuclei play an important role in nuclear astropyhsics, in particular during core-collapse supernova evolution and nucleosynthesis Bru.91 (); Jan.07 (). Elastic neutrino scattering on nuclei and nucleons determines the neutrino trapping and the diffusion time scale of the outwards streaming neutrinos Pin.06 (). A variety of processes contribute to the energy loss in stellar interiors, e.g. pair, photo-, plasma, nucleon-nucleon bremsstrahlung, and recombination neutrino processes Ito.96 (). In addition, inelastic neutrino-electron scattering mainly thermalizes the neutrino spectra Mez.93 () and inelastic neutral-current neutrino-nucleus scattering has a large contribution to the opacity Ger.75 (); Ger.75b (), during the core collapse and subsequent explosion phase. Additionally, in the environment of exploding massive stars, which has long been explored as a possible site for the r-process nucleosynthesis, charged-current neutrino-nucleus reactions play an important role in the production of chemical elements. As pointed out in Ref. Hek.00 (), there are many interesting effects of neutrino-induced reactions before, during, and after the r-process. To explore their role in nucleosynthesis reliable neutrino-nucleus reaction rate compilations involving also nuclei with large neutron excess are required. Available data on neutrino-nucleus cross sections are limited to deuterium Wil.80 (), and C and Fe target nuclei, obtained by the LSND Ath.97 () and KARMEN Bod.94 (); Mas.98 () collaborations, and at LAMPF Kra.92 (). Therefore, only theoretical approaches can provide cross sections for a large number of target nuclei that are involved in various applications of neutrino physics and astrophysics. Modeling neutrino-nucleus interactions is also important in view of the current research and development of neutrino detectors, e.g., for supernova and solar neutrinos, neutrinos produced in laboratories, and geoneutrinos. The ongoing and planned neutrino detector facilities involve a variety of target materials, induced reactions and scientific objectives, e.g., MOON Eji.08 (), MiniBooNE Agu.09 (), MINOS Ada.06 (), SNO+ Chr.09 (), OPERA Aga.10 (), LVD (Large Volume Detector)Aga.07 (), ORLaND experiment proposal at the Spallation Neutron Source (SNS) Efr.05 (), NOvA neutrino experiment Har.05 (), Daya Bay reactor neutrino experiment An.12 (), etc. There is also interesting concept of beta-beams for the production of neutrinos by using -decay of boosted radioactive ions Zuc.02 (); Vol.04 (). The proposal to establish a beta-beam facility that could produce low-energy neutrino beams in the 100 MeV energy range would allow direct insight into the neutrino-induced reactions in nuclei as well as its underlying structure involving a non-trivial combination of nuclear allowed and forbidden transitions Vol.04 ().
At low neutrino energies, the neutrino-induced reactions are sensitive to the properties of nuclei involved, i.e. their initial and excited states. Therefore it is necessary to employ a microscopic framework providing reasonable description of nuclear structure properties. Over the past years, a variety of advanced microscopic models have been developed and employed in studies of charged-current neutrino-induced reactions at low energies, also including various particle decay channels. In particular, these include the nuclear shell model Hax.87 (); Eng.96 (); Hay.00 (); Sam.02 (); Yos.08 (); Suz.09 (); Suz.11 (), random phase approximation (RPA) Hek.00 (); Aue.97 (); Sin.98 (); Vol.00 (); Vol.02 (); Eng.03 (); Ath.06 (), continuum RPA (CRPA) Kol.92 (); Kol.95 (); Jac.99 (); Jac.02 (); Bot.05 (), and the hybrid model which combines the shell model for allowed transitions, with the RPA to account for the forbidden transitions Kol.99 (); Toi.01 (); Kol.03 (); Suz.09 (). Frameworks based on quasiparticle RPA (QRPA) have also recently been developed, based on Skyrme functionals Vol.00 (); Laz.07 (), Brueckner G matrix employed for two-body interaction by solving the Bethe-Salpeter equation using Bonn-CD potential Che.10 (); Che2.10 (); Che.11 (), and projected QRPA Krm.05 (); Sam.08 (); Sam.11 (). The Fermi gas model Wal.75 (); Gai.86 (); Kur.90 () and Fermi liquid theory (FLT) Lyk.07 () have also been employed in studies of low-energy charged-current neutrino-nucleus reactions. At finite temperature in stellar environments, thermal population of the excited states may enhance the weak interaction rates and cross sections at low neutrino energies Dzh.10 (); Dzh.11 ().
Despite considerable progress in the development of advanced theoretical frameworks, up to this time only a limited number of microscopic models have been employed in large-scale calculations of neutrino-induced reactions and their implementation in supernova simulations. Because the complete modeling of neutrino-induced reactions necessitates the inclusion not only of Gamow-Teller transitions but also contributions from forbidden transitions and other higher multipoles, covering a large pool of target nuclei represents a computationally very demanding problem. Although the shell model provides a very accurate description of ground state wave functions, the description of high-lying excitations necessitates the use of large model spaces which often leads to computational difficulties, making the approach applicable essentially only to allowed transitions in light and medium-mass nuclei. For systematic studies of neutrino-nucleus cross sections throughout the nuclide chart including the heavy nuclei, microscopic calculations must therefore be performed using models based on the QRPA. The first global calculations of electron neutrino-nucleus cross sections have been conducted using the gross theory of -decay Ful.95 (); McL.95 (). Two microscopic frameworks have been employed in large-scale calculations of neutrino-induced reactions i) extended Thomas-Fermi plus Strutinsky integral (ETFSI) and continuum quasiparticle random phase approximation (CQRPA) Bor.00 (), and ii) RPA with Landau-Migdal (LM) force, using Wood-Saxon potential (WS) to determine the single-particle basis of target nuclei Lan.01 (). There are no up-to-date systematic calculations conducted in the framework of energy density functional.
Recently the framework based on relativistic nuclear energy density functional (RNEDF) has been introduced in modeling the charged-current neutrino-nucleus cross sections Paa.08 (), and it has also been extended for applications in modeling neutral-current neutrino-induced reactions Djapo.12 (). In the case of iron group nuclei, comparison of the results for charge exchange reactions obtained using Skyrme functionals and the RNEDF, as well as with the shell model, showed reasonable theoretical uncertainty inherent in modeling neutrino-nucleus cross sections. However, in view of applications in supernova and r-process simulations, which also involve unstable nuclei far from the valley of stability, where no experimental data are available, it is necessary to provide independent insight into relevant neutrino-induced processes from various models and effective interactions. The main objective of this work is conducting large-scale calculations of charged-current neutrino-nucleus cross sections in a large pool of nuclei from oxygen towards lead, based on the RNEDF. In addition to the overall cross sections covering the range of neutrino energies up to 100 MeV, calculations also include the cross sections averaged over neutrino fluxes for the range of temperatures characteristic of various stages of stellar evolution. In addition to the Fermi-Dirac neutrino spectra we also apply supernova spectra from a recent simulation that includes three-flavor Boltzmann neutrino transport. Particular aim of this work is to emphasize the role of forbidden transitions in modeling neutrino-nucleus cross sections in large pool of nuclei, including both stable and nuclei away from the valley of stability. The relevance of forbidden transitions has already been discussed in several studies, e.g. Refs. Sur.98 (); Vol.04 (); Laz.07 (). In the study of neutrino capture by r-process waiting point nuclei, first forbidden strength, together with the low-lying Gamow-Teller transitions, increased the rate of neutrino scattering from very neutron-rich nuclei by a factor of at least 2 and in some instances even by a factor of 5 Sur.98 (). As pointed out in microscopic calculations based on Skyrme functionals in Ref. Laz.07 (), the properties of forbidden states are closely related to the neutrino-nucleus cross sections, and could be extracted by using neutrinos from low-energy beta-beams.
The paper is organized as follows. Section II contains the basic theoretical background for the neutrino-nucleus cross sections in the charged-current channel based on weak Hamiltonian and the RNEDF. The results of large-scale calculations of the neutrino-nucleus cross sections are illustrated and discussed in Sec. III. The conclusions of the present work are summarized in Sec. IV.
Ii Theoretical background
In the present work we explore the charged-current neutrino-nucleus process,
(1) |
where the incoming electron neutrino induces charge-exchange reaction in target nucleus . The general formalism of the neutrino-nucleus cross sections, derived assuming the Hamiltonian of the weak interaction in the current-current form, is given in details in Refs. Con.72 (); Wal.75 (). The cross sections include the transition matrix elements between the initial and final nuclear state , for the charge , longitudinal , transverse electric , and transverse magnetic multipole operators Con.72 (). In the present work the RNEDF is employed in calculations of the transition matrix elements contributing to the neutrino-nucleus cross sections. The RNEDF has already been successfully employed in studies of giant resonances and exotic modes of excitation Vre.99 (); Vre.04 (); PVKC.07 (); PPP.05 (); Paa2.09 (); Kha.11 (), -decay rates of r-process nuclei Nik.05 (), muon capture Mar.09 () and stellar electron capture rates Niu.11 (), and in constraining the neutron skin in nuclei Vre.03 (); Kli.07 (). More details about the implementation of the RNEDF in modeling charged-current neutrino-nucleus reactions are given in Refs. Paa.08 (); Paa.11 ().
The RNEDF based framework employs the self-consistent mean field for nucleons and minimal set of meson fields; isoscalar scalar -meson , isoscalar vector -meson and the isovector vector -meson , supplemented with the electromagnetic field. The meson-nucleon interaction is included with a minimal set of the interaction terms, where the vertex functionals include explicit dependence on the vector density. The nuclear ground state properties are described using the relativistic Hartree-Bogoliubov model (RHB), and relevant transitions induced by neutrinos are calculated in the relativistic quasiparticle random phase approximation (RQRPA). More details on the RHB model based on effective density-dependent interactions are given in Ref. Nik.02 (). The RQRPA is formulated in the canonical single-nucleon basis of the RHB model Paa.03 (); NVR.02 (). In modeling the neutrino-nucleus cross sections, important advantage of this framework is that it is fully consistent in view of the effective interactions employed. In the particle-hole and pairing channels, the same interactions are used in the RHB equations that determine the canonical quasiparticle basis, and in the matrix equations of the RQRPA. In this way, one can employ the same RNEDF in description of the weak processes throughout the nuclide map without any additional adjustments to the properties of specific target nuclei under consideration. For the model parameters that determine the density-dependent couplings and the meson masses, in this work DD-ME2 parameterization is used Lal.05 (). The pairing correlations in open shell nuclei are described by the finite range Gogny interaction, with parameterization D1S Ber.91 ().
Complete calculation of inelastic neutrino-nucleus reactions spanning the range of neutrino energies up to 100 MeV necessitates the inclusion of a number of transitions with various multipoles Paa.08 (). Although higher-order multipoles have rather small contributions at low incoming neutrino energies, these can not be neglected at energies about tens of MeV. In the present study, multipoles up to have been taken into account. The large scale calculations of neutrino-nucleus cross sections, involving more than 1000 nuclei and a complete set of all multipoles up to with both parities, necessitate considerable computational effort. Therefore, for the purposes of the present work, computational framework has been developed using parallel computing methods based on Message Passing Interface (MPI) Gro.94 () for the implementation on cluster and grid computer systems.
Iii Results and discussion
By employing the model outlined in the previous section, we have conducted large-scale calculations of reactions in the OPb pool of even-even target nuclei spanning the range from oxygen toward lead (Z=8-82) with neutron number N=8-182. The calculations include excitations of all multipoles up to and both parities. In the following we explore the systematic behavior of the overall cross sections throughout the pool under consideration, including the cross sections averaged over the Michel neutrino flux obtained from the decay at rest (DAR) of Kra.92 ()
(2) |
To simulate the supernova neutrino spectrum, at first instance we use the Fermi-Dirac distribution
(3) |
In Fig. 1 the inclusive neutrino-nucleus cross sections, averaged over the Michel neutrino flux, are shown for the OPb pool of nuclei. Stable nuclei in the pool are denoted by filled circles. In addition, the calculated cross sections for C and Fe are especially emphasized in comparison to the KARMEN experimental data Bod.94 (); Mas.98 (). One can observe that the pool of nuclei calculated with the same energy density functional without any adjustments of the model parameters fits nicely into two experimental data points. The cross sections increase systematically with increasing neutron number, resulting in values larger up to a factor in comparison to those in the valley of stability. We note that sharp edge of the data at large denotes boundary values for the pool under consideration, obtained for Pb isotope chain.
Figure 2 shows the cross sections in the OPb pool of nuclei, averaged over the supernova neutrino flux given by Eq. (3) in the case of MeV and . The results can be compared to those of the stable target nuclei. While the cross sections display a similar pattern as in the previous case with the Michel spectrum, the details of the data, however, depend on of the neutrino spectrum. The cross sections for various groups of target nuclei are separately displayed in Fig. 2: stable nuclei, proton-rich nuclei constrained by , and neutron-rich nuclei with with . The cross sections exhibit a systematic behavior throughout the nuclide map. In the case of neutron rich nuclei, the cross sections are larger in comparison to the stable nuclei, due to the increased number of neutrons that participate in charge-exchange neutrino-induced reactions. For proton-rich nuclei, the reaction pattern is opposite due to þthe smaller number of neutron-proton configurations, i.e., the cross sections are considerably reduced in comparison to those of stable nuclei. This point is further illustrated in Fig. 3 where we plot the same cross sections versus the difference between the neutron and the proton numbers. For nuclei with the cross sections increase significantly with increasing number of excess neutrons. Although some scattering is apparent in the data, most of the results cluster along an almost linear function of the number of excess neutrons. Temperature dependence of neutrino-induced reactions within the OPb pool is illustrated in Fig. 4, where the cross sections averaged over the Fermi-Dirac distribution () are shown for stable target nuclei in the range MeV. At low temperatures, the cross sections appear rather scattered due to the stronger dependence on low-energy excitations in nuclei. However, for higher T, the cross sections are rather smooth and considerably larger due to the inclusion of a number of multipoles contributing at higher neutrino energies. The full set of calculations has been completed for the OPb pool in the range of temperatures MeV and complete tables are available on request.
In view of applications in astrophysical models and in predicting the detector response to neutrinos involving various target nuclei, it is crucial to assess the sensitivity of the neutrino-nucleus cross sections on the theoretical frameworks and effective interactions employed. In the case of iron group nuclei, it has been shown that by employing different microscopic models, one can estimate reasonable theoretical uncertainty in neutrino-nucleus cross sections averaged over the Michel spectrum, i.e., for Fe cm Paa.11 (). Since for nuclei far from the valley of stability various models result in larger differences in nuclear structure properties, one could also expect larger sensitivity in the results of modeling neutrino-induced reactions. In the present analysis theoretical uncertainties are assessed on the basis of the OPb pool of nuclei. Figures 5, 6 and 7 show the dependence of the flux-averaged neutrino-nucleus cross sections per nucleon ( MeV, ) on the neutron number in the cases of Ni, Sn, and Pb isotopic chains, respectively. In these figures we compare our results to the available sets of cross sections based on two other theoretical frameworks: ETFSI + CQRPA Bor.00 () and RPA (WS+LM) Lan.01 (), where the former contains only contributions from the Gamow-Teller and Fermi transitions, while the latter includes all transitions up to . The results of the present RQRPA analysis for are shown separately for the full calculation including all transitions up to , and partial cross sections obtained only for the isobaric analog and Gamow-Teller transitions.
The RQRPA results for nuclei, obtained taking into account only the Fermi and the Gamow-Teller transitions, are consistently higher than the corresponding values calculated using the ETFSI + CQRPA framework for all three isotopic chains. The deviation is rather small in the region around the valley of stability in Ni and Sn isotopes, but increases with additional neutrons. In the case of the Pb isotopic chain, the difference is significant even for the lightest isotopes considered in this study. For nuclei with the proton-to-neutron ratio , the ETFSI + CQRPA framework predicts a change of the general trend and an anomalous increase of the cross sections with a corresponding decrease of the number of neutrons. The agrement between results that take into accout forbidden transitions, RQRPA and RPA(WS+LM), is much better in all three isotopic chains, although for very neutron-rich Ni and Sn isotopes the RQRPA predicts higher values of the flux-averaged cross section. In the case of the Pb isotopic chain, the agreement between the RQRPA and RPA(WS+LM) cross sections is excellent for all the isotopes under consideration. There are several reasons for the shown deviations. Each model employs different effective interactions that result in variations of the excitation pattern contributing to the cross sections, and the present calculations are based on a fully self-consistent approach to the neutrino-induced reactions. The three models agree best for nuclei around the valley of stability, where most of the experimental data on nuclear structure properties are available. In the very neutron rich region the differences between various models increase and grow relatively large. However, in all three isotopic chains the agreement of the RQRPA results is much better with the values obtained using the RPA (WS + LM) model than with the ETFSI + CQRPA, partly indicating the importance of the forbidden transitions.
In the case of cross sections averaged over the supernova neutrino flux for (, ), the impact of the self-consistent and complete calculations including all relevant multipoles is explored in view of previous knowledge on neutrino induced reactions on a large-scale basis within the OPb pool. Figure 8 shows the ratio of the cross sections averaged over supernova neutrino flux of the present work and ETFSI + CQRPA Bor.00 () model. In order to identify nuclei with pronounced discrepancies, the following groups of nuclei are separately denoted in figure: stable nuclei, proton-rich nuclei defined by and neutron-rich nuclei with . For target nuclei where the cross sections are available for both models, the present large-scale calculations result mainly in systematically larger cross sections, for most of nuclei up to the factor of 4, and for smaller number of nuclei up to the factor of 7. However, for a limited set of medium mass nuclei the cross sections of this work are smaller than the ETFSI+CQRPA ones, and these are mainly limited to neutron deficient nuclei in the range . The reason is the anomalous increase of the cross sections in the ETFSI + CQRPA model, when moving along a particular isotope chain from nucleus toward proton rich nuclei (see Figs. 5 and 6). This behavior has not been observed in the present analysis. Apart from this anomaly, one can observe systematic linear increase of the ratio of cross sections with the nuclear mass. As already discussed in the case of isotopic chains, in addition to the self-consistent implementation of the RNEDF, the model employed in the present work includes a complete set of transition operators up to multipoles which enhance the overall cross sections. As shown in Fig. 8 the linear trend is particularly apparent in the case of neutron-rich nuclei (), where the absolute values of the cross sections are relatively large in comparison to other nuclei under consideration.
To further explore the role of forbidden transitions in the OPb pool, in Fig. 9 we plot the ratio of cross sections taking into account all multipoles up to and including with cross sections that only include the contributions from the Fermi and Gamow-Teller transitions. A trend very similar to the one in the previous figure appears, which indicates that the forbidden transitions contribute a larger fraction of the total cross sections in nuclei with higher mass. In heavier nuclei, and neutron-rich nuclei in particular, the difference between the proton and neutron numbers grows large, to the extent that in the very heavy nuclei neutrons occupy a full shell more above the protons. However, even in cases where there is not enough neutrons to occupy a full shell, additional neutrons move closer in energy to the orbits of the next shell, and therefore make the forbidden transitions less suppressed. The present analysis provides quantitative estimates for the mass dependence of the underlying structure in the neutrino-nucleus cross sections: in the mass region below the forbidden transitions contribute less than 10% of the total flux-averaged cross sections, but their contribution grows practically linear with mass and can provide up to 50% of the total cross section. Therefore, in order to provide a realistic description of neutrino-nucleus reactions, as well as in the corresponding astrophysical applications, models must take into account more than the simplest Fermi and Gamow-Teller terms.
In Fig. 10 the ratio of the averaged cross sections of this work and RPA (WS + LM) Lan.01 () is shown. Qualitative agreement between the two models is obtained and differences are within a factor of 2. The ratio exhibits a mild mass dependence, i.e. the model based on RQRPA provides larger cross sections mainly for nuclei up to , while for heavier systems the cross sections are up to 50% smaller when compared to the ones previously reported with the RPA (WS + LM) model. These discrepancies occur for particular isotopic chains, however, there is no global systematic over- or underestimation of the cross-sections with respect to the RPA (WS + LM) approach, indicating that the deviations mainly originate from the different nuclear ground state and the residual interaction. The lack of a linear trend for neutron-rich nuclei (in comparison to Figs. 9 and 8), along with the fact that RPA (WS+LM) model includes all transitions up to and including points to the importance of the lowest forbidden transitions in the quantitative description of semileptonic processes in nuclei. The largest deviations appear for low mass proton-rich nuclei, where the reactions become suppressed and the absolute values of the cross sections are actually very low, making them very sensitive to the structural details of the ground state. In comparison to the RPA (WS + LM) model, the main advantage of the present work is implementation of fully self-consistent theory framework in description of nuclear sector in neutrino induced reactions. In this way predictions for the cross sections become feasible not only in the region of stable nuclei and those that are covered by the RPA (WS+LM) model , but also in the unknown regions of the nuclide map of relevance for the stellar processes and nucleosynthesis.
In addition to the Fermi-Dirac spectrum Eq. (3), we also apply the distribution obtained from a core-collapse supernova simulation Fischer2012 (). It can be described by the -fit as follows
(4) |
where the parameter is determined from the mean neutrino energy and the root-mean-square of the neutrino energy as Keil2003 ()
(5) |
These data are taken from the core-collapse supernova explosion simulations of Ref. Fischer:2010 (). The supernova model is based on general relativistic radiation hydrodynamics and three-flavor Boltzmann neutrino transport in spherical symmetry. It also employs a nuclear equation of state Shen:1998 () and a modern set of neutrino opacities (for details, see Table 1 in Ref. Fischer2012 () and the references therein). The simulations were launched from massive progenitor stars, the 8.8 M O-Ne-Mg-core and the more massive Fe-core progenitors of 10.8 and 18 M. Note that while for the low-mass star neutrino driven explosions can be obtained in spherically symmetric models, neutrino heating/cooling was enhanced artificially in order to trigger the onset of explosion for the more massive stars (for details, see Ref. Fischer:2010 ()). For the current discussions, we select the data from the 18 M explosion model. Nevertheless, all simulations exhibit the same intrinsic feature of the continuously decreasing neutrino luminosities and average energies (see, e.g., Figure 14 in Ref. Fischer:2010 ()) after the supernova explosion has been launched. This important aspect is related to the continuous emission of neutrinos of all flavors, which deleptonizes the central protoneutron star on timescales of tens of seconds. Note that at the onset of the supernova explosion, the protoneutron star is hot (10–40 MeV) and still lepton rich, in which sense it differs from the final supernova remnant neutron star. In order to characterize the dynamical evolution of the neutrino spectra during the protoneutron star deleptonization, we select three different post bounce times after the explosion has been launched, , 5, and 20 seconds (see Table 1 for the corresponding neutrino energies).
Once the neutrinos have decoupled from matter at the sphere of last scattering, which depends on neutrino flavor and energy, they can still contribute to charge and neutral current processes involving heavy nuclei at large distance ( km) relevant for the nucleosynthesis, e.g., the process Froehlich:2006 () and the process Qian:2007 (). At conditions where these processes can occur, the zero-temperature approach employed in the current paper for the calculation of the cross sections is a valid approximation. The protoneutron star deleptonization with the continuously reducing average neutrino energies has important consequences for the neutrino flux-averaged cross sections, as shown in Fig. 11. In comparison to the Fermi-Dirac spectrum with MeV shown in Fig. 2, the supernova simulation spectrum results in significantly lower flux-averaged cross sections by more than a factor of two for all nuclei under investigation, even shortly ( s) after the onset of explosion. Note that the supernova neutrino temperature, which can be defined as the matter temperature at the neutrinosphere, is slightly higher with MeV than the assumed MeV of the Fermi-Dirac spectrum. However, the Fermi-Dirac neutrino spectrum is significantly broader than those taken from the supernova simulations, which can be seen by comparing mean neutrino energies and the mean-square energies of and . For comparison, we also plot the spectra in Fig. 12. The above mentioned difference between Fermi-Dirac spectra and supernova simulation spectra even increases during the ongoing explosion at late times, i.e. the supernova spectra shift increasingly towards Maxwell-Boltzmann like spectra. These findings are summarized in Table 1, comparing the Fermi-Dirac spectra () and the supernova simulation spectra characterized by the -fit () at different post-bounce times after the explosion has been launched.
^{1}^{1}1Post bounce simulation time for the supernova spectra | [MeV] | [MeV] | |
---|---|---|---|
12.60 | 207.03 | ||
1 | 9.32 | 108.57 | |
5 | 8.68 | 95.29 | |
20 | 6.73 | 58.32 |
Iv Conclusion
We have employed a fully self-consistent model based on the RNEDF in large-scale calculations of charged-current neutrino-nucleus cross sections in the OPb pool of target nuclei, spanning the range and . The two main advantages of this approach are: (i) self-consistent modeling of all relevant transition matrix elements involving open-shell nuclei, without any additional adjustments of the model parameters to the nuclear target under consideration, and (ii) treating the allowed and the forbidden transitions on an equal footing, i.e., transitions of all relevant multipoles are taken into account. For the purpose of the present work, advanced computational framework has been developed, based on a parallel computing scheme using the Message Passing Interface for the implementation on cluster and grid computer systems. Reasonable agreement of the cross sections with the only available experimental data, for C and Fe, support the course of systematic calculations throughout the nuclide map. The results include a complete set of the cross sections for the OPb pool, calculated for the range of neutrino energies MeV. The cross sections were averaged over the experimental Michel spectrum, the neutrino spectrum described by the Fermi-Dirac distribution for the range of temperatures MeV and chemical potential . In addition, the averaged cross sections have been calculated using -fit neutrino fluxes, taken from a recent core-collapse supernova explosion model that is based on three-flavor Boltzmann neutrino transport. Comparing Fermi-Dirac and supernova simulation neutrino spectra, we found that the latter result in significantly smaller flux-averaged cross sections due to the more pinched supernova neutrino spectra. This aspect even increases during the ongoing deleptonization on the order of 10–20 s after which the cross sections become negligible and charged-current processes with nuclei have ceased.
Note that the neutrino spectra form the supernova simulations do not take into account possible collective neutrino-flavor oscillations. These phenomena can take place after neutrino decoupling from matter and result in complete spectral swaps above a certain neutrino energy for inverted neutrino-mass hierarchy Dasgupta:2009mg (). Because heavy-lepton flavor neutrinos have higher average energies, being less strongly bound to matter in the absence of charged-current weak processes, it enhances the high-energy tail of the electron-flavor neutrino spectra. This in turn may impact neutrino-induced nucleosynthesis of heavy elements at large distance Martinez:11 (); Duan:11 ().
Furthermore, it has been realized that medium modifications for the charged-current weak processes with neutrons and protons must be taken into account when computing neutrino transport and neutrino decoupling from matter at high densities using a nuclear equation of state Reddy:98 (). These modifications have been explored at the mean-field level and shown to increase spectral differences between electron neutrinos and antineutrinos Martinez:12 (); Roberts:12 (); Horowitz:12 (). Improved supernova simulations that take these effects into account may ultimately alter the results discussed in the current paper.
Several key features of the neutrino-nucleus cross sections have been illustrated for various supernova neutrino fluxes, indicating systematic increase of the averaged cross sections with increase of the number of neutrons in target nuclei as well as with temperature of the neutrino distribution. When going away from the valley of stability toward neutron-rich nuclei, the cross sections become considerably enhanced, while on the proton-rich side they are suppressed due to the blocking of orbitals available for neutrino induced transitions.
The cross sections from the present analysis have been discussed in view of previously reported large scale calculations. Current results are consistently higher than the results obtained using the ETFSI + CQRPA model which includes only the Fermi and the Gamow-Teller transitions. For lighter nuclei the RQRPA cross sections are approximately 50% larger, while for the heaviest nuclei studied they are up to 3 times larger. This enhancement originates in part from the forbidden transitions that are fully taken into account in the present work. We provide important quantitative insight into the underlying structure of the neutrino-nucleus cross sections throughout the nuclide map: in the mass region below the forbidden transitions contribute less than 10% of the total flux-averaged cross sections, but their contribution increases with mass and can provide up to 50% of the total cross section. For the limited set of nuclei, the differences between the present results and those of the RPA (LM + WS) model are constrained within a factor of three. There is no general trend in the cross sections ratio between the two models, and the differences arise mainly due to the different approach in the calculations of the nuclear ground state and the residual interaction in (Q)RPA.
We have shown that the self-consistent calculations based on the RNEDF, which take into account transitions of all relevant multipoles, result in differences in comparison to previously reported cross sections, and for a considerable number of target nuclei the cross sections are larger. Therefore, revised calculations in modeling stellar evolution and r-process nucleosynthesis based on self-consistent descriptions of neutrino-induced reactions may allow an updated insight into the origin of elements in the Universe. On the other hand, by employing various theory frameworks for neutrino-nucleus cross sections, one could estimate the uncertainties in the calculated element abundance patterns. The tables with the cross sections vs. neutrino energies of this work are available from the authors on request.
ACKNOWLEDGMENTS
We gratefully acknowledge S. Goriely for providing data for the ETFSI+CQRPA calculations from Ref. Bor.00 (). This work is supported by MZOS - project 1191005-1010, the Croatian Science Foundation and by the Helmholtz International Center for FAIR within the framework of the LOEWE program launched by the State of Hesse. The computational resources have been provided in part by the Croatian National Grid Infrastructure (CRO-NGI). T.F. is supported by the Swiss National Science Foundation under project no. PBBSP2-133378.
References
- (1) S.W. Bruenn and W.C. Haxton, Astrophys. J. 376, 678 (1991).
- (2) H.-Th. Janka, K. Langanke, A. Marek, G. Martínez-Pinedo, and B. Müller, Phys. Rep. 442, 38 (2007).
- (3) G. Martínez-Pinedo, M. Liebendörfer, and D. Frekers, Nucl. Phys. A 777, 395 (2006).
- (4) N. Itoh, H. Hayashi, A. Nishikawa, and Y. Kohyama, Astrophys. J. Supp. 102, 411 (1996).
- (5) A. Mezzacappa and S. W. Bruenn, Astrophys. J. Lett. 410, 740 (1993).
- (6) S. S. Gershtein, V. N. Folomeschkin, M.Yu.Khlopov, and R.A.Eramzhyan, J. Exp. Theor. Phys. 42, 571 (1975).
- (7) S. S. Gershtein, V. S. Imshennik, V. N. Folomeschkin, M. Yu. Khlopov, V. M. Chechetkin, and R. A. Eramzhyan, J. Exp. Theor. Phys. 42, 751 (1975).
- (8) A. Hektor, E. Kolbe, K. Langanke, and J. Toivanen, Phys. Rev. C 61, 055803 (2000).
- (9) S. E. Willis et al., Phys. Rev. Lett. 44, 522 (1980).
- (10) C. Athanassopoulos et al., Phys. Rev. C 55, 2078 (1997).
- (11) B. E. Bodmann et al., Phys. Lett. B 332, 251 (1994).
- (12) R. Maschuw, Prog. Part. Nucl. Phys. 40, 183 (1998).
- (13) D. A. Krakauer et al., Phys. Rev. C 45, 2450 (1992).
- (14) H. Ejiri et al., Eur. Phys. J. Special Topics 162, 239 (2008).
- (15) A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Nucl. Instrum. Methods A 599, 28 (2009).
- (16) P. Adamson et al., Phys. Rev. D 73, 072002 (2006).
- (17) C. Kraus (SNO+ collaboration), Prog. Part. Nucl. Phys. 64, 273 (2010).
- (18) N. Agafonova et al., Phys. Lett. B 691, 138 (2010).
- (19) N. Yu Agafonova et al., Astron. Phys. 27, 254 (2007).
- (20) Y. Efremenko, Nucl. Phys. B, Proc. Suppl. 138, 343 (2005); F. T. Avignone III and Y. V. Efremenko, J. Phys. G 29, 2615 (2003).
- (21) D. A. Harris (MINOS NOvA Collaborations), Nucl. Phys. B, Proc. Suppl. 149, 150 (2005).
- (22) F. P. An et al., Phys. Rev. Lett. 108, 171803 (2012).
- (23) P. Zucchelli, Phys. Lett. B 532, 166 (2002).
- (24) C. Volpe, J. Phys. G 30, L1 (2004).
- (25) W. C. Haxton, Phys. Rev. D 36, 2283 (1987).
- (26) J. Engel, E. Kolbe, K. Langanke, and P. Vogel, Phys. Rev. C 54, 2740 (1996).
- (27) A. C. Hayes and I. S. Towner, Phys. Rev. C 61, 044603 (2000).
- (28) J. M. Sampaio, K. Langanke, G. Martinez-Pinedo, and D. J. Dean, Phys. Lett. B 529, 19 (2002).
- (29) T. Yoshida, T. Suzuki, S. Chiba, T. Kajino, H. Yokomakura, K. Kimura, A. Takamura, and D. H. Hartmann, Astrophys. J. 686, 448 (2008).
- (30) T. Suzuki, M. Honma, K. Higashiyama, T. Yoshida, T. Kajino, T. Otsuka, H. Umeda, and K Nomoto, Phys. Rev. C 79, 061603(R) (2009).
- (31) T. Suzuki, M. Honma, T. Yoshida, H. Mao, T. Kajino, and T. Otsuka, Prog. Part. Nucl. Phys. 66, 385 (2011).
- (32) N. Auerbach, N. Van Giai, and O. K. Vorov, Phys. Rev. C 56, R2368 (1997).
- (33) S. K. Singh, N. C. Mukhopadhyay, and E. Oset, Phys. Rev. C 57, 2687 (1998).
- (34) C. Volpe, N. Auerbach, G. Colò, T. Suzuki, and N. Van Giai, Phys. Rev. C 62, 015501 (2000).
- (35) C. Volpe, N. Auerbach, G. Colò, and N. Van Giai, Phys. Rev. C 65, 044603 (2002).
- (36) J. Engel, G. C. McLaughlin, and C. Volpe, Phys. Rev. D 67, 013005 (2003).
- (37) M. S. Athar, S. Ahmad, and S. K. Singh, Nucl. Phys. A 764, 551 (2006).
- (38) E. Kolbe, K. Langanke, S. Krewald, and F.-K. Thielemann, Nucl. Phys. A 540, 599 (1992).
- (39) E. Kolbe, K. Langanke, F.-K. Thielemann, and P. Vogel, Phys. Rev. C 52, 3437 (1995).
- (40) N. Jachowicz, S. Rombouts, K. Heyde, and J. Ryckebusch, Phys. Rev. C 59, 3246 (1999).
- (41) N. Jachowicz, K. Heyde, J. Ryckebusch, and S. Rombouts, Phys. Rev. C 65, 025501 (2002).
- (42) A. Botrugno and G. Co’, Nucl. Phys. A 761, 200 (2005).
- (43) E. Kolbe, K. Langanke, and P. Vogel, Nucl. Phys. A 652, 91 (1999).
- (44) J. Toivanen, E. Kolbe, K. Langanke, G. Martinez-Pinedo, and P. Vogel, Nucl. Phys. A 694, 395 (2001).
- (45) E. Kolbe, K. Langanke, G. Martinez-Pinedo, and P. Vogel, J. Phys. G 29, 2569 (2003).
- (46) R. Lazauskas and C. Volpe, Nucl. Phys. A 792, 219 (2007).
- (47) Myung-Ki Cheoun, Eunja Ha, K. S. Kim, and Toshitaka Kajino, J. Phys. G 37, 055101 (2010).
- (48) Myung-Ki Cheoun, Eunja Ha, T. Hayakawa, T. Kajino, and S. Chiba, Phys. Rev. C 82, 035504 (2010).
- (49) Myung-Ki Cheoun, Eunja Ha, and T. Kajino, Phys. Rev. C 83, 028801 (2011).
- (50) F. Krmpotic, A. Samana, and A. Mariano, Phys. Rev. C 71, 044319 (2005).
- (51) A. R. Samana and C. A. Bertulani, Phys. Rev. C 78, 024312 (2008).
- (52) A. R. Samana, F. Krmpotic, N. Paar, and C. A. Bertulani, Phys. Rev. C 83, 024303 (2011).
- (53) J. D. Walecka Muon Physics, edited by V. M. Hughes and C. S. Wu, Academic, New York (1975).
- (54) J. S. O’Connel, T. W. Donnelly, and J. D. Walecka, Phys. Rev. C 6, 719 (1972).
- (55) T. K. Gaisser and J. S. O’Connell, Phys. Rev. D 34, 822 (1986).
- (56) T. Kuramoto, M. Fukugita, Y. Kohyama, and K. Kubodera, Nucl. Phys. A 512, 711 (1990).
- (57) G. I. Lykasov and V. A. Bednyakov, Phys. Rev. C 76, 014622 (2007).
- (58) A. A. Dzhioev, A. I. Vdovin, V. Yu. Ponomarev, J. Wambach, K. Langanke, and G. Martinez-Pinedo, Phys. Rev. C 81, 015804 (2010).
- (59) A. A. Dzhioev, A. I. Vdovin, V. Yu. Ponomarev, and J. Wambach, Phys. At. Nucl. 74, 1162 (2011).
- (60) G. M. Fuller and B. S. Meyer, Astrophys. J. 453, 792 (1995).
- (61) G. McLaughlin and G. M. Fuller, Astrophys. J. 455, 202 (1996).
- (62) I. N. Borzov and S. Goriely, Phys. Rev. C 62, 035501 (2000).
- (63) K. Langanke and E. Kolbe, At. Data. Nucl. Data Tables 79, 293 (2001).
- (64) N. Paar, D. Vretenar, T. Marketin, and P. Ring, Phys. Rev. C 77, 024608 (2008).
- (65) H. Đapo, N. Paar, Phys. Rev. C 86, 035804 (2012).
- (66) R. Surman and J. Engel, Phys. Rev. C 58, 2526 (1998).
- (67) N. Paar, T. Suzuki, M. Honma, T. Marketin, and D. Vretenar, Phys. Rev. C 84, 047305 (2011).
- (68) D. Vretenar, P. Ring, G. A. Lalazissis, and N. Paar, Nucl. Phys. A 649, 29c (1999).
- (69) D. Vretenar, T. Nikšić, N. Paar, and P. Ring, Nucl. Phys. A 731, 281 (2004).
- (70) N. Paar, D. Vretenar, E. Khan, and G. Colo, Rep. Prog. Phys. 70, 691 (2007).
- (71) N. Paar, P. Papakonstantinou, V. Y. Ponomarev, and J. Wambach, Phys. Lett. B 624, 195 (2005).
- (72) N. Paar, Y. F. Niu, and D. Vretenar, J. Meng, Phys. Rev. Lett. 103, 032502 (2009).
- (73) E. Khan, N. Paar, and D. Vretenar, Phys. Rev. C 84, 051301 (2011).
- (74) T. Nikšić, T. Marketin, D. Vretenar, N. Paar, and P. Ring, Phys. Rev. C 71, 014308 (2005).
- (75) T. Marketin, N. Paar, T. Nikšić, and D. Vretenar, Phys. Rev. C 79, 054323 (2009).
- (76) Y. F. Niu, N. Paar, D. Vretenar, and J. Meng, Phys. Rev. C 83, 045807 (2011).
- (77) D. Vretenar, N. Paar, T. Nikšić, and P. Ring, Phys. Rev. Lett. 91, 262502 (2003).
- (78) A. Klimkiewicz et al., Phys. Rev. C 76, 051603 (2007).
- (79) T. Nikšić, D. Vretenar, P. Finelli, and P. Ring, Phys. Rev. C 66, 024306 (2002).
- (80) N. Paar, P. Ring, T. Nikšić, and D. Vretenar, Phys. Rev. C 67, 034312 (2003).
- (81) T. Nikšić, D. Vretenar, and P. Ring, Phys. Rev. C 66, 064302 (2002).
- (82) G. A. Lalazissis, T. Nikšić, D. Vretenar, and P. Ring, Phys. Rev. C 71, 024312 (2005).
- (83) J. Berger, M. Girod, and D. Gogny, Comput. Phys. Commun. 63, 365 (1991).
- (84) W. Gropp, E. Lusk, and A. Skjellum, Using MPI: portable parallel programming with the message-passing interface. Cambridge, MA, USA: MIT Press Scientific And Engineering Computation Series. ISBN 0-262-57104-8 (1994).
- (85) T. Fischer, G. Martínez-Pinedo, M. Hempel, and M. Liebendörfer, Phys. Rev. D 85, 083003 (2012).
- (86) M. T. Keil, G. G. Raffelt, and H.-T. Janka, Astrophys. J. 590, 971 (2003).
- (87) T. Fischer, S. C. Whitehouse, A. Mezzacappa, F.-K. Thilelemann, and M. Liebendörfer, Astronom. Astrophys. 517, A80 (2010).
- (88) H. Shen, H. Toki, K Oyamatsu, and K. Sumiyoshi, Nucl. Phys. A637, 435 (1998).
- (89) C. Fröhlich, G. Martínez-Pinedo, M. Liebendörfer, F.-K. Thielemann, E. Bravo, W. R. Hix, K. Langanke, and N. T. Zinner, Phys. Rev. Lett. 96, 142502 (2006).
- (90) Y.-Z. Qian and G. J. Wasserberg, Phys. Rep. 442, 1 (2007).
- (91) B. Dasgupta, A. Dighe, G. G. Raffelt and A. Y. Smirnov, Phys. Rev. Lett. 103, 051105 (2009)
- (92) G. Martínez-Pinedo, B. Ziebarth, T. Fischer and K. Langanke, European Physical Journal A 47, 98, (2011)
- (93) H Duan, A. Friedland, G. G. McLaughlin and R. Surman, Phys. G Nucl. Part. Phys. 38, 035201 (2011)
- (94) S. Reddy, M. Prakash, and J. M. Lattimer, Phys. Rev. D 58, 013009 (1998)
- (95) G. Martínez-Pinedo, T. Fischer, A. Lohs, and L. Huther, Phys. Rev. Lett. 109, 251104 (2012)
- (96) L. Roberts, S. Reddy, and G. Shen, eprints, astro-ph.HE/1205.4066 (2012)
- (97) C. J. Horowitz, G. Shen, E. O’Connor, and C. D. Ott, eprints, astro-ph.HE/1209.3173 (2012)