Reliable predictions of weak rapid neutron capture (r-process) abundances require a systematic treatment of nuclear physics uncertainties, especially neutron capture rates far from stability. We employ new neutron capture rates from cross sections calculated with Yet Another Hauser-Feshbach Code () using an uncertainty-quantified Koning-Delaroche potential modified for use on neutron-rich systems. Using these rates as a baseline, we perform Monte-Carlo studies with independent rate variations (uncorrelated Monte Carlo) and find correlations between specific neutron capture rates and the resulting elemental abundances for the three weak r-process scenarios: two separate simulations of neutron star merger remnant accretion disks and a simulation of a magnetorotational supernova. We discuss the underlying nuclear dynamics that give rise to these correlations and the role of astrophysical conditions in them. We demonstrate how reducing the uncertainty in these rates would improve the prospects for conducting precision r-process studies in the future. We additionally present a correlated Monte Carlo study, which incorporates the full covariance matrix that describes the relationships between individual neutron capture rates that arise from an uncertainty-quantified optical potential. We find that the magnitude of the uncertainty in the abundance pattern is similar to that produced by an uncorrelated Monte Carlo that employs only the on-diagonal components of the covariance matrix. We show how correlations restructure how the abundances co-vary, but do not necessarily decrease the overall uncertainty envelope.

Radiative strength functions (RSFs) model the bulk electromagnetic response of highly excited nuclei and are critical inputs for statistical reaction codes. In this paper, we present a definition of the RSF that is consistent with Hauser-Feshbach reaction codes and that can be efficiently computed with the shell model using the Lanczos strength-function (LSF) method. We introduce a variant of the shell-model LSF method that exploits the energy-localized Brink-Axel hypothesis, which makes it possible to compute both electric and magnetic RSFs across all energies relevant to capture reactions. We verify agreement with the conventional definition of RSFs with benchmark calculations of 24Mg and then present novel results for 56Fe. For 56Fe we find that (i) the M1 RSF shape evolves smoothly with excitation energy, consistent with the energy-localized Brink-Axel hypothesis; (ii) both M1 and E1 transitions contribute significantly to the radiative strength below the photoabsorption threshold; and (iii) within the sdpf model space, the strength below 3 MeV observed in Oslo-type experiments cannot be fully reproduced. These results pave the way for a coherent microscopic description of RSFs and further motivate the use of energy-dependent RSFs in modern reaction codes.

The surrogate reaction method is an alternative to direct measurements of compound nuclear reaction cross sections. We introduce theory tools for extracting capture cross sections from experiments that use proton inelastic scattering as a surrogate reaction mechanism. This makes it possible to constrain compound nucleus decay models which are typically the largest source of uncertainty in capture cross section calculations. This letter describes the theory developments that were used to simultaneously infer 89Y(p, γ) and 89Zr(n, γ) cross sections from 90Zr(p, p′γ) surrogate measurements.

In this review, we discuss recent applications of the ab initio symmetry-adapted no-core shell-model (SA-NCSM) theory for study and prediction of structure and reactions of stable and unstable nuclei from light to medium mass range. We explore structure properties of neutron-rich He, Li, and Mg isotopes, with a focus on nuclear collectivity, clustering, and spectroscopic factors, as well as multi-particle excitations of utmost significance in the proximity of the drip lines. In addition, we present extensions of the SA-NCSM with continuum for determining the microscopic structure of reaction fragments, which enables calculations of reaction cross sections for targets from the lightest 4,6He to 40Ca, rooted in first principles. We illustrate this for neutron and proton elastic scattering, deuteron and alpha capture reactions, and alpha knock-out reactions. Furthermore, we discuss microscopic optical potentials with uncertainty quantification, a critical ingredient in many reaction models. We also discuss the impact of alpha clustering on reactions of significance to nuclear astrophysics, as well as on beta decays and beyond-the-standard-model physics.

The pygmy dipole resonance (PDR), marked by enhanced electric dipole strength near particle emission energies, offers a unique perspective on the collective dynamics of nuclear structure. Its precise nature, particularly its degree of collectivity, remains a topic of debate. In this study, we investigate low-energy dipole excitations in spherical Mo isotopes (^82Mo to ^98Mo) using a fully consistent Hartree-Fock-Bogoliubov (HFB) and quasiparticle random phase approximation (QRPA) framework. We observe that an enhancement in dipole strength near particle emission energies is closely correlated with the development of either neutron or proton skins. To further understand the nature of this enhancement, we examine the behavior of proton and neutron transition densities. Our analysis shows that these (low-lying dipole) states exhibit distinct characteristics involving in-phase oscillations within the nucleus and neutron- or proton-dominated oscillations at the surface, while the primary contributor to this enhancement displays an intricate underlying structure. We also investigate the collectivity of these excitations by analyzing two-quasiparticle fragmentations and relative energy shifts. Our findings reveal that skin oscillation states exhibit moderate collectivity, as indicated by substantial configuration mixing, but limited coherence, whereas the GDR states exhibit strong coherence and large energy shifts characteristic of fully developed collective motion. This study paves the way for future investigations into the collective nature of low-energy dipole states in the enhancement region, particularly in deformed nuclei, where nuclear shape effects may play a crucial role in their excitation dynamics.

Phenomenological optical-model potentials (OMPs) are a key ingredient for nuclear cross section libraries that enter astrophysical nucleosynthesis simulations. While existing OMPs can reliably reproduce direct reaction cross sections on β-stable targets, the lack of scattering data on β-unstable targets limits the credibility of OMPs extrapolated to the extremely neutron-rich regime reached during explosive nucleosynthesis. Recent work with fully non-local dispersive OMPs indicates that even in regions where scattering data are unavailable, bound-state quantities, such as the proton and neutron number and binding energy, can serve as powerful constraints on the OMP. In this proceeding, we describe first steps toward a global, non-local, uncertainty-quantified, and fully dispersive OMP capable of leveraging both scattering and bound-state observables from stability to the driplines. As an example, we show how single-nucleon scattering data on traditional OMP training nuclei 40,48Ca, 90Zr, 112,124Sn, 208Pb can be combined with structural information to improve neutron capture cross sections at astrophysical energies relevant for the i-process and weak r-process.

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly ^96Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

Cross sections for compound-nuclear (CN) reactions are important for nuclear astrophysics and other applications. Direct measurements are not always possible for the reactions of interest and calculations without experimental constraints can be quite uncertain. Thus indirect approaches, such as the surrogate reaction method (SRM), are being developed to fill the gaps. The SRM, which uses a (direct) inelastic scattering or transfer reaction to obtain information on the decay of a specific compound nucleus, has a long history of providing probabilities for fission, γ, and particle emission. While earlier implementations of the method used minimal theory to provide approximate cross sections for (n,f) reactions, better theoretical descriptions of the underlying reaction mechanisms have made it possible to also obtain (n,γ), (n,n’), and (n,2n) cross sections that agree well with benchmarks. I discuss multiple applications of the modern implementation of the SRM, high-light theory advances that enable them, and comment on opportunities offered at new experimental facilities.

The past couple of decades have seen tremendous advances in nuclear structure and reaction theory. Innovative theory frameworks for describing the nuclear many-body system, increasingly powerful computers, and opportunities for confronting theory predictions with data on unstable nuclei, have been driving the field. An important goal is to move from phenomenological ingredients in reaction calculations to predictive theories based on microscopic frameworks. We discuss ongoing efforts aimed at integrating microscopic descriptions of nuclear structure into reaction predictions for medium-mass and heavy nuclei. This contribution highlights areas where Eric Bauge, a champion for building bridges, has made important contributions by encouraging and enabling collaborations between communities with complementary expertise.

The zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of-the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the quasiparticle random phase approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K pi = 0+ QRPA components and demonstrate an analogous dipole-octupole coupling through the K pi = 0-and K pi = 1-components. Intrinsic transition densities and associated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e., we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

Optical-model potentials (OMPs) are critical ingredients for basic and applied nuclear physics. Present-day computational capabilities allow us to generate data-driven nucleon-nucleus OMPs that are non-local and exactly dispersive (as theoretically required to be), include statistically-sound uncertainty quantification, and are trained on both scattering and bound-state data from a wide area of the nuclear chart. Combined together, these features allow for significant improvement in fidelity and extrapolative power of the model. Here, we present preliminary work toward the development and training of such an OMP. The capability of the model to describe data at this first stage is encouraging.

This white paper is the result of a collaboration by those that attended a workshop at the Facility for Rare Isotope Beams (FRIB), organized by the FRIB Theory Alliance (FRIB-TA), on Theoretical Justifications and Motivations for Early High-Profile FRIB Experiments. It covers a wide range of topics related to the science that will be explored at FRIB. After a brief introduction, the sections address: (II) Overview of theoretical methods, (III) Experimental capabilities, (IV) Structure, (V) Near-threshold Physics, (VI) Reaction mechanisms, (VII) Nuclear equations of state, (VIII) Nuclear astrophysics, (IX) Fundamental symmetries, and (X) Experimental design and uncertainty quantification.

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the effective interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial, and are an important current line of research. In this work, we present a new approach for the systematic derivation of optical potentials using input from nuclear structure models. We provide an application addressing the n+24Mg elastic scattering reaction based on the valence shell model.

Porter-Thomas fluctuations of neutron widths skew compound nuclear decay probabilities from their statistical Hauser-Feshbach values. We present a straightforward method to correct Hauser-Feshbach calculations for these fluctuations, useful for modeling near-threshold competition between gamma and neutron emission following beta decay or when standard width fluctuation corrections are inadequate.

We report on the first ab initio informed a knock-out reaction in the intermediate-mass region, with the aim to probe the underlying chiral potential and its impact on the emergence of alpha clustering in this mass region. The theoretical predictions of the a+16O clustering in the 20Ne ground state, based on the ab initio symmetry-adapted no-core shell model with continuum, yield a triple differential cross section for 20Ne(p, pa)16O that is in a remarkable agreement with the data. This allows us to examine predictions of surface and in-medium acluster features that emerge from the underlying realistic nucleon-nucleon interaction with no parameters fitted to nuclear data beyond the two-body system, and to compare these to the successful antisymmetrized molecular dynamics approach.

We investigated decays of ^{51,52,53}K at the ISOLDE Decay Station at CERN in order to understand the mechanism of the β-delayed neutron-emission (βn) process. The experiment quantified neutron and γ-ray emission paths for each precursor. We used this information to test the hypothesis, first formulated by Bohr in 1939, that neutrons in the βn process originate from the structureless "compound nucleus." The data are consistent with this postulate for most of the observed decay paths. The agreement, however, is surprising because the compound-nucleus stage should not be achieved in the studied β decay due to insufficient excitation energy and level densities in the neutron emitter. In the ^{53}K βn decay, we found a preferential population of the first excited state in ^{52}Ca that contradicted Bohr's hypothesis. The latter was interpreted as evidence for direct neutron emission sensitive to the structure of the neutron-unbound state. We propose that the observed nonstatistical neutron emission proceeds through the coupling with nearby doorway states that have large neutron-emission probabilities. The appearance of "compound-nucleus" decay is caused by the aggregated small contributions of multiple doorway states at higher excitation energy.

Developing a predictive capability for inelastic scattering will find applications in multiple areas. Experimental data for neutron-nucleus inelastic scattering is limited and thus one needs a robust theoretical framework to complement it. Charged-particle inelastic scattering can be used as a surrogate for $(n, \gamma)$ reactions to predict capture cross sections for unstable nuclei. Our work uses microscopic nuclear structure calculations for spherical nuclei to obtain nucleon-nucleus scattering potentials and calculate cross sections for these processes. We implement the Jeukenne, Lejeune, Mahaux (JLM) semi-microscopic folding approach, where the medium effects on nuclear interaction are parameterized in nuclear matter to obtain the nucleon-nucleon $(NN)$ interaction in a medium at positive energies. We solve for the nuclear ground state using the Hartree-Fock-Bogliubov (HFB) many-body method, assuming the nucleons within the nucleus interact via the Gogny-D1M potential. The vibrational excited states of the target nucleus are calculated using the quasi-particle random phase approximation (QRPA). We demonstrate our approach for spherical nuclei in the medium-mass region, showing scattering results for the $^{90}$Zr nucleus.

Nuclear reaction data required for astrophysics and applications is incomplete, as not all nuclear reactions can be measured or reliably predicted. Neutron-induced reactions involving unstable targets are particularly challenging, but often critical for simulations. In response to this need, indirect approaches, such as the surrogate reaction method, have been developed. Nuclear theory is key to extract reliable cross sections from such indirect measurements. We describe ongoing efforts to expand the theoretical capabilities that enable surrogate reaction measurements. We focus on microscopic predictions for charged-particle inelastic scattering, uncertainty -quantified optical nucleon-nucleus models, and neural-network enhanced parameter inference.

Reaction measurements on fission products are being planned at both Argonne National Lab and at the Facility for Rare Isotope Beams. These indirect experiments produce specific short-lived nuclei via beta decay, and the subsequent neutron and gamma emission are studied. Some initial experiments found a surprising overabundance of gamma emission, which theory has yet to explain. To remedy this, we are developing an integrated nuclear data workflow that connects advanced nuclear shell model codes for describing the beta decay with a contemporary nuclear reaction model code.