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 around the peak brightness 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.

We develop a method to compute thermally mediated transition rates between the ground state and long-lived isomers in nuclei. We also establish criteria delimiting a thermalization temperature above which a nucleus may be considered a single species and below which it must be treated as two separate species: a ground-state species and an astrophysical isomer (“astromer”) species. Below the thermalization temperature, the destruction rates dominate the internal transition rates between the ground state and the isomer. If the destruction rates also differ greatly from one another, the nuclear levels fall out of or fail to reach thermal equilibrium. Without thermal equilibrium, there may not be a safe assumption about the distribution of occupation probability among the nuclear levels when computing nuclear reaction rates. In these conditions, the isomer has astrophysical consequences and should be treated as a separate astromer species which evolves separately from the ground state in a nucleosynthesis network. We apply our transition-rate methods and perform sensitivity studies on a few well-known astromers. We also study transitions in several other isomers of likely astrophysical interest.

Certain nuclear isomers are well known to affect nucleosynthesis with important observable consequences (e.g., ²⁶Al and ¹⁸⁰Ta). We study the impact of nuclear isomers in the context of rapid neutron capture process (r-process) nucleosynthesis. We demonstrate that nuclear isomers are dynamically populated in the r process and that some are populated far from thermal equilibrium; this makes them astrophysical isomers, or “astromers.” We compute thermally mediated transition rates between long-lived isomers and the corresponding ground states in neutron-rich nuclei. We calculate the temperature-dependent β-decay feeding factors, which represent the fraction of material going to each of the isomer and ground state daughter species from the β-decay parent species. We simulate nucleosynthesis following the decay of a solar-like r-process composition and include as separate species nuclear excited states with measured terrestrial half-lives greater than 100 μs. We introduce a new metric to identify those astromers most likely to be influential and summarize them in a table. Notable entries include many second peak nuclei (e.g., the Te isotopic chain) and previously overlooked isomers in stable nuclei (e.g., ¹¹⁹Sn, ¹³¹Xe, and ¹⁹⁵Pt). Finally, we comment on the capacity of isomer production to alter radioactive heating in an r-process environment.

The excited-state structure of atomic nuclei can modify nuclear processes in stellar environments. In this Letter, we study the influence of nuclear excitations on Urca cooling (repeated back-and-forth $\beta$ decay and electron capture in a pair of nuclear isotopes) in the crust and ocean of neutron stars. We provide for the first time an expression for Urca process neutrino luminosity which accounts for a thermal Boltzmann distribution of excited states in both members of an Urca pair. We use our new formula with state-of-the-art nuclear structure inputs to compute neutrino luminosities of candidate Urca cooling pairs. Our nuclear inputs consist of the latest experimental data supplemented with calculations using the projected shell model. We show that, in contrast to previous results that only consider the ground states of both nuclei in the pair, our calculated neutrino luminosities for different Urca pairs vary sensitively with the environment temperature and can be radically different from those obtained in the one-transition approximation. We find that nuclear excitations can lead to an enhancement in total Urca neutrino luminosities in the accreted neutron star crust by about 5 times as compared with the previous Urca results, which is expected to cause significant observational effects.

Nuclear isomers are populated in the rapid neutron capture process (r process) of nucleosynthesis. The r process may cover a wide range of temperatures, potentially starting from several tens of GK (several MeV) and then cooling as material is ejected from the event. As the r-process environment cools, isomers can freeze out of thermal equilibrium or be directly populated as astrophysically metastable isomers (astromers). Astromers can undergo reactions and decays at rates very different from the ground state, so they may need to be treated independently in nucleosythesis simulations. Two key behaviors of astromers—ground state ↔ isomer transition rates and thermalization temperatures—are determined by direct transition rates between pairs of nuclear states. We perform a sensitivity study to constrain the effects of unknown transitions on astromer behavior. Detailed balance ensures that ground → isomer and isomer → ground transitions are symmetric, so unknown transitions are equally impactful in both directions. We also introduce a categorization of astromers that describes their potential effects in hot environments. We provide a table of neutron-rich isomers that includes the astromer type, thermalization temperature, and key unmeasured transition rates.