The $ u$-process with fully time-dependent supernova neutrino emission spectra

The neutrino process that occurs in the outer stellar shells during a supernova explosion and involves neutrino-nucleus reactions produces a range of rare, stable and radioactive isotopes. We improve previous $ u$-process studies by using, for the first time, the time-dependent neutrino emission spectra, as predicted from supernova simulations, rather than a simplified parametric description modeled after the neutron-star cooling phase. In particular, our calculations use time-dependent neutrino spectra for all neutrino species, consider their deviation from a Fermi-Dirac distribution and account for the neutrino emission from the neutrino burst and accretion phases. We find that the time-dependent treatment of the neutrino emission spectra results in higher yields for the selected nuclei produced by the $ u$~process as compared to previous studies and also compared to the approximation of assuming constant neutrino energies corresponding to the time-averaged mean energy radiated in each species. The effect is largest for nuclides produced by charged-current reactions. Our results reflect the dynamical competition between neutrino-induced reactions and the effect of the shock passage through the star. By varying the neutrino burst luminosity and the duration of the accretion phase, we study the impact of these early emission phases and their uncertainties on the $ u$-process nucleosynthesis. We find that the deviation of the neutrino spectra from a Fermi-Dirac distribution calculated in supernova simulations has a negligible effect on the $ u$-process yields.

Neutrino nucleosynthesis: An overview

Neutrinos produced during a supernova explosion induce reactions on abundant nuclei in the outer stellar shells and contribute in this way to the synthesis of the elements in the Universe. This neutrino nucleosynthesis process has been identified as an important contributor to the origin of $^7$Li, $^{11}$B,$^{19}$F, $^{138}$La, and $^{180}$Ta, but also to the long-lived radionuclides $^{22}$Na and $^{26}$Al, which are both key isotopes for $gamma$-ray astronomy. The manuscript summarizes the recent progress achieved in simulations of neutrino nucleosynthesis.

On the nuclear robustness of the r process in neutron-star mergers

We have performed r-process calculations for matter ejected dynamically in neutron star mergers based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic simulation. Our calculations consider an extended nuclear network, including spontaneous, $beta$- and neutron-induced fission and adopting fission yield distributions from the ABLA code. We have studied the sensitivity of the r-process abundances to nuclear masses by using different models. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the abundance yields at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number $N=130$. The remaining trajectories, which contribute 10% by mass to the total integrated abundances, follow such a fast expansion that the r process does not use all the neutrons. This also leads to a larger variation of abundances among trajectories as fission does not dominate the r-process dynamics. The total integrated abundances are dominated by contributions from the slow abundances and hence reproduce the general features of the observed r-process abundances. We find that at timescales of weeks relevant for kilonova light curve calculations, the abundance of actinides is larger than the one of lanthanides. Hence actinides can be even more important than lanthanides to determine the photon opacities under kilonova conditions. (Abridged)

Neutrino-pair emission from nuclear de-excitation in core-collapse supernova simulations

We study the impact of neutrino-pair production from the de-excitation of highly excited heavy nuclei on core-collapse supernova simulations, following the evolution up to several 100 ms after core bounce. Our study is based on the AGILE-Boltztran supernova code, which features general relativistic radiation hydrodynamics and accurate three-flavor Boltzmann neutrino transport in spherical symmetry. In our simulations the nuclear de-excitation process is described in two different ways. At first we follow the approach proposed by Fuller and Meyer [Astrophys. J. 376,701 (1991)], which is based on strength functions derived in the framework of the nuclear Fermi-gas model of non-interacting nucleons. Secondly, we parametrize the allowed and forbidden strength distributions in accordance with measurements for selected nuclear ground states. We determine the de-excitation strength by applying the Brink hypothesis and detailed balance. For both approaches, we find that nuclear de-excitation has no effect on the supernova dynamics. However, we find that nuclear de-excitation is the leading source for the production of electron antineutrinos as well as heavy-lepton flavor (anti)neutrinos during the collapse phase. At sufficiently high densities, the associated neutrino spectra are influenced by interactions with the surrounding matter, making proper simulations of neutrino transport important for the determination of the neutrino-energy loss rate. We find that even including nuclear de-excitations, the energy loss during the collapse phase is overwhelmingly dominated by electron neutrinos produced by electron captures.

Microscopic Nuclear Structure and Reaction Calculations in the FMD Approach

We present here a first application of the Fermionic Molecular Dynamics (FMD) approach to low-energy nuclear reactions, namely the $^3$He($alpha$,$gamma$)$^7$Be radiative capture reaction. We divide the Hilbert space into an external region where the system is described as $^3$He and $^4$He clusters interacting only via the Coulomb interaction and an internal region where the nuclear interaction will polarize the clusters. Polarized configurations are obtained by a variation after parity and angular momentum projection procedure with respect to the parameters of all single particle states. A constraint on the radius of the intrinsic many-body state is employed to obtain polarized clusters at desired distances. The boundary conditions for bound and scattering states are implemented using the Bloch operator. The FMD calculations reproduce the correct energy for the centroid of the $3/2^-$ and $1/2^-$ bound states in $^7$Be. The charge radius of the ground state is in good agreement with recent experimental results. The FMD calculations also describe well the experimental phase shift data in the $1/2^+$, $3/2^+$ and $5/2^+$ channels that are important for the capture reaction at low energies. Using the bound and scattering many-body wave functions we calculate the radiative capture cross section. The calculated $S$ factor agrees very well, both in absolute normalization and energy dependence, with the recent experimental data from the Weizmann, LUNA, Seattle and ERNA experiments.

Towards Microscopic Ab Initio Calculations of Astrophysical S-Factors

Low energy capture cross sections are calculated within a microscopic many-body approach using an effective Hamiltonian derived from the Argonne V18 potential. The dynamics is treated within Fermionic Molecular Dynamics (FMD) which uses a Gaussian wave-packet basis to represent the many-body states. A phase-shift equivalent effective interaction derived within the Unitary Correlation Operator Method (UCOM) that treats explicitly short-range central and tensor correlations is employed. As a first application the 3He(alpha,gamma)7Be reaction is presented. Within the FMD approach the microscopic many-body wave functions of the 3/2- and 1/2- bound states in 7Be as well as the many-body scattering states in the 1/2+, 3/2+ and 5/2+ channels are calculated as eigenstates of the same microscopic effective Hamiltonian. Finally the S-factor is calculated from E1 transition matrix elements between the many-body scattering and bound states. For 3He(alpha,gamma)7Be the S-factor agrees very well, both in absolute normalization and energy dependence, with the recent experimental data from the Weizmann, LUNA, Seattle and ERNA experiments. For the 3H(alpha,gamma)7Li reaction the calculated S-factor is about 15% above the data.

Complete inclusion of parity-dependent level densities in the statistical description of astrophysical reaction rates

Microscopic calculations show a strong parity dependence of the nuclear level density at low excitation energy of a nucleus. Previously, this dependence has either been neglected or only implemented in the initial and final channels of Hauser-Feshbach calculations. We present an indirect way to account for a full parity dependence in all steps of a reaction, including the one of the compound nucleus formed in a reaction. To illustrate the impact on astrophysical reaction rates, we present rates for neutron captures in isotopic chains of Ni and Sn. Comparing with the standard assumption of equipartition of both parities, we find noticeable differences in the energy regime of astrophysical interest caused by the parity dependence of the nuclear level density found in the compound nucleus even at sizeable excitation energies.