We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We
have studied the material ejected from four NS-NS merger systems based on the hydrodynamical simulations of Ardevol-Pulpillo et al. (2019) which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r-process elements with mass numbers $A gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models.
Carbon-enhanced metal-poor (CEMP) r/s-stars show surface-abundance distributions characteristic of the so-called intermediate neutron capture process (i-process) of nucleosynthesis. We previously showed that the ingestion of protons in the convective
helium-burning region of a low-mass low-metallicity star can explain the surface abundance distribution observed in CEMP r/s stars relatively well. Such an i-process requires detailed reaction network calculations involving hundreds of nuclei for which reaction rates have not yet been determined experimentally. We investigate the nuclear physics uncertainties affecting the i-process during the AGB phase of low-metallicity low-mass stars by propagating the theoretical uncertainties in the radiative neutron capture cross sections, as well as the 13C(a,n)16O reaction rate, and estimating their impact on the surface-abundance distribution. It is found that considering systematic uncertainties on the various nuclear ingredients affecting the radiative neutron capture rates, surface elemental abundances are typically predicted within +/-0.4 dex. The 56 < Z < 59 region of the spectroscopically relevant heavy-s elements of Ba-La-Ce-Pr as well as the r-dominated Eu element remain relatively unaffected by nuclear uncertainties. In contrast, the inclusion of the direct capture contribution impacts the rates in the neutron-rich A~45, 100, 160, and 200 regions, and the i-process production of the Z~45 and 65-70 elements. Uncertainties in the photon strength function also impact the overabundance factors by typically 0.2-0.4 dex. Nuclear level densities tend to affect abundance predictions mainly in the Z=74-79 regions. The uncertainties associated with the neutron-producing reaction 13C(a,n)16O and the unknown beta-decay rates are found to have a low impact on the overall surface enrichment
Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A >140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncer
tainties. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-rich fissioning nuclei of relevance in r-process calculations. These improvements include a phenomenological modification of the scission distance and a smoothing procedure of the distribution. Such corrections lead to a much better agreement with experimental fission yields. Those yields are also used to estimate the number of neutrons emitted by the excited fragments on the basis of different neutron evaporation models. Our new fission yields are extensively compared to those predicted by the so-called GEF model. The impact of fission on the r-process nucleosynthesis in binary neutron mergers is also reanalyzed. Two scenarios are considered, the first one with low initial electron fraction is subject to intense fission recycling, in contrast to the second one which includes weak interactions on nucleons. The various regions of the nuclear chart responsible for fission recycling during the neutron irradiation as well as after freeze-out are discussed. The contribution fission processes may have to the final abundance distribution is also studied in detail in the light of newly defined quantitative indicators describing the fission recycling, the fission seeds and the fission progenitors. In particular, those allow us to estimate the contribution of fission to the final abundance distribution stemming from specific heavy nuclei. Calculations obtained with SPY and GEF fission fragment distributions are compared for both r-process scenarios.
The radiative neutron capture rates for isotopes of astrophysical interest are commonly calculated within the statistical Hauser-Feshbach reaction model. Such an approach, assuming a high level density in the compound system, can be questioned in lig
ht and neutron-rich nuclei for which only a few or no resonant states are available. Therefore, in this work we focus on the direct neutron-capture process. We employ a shell-model approach in several model spaces with well-established effective interactions to calculate spectra and spectroscopic factors in a set of 50 neutron-rich target nuclei in different mass regions, including doubly-, semi-magic and deformed ones. Those theoretical energies and spectroscopic factors are used to evaluate direct neutron capture rates and to test global theoretical models using average spectroscopic factors and level densities based on the Hartree-Fock-Bogoliubov plus combinatorial method. The comparison of shell-model and global model results reveals several discrepancies that can be related to problems in level densities. All the results show however that the direct capture is non-negligible with respect to the by-default Hauser-Feshbach predictions and can be even 100 times more important for the most neutron-rich nuclei close to the neutron drip line.
A century ago, nuclear physics entered astrophysics, giving birth to a new field of science referred to as Nuclear Astrophysics. With time, it developed at an impressive pace into a vastly inter- and multidisciplinary discipline bringing into its wak
e not only astronomy and cosmology, but also many other sub-fields of physics, especially particle, solid-state and computational physics, as well as chemistry, geology and even biology. The present Astronuclear Physics review focusses primarily on the facets of nuclear physics that are of relevance to astronomy and astrophysics, the theoretical aspects being of special concern here.
Photon strength functions describing the average response of the nucleus to an electromagnetic probe are key input information in the theoretical modelling of nuclear reactions. Consequently they are important for a wide range of fields such as nucle
ar structure, nuclear astrophysics, medical isotope production, fission and fusion reactor technologies. They are also sources of information for widely used reaction libraries such as the IAEA Reference Input Parameter Library and evaluated data files such as EGAF. In the past two decades, the amount of reaction gamma-ray data measured to determine photon strength functions has grown rapidly. Different experimental techniques have led to discrepant results and users are faced with the dilemma which (if any) of the divergent data to adopt. We report on a coordinated effort to compile and assess the existing experimental data on photon strength functions from the giant dipole resonance region to energies below the neutron separation energy. The assessment of the discrepant data at energies around or below the neutron separation energy has been possible only in a few cases where adequate information on the model-dependent analysis and estimation of uncertainties was available. In the giant dipole resonance region, we adopt the recommendations of the new IAEA photonuclear data library. We also present global empirical and semi-microscopic models that describe the photon strength functions in the entire energy region and reproduce reasonably well most of the experimental data. The compiled experimental photon strengths and recommended model calculations are available from the PSF database hosted at the IAEA (URL:www-nds.iaea.org/PSFdatabase).
Nuclear level densities (NLDs) and $gamma$-ray strength functions ($gamma$SFs) have been extracted from particle-$gamma$ coincidences of the $^{92}$Zr($p,p gamma$)$^{92}$Zr and $^{92}$Zr($p,d gamma$)$^{91}$Zr reactions using the Oslo method. The new
$^{91,92}$Zr $gamma$SF data, combined with photonuclear cross sections, cover the whole energy range from $E_{gamma} approx 1.5$~MeV up to the giant dipole resonance at $E_{gamma} approx 17$~MeV. The wide-range $gamma$SF data display structures at $E_{gamma} approx 9.5$~MeV, compatible with a superposition of the spin-flip $M1$ resonance and a pygmy $E1$ resonance. Furthermore, the $gamma$SF shows a minimum at $E_{gamma} approx 2-3$~MeV and an increase at lower $gamma$-ray energies. The experimentally constrained NLDs and $gamma$SFs are shown to reproduce known ($n, gamma$) and Maxwellian-averaged cross sections for $^{91,92}$Zr using the {sf TALYS} reaction code, thus serving as a benchmark for this indirect method of estimating ($n, gamma$) cross sections for Zr isotopes.
Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment
of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-capture process (r-process). Supernovae and merging neutron-star (NS) or NS-black hole binaries are considered as most likely sources of the r-nuclei. But arguments in favour of one or the other or both are indirect and make use of assumptions; they are based on theoretical models with remaining simplifications and shortcomings. An unambiguous observational determination of a production event is still missing. In order to facilitate searches in this direction, e.g. by looking for radioactive tracers in stellar envelopes, the interstellar medium or terrestrial reservoirs, we provide improved theoretical estimates and corresponding uncertainty ranges for the actinide production ($^{232}$Th, $^{235,236,238}$U, $^{237}$Np, $^{244}$Pu, and $^{247}$Cm) in neutrino-driven winds of core-collapse supernovae. Since state-of-the-art supernova models do not yield r-process viable conditions ---but still lack, for example, the effects of strong magnetic fields--- we base our investigation on a simple analytical, Newtonian, adiabatic and steady-state wind model and consider the superposition of a large number of contributing components, whose nucleosynthesis-relevant parameters (mass weight, entropy, expansion time scale, and neutron excess) are constrained by the assumption that the integrated wind nucleosynthesis closely reproduces the solar system distribution of r-process elements. We also test the influence of uncertain nuclear physics.
Partial cross sections of the $^{89}$Y(p,$gamma$)$^{90}$Zr reaction have been measured to investigate the $gamma$-ray strength function in the neutron-magic nucleus $^{90}$Zr. For five proton energies between $E_p=3.65$ MeV and $E_p=4.70$ MeV, partia
l cross sections for the population of seven discrete states in $^{90}$Zr have been determined by means of in-beam $gamma$-ray spectroscopy. Since these $gamma$-ray transitions are dominantly of $E1$ character, the present measurement allows an access to the low-lying dipole strength in $^{90}$Zr. A $gamma$-ray strength function based on the experimental data could be extracted, which is used to describe the total and partial cross sections of this reaction by Hauser-Feshbach calculations successfully. Significant differences with respect to previously measured strength functions from photoabsorption data point towards deviations from the Brink-Axel hypothesis relating the photo-excitation and de-excitation strength functions.
The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to com
pare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic process. We compute the early cosmic evolution of a typical r process element, Europium. Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated evolution of Eu. We find that the bulk of Eu observations at high iron abundance can be rather well fitted by either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution favors NSM as the main astrophysical site for the r process. A comparison between our calculations and spectroscopic observations at very low metallicities allows to constrain the coalescence timescale in the NSM scenario to about 0.1 to 0.2 Gyr. These values are in agreement with the coalescence timescales of some observed binary pulsars. Finally, the cosmic evolution of Eu is used to put constraints on the NSM rate, the merger rate in the horizon of the gravitational wave detectors advanced Virgo/ad LIGO, as well as the expected rate of electromagnetic counterparts to mergers (kilonovae) in large near-infrared surveys.
