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.

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.

Differing reaction rates of long-lived nuclear states can force the level occupations out of thermal equilibrium, causing calculations of overall rates which rely on thermal equilibrium to be inaccurate. Therefore, nucleosynthesis calculations which include nuclei with isomers must use techniques that do not assume thermal equilibrium, and it is imperative that such techniques appropriately account for transitions between the ground and isomeric states via higher-lying levels. We develop a formalism to compute the steady state occupations of nuclear levels and apply it to the examples 26Al, 34Cl, and 85Kr. We show that this approach is useful both for assessing the required number of nuclear levels and for determining the temperature above which thermal equilibrium rates are appropriate.

Isotopes with low-lying long-lived isomers can behave very differently from other isotopes in astrophysical environments. In particular, the assumption of thermal equilibrium in computing the temperature-dependent $\beta$-decay rates of such isotopes can fail below certain temperatures. We focus on the $\beta$-decay of $^{26}$Al since it is one of the most important isotopes in observational astrophysics and has a low-lying isomeric state; we compare and contrast these results with $^{34}$Cl. We rule out recently reported $^{26}$Al effective $\beta$-decay rates that showed large differences from previous calculations, finding that we agree with the earlier results. We conclude that in general, effective $\beta$-decay rates should be defined separately for the ground and isomeric states at temperatures where thermal equilibrium cannot be achieved.

We present shell model calculations of nuclear neutrino energy spectra for 70 $sd$-shell nuclei over the mass number range $A=21-35$. Our calculations include nuclear excited states as appropriate for the hot and dense conditions characteristic of pre-collapse massive stars. We consider neutrinos produced by charged lepton captures and decays and, for the first time in tabular form, neutral current nuclear deexcitation, providing neutrino energy spectra on the Fuller-Fowler-Newman temperature-density grid for these interaction channels for each nucleus. We use the full $sd$-shell model space to compute initial nuclear states up to 20 MeV excitation with transitions to final states up to 35-40 MeV, employing a modification of the Brink-Axel hypothesis to handle high temperature population factors and the nuclear partition functions.

Detailed balance is often invoked in discussions of nuclear weak transitions in astrophysical environments. Satisfaction of detailed balance is rightly touted as a virtue of some methods of computing nuclear transition strengths, but we argue that it need not necessarily be strictly obeyed, especially when the system is far from weak equilibrium. We present the results of shell model calculations of nuclear weak strengths in both charged current and neutral current channels at astrophysical temperatures. Using these strengths to compute some reaction rates, we find that, despite some violation of detailed balance, our method is robust up to high temperature, and we comment on the relationship between detailed balance and weak equilibrium in astrophysical conditions.

Astrophysical environments that reach temperatures greater than $\sim$ 100 keV can have significant neutrino energy loss via both plasma processes and nuclear weak interactions. We find that nuclear processes likely produce the highest-energy neutrinos. Among the important weak nuclear interactions are both charged current channels (electron capture/emission and positron capture/emission) and neutral current channels (de-excitation of nuclei via neutrino pair emission). We show that in order to make a realistic prediction of the nuclear neutrino spectrum, one must take nuclear structure into account; in some cases, the most important transitions may involve excited states, possibly in both parent and daughter nuclei. We find that the standard technique of producing a neutrino energy spectrum by using a single transition with a Q-value and matrix element chosen to fit published neutrino production rates and energy losses will not accurately capture important spectral features.

We present shell model calculations of electron capture strength distributions in A=28 nuclei and computations of the corresponding capture rates in supernova core conditions. We find that in these nuclei the Brink-Axel hypothesis for the distribution of Gamow-Teller strength fails at low and moderate initial excitation energy, but may be a valid tool at high excitation. The redistribution of GT strength at high initial excitation may affect capture rates during collapse. If these trends which we have found in lighter nuclei also apply for the heavier nuclei which provide the principal channels for neutronization during stellar collapse, then there could be two implications for supernova core electron capture physics. First, a modified Brink-Axel hypothesis could be a valid approximation for use in collapse codes. Second, the electron capture strength may be moved down significantly in transition energy, which would likely have the effect of increasing the overall electron capture rate during stellar collapse.

We present shell-model calculations showing that residual interaction-induced configuration mixing enhances the rate of neutral current de-excitation of thermally excited nuclei into neutrino-antineutrino pairs. Though our calculations reinforce the conclusions of previous studies that this process is the dominant source of neutrino pairs near the onset of neutrino trapping during stellar collapse, our shell-model result has the effect of increasing the energy of these pairs, possibly altering their role in entropy transport in supernovae.