We apply for the first time fission yields determined across the chart of nuclides from the macroscopic-microscopic theory of the Finite Range Liquid Drop Model to simulations of rapid neutron capture ($r$-process) nucleosynthesis. With the fission rates and yields derived within the same theoretical framework utilized for other relevant nuclear data, our results represent an important step toward self-consistent applications of macroscopic-microscopic models in $r$-process calculations. The yields from this model are wide for nuclei with extreme neutron excess. We show that these wide distributions of neutron-rich nuclei, and particularly the asymmetric yields for key species that fission at late times in the $r$ process, can contribute significantly to the abundances of the lighter heavy elements, specifically the light precious metals palladium and silver. Since these asymmetric yields correspondingly also deposit into the lanthanide region, we consider the possible evidence for co-production by comparing our nucleosynthesis results directly with the trends in the elemental ratios of metal-poor stars rich in $r$-process material. We show that for $r$-process enhanced stars palladium over europium and silver over europium display mostly flat trends suggestive of co-production and compare to the lanthanum over europium trend which is often used to justify robustness arguments in the lanthanide region. We find that such robustness arguments may be extendable down to palladium and heavier and demonstrate that fission deposition is a mechanism by which such a universality or robustness can be achieved.

We review recent work examining the influence of fission in rapid neutron capture ($r$-process) nucleosynthesis which can take place in astrophysical environments. We briefly discuss the impact of uncertain fission barriers and fission rates on the population of heavy actinide species. We demonstrate the influence of the fission fragment distributions for neutron-rich nuclei and discuss currently available treatments, including recent macroscopic-microscopic calculations. We conclude by comparing our nucleosynthesis results directly with stellar data for metal-poor stars rich in $r$-process elements to consider whether fission plays a role in the so-called `universality' of $r$-process abundances observed from star to star.

Meteoritic analysis demonstrates that radioactive nuclei heavier than iron were present in the early Solar System. Among them, $^{129}$I and $^{247}$Cm both have a rapid neutron-capture process ($r$ process) origin and decay on the same timescale ($\simeq$ 15.6 Myr). We show that the $^{129}$I/$^{247}$Cm abundance ratio in the early Solar System (438$\pm$184) is immune to galactic evolution uncertainties and represents the first direct observational constraint for the properties of the last $r$-process event that polluted the pre-solar nebula. We investigate the physical conditions of this event using nucleosynthesis calculations and demonstrate that moderately neutron-rich ejecta can produce the observed ratio. We conclude that a dominant contribution by exceedingly neutron-rich ejecta is highly disfavoured.

Lanthanide element signatures are key to understanding many astrophysical observables from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis which enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the $r$-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multi-dimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of $r$-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the $r$-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. We find our mass predictions given outflows which undergo an extended (n,$\gamma$)$\rightleftarrows$($\gamma$,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

Neutron star mergers (NSMs) are rapid neutron capture ($\textit{r}$-process) nucleosynthesis sites that expel matter at high velocities, from $0.1c$ to as high as $0.6c$. Nuclei ejected at these speeds are sufficiently energetic to initiate spallation nuclear reactions with interstellar medium particles. We adopt a thick-target model for the propagation of high-speed heavy nuclei in the interstellar medium, similar to the transport of cosmic rays. We find that spallation may create observable perturbations to NSM isotopic abundances, particularly around the low-mass edges of the $\textit{r}$-process peaks where neighboring nuclei have very different abundances. The extent to which spallation modifies the final NSM isotopic yields depends on: (1) the ejected abundances, which are determined by the NSM astrophysical conditions and the properties of nuclei far from stability, (2) the ejecta velocity distribution and propagation in interstellar matter, and (3) the spallation cross-sections. Observed solar and stellar $\textit{r}$-process yields could thus constrain the velocity distribution of ejected neutron star matter, assuming NSMs are the dominant $\textit{r}$-process source. We suggest avenues for future work, including measurement of relevant cross sections.

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of $r$-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. We show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining", which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

Neutron star mergers (NSMs) are rapid neutron capture (r-process) nucleosynthesis sites, which eject materials at high velocities, from 0.1c to as high as 0.6c. Thus the r-process nuclei ejected from a NSM event are sufficiently energetic to initiate spallation reactions with the interstellar medium (ISM) particles. With a thick-target model for the propagation of high-speed heavy nuclei in the ISM, we find that spallation reactions may shift the r-process abundance patterns towards solar data, particularly around the low-mass edges of the r-process peaks where neighboring nuclei have very different abundances. The spallation effects depend both on the astrophysical conditions of the r-process nuclei and nuclear physics inputs for the nucleosynthesis calculations and the propagation process. This work extends that of [Wang et al.(2019)] by focusing on the influence of nuclear physics variations on spallation effects.

Of the variations in the elemental abundance patterns of stars enhanced with $r$-process elements, the variation in the relative actinide-to-lanthanide ratio is among the most significant. We investigate the source of these actinide differences in order to determine whether these variations are due to natural differences in astrophysical sites, or due to the uncertain nuclear properties that are accessed in $r$-process sites. We find that variations between relative stellar actinide abundances is most likely astrophysical in nature, owing to how neutron-rich the ejecta from an $r$-process event may be. Furthermore, if an $r$-process site is capable of generating variations in the neutron-richness of its ejected material, then only one type of $r$-process site is needed to explain all levels of observed relative actinide enhancements.

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above $Y_e \sim 0.3$ and preventing third peak r-process material from being synthesized. We compare our model to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

We have calculated a complete set of primary fission fragment mass yields, $Y(A)$, for heavy nuclei across the chart of nuclides, including those of particular relevance to the rapid neutron capture process ($r$ process) of nucleosynthesis. We assume that the nuclear shape dynamics are strongly damped which allows for a description of the fission process via Brownian shape motion across nuclear potential-energy surfaces. The macroscopic energy of the potential was obtained with the Finite-Range Liquid-Drop Model (FRLDM), while the microscopic terms were extracted from the single-particle level spectra in the fissioning system by the Strutinsky procedure for the shell energies and the BCS treatment for the pairing energies. For each nucleus considered, the fission fragment mass yield, $Y(A)$, is obtained from 50,000 -- 500,000 random walks on the appropriate potential-energy surface. The full mass and charge yield, $Y(Z,A)$, is then calculated by invoking the Wahl systematics. With this method, we have calculated a comprehensive set of fission-fragment yields from over 3,800 nuclides bounded by $80\leq Z \leq 130$ and $A\leq330$; these yields are provided as an ASCII formatted database in the supplemental material. We compare our yields to known data and discuss general trends that emerge in low-energy fission yields across the chart of nuclides.

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution.

We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, $10^{-4}$, or $10^{-3}$ and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B, V, and I time evolution. However, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source.

The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as $Y_e\sim 0.4$ in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances.

The astrophysical production site of the heaviest elements in the universe remains a mystery. Incorporating heavy element signatures of metal-poor, $r$-process enhanced stars into theoretical studies of $r$-process production can offer crucial constraints on the origin of heavy elements. In this study, we introduce and apply the "Actinide-Dilution with Matching" model to a variety of stellar groups ranging from actinide-deficient to actinide-enhanced to empirically characterize $r$-process ejecta mass as a function of electron fraction. We find that actinide-boost stars do not indicate the need for a unique and separate $r$-process progenitor. Rather, small variations of neutron richness within the same type of $r$-process event can account for all observed levels of actinide enhancements. The very low-$Y_e$, fission-cycling ejecta of an $r$-process event need only constitute 10-30% of the total ejecta mass to accommodate most actinide abundances of metal-poor stars. We find that our empirical $Y_e$ distributions of ejecta are similar to those inferred from studies of GW170817 mass ejecta ratios, which is consistent with neutron-star mergers being a source of the heavy elements in metal-poor, $r$-process enhanced stars.

We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,\gamma)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,\gamma)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $\gamma$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

The rapid-neutron-capture ("r") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far fail to account for the observed over-enhancement of actinides, present in about 30% of r-process-enhanced stars. In this work, we investigate actinide production in the dynamical ejecta of a neutron star merger and explore if varying levels of neutron richness can reproduce the actinide boost. We also investigate the sensitivity of actinide production on nuclear physics properties: fission distribution, beta-decay, and mass model. For most cases, the actinides are over-produced in our models if the initial conditions are sufficiently neutron-rich for fission cycling. We find that actinide production can be so robust in the dynamical ejecta that an additional lanthanide-rich, actinide-poor component is necessary in order to match observations of actinide-boost stars. We present a simple actinide-dilution model that folds in estimated contributions from two nucleosynthetic sites within a merger event. Our study suggests that while the dynamical ejecta of a neutron star merger is a likely production site for the formation of actinides, a significant contribution from another site or sites (e.g., the neutron star merger accretion disk wind) is required to explain abundances of r-process-enhanced, metal-poor stars.

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.

Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.

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.

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.

We present $\beta$-delayed neutron emission and $\beta$-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 $\beta$-decay, we follow the statistical decay taking into account competition between neutrons, $\gamma$-rays, and fission. We find a region of the chart of nuclides where the probability of $\beta$-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 $\beta$-delayed fission, a new decay mode for extremely neutron-rich heavy nuclei. In this new decay mode, $\beta$-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 $\beta$-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.

Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated Gamma-ray burst and analyze their implication on the origin of r-process elements. The range of NS-NS merger rate densities of 320-4740 Gpc$^{-3}$ yr$^{-1}$ provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1-5 Earth masses of Eu, and 3-13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site.

Radiative neutron capture is an important nuclear reaction whose accurate description is needed for many applications ranging from nuclear technology to nuclear astrophysics. The description of such a process relies on the Hauser-Feshbach theory which requires the nuclear optical potential, level density and $\gamma$-strength function as model inputs. It has recently been suggested that the M1 scissors mode may explain discrepancies between theoretical calculations and evaluated data. We explore statistical model calculations with the strength of the M1 scissors mode estimated to be dependent on the nuclear deformation of the compound system. We show that the form of the M1 scissors mode improves the theoretical description of evaluated data and the match to experiment in both the fission product and actinide regions. Since the scissors mode occurs in the range of a few keV $\sim$ a few MeV, it may also impact the neutron capture cross sections of neutron-rich nuclei that participate in the rapid neutron capture process of nucleosynthesis. We comment on the possible impact to nucleosynthesis by evaluating neutron capture rates for neutron-rich nuclei with the M1 scissors mode active.

The bulk of the rare earth elements are believed to be synthesized in the rapid neutron capture process or $r$ process of nucleosynthesis. The solar $r$-process residuals show a small peak in the rare earths around $A\sim 160$, which is proposed to be formed dynamically during the end phase of the $r$ process by a pileup of material. This abundance feature is of particular importance as it is sensitive to both the nuclear physics inputs and the astrophysical conditions of the main $r$ process. We explore the formation of the rare earth peak from the perspective of an inverse problem, using Monte Carlo studies of nuclear masses to investigate the unknown nuclear properties required to best match rare earth abundance sector of the solar isotopic residuals. When nuclear masses are changed, we recalculate the relevant $\beta$-decay properties and neutron capture rates in the rare earth region. The feedback provided by this observational constraint allows for the reverse engineering of nuclear properties far from stability where no experimental information exists. We investigate a range of astrophysical conditions with this method and show how these lead to different predictions in the nuclear properties influential to the formation of the rare earth peak. We conclude that targeted experimental campaigns in this region will help to resolve the type of conditions responsible for the production of the rare earth nuclei, and will provide new insights into the longstanding problem of the astrophysical site(s) of the $r$ process.

The primary astrophysical source of the rare earth elements is the rapid neutron capture process ($r$ process). The rare earth peak that is seen in the solar $r$-process residuals has been proposed to originate as a pile-up of nuclei during the end of the $r$ process. We introduce a new method utilizing Monte Carlo studies of nuclear masses in the rare earth region, that includes self-consistently adjusting $\beta$-decay rates and neutron capture rates, to find the mass surfaces necessary for the formation of the rare earth peak. We demonstrate our method with two types of astrophysical scenarios, one corresponding conditions typical of core-collapse supernova winds and one corresponding to conditions typical of the ejection of the material from the tidal tails of neutron star mergers. In each type of astrophysical conditions, this method successfully locates a region of enhanced stability in the mass surface that is responsible for the rare earth peak. For each scenario, we find that the change in the mass surface has qualitatively different features, thus future measurements can shed light on the type of environment in which the $r$ process occurred.

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.

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.

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.

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.

We study the sensitivity of the r-process abundance pattern to neutron capture rates along the rare earth region (A~150 to A~180). We introduce the concepts of large nuclear flow and flow saturation which determine the neutron capture rates that are influential in setting the rare earth abundances. We illustrate the value of the two concepts by considering high entropy conditions favorable for rare earth peak production and identifying important neutron capture rates among the rare earth isotopes. We also show how these rates influence nuclear flow and specific sections of the abundance pattern.

The astrophysical site or sites responsible for the r-process of nucleosynthesis still remains an enigma. Since the rare earth region is formed in the latter stages of the r-process it provides a unique probe of the astrophysical conditions during which the r-process takes place. We use features of a successful rare earth region in the context of a high entropy r-process (S>100k_B) and discuss the types of astrophysical conditions that produce abundance patterns that best match meteoritic and observational data. Despite uncertainties in nuclear physics input, this method effectively constrains astrophysical conditions.

We study the formation and final structure of the rare earth peak ($A\sim160$) of the $r$-process nucleosynthesis. The rare earth peak forms at late times in the $r$-process after neutron exhaustion (neutron-to-seed ratio unity or R=1) as matter decays back to stability. Since rare earth peak formation does not occur during \nggn \ equilibrium it is sensitive to the strong interplay between late time thermodynamic evolution and nuclear physics input. Depending on the conditions the peak forms either because of the pattern of the neutron capture rates or because of the pattern of the separation energies. We analyze three mass models under different thermodynamic conditions. We find that the subtleties of each mass model, including separation energies and neutron capture rates, influence not only the final shape of the peak but also when it forms. We identify the range of nuclei which are influential in rare earth peak formation.