$\textbf{Background}$ More than half of all the elements heavier than iron are made by the rapid neutron capture process (or r process). For very neutron-rich astrophysical conditions, such at those found in the tidal ejecta of neutron stars, nuclear fission determines the r-process endpoint, and the fission fragment yields shape the final abundances of $110\le A \le 170$ nuclei. The knowledge of fission fragment yields of hundreds of nuclei inhabiting very neutron-rich regions of the nuclear landscape is thus crucial for the modeling of heavy-element nucleosynthesis. $\textbf{Purpose}$ In this study, we propose a model for the fast calculation of fission fragment yields based on the concept of shell-stabilized prefragments defined with help of the nucleonic localization functions. $\textbf{Methods}$ To generate realistic potential energy surfaces and nucleonic localizations, we apply Skyrme Density Functional Theory. The distribution of the neck nucleons among the two prefragments is obtained by means of a statistical model. $\textbf{Results}$ We benchmark the method by studying the fission yields of $^{178}$Pt, $^{240}$Pu, $^{254}$Cf, and $^{254,256,258}$Fm and show that it satisfactorily explains the experimental data. We then make predictions for $^{254}$Pu and $^{290}$Fm as two representative cases of fissioning nuclei that are expected to significantly contribute during the r-process nucleosynthesis occurring in neutron star mergers. $\textbf{Conclusions}$ The proposed framework provides an efficient alternative to microscopic approaches based on the evolution of the system in a space of collective coordinates all the way to scission. It can be used to carry out global calculations of fission fragment distributions across the r-process region.

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.

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.

Background: Nuclear pasta phases, present in the inner crust of neutron stars, are associated with nucleonic matter at sub-saturation densities arranged in regular shapes. Those complex phases, residing in a layer which is approximately 100 m thick, impact many features of neutron stars. Theoretical quantum-mechanical simulations of nuclear pasta are usually carried out in finite 3D boxes assuming periodic boundary conditions (PBC). The resulting solutions are affected by spurious finite-size effects. Purpose: In order to remove spurious finite-size effects, it is convenient to employ twist-averaged boundary conditions (TABC) used in condensed matter, nuclear matter, and lattice QCD applications. In this work, we study the effectiveness of TABC in the context of pasta phases simulations within nuclear density functional theory. Methods: We perform Skyrme-Hartree-Fock calculations in three dimensions by implementing Bloch boundary conditions. The TABC averages are obtained by means of Gauss-Legendre integration over twist angles. Results: We benchmark the TABC for a free nucleonic gas and apply it to simple cases such as the rod and slab phases, as well as to more elaborate P-surface and gyroidal phases. Conclusions: We demonstrate that by applying TABC reliable results can be obtained from calculations performed in relatively small volumes. By studying various contributions to the total energy, we gain insights into pasta phases in mid-density range.

Nuclear masses play a fundamental role in understanding how the heaviest elements in the Universe are created in the $r$-process. We predict $r$-process nucleosynthesis yields using neutron capture and photodissociation rates that are based on nuclear density functional theory. Using six Skyrme energy density functionals based on different optimization protocols, we determine for the first time systematic uncertainty bands -- related to mass modeling -- for $r$-process abundances in realistic astrophysical scenarios. We find that features of the underlying microphysics make an imprint on abundances especially in the vicinity of neutron shell closures: abundance peaks and troughs are reflected in trends of neutron separation energy. Further advances in nuclear theory and experiments, when linked to observations, will help in the understanding of astrophysical conditions in extreme $r$-process sites.

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

Background: An asymptotic normalization coefficient (ANC) characterizes the asymptotic form of a one-nucleon overlap integral required for description of nucleon-removal reactions. Purpose: We investigate the impact of the particle continuum on proton and neutron ANCs for mirror systems from $p$- and $sd$-shell regions. Method: We use the real-energy and complex-energy continuum shell model approaches. Results: We studied the general structure of the single-particle ANCs as a function of the binding energy and orbital angular momentum. We computed ANCs in mirror nuclei for different physical situations, including capture reactions to weakly-bound and unbound states. Conclusions: We demonstrated that the single-particle ANCs exhibit generic behavior that is different for charged and neutral particles. We verified the previously proposed relation [N.K. Timofeyuk, R.C. Johnson, and A.M. Mukhamedzhanov, Phys. Rev. Lett. 91, 232501 (2003); Phys. Rev. Lett. 97, 069904(E) (2006); N. K. Timofeyuk and P. Descouvemont, Phys. Rev. C 72, 064324 (2005)] between proton and neutron mirror ANCs. We find minor modifications if the spectroscopic strength is either localized in a single state or broadly distributed. For cases when several states couple strongly to the decay channel, these modifications may reach 30%.

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.