The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip line, where presently only limited experimental information is available, and one has to rely strongly on theoretical predictions for nuclear properties. (b) While for many years the occurrence of the r-process has been associated with supernovae, more recent studies have cast substantial doubts on this environment. Alternative scenarios include the mergers of neutron stars, neutron-star black hole mergers, but possibly also rare classes of supernovae as well as hypernovae/collapsars with polar jet ejecta and also accretion disk outflows related to the collapse of fast rotating massive stars with high magnetic fields. Stellar r-process abundance observations, have provided insights into, and constraints on the frequency of and conditions in the responsible stellar production sites. One of them, neutron star mergers, was just identified and related to the Gravitational Wave event GW170817. High resolution observations, increasingly more precise due to improved experimental atomic data, have been particularly important in defining the heavy element abundance patterns of the old halo stars, and thus determining the extent, and nature, of the earliest nucleosynthesis in our Galaxy. Combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science, this review attempts to provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged)

Beta-delayed-neutron ($\beta$n) spectroscopy was performed using the Beta-decay Paul Trap and an array of radiation detectors. The $\beta$n branching ratios and energy spectra for $^{135,136}$Sb and $^{140}$I were obtained by measuring the time of flight of recoil ions emerging from the trapped ion cloud. These nuclei are located at the edge of an isotopic region identified as having $\beta$n branching ratios that impact the r-process abundance pattern around the A~130 peak. For $^{135,136}$Sb and $^{140}$I, $\beta$n branching ratios of 14.6(11)%, 17.6(28)%, and 7.6(28)% were determined, respectively. The $\beta$n energy spectra obtained for $^{135}$Sb and $^{140}$I are compared with results from direct neutron measurements, and the $\beta$n energy spectrum for $^{136}$Sb has been measured for the first time.

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.

The rare-earth peak in the $r$-process abundance pattern depends sensitively on both the astrophysical conditions and subtle changes in nuclear structure in the region. This work takes an important step elucidating the nuclear structure and reducing the uncertainties in $r$-process calculations via precise atomic mass measurements at the JYFLTRAP double Penning trap. $^{158}$Nd, $^{160}$Pm, $^{162}$Sm, and $^{164-166}$Gd have been measured for the first time and the precisions for $^{156}$Nd, $^{158}$Pm, $^{162,163}$Eu, $^{163}$Gd, and $^{164}$Tb have been improved considerably. Nuclear structure has been probed via two-neutron separation energies $S_{2n}$ and neutron pairing energy metrics $D_n$. The data do not support the existence of a subshell closure at $N=100$. Neutron pairing has been found to be weaker than predicted by theoretical mass models. The impact on the calculated $r$-process abundances has been studied. Substantial changes resulting in a smoother abundance distribution and a better agreement with the solar $r$-process abundances are observed.

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.

We have measured the ionization efficiency of silicon nuclear recoils with kinetic energy between 1.8 and 20 keV. We bombarded a silicon-drift diode with a neutron beam to perform an elastic-scattering experiment. A broad-energy neutron spectrum was used and the nuclear recoil energy was reconstructed using a measurement of the time of flight and scattering angle of the scattered neutron. The overall trend of the results of this work is well described by the theory of Lindhard et al. above 4 keV of recoil energy. Below this energy, the presented data shows a deviation from the model. The data indicates a faster drop than the theory prediction at low energies.

The low-lying energy levels of proton-rich $^{56}$Cu have been extracted using in-beam $\gamma$-ray spectroscopy with the state-of-the-art $\gamma$-ray tracking array GRETINA in conjunction with the S800 spectrograph at the National Superconducting Cyclotron Laboratory at Michigan State University. Excited states in $^{56}$Cu serve as resonances in the $^{55}$Ni(p,$\gamma$)$^{56}$Cu reaction, which is a part of the rp-process in type I x-ray bursts. To resolve existing ambiguities in the reaction Q-value, a more localized IMME mass fit is used resulting in $Q=639\pm82$~keV. We derive the first experimentally-constrained thermonuclear reaction rate for $^{55}$Ni(p,$\gamma$)$^{56}$Cu. We find that, with this new rate, the rp-process may bypass the $^{56}$Ni waiting point via the $^{55}$Ni(p,$\gamma$) reaction for typical x-ray burst conditions with a branching of up to $\sim$40$\%$. We also identify additional nuclear physics uncertainties that need to be addressed before drawing final conclusions about the rp-process reaction flow in the $^{56}$Ni region.

The astrophysical rapid neutron capture process or `$r$ process' of nucleosynthesis is believed to be responsible for the production of approximately half the heavy element abundances found in nature. This multifaceted problem remains one of the greatest open challenges in all of physics. Knowledge of nuclear physics properties such as masses, $\beta$-decay and neutron capture rates, as well as $\beta$-delayed neutron emission probabilities are critical inputs that go into calculations of $r$-process nucleosynthesis. While properties of nuclei near stability have been established, much still remains unknown regarding neutron-rich nuclei far from stability that may participate in the $r$ process. Sensitivity studies gauge the astrophysical response of a change in nuclear physics input(s) which allows for the isolation of the most important nuclear properties that shape the final abundances observed in nature. This review summarizes the extent of recent sensitivity studies and highlights how these studies play a key role in facilitating new insight into the $r$ process. The development of these tools promotes a focused effort for state-of-the-art measurements, motivates construction of new facilities and will ultimately move the community towards addressing the grand challenge of `How were the elements from iron to uranium made?'.

We have performed for the first time a comprehensive study of the sensitivity of $r$-process nucleosynthesis to individual nuclear masses across the chart of nuclides. Using the latest version (2012) of the Finite-Range Droplet Model, we consider mass variations of $\pm0.5$ MeV and propagate each mass change to all affected quantities, including $Q$-values, reaction rates, and branching ratios. We find such mass variations can result in up to an order of magnitude local change in the final abundance pattern produced in an $r$-process simulation. We identify key nuclei whose masses have a substantial impact on abundance predictions for hot, cold, and neutron star merger $r$-process scenarios and could be measured at future radioactive beam facilities.

Nuclear masses are one of the key ingredients of nuclear physics that go into astrophysical simulations of the $r$ process. Nuclear masses effect $r$-process abundances by entering into calculations of Q-values, neutron capture rates, photo-dissociation rates, beta-decay rates, branching ratios and the properties of fission. Most of the thousands of short-lived neutron-rich nuclei which are believed to participate in the $r$ process lack any experimental verification, thus the identification of the most influential nuclei is of paramount importance. We have conducted mass sensitivity studies near the $N=82$ closed shell in the context of a main $r$-process. Our studies take into account how an uncertainty in a single nuclear mass propagates to influence the relevant quantities of neighboring nuclei and finally to $r$-process abundances. We identify influential nuclei in various astrophysical conditions using the FRDM mass model. We show that our conclusions regarding these key nuclei are still retained when a superposition of astrophysical trajectories is considered.

Rapid neutron capture or `$r$-process' nucleosynthesis may be responsible for half the production of heavy elements above iron on the periodic table. Masses are one of the most important nuclear physics ingredients that go into calculations of $r$-process nucleosynthesis as they enter into the calculations of reaction rates, decay rates, branching ratios and Q-values. We explore the impact of uncertainties in three nuclear mass models on $r$-process abundances by performing global monte carlo simulations. We show that root-mean-square (rms) errors of current mass models are large so that current $r$-process predictions are insufficient in predicting features found in solar residuals and in $r$-process enhanced metal poor stars. We conclude that the reduction of global rms errors below $100$ keV will allow for more robust $r$-process predictions.

We have performed for the first time a complete $r$-process mass sensitivity study in the $N=82$ region. We take into account how an uncertainty in a single nuclear mass propagates to influence important quantities of neighboring nuclei, including Q-values and reaction rates. We demonstrate that nuclear mass uncertainties of $\pm0.5$ MeV in the $N=82$ region result in up to an order of magnitude local change in $r$-process abundance predictions. We identify key nuclei in the study whose mass has a substantial impact on final $r$-process abundances and could be measured at future radioactive beam facilities.

A JINA/VISTARS r-process campaign was completed at the A1900 Fragment Separator of the National Superconducting Cyclotron Laboratory in the fall of 2005. The purpose of the campaign was the measurement of the beta-decay half-lives and beta-delayed neutron-emission probabilities of different unknown neutron-rich nuclei participating in the r-process. Details of this campaign will be presented.

Isomeric states in the nuclei along the rapid proton capture process path are studied by the projected shell model. Emphasis is given to two waiting point nuclei Se-68 and Kr-72 that are characterized by shape coexistence. Energy surface calculations indicate that the ground state of these nuclei corresponds to an oblate-deformed minimum, while the lowest state at the prolate-deformed minimum can be considered as a shape isomer. Due to occupation of the orbitals with large K-components, states built upon two-quasiparticle excitations at the oblate-deformed minimum may form high $K$-isomers. The impact of the isomer states on isotopic abundance in X-ray bursts is studied in a multi-mass-zone X-ray burst model by assuming an upper-lower limit approach.

Nuclei with magic numbers serve as important benchmarks in nuclear theory. In addition, neutron-rich nuclei play an important role in the astrophysical rapid neutron-capture process (r-process). 78Ni is the only doubly-magic nucleus that is also an important waiting point in the r-process, and serves as a major bottleneck in the synthesis of heavier elements. The half-life of 78Ni has been experimentally deduced for the first time at the Coupled Cyclotron Facility of the National Superconducting Cyclotron Laboratory at Michigan State University, and was found to be 110 (+100 -60) ms. In the same experiment, a first half-life was deduced for 77Ni of 128 (+27 -33) ms, and more precise half-lives were deduced for 75Ni and 76Ni of 344 (+20 -24) ms and 238 (+15 -18) ms respectively.

Isomeric states in the nuclei along the rapid proton capture process path are studied by the projected shell model. Emphasis is given to two waiting point nuclei 68Se and 72Kr that are characterized by shape coexistence. Energy surface calculations indicate that the ground state of these nuclei corresponds to an oblate-deformed minimum, while the lowest state at the prolate-deformed minimum can be considered as a shape isomer. The impact of these isomer states on isotopic abundance in x-ray bursts is studied in a multi-mass-zone x-ray burst model by assuming an upper-lower limit approach.

The structure of collective states and quasi-particle excitations in $^{68}$Se is analyzed by using the projected shell model. It is found that the $g_{9/2}$ orbitals play a decisive role for the structure of this nucleus. At the prolate minimum, alignment of the low-K states causes the observed first backbending in moment of inertia and we predict a second backbending that has not yet been reported. At the oblate minimum, quasi-particles having high-K can form long-lived isomeric states. The structure of $^{68}$Se is discussed with respect to the isomers and their possible impact on the determination of its nuclear mass.

We calculate the rapid proton (rp) capture process of hydrogen burning on the surface of an accreting neutron star with an updated reaction network that extends up to Xe, far beyond previous work. In both steady-state nuclear burning appropriate for rapidly accreting neutron stars (such as the magnetic polar caps of accreting X-ray pulsars) and unstable burning of Type I X-ray bursts, we find that the rp process ends in a closed SnSbTe cycle. This prevents the synthesis of elements heavier than Te and has important consequences for X-ray burst profiles, the composition of accreting neutron stars, and potentially galactic nucleosynthesis of light p nuclei.