Search

Additional file 1. Supplementary analyses.

Multimessenger observations of the neutron star merger event GW170817 have re-energized the debate over the astrophysical origins of the most massive elements via the r-process nucleosynthesis. A key aspect of such studies is comparing astronomical observations to theoretical nucleosynthesis yields in a meaningful way. To perform realistic nucleosynthesis calculations, understanding the uncertainty in microphysics details such as nuclear reaction rates is as essential as understanding uncertainties in modeling the astrophysical environment. We present an investigation of neutron capture rate calculations' uncertainty away from stability using the Hauser-Feshbach model. We provide a quantitative measure of the calculations' dependability when we extrapolate models of statistical properties to nuclei in an r-process network. We select several level density and gamma-ray strength models appropriate for neutron-capture and use them to calculate the reaction rate for each nucleus in the network. We observe how statistical properties affect the theoretical reaction rates. The rates are then sampled with the Monte Carlo technique and used in network calculations to map the range of possible r-process abundances. The results show that neutron capture rates can vary by a couple of orders of magnitude between calculations. Phenomenological models provide smoother results than semi-microscopic. They cannot, however, reproduce nuclear structure changes such as shell closures. While semi-microscopic models predict nuclear structure effects away from stability, it is not clear that these results are quantitatively accurate. The effect of the uncertainty on r-process yields is large enough to impede comparisons between observation and calculations. Progress in developing better microscopic models of gamma strengths and level densities is urgently needed to improve the fidelity of r-process models., Comment: Submitted to Phys. Rev. C (abstract adapted to arXiV character limit and style)

We present $��$-delayed neutron emission and $��$-delayed fission calculations for heavy, neutron-rich nuclei using the coupled Quasi-Particle Random Phase Approximation plus Hauser-Feshbach (QRPA+HF) approach. From the initial population of a compound nucleus after $��$-decay, we follow the statistical decay taking into account competition between neutrons, $��$-rays, and fission. We find a region of the chart of nuclides where the probability of $��$-delayed fission is $\sim100$%, that likely prevents the production of superheavy elements in nature. For a subset of nuclei near the neutron dripline, neutron multiplicity and the probability of fission are both large, leading to the intriguing possibility of multi-chance $��$-delayed fission, a new decay mode for extremely neutron-rich heavy nuclei. In this new decay mode, $��$-decay can be followed by multiple neutron emission leading to subsequent daughter generations which each have a probability to fission. We explore the impact of $��$-delayed fission in rapid neutron capture process ($r$-process) nucleosynthesis in the tidal ejecta of a neutron star--neutron star merger and show that it is a key fission channel that shapes the final abundances near the second $r$-process peak., 9 pages, 5 figures, submitted

After identifying the nuclei in the regions near A =80 and A = 160 for which beta-decay rates have the greatest effect on weak and main r-process abundance patterns, we apply the finite-amplitude method (FAM) with Skyrme energy-density functionals (EDFs) to calculate beta-decay half-lives of those nuclei in the quasiparticle random-phase approximation (QRPA). We use the equal filling approximation to extend our implementation of the charge-changing FAM, which incorporates pairing correlations and allows axially symmetric deformation, to odd-A and odd-odd nuclei. Within this framework we find differences of up to a factor of seven between our calculated beta-decay half-lives and those of previous efforts. Repeated calculations with nuclei near A = 160 and multiple EDFs show a spread of two to four in beta-decay half-lives, with differences in calculated Q values playing an important role. We investigate the implications of these results for r-process simulations., Comment: 18 pages, 16 figures

Beta-decay rates for spherical neutron-rich r-process waiting-point nuclei are calculated within a fully self-consistent Quasiparticle Random-Phase Approximation, formulated in the Hartree-Fock-Bogolyubov canonical single-particle basis. The same Skyrme force is used everywhere in the calculation except in the proton-neutron particle-particle channel, where a finite-range force is consistently employed. In all but the heaviest nuclei, the resulting half-lives are usually shorter by factors of 2 to 5 than those of calculations that ignore the proton-neutron particle-particle interaction. The shorter half-lives alter predictions for the abundance distribution of r-process elements and for the time it takes to synthesize them., Comment: 14 pages RevTex, 10 eps figures, submitted to Phys. Rev. C

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event., Comment: Proceedings for the FRIB Theory Alliance workshop "FRIB and the GW170817 kilonova", held 16-27 July 2018 at the Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI USA; 44 pages, 14 figures