The $^{24}$Mg($\alpha,\gamma$)$^{28}$Si reaction influences the production of magnesium and silicon isotopes during carbon burning and is one of eight reaction rates found to significantly impact the shape of calculated X-ray burst light curves. The reaction rate is based on measured resonance strengths and known properties of levels in $^{28}$Si. The $^{24}$Mg($\alpha,\gamma$)$^{28}$Si reaction rate has been re-evaluated including recent additional indirect data. The reaction rate is substantially unchanged from previously calculated rates, especially at astrophysically important temperatures. Increases in the reaction rate could occur at lower temperatures due to as-yet unmeasured resonances but these increases have little astrophysical impact. The $^{24}$Mg($\alpha,\gamma$)$^{28}$Si reaction rate at temperatures relevant to carbon burning and Type I X-ray bursts is well constrained by the available experimental data. This removes one reaction from the list of eight previously found to be important for X-ray burst light curve model-observation comparisons.

Recent observational advances have enabled high resolution mapping of ${^{44}}$Ti in core-collapse supernova (CCSN) remnants. Comparisons between observations and models provide stringent constraints on the CCSN mechanism. However, past work has identified several uncertain nuclear reaction rates that influence ${^{44}}$Ti and ${^{56}}$Ni production in post-processing model calculations. We evolved one dimensional models of $15~M_{\odot}$, $18~M_{\odot}$, $22~M_{\odot}$ and $25~M_{\odot}$ stars from zero-age main sequence through CCSN using {\tt MESA} (Modules for Experiments in Stellar Astrophysics) and investigated the previously identified reaction rate sensitivities of ${^{44}}$Ti and ${^{56}}$Ni production. We tested the robustness of our results by making various assumptions about the CCSN explosion energy and mass-cut. We found a number of reactions that have a significant impact on the nucleosynthesis of ${^{44}}$Ti and ${^{56}}$Ni, particularly for lower progenitor masses. Notably, the reaction rates $^{13}{\rm N}(\alpha,p)^{16}{\rm O}$, $^{17}{\rm F}(\alpha,p)^{20}{\rm Ne}$, $^{52}{\rm Fe}(\alpha,p)^{55}{\rm Co}$, $^{56}{\rm Ni}(\alpha,p)^{59}{\rm Cu}$, $^{57}{\rm Ni}(n,p)^{57}{\rm Co}$, $^{56}{\rm Co}(p,n)^{56}{\rm Ni}$, $^{39}{\rm K}(p,\gamma)^{40}{\rm Ca}$, $^{47}{\rm V}(p,\gamma)^{48}{\rm Cr}$, $^{52}{\rm Mn}(p,\gamma)^{53}{\rm Fe}$, $^{57}{\rm Co}(p,\gamma)^{58}{\rm Ni}$, and $^{39}{\rm K}(p,\alpha)^{36}{\rm Ar}$ are influential for a large number of model conditions. Furthermore, we found the list of influential reactions identified by previous post-processing studies of CCSN shock-driven nucleosynthesis is likely incomplete, motivating future larger-scale sensitivity studies.

Dependable predictions of the X-ray burst ashes and light curves require a stringent constraint on the composition of the accreted material as an input parameter. Lower metallicity models are generally based on a metal deficient donor and all metals are summed up in CNO abundances or solar metal distribution is assumed. In this work, we study the alteration of accreted composition due to spallation in the atmosphere of accreting neutron stars considering a cascading destruction process. We find that the inclusion of the cascading process brings the replenishment of CNO elements and overall survival probability is higher compared to isolated destruction of CNO elements. Spallation model provides the distribution of metals as a function of mass accretion rate. Multi-zone X-ray burst models calculated with reduced metallicities have enhanced abundances for high-mass nuclei in X-ray burst ashes. The increased metallicity due to the replenishment of CNO elements changes the composition of burst ashes compared to lower metallicity conditions. This will modify the thermal and compositional structure of accreted neutron star crusts.

The $^{23}$Al($p,\gamma$)$^{24}$Si reaction is among the most important reactions driving the energy generation in Type-I X-ray bursts. However, the present reaction-rate uncertainty limits constraints on neutron star properties that can be achieved with burst model-observation comparisons. Here, we present a novel technique for constraining this important reaction by combining the GRETINA array with the neutron detector LENDA coupled to the S800 spectrograph at the National Superconducting Cyclotron Laboratory. The $^{23}$Al($d,n$) reaction was used to populate the astrophysically important states in $^{24}$Si. This enables a measurement in complete kinematics for extracting all relevant inputs necessary to calculate the reaction rate. For the first time, a predicted close-lying doublet of a 2$_2^+$ and (4$_1^+$,0$_2^+$) state in $^{24}$Si was disentangled, finally resolving conflicting results from two previous measurements. Moreover, it was possible to extract spectroscopic factors using GRETINA and LENDA simultaneously. This new technique may be used to constrain other important reaction rates for various astrophysical scenarios.

Neutron stars cooling after sustained accretion outbursts provide unique information about the neutron star crust and underlying dense matter. Comparisons between astronomical observations of these cooling transients and model calculations of neutron star crust cooling have frequently been used to constrain neutron star properties such as the mass, radius, crust composition, and presence of nuclear pasta. These calculations often use a fixed pressure at which the crust-core transition happens, though this quantity depends on the dense matter equation of state. We demonstrate that varying the crust-core transition pressure in a manner consistent with adopting various equations of state results in modest changes to the crust cooling light curve. This validates the approach adopted in most crust cooling studies to date, where the neutron star mass and radius are varied while leaving the crust-core transition pressure constant.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Type-I X-ray bursts can be used to determine properties of accreting neutron stars via comparisons between model calculations and astronomical observations, exploiting the sensitivity of models to astrophysical conditions. However, the sensitivity of models to nuclear physics uncertainties calls into question the fidelity of constraints derived in this way. Using X-ray burst model calculations performed with the code {\tt MESA}, we investigate the impact of uncertainties for nuclear reaction rates previously identified as influential and compare them to the impact of changes in astrophysical conditions, using the conditions that are thought to best reproduce the source GS 1826-24 as a baseline. We find that reaction rate uncertainties are unlikely to significantly change conclusions about the properties of accretion onto the neutron star surface for this source. However, we find that reaction rate uncertainties significantly hinder the possibility of extracting the neutron star mass-radius ratio by matching the modeled and observed light curves due to the influence of reaction rates on the modeled light curve shape. Particularly influential nuclear reaction rates are $^{15}{\rm O}(\alpha,\gamma)$ and $^{23}{\rm Al}(p,\gamma)$, though other notable impacts arise from $^{14}{\rm O}(\alpha,p)$, $^{18}{\rm Ne}(\alpha,p)$, $^{22}{\rm Mg}(\alpha,p)$, $^{24}{\rm Mg}(\alpha,\gamma)$, $^{59}{\rm Cu}(p,\gamma)$, and $^{61}{\rm Ga}(p,\gamma)$. Furthermore, we find that varying some nuclear reaction rates within their uncertainties has an impact on the neutron star crust composition and thermal structure that is comparable to relatively significant changes accretion conditions.

Now 50 years since the existence of the neutron star crust was proposed, we review the current understanding of the nuclear physics of the outer layers of accreting neutron stars. Nuclei produced during nuclear burning replace the nascent composition of the neutron star ocean and crust. Non-equilibrium nuclear reactions driven by compression alter the outer thermal structure and chemical composition, leaving observable imprints on astronomical phenomena. As observations of bursting neutron stars and cooling neutron stars have increased, the recent volume of astronomical data allows new insights into the microphysics of the neutron star interior and the possibility to test nuclear physics input in model calculations. Despite numerous advances in our understanding of neutron star interiors and observed neutron star phenomena, many challenges remain in the astrophysics theory of accreting neutron stars, the nuclear theory of neutron-rich nuclei, and the reach and precision of terrestrial nuclear physics experiments.

Type-I X-ray burst light curves encode unique information about the structure of accreting neutron stars and the nuclear reaction rates of the rp-process that powers bursts. Using the first model calculations of hydrogen/helium burning bursts for a large range of astrophysical conditions performed with the code MESA, this work shows that simultaneous model-observation comparisons for bursts from several accretion rates $\dot{M}$ are required to remove degeneracies in astrophysical conditions that otherwise reproduce bursts for a single $\dot{M}$ and that such consistent multi-epoch modeling could possibly limit the $^{15}\rm{O}(\alpha,\gamma)^{19}\rm{Ne}$ reaction rate. Comparisons to the year 1998, 2000, and 2007 bursting epochs of the neutron star GS 1826-24 show that $\dot{M}$ must be larger than previously inferred and that the shallow heating in this source must be below 0.5 MeV/u, providing a new method to constrain the shallow heating mechanism in the outer layers of accreting neutron stars. Features of the light curve rise are used to demonstrate that a lower-limit could likely be placed on the $^{15}\rm{O}(\alpha,\gamma)$ reaction rate, demonstrating the possibility of constraining nuclear reaction rates with X-ray burst light curves.

Accreting neutron stars host a number of astronomical observables which can be used to infer the properties of the underlying dense matter. These observables are sensitive to the heating and cooling processes taking place in the accreted neutron star (NS) crust. Within the past few years it has become apparent that electron-capture/beta-decay (urca) cycles can operate within the NS crust at high temperatures. Layers of nuclei undergoing urca cycling can create a thermal barrier, or Great Wall, between heating occurring deep in the crust and the regions above the urca layers. This paper briefly reviews the urca process and the implications for observables from accreting neutron stars.

The low-lying energy levels of proton-rich $^{56}$Cu have been extracted using in-beam $\gamma$-ray spectroscopy with the state-of-the-art $\gamma$-ray tracking array GRETINA in conjunction with the S800 spectrograph at the National Superconducting Cyclotron Laboratory at Michigan State University. Excited states in $^{56}$Cu serve as resonances in the $^{55}$Ni(p,$\gamma$)$^{56}$Cu reaction, which is a part of the rp-process in type I x-ray bursts. To resolve existing ambiguities in the reaction Q-value, a more localized IMME mass fit is used resulting in $Q=639\pm82$~keV. We derive the first experimentally-constrained thermonuclear reaction rate for $^{55}$Ni(p,$\gamma$)$^{56}$Cu. We find that, with this new rate, the rp-process may bypass the $^{56}$Ni waiting point via the $^{55}$Ni(p,$\gamma$) reaction for typical x-ray burst conditions with a branching of up to $\sim$40$\%$. We also identify additional nuclear physics uncertainties that need to be addressed before drawing final conclusions about the rp-process reaction flow in the $^{56}$Ni region.

Nuclear burning near the surface of an accreting neutron star produces ashes that, when compressed deeper by further accretion, alter the star's thermal and compositional structure. Bygone nucleosynthesis can be constrained by the impact of compressed ashes on the thermal relaxation of quiescent neutron star transients. In particular, Urca cooling nuclei pairs in nuclear burning ashes, which cool the neutron star crust via neutrino emission from electron-capture/beta-decay cycles, provide signatures of prior nuclear burning over the ~century timescales it takes to accrete to the electron-capture depth of the strongest cooling pairs. Using crust cooling models of the accreting neutron star transient MAXI J0556-332, we show that this source likely lacked Type I X-ray bursts and superbursts >120 years ago. Reduced nuclear physics uncertainties in rp-process reaction rates and electron-capture ft-values for low-lying transitions will improve nucleosynthesis constraints using this technique.

An accretion outburst onto a neutron star deposits hydrogen-rich and/or helium-rich material into the neutron star's envelope. Thermonuclear burning of accreted material robustly produces Urca pairs\textrm{---}pairs of nuclei that undergo cycles of $e^-$-capture and $\beta ^-$-decay. The strong $T^5$ dependence of the Urca cooling neutrino luminosity means that Urca pairs in the neutron star interior potentially remove heat from accretion-driven nuclear reactions. In this study, we identify Urca pairs in the neutron star's ocean \textrm{---} a plasma of ions and electrons overlaying the neutron star crust \textrm{---} and demonstrate that Urca cooling occurs at all depths in the ocean. We find that Urca pairs in the ocean and crust lower the ocean's steady-state temperature during an accretion outburst and that unstable carbon ignition, which is thought to trigger superbursts, occurs deeper than it would otherwise. Cooling superburst light curves, however, are only marginally impacted by cooling from Urca pairs because the superburst peak radiative luminosity $L_{\rm peak}$ is always much greater than the Urca pair neutrino luminosity $L_{\nu}$ in the hot post-superburst ocean.

X-ray bursts are thermonuclear flashes on the surface of accreting neutron stars and reliable burst models are needed to interpret observations in terms of properties of the neutron star and the binary system. We investigate the dependence of X-ray burst models on uncertainties in (p,$\gamma$), ($\alpha$,$\gamma$), and ($\alpha$,p) nuclear reaction rates using fully self-consistent burst models that account for the feedbacks between changes in nuclear energy generation and changes in astrophysical conditions. A two-step approach first identified sensitive nuclear reaction rates in a single-zone model with ignition conditions chosen to match calculations with a state-of-the-art 1D multi-zone model based on the {\Kepler} stellar evolution code. All relevant reaction rates on neutron deficient isotopes up to mass 106 were individually varied by a factor of 100 up and down. Calculations of the 84 highest impact reaction rate changes were then repeated in the 1D multi-zone model. We find a number of uncertain reaction rates that affect predictions of light curves and burst ashes significantly. The results provide insights into the nuclear processes that shape X-ray burst observables and guidance for future nuclear physics work to reduce nuclear uncertainties in X-ray burst models.