We measured carbon abundances and the $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratio in 31 giant branch stars with previous CN and CH band measurements that span -2.33 $<$ M$_{\rm V}$ $<$ 0.18 in the globular cluster M10 (NGC 6254). Abundances were determined by comparing CO features at $\sim 2.3\, \mu \mathrm{m}$ and specifically the $^{13}$CO bandhead at $2.37\, \mu \mathrm{m}$, to synthetic spectra generated with MOOG. The observed spectra were obtained with GNIRS on Gemini North with a resolution of R $\approx 3500$. The carbon abundances derived from the IR spectra agree with previous [C/Fe] measurements found using CN and CH features at the near-UV/blue wavelength range. We found an average carbon isotope ratio of $^{12}\mathrm{C}/^{13}\mathrm{C}$ = 5.10$_{-0.17}^{+0.18}$ for first generation stars (CN-normal; 13 stars total) and $^{12}\mathrm{C}/^{13}\mathrm{C}$ = 4.84$_{-0.22}^{+0.27}$ for second generation stars (CN-enhanced; 15 stars). We therefore find no statistically significant difference in $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratio between stars in either population for the observed magnitude range. Finally, we created models of the expected carbon, nitrogen, and $^{12}\mathrm{C}/^{13}\mathrm{C}$ surface abundance evolution on the red giant branch due to thermohaline mixing using the MESA stellar evolution code. The efficiency of the thermohaline mixing must be increased to a factor of $\approx 60$ to match [C/Fe] abundances, and by a factor of $\approx 666$ to match $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratios. We could not simultaneously fit the evolution of both carbon and the $^{12}\mathrm{C}/^{13}\mathrm{C}$ ratio with models using the same thermohaline efficiency parameter.

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.

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

The transient neutron star (NS) low-mass X-ray binary MAXI J0556$-$332 provides a rare opportunity to study NS crust heating and subsequent cooling for multiple outbursts of the same source. We examine {\it MAXI}, {\it Swift}, {\it Chandra}, and {\it XMM-Newton} data of MAXI J0556$-$332 obtained during and after three accretion outbursts of different durations and brightness. We report on new data obtained after outburst III. The source has been tracked up to $\sim$1800 d after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from $\sim$333 eV to $\sim$146 eV was observed during the first $\sim$1200 d. Outburst III reheated the crust up to $\sim$167 eV, after which the crust cooled again to $\sim$131 eV in $\sim$350 d. We model the thermal evolution of the crust and find that this source required a different strength and depth of shallow heating during each of the three outbursts. The shallow heating released during outburst I was $\sim$17 MeV nucleon$^{-1}$ and outburst III required $\sim$0.3 MeV nucleon$^{-1}$. These cooling observations could not be explained without shallow heating. The shallow heating for outburst II was not well constrained and could vary from $\sim$0--2.2 MeV nucleon$^{-1}$, i.e., this outburst could in principle be explained without invoking shallow heating. We discuss the nature of the shallow heating and why it may occur at different strengths and depths during different outbursts.

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.

In AM CVn binaries, a white dwarf primary accretes material from a helium-rich white dwarf or stellar companion. The unstable ignition of nuclear burning via the $3\alpha$ reaction in an accumulated helium layer powers a thermonuclear runaway near accretion rates $\dot{M} \lesssim 10^{-6} \, \mathrm{M_{\odot} \ yr^{-1}}$ that may be observed as helium nova or .Ia supernova. Helium burning in the primary's envelope at temperatures $T \lesssim 10^{8} \, \mathrm{K}$ may proceed via the direct three-body fusion of $\alpha$-particles. Here we show that the direct three-body rate by Nguyen et al. (2012) -- which is reduced relative to the extrapolated resonant rate at temperatures $T \gtrsim 5 \times 10^{7} \, \mathrm{K}$ -- results in novae with longer recurrence times and larger ignition masses. By contrast, we find that the enhancement in the direct three-body rate at temperatures below $T \lesssim 5 \times 10^{7} \, \mathrm{K}$ does not result in significant differences in nova outburst properties. The most massive envelopes in our models are near the density threshold for detonation of the helium layer, where an increase in the density at ignition due to the $3\alpha$ rate may be important.

The neutron star ocean is a plasma of ions and electrons that extends from the base of the neutron star's envelope to a depth where the plasma crystallizes into a solid crust. During an accretion outburst in an X-ray transient, material accumulates in the envelope of the neutron star primary. This accumulation compresses the neutron star's outer layers and induces nuclear reactions in the ocean and crust. Accretion-driven heating raises the ocean's temperature and increases the frequencies of g-modes in the ocean; when accretion halts, the ocean cools and ocean g-mode frequencies decrease. If the observed low-frequency quasi-periodic oscillations on accreting neutron stars are g-modes in the ocean, the observed quasi-periodic oscillation frequencies will increase during outburst---reaching a maximum when the ocean temperature reaches steady state --- and subsequently decrease during quiescence. For time-averaged accretion rates during outbursts between $\langle \dot{M} \rangle = 0.1 \textrm{--} 1.0\, \dot{M}_{\rm Edd}$ the predicted g-mode fundamental $n=1$ $l=2$ frequency is between $\approx 3 \textrm{--} 7 \, \mathrm{Hz}$ for slowly rotating neutron stars. Accreting neutron stars that require extra shallow heating, such as the Z-sources MAXI J0556-332, MXB 1659-29, and XTE J1701-462, have predicted g-mode fundamental frequencies between $\approx 3 \textrm{--} 16 \, \mathrm{Hz}$. Therefore, observations of low-frequency quasi-periodic oscillations between $\approx 8 \textrm{--} 16\, \mathrm{Hz}$ in these sources, or in other transients that require shallow heating, will support a g-mode origin for the observed quasi-periodic oscillations.

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.

Crustal cooling of accretion-heated neutron stars provides insight into the stellar interior of neutron stars. The neutron star X-ray transient, KS~1731$-$260, was in outburst for 12.5 years before returning to quiescence in 2001. We have monitored the cooling of this source since then through {\it Chandra} and {\it XMM-Newton} observations. Here, we present a 150 ks {\it Chandra} observation of KS~1731$-$260 taken in August 2015, about 14.5 years into quiescence, and 6 years after the previous observation. We find that the neutron star surface temperature is consistent with the previous observation, suggesting that crustal cooling has likely stopped and the crust has reached thermal equilibrium with the core. Using a theoretical crust thermal evolution code, we fit the observed cooling curves and constrain the core temperature (T$_c = 9.35\pm0.25\times10^7$ K), composition (Q$_{imp} = 4.4^{+2.2}_{-0.5}$) and level of extra shallow heating required (Q$_{sh} = 1.36\pm0.18$ MeV/nucleon). We find that the presence of a low thermal conductivity layer, as expected from nuclear pasta, is not required to fit the cooling curve well, but cannot be excluded either.

An accretion outburst onto a neutron star transient heats the neutron star's crust out of thermal equilibrium with the core. After the outburst the crust thermally relaxes toward equilibrium with the neutron star core and the surface thermal emission powers the quiescent X-ray light curve. Crust cooling models predict that thermal equilibrium of the crust will be established $\approx 1000 \, \mathrm{d}$ into quiescence. Recent observations of the cooling neutron star transient MXB 1659-29, however, suggest that the crust did not reach thermal equilibrium with the core on the predicted timescale and continued to cool after $\approx 2500 \, \mathrm{d}$ into quiescence. Because the quiescent light curve reveals successively deeper layers of the crust, the observed late time cooling of MXB 1659-29 depends on the thermal transport in the inner crust. In particular, the observed late time cooling is consistent with a low thermal conductivity layer near the depth predicted for nuclear pasta that maintains a temperature gradient between the neutron star's inner crust and core for thousands of days into quiescence. As a result, the temperature near the crust-core boundary remains above the critical temperature for neutron superfluidity and a layer of normal neutrons forms in the inner crust. We find that the late time cooling of MXB 1659-29 is consistent with heat release from a normal neutron layer near the crust-core boundary with a long thermal time.

An accretion outburst in an X-ray transient deposits material onto the neutron star primary; this accumulation of matter induces reactions in the neutron star's crust. During the accretion outburst these reactions heat the crust out of thermal equilibrium with the core. When accretion halts, the crust cools to its long-term equilibrium temperature on observable timescales. Here we examine the accreting neutron star transient MAXI J0556-332, which is the hottest transient, at the start of quiescence, observed to date. Models of the quiescent light curve require a large deposition of heat in the shallow outer crust from an unknown source. The additional heat injected is $\approx 4\textrm{-}10\,\mathrm{MeV}$ per accreted nucleon; when the observed decline in accretion rate at the end of the outburst is accounted for, the required heating increases to $\approx 6\textrm{-}16\,\mathrm{MeV}$. This shallow heating is still required to fit the lightcurve even after taking into account a second accretion episode, uncertainties in distance, and different surface gravities. The amount of shallow heating is larger than that inferred for other neutron star transients and is larger than can be supplied by nuclear reactions or compositionally driven convection; but it is consistent with stored mechanical energy in the accretion disk. The high crust temperature ($T_b \gtrsim 10^{9} \, {\rm K}$) makes its cooling behavior in quiescence largely independent of the crust composition and envelope properties, so that future observations will probe the gravity of the source. Fits to the lightcurve disfavor the presence of Urca cooling pairs in the crust.

If the observed quasi-periodic oscillations in magnetar flares are partially confined to the crust, then the oscillation frequencies are unique probes of the nuclear physics of the neutron star crust. We study crustal oscillations in magnetars including corrections for a finite Alfven velocity. Our crust model uses a new nuclear mass formula that predicts nuclear masses with an accuracy very close to that of the finite range droplet model. This mass model for equilibrium nuclei also includes shell corrections and an updated neutron-drip line. We perturb our crust model to predict axial crust modes and assign them to observed giant flare quasi-periodic oscillation frequencies from the soft gamma-ray repeater SGR 1806-20. We find magnetar crusts that match observations for various magnetic field strengths, entrainment of the free neutron gas in the inner crust, and crust-core transition densities. We find that observations can be reconciled with smaller values of the symmetry energy slope parameter, L, if there is a significant amount of entrainment of the neutrons by the superfluid or if the crust-core transition density is large. We also find neutron star masses and radii which are in agreement with expectations from what is known about low-density matter from nuclear experiment. Matching observations with a field-free model we obtain the approximate values of M =1.35 Msun and R = 11.9 km. Matching observations using a model with the surface dipole field of SGR 1806-20 (B=2.4x10^15 G) we obtain the approximate values of M = 1.25 Msun and R = 12.4 km.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

We survey the properties of all orbit families in the rotating frame of a family of realistic triaxial potentials with central supermassive black holes (SMBHs). In such galaxies, most regular box orbits (vital for maintaining triaxiality) are associated with resonances which occupy two-dimensional surfaces in configuration space. For slow figure rotation all orbit families are largely stable. At intermediate pattern speeds a significant fraction of the resonant box orbits as well as inner long-axis tubes are destabilized by the "envelope doubling" that arises from the Coriolis forces and are driven into the destabilizing center. Thus, for pattern rotation periods .2 Gyr < Tp < 5 Gyr, the two orbit families that are most important for maintaining triaxiality are highly chaotic. As pattern speed increases there is also a sharp decrease in the overall fraction of prograde short-axis tubes and a corresponding increase in the retrograde variety. At the highest pattern speeds (close to that of triaxial bars), box-like orbits undergo a sudden transition to a new family of stable retrograde loop-like orbits, which resemble orbits in three-dimensional bars, and circulate about the short axis. Our analysis implies that triaxial systems (with central cusps and SMBHs) can either have high pattern speeds like fast bars or low patten speeds like triaxial elliptical galaxies or dark matter halos found in N-body simulations. Intermediate pattern speeds produce a high level of stochasticity in both the box and inner long-axis tube orbit families implying that stable triaxial systems are unlikely to have such pattern speeds.