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.

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.

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.

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.

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.

We examine the physics of the early universe when Majorana neutrinos (electron neutrino, muon neutrino, tau neutrino) possess transition magnetic moments. These extra couplings beyond the usual weak interaction couplings alter the way neutrinos decouple from the plasma of electrons/positrons and photons. We calculate how transition magnetic moment couplings modify neutrino decoupling temperatures, and then use a full weak, strong, and electromagnetic reaction network to compute corresponding changes in Big Bang Nucleosynthesis abundance yields. We find that light element abundances and other cosmological parameters are sensitive to magnetic couplings on the order of 10^{-10} Bohr magnetons. Given the recent analysis of sub-MeV Borexino data which constrains Majorana moments to the order of 10^{-11} Bohr magnetons or less, we find that changes in cosmological parameters from magnetic contributions to neutrino decoupling temperatures are below the level of upcoming precision observations.