Transient LMXBs that host neutron stars (NSs) provide excellent laboratories for probing the dense matter physics present in NS crusts. During accretion outbursts in LMXBs, exothermic reactions may heat the NS crust, disrupting the crust-core equilibrium. When the outburst ceases, the crust cools to restore thermal equilibrium with the core. Monitoring this evolution allows us to probe the dense matter physics in the crust. Properties of the deeper crustal layers can be probed at later times after the end of the outburst. We report on the unexpected late-time temperature evolution (>2000 days after the end of their outbursts) of two NSs in LMXBs, XTE J1701-462 and EXO 0748-676. Although both these sources exhibited very different outbursts (in terms of duration and the average accretion rate), they exhibit an unusually steep decay of ~7 eV in the observed effective temperature (occurring in a time span of ~700 days) around ~2000 days after the end of their outbursts. Furthermore, they both showed an even more unexpected rise of ~3 eV in temperature (over a time period of ~500-2000 days) after this steep decay. This rise was significant at the 2.4{\sigma} and 8.5{\sigma} level for XTE J1701-462 and EXO 0748-676, respectively. The physical explanation for such behaviour is unknown and cannot be straightforwardly be explained within the cooling hypothesis. In addition, this observed evolution cannot be well explained by low-level accretion either without invoking many assumptions. We investigate the potential pathways in the theoretical heating and cooling models that could reproduce this unusual behaviour, which so far has been observed in two crust-cooling sources. Such a temperature increase has not been observed in the other NS crust-cooling sources at similarly late times, although it cannot be excluded that this might be a result of the inadequate sampling obtained at such late times.

We present the largest sample of type-I (thermonuclear) X-ray bursts yet assembled, comprising 7083 bursts from 85 bursting sources. The sample is drawn from observations with Xenon-filled proportional counters on the long-duration satellites RXTE, BeppoSAX and INTEGRAL, between 1996 February 8, and 2012 May 3. The burst sources were drawn from a comprehensive catalog of 115 burst sources, assembled from earlier catalogs and the literature. We carried out a consistent analysis for each burst lightcurve (normalised to the relative instrumental effective area), and provide measurements of rise time, peak intensity, burst timescale, and fluence. For bursts observed with the RXTE/PCA and BeppoSAX/WFC we also provide time-resolved spectroscopy, including estimates of bolometric peak flux and fluence, and spectral parameters at the peak of the burst. For 950 bursts observed with the PCA from sources with previously detected burst oscillations, we include an analysis of the high-time resolution data, providing information on the detectability and amplitude of the oscillations, as well as where in the burst they are found. We also present analysis of 117883 observations of the burst sources within the sample timeframe. We extracted 3-25 keV X-ray spectra from most observations, and (for observations meeting our signal-to-noise criterion), we provide measurements of the flux, spectral colours, and for selected sources, the position on the colour-colour diagram, for the best-fit spectral model. We present a description of the sample, a summary of the science investigations completed to date, and suggestions for further studies.

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.

Transient neutron star (NS) LMXBs undergo episodes of accretion, alternated with quiescent periods. During an accretion outburst, the NS heats up due to exothermic accretion-induced processes taking place in the crust. Besides the long-known deep crustal heating of nuclear origin, a likely non-nuclear source of heat, dubbed 'shallow heating', is present at lower densities. Most of the accretion-induced heat slowly diffuses into the core on a timescale of years. Over many outburst cycles, a state of equilibrium is reached when the core temperature is high enough that the heating and cooling (photon and neutrino emission) processes are in balance. We investigate how stellar characteristics and outburst properties affect the long-term temperature evolution of a transiently accreting NS. For the first time the effects of crustal properties are considered, particularly that of shallow heating. Using our code NSCool, we tracked the thermal evolution of a NS undergoing outbursts over a period of $10^5$ yr. The outburst sequence is based on the regular outbursts observed from Aql X-1. For each model, we calculated the timescale over which equilibrium was reached and we present these timescales along with the temperature and luminosity parameters of the equilibrium state. We find that shallow heating significantly contributes to the equilibrium state. Increasing its strength raises the equilibrium core temperature. We find that if deep crustal heating is replaced by shallow heating alone, the core would still heat up, reaching only a 2% lower equilibrium core temperature. Deep crustal heating may therefore not be vital to the heating of the core. Additionally, shallow heating can increase the quiescent luminosity to values higher than previously expected.

The structure and composition of the crust of neutron stars plays an important role in their thermal and magnetic evolution, hence in setting their observational properties. One way to study the crust properties is to measure how it cools after it has been heated during an accretion outburst in a low-mass X-ray binary (LMXB). Such studies have shown that there is a tantalizing source of heat, of currently unknown origin, that is located in the outer layers of the crust and has a strength that varies between different sources and different outbursts. With the aim of understanding the mechanism behind this "shallow heating", we present Chandra and Swift observations of the neutron star LMXB Aql X-1, obtained after its bright 2016 outburst. We find that the neutron star temperature was initially much lower, and started to decrease at much later time, than observed after the 2013 outburst of the source, despite the fact that the properties of the two outbursts were very similar. Comparing our data to thermal evolution simulations, we infer that the depth and magnitude of shallow heating must have been much larger during the 2016 outburst than during the 2013 one. This implies that basic neutron star parameters that do not change between outbursts, do not play a strong role in shallow heating. Furthermore, it suggests that outbursts with a similar accretion morphology can give rise to very different shallow heating. We also discuss alternative explanations for the difference in quiescent evolution after the 2016 outburst.

In 2017, the Be/X-ray transient 4U 0115+63 exhibited a new type-II outburst that was two times fainter than its 2015 giant outburst (in the Swift/BAT count rates). Despite this difference between the two bright events, the source displayed similar X-ray behaviour after these periods. Once the outbursts ceased, the source did not transit towards quiescence directly, but was detected about a factor of 10 above its known quiescent level. It eventually decayed back to quiescence over time scales of months. In this paper we present the results of our Swift monitoring campaign, and an XMM-Newton observation of 4U 0115+63 during the decay of the 2017 type-II outburst, and its subsequent low-luminosity behaviour. We discuss the possible origin of the decaying source emission at this low-level luminosity, which has now been shown as a recurrent phenomenon, in the framework of the two proposed scenarios to explain this faint state: cooling from an accretion-heated neutron-star crust or continuous low-level accretion. In addition, we compare the outcome of our study with the results we obtained from the 2015/2016 monitoring campaign on this source.

The Be/X-ray transient GRO J1750-27 exhibited a type-II (giant) outburst in 2015. After the source transited to quiescence, we triggered our multi-year Chandra monitoring programme to study its quiescent behaviour. The programme was designed to follow the cooling of a potentially heated neutron-star crust due to accretion of matter during the preceding outburst, similar to what we potentially have observed before in two other Be/X-ray transients, namely 4U 0115+63 and V 0332+53. However, unlike for these other two systems, we do not find any strong evidence that the neutron-star crust in GRO J1750-27 was indeed heated during the accretion phase. We detected the source at a rather low X-ray luminosity (~10^33 erg/s) during only three of our five observations. When the source was not detected it had very low-luminosity upper limits (<10^32 erg/s; depending on assumed spectral model). We interpret these detections and the variability observed as emission likely due to very low-level accretion onto the neutron star. We also discuss why the neutron-star crust in GRO J1750-27 might not have been heated while the ones in 4U 0115+63 and V 0332+53 possibly were.

We present a new Chandra observation (performed in July 2016) of the neutron star X-ray transient IGR J17480-2446, located in the globular cluster Terzan 5. We study the continued cooling of the neutron star crust in this system that was heated during the 2010 outburst of the source. This new observation was performed two years after the last observation of IGR J17480-2446, hence, significantly extending the cooling baseline. We reanalysed all available Chandra observations of the source (but excluding observations during which one of the known transients in Terzan 5 was in outburst) and fitted the obtained cooling curve with our cooling code NSCool, which allows for much improved modelling than what was previously performed for the source. The data and our fit models indicate that the crust was still cooling ~5.5 years after the outburst ended. The neutron star crust has likely not reached crust-core thermal equilibrium yet, and further cooling is predicted (which can be confirmed with additional Chandra observations in >5 years). Intriguingly, we find indications that the thermal conductivity might be relatively low in part of the crust compared to what has been inferred for other crust-cooling sources and tentatively suggest that this layer might be located around the neutron drip. The reason for this difference is unclear, but might be related to the fact that IGR J17480-2446 harbours a relatively slowly rotating neutron star (with a spin of 11 Hz) that has a relatively strong inferred surface magnetic field ($10^{9-10}$ Gauss) compared to what is known or typically assumed for other cooling sources.

Burst oscillations are brightness asymmetries that develop in the burning ocean during thermonuclear bursts on accreting neutron stars. They have been observed during H/He-triggered (Type I) bursts and Carbon-triggered superbursts. The mechanism responsible is not unknown, but the dominant burst oscillation frequency is typically within a few Hz of the spin frequency, where this is independently known. One of the best-studied burst oscillation sources, 4U 1636-536, has oscillations at $581\,\text{Hz}$ in both its regular Type I bursts and in one superburst. Recently however, Strohmayer \& Mahmoodifar reported the discovery of an additional signal at a higher frequency, $835\,\text{Hz}$, during the superburst. This higher frequency is consistent with the predictions for several types of global ocean mode, one of the possible burst oscillation mechanisms. If this is the case then the same physical mechanism may operate in the normal Type I bursts of this source. In this paper we report a stacked search for periodic signals in the regular Type I bursts: we found no significant signal at the higher frequency, with upper limits for the single trial root mean square (rms) fractional amplitude of 0.57(6)\%. Our analysis did however reveal that the dominant $581\,\text{Hz}$ burst oscillation signal is present at a weak level even in the sample of bursts where it cannot be detected in individual bursts. This indicates that any cutoff in the burst oscillation mechanism occurs below the detection threshold of existing X-ray telescopes.

Monitoring the cooling of neutron-star crusts heated during accretion outbursts allows us to infer the physics of the dense matter present in the crust. We examine the crust cooling evolution of the low-mass X-ray binary MXB 1659-29 up to ~505 days after the end of its 2015 outburst (hereafter outburst II) and compare it with what we observed after its previous 1999 outburst (hereafter outburst I) using data obtained from the Swift, XMM-Newton, and Chandra observatories. The observed effective surface temperature of the neutron star in MXB 1659-29 dropped from ~92 eV to ~56 eV from ~12 days to ~505 days after the end of outburst II. The most recently performed observation after outburst II suggests that the crust is close to returning to thermal equilibrium with the core. We model the crust heating and cooling for both its outbursts collectively to understand the effect of parameters that may change for every outburst (e.g., the average accretion rate, the length of outburst, the envelope composition of the neutron star at the end of the outburst) and those which can be assumed to remain the same during these two outbursts (e.g., the neutron star mass, its radius). Our modelling indicates that all parameters were consistent between the two outbursts with no need for any significant changes. In particular, the strength and the depth of the shallow heating mechanism at work (in the crust) were inferred to be the same during both outbursts, contrary to what has been found when modelling the cooling curves after multiple outburst of another source, MAXI J0556-332. This difference in source behaviour is not understood. We discuss our results in the context of our current understanding of cooling of accretion-heated neutron-star crusts, and in particular with respect to the unexplained shallow heating mechanism.

Runaway thermonuclear burning of a layer of accumulated fuel on the surface of a compact star provides a brief but intense display of stellar nuclear processes. For neutron stars accreting from a binary companion, these events manifest as thermonuclear (type-I) X-ray bursts, and recur on typical timescales of hours to days. We measured the burst rate as a function of accretion rate, from seven neutron stars with known spin rates, using a burst sample accumulated over several decades. At the highest accretion rates, the burst rate is lower for faster spinning stars. The observations imply that fast (> 400 Hz) rotation encourages stabilization of nuclear burning, suggesting a dynamical dependence of nuclear ignition on the spin rate. This dependence is unexpected, because faster rotation entails less shear between the surrounding accretion disk and the star. Large-scale circulation in the fuel layer, leading to enhanced mixing of the burst ashes into the fuel layer, may explain this behavior; further numerical simulations are required to confirm.

With our neutron star crust cooling code {\tt NSCool} we track the thermal evolution of the neutron star in Aql X-1 over the full accretion outburst history from 1996 until 2015. For the first time, we model many outbursts (23 outbursts were detected) collectively and in great detail. This allows us to investigate the influence of previous outbursts on the internal temperature evolution and to test different neutron star crust cooling scenarios. Aql X-1 is an ideal test source for this purpose, because it shows frequent, short outbursts and thermally dominated quiescence spectra. The source goes into outburst roughly once a year for a few months. Assuming that the quiescent {\it Swift}/XRT observations of Aql X-1 can be explained within the crust cooling scenario (Waterhouse et al. 2016), we find three main conclusions. Firstly, the data are well reproduced by our model if the envelope composition and shallow heating parameters are allowed to change between outbursts. This is not the case if both shallow heating parameters (strength and depth) are tied throughout all accretion episodes, supporting earlier results that the properties of the shallow heating mechanism are not constant between outbursts. Second, from our models shallow heating could not be connected to one specific spectral state during outburst. Third, and most importantly, we find that the neutron star in Aql X-1 does not have enough time between outbursts to cool down to crust-core equilibrium and that heating during one outburst influences the cooling curves of the next.

We studied the transient neutron-star low-mass X-ray binary GRS 1747-312, located in the globular cluster Terzan 6, in its quiescent state after its outburst in August 2004, using an archival XMM-Newton observation. A source was detected in this cluster and its X-ray spectrum can be fitted with the combination of a soft, neutron-star atmosphere model and a hard, power-law model. Both contributed roughly equally to the observed 0.5-10 keV luminosity (~4.8X10^33 erg/s). This type of X-ray spectrum is typically observed for quiescent neutron-star X-ray transients that are perhaps accreting in quiescence at very low rates. Therefore, if this X-ray source is the quiescent counterpart of GRS 1747-312, then this source is also accreting at low levels in-between outbursts. Since source confusion a likely problem in globular clusters, it is quite possible that part, if not all, of the emission we observed is not related to GRS 1747-312, and is instead associated with another source or conglomeration of sources in the cluster. Currently, it is not possible to determine exactly which part of the emission truly originates from GRS 1747-312, and a Chandra observation (when no source is in outburst in Terzan 6) is needed to be conclusive. Assuming that the detected emission is due to GRS 1747-312, we discuss the observed results in the context of what is known about other quiescent systems. We also investigated the thermal evolution of the neutron star in GRS 1747-312, and inferred that GRS 1747-312 can be considered a typical quiescent system under our assumptions.

We report on two new quiescent {\it XMM-Newton} observations (in addition to the earlier {\it Swift}/XRT and {\it XMM-Newton} coverage) of the cooling neutron star crust in the low-mass X-ray binary 1RXS J180408.9$-$342058. Its crust was heated during the $\sim$4.5 month accretion outburst of the source. From our quiescent observations, fitting the spectra with a neutron star atmosphere model, we found that the crust had cooled from $\sim$ 100 eV to $\sim$73 eV from $\sim$8 days to $\sim$479 days after the end of its outburst. However, during the most recent observation, taken $\sim$860 days after the end of the outburst, we found that the crust appeared not to have cooled further. This suggested that the crust had returned to thermal equilibrium with the neutron star core. We model the quiescent thermal evolution with the theoretical crustal cooling code NSCool and find that the source requires a shallow heat source, in addition to the standard deep crustal heating processes, contributing $\sim$0.9 MeV per accreted nucleon during outburst to explain its observed temperature decay. Our high quality {\it XMM-Newton} data required an additional hard component to adequately fit the spectra. This slightly complicates our interpretation of the quiescent data of 1RXS J180408.9$-$342058. The origin of this component is not fully understood.

In accreting neutron star X-ray transients, the neutron star crust can be substantially heated out of thermal equilibrium with the core during an accretion outburst. The observed subsequent cooling in quiescence (when accretion has halted) offers a unique opportunity to study the structure and thermal properties of the crust. Initially crust cooling modelling studies focussed on transient X-ray binaries with prolonged accretion outbursts (> 1 year) such that the crust would be significantly heated for the cooling to be detectable. Here we present the results of applying a theoretical model to the observed cooling curve after a short accretion outburst of only ~10 weeks. In our study we use the 2010 outburst of the transiently accreting 11 Hz X-ray pulsar in the globular cluster Terzan 5. Observationally it was found that the crust in this source was still hot more than 4 years after the end of its short accretion outburst. From our modelling we found that such a long-lived hot crust implies some unusual crustal properties such as a very low thermal conductivity (> 10 times lower than determined for the other crust cooling sources). In addition, we present our preliminary results of the modelling of the ongoing cooling of the neutron star in MXB 1659-298. This transient X-ray source went back into quiescence in March 2017 after an accretion phase of ~1.8 years. We compare our predictions for the cooling curve after this outburst with the cooling curve of the same source obtained after its previous outburst which ended in 2001.

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.

We have monitored the transient neutron star low-mass X-ray binary 1RXS J180408.9-342058 in quiescence after its ~4.5 month outburst in 2015. The source has been observed using Swift and XMM-Newton. Its X-ray spectra were dominated by a thermal component. The thermal evolution showed a gradual X-ray luminosity decay from ~18 x 10^32 to ~4 x 10^32 (D/5.8 kpc)^2 erg s^{-1} between ~8 and ~379 d in quiescence, and the inferred neutron star surface temperature (for an observer at infinity; using a neutron star atmosphere model) decreased from ~100 to ~71 eV. This can be interpreted as cooling of an accretion-heated neutron star crust. Modelling the observed temperature curve (using NSCOOL) indicated that the source required ~1.9 MeV per accreted nucleon of shallow heating in addition to the standard deep crustal heating to explain its thermal evolution. Alternatively, the decay could also be modelled without the presence of deep crustal heating, only having a shallow heat source (again ~1.9 MeV per accreted nucleon was required). However, the XMM-Newton data statistically required an additional power-law component. This component contributed ~30 per cent of the total unabsorbed flux in 0.5-10 keV energy range. The physical origin of this component is unknown. One possibility is that it arises from low-level accretion. The presence of this component in the spectrum complicates our cooling crust interpretation because it might indicate that the smooth luminosity and temperature decay curves we observed may not be due to crust cooling but due to some other process.

Neutron stars in low mass X-ray binaries exhibit oscillations during thermonuclear bursts, attributed to asymmetric brightness patterns on the burning surfaces. All models that have been proposed to explain the origin of these asymmetries (spreading hotspots, surface waves, and cooling wakes) depend on the accretion rate. By analysis of archival RXTE data of six oscillation sources, we investigate the accretion rate dependence of the amplitude of burst oscillations. This more than doubles the size of the sample analysed previously by Muno et al. (2004), who found indications for a relationship between accretion rate and oscillation amplitudes. We find that burst oscillation signals can be detected at all observed accretion rates. Moreover, oscillations at low accretion rates are found to have relatively small amplitudes ($A_\text{rms}\leq0.10$) while oscillations detected in bursts observed at high accretion rates cover a broad spread in amplitudes ($0.05\leq A_\text{rms}\leq0.20$). In this paper we present the results of our analysis and discuss these in the light of current burst oscillation models. Additionally, we investigate the bursts of two sources without previously detected oscillations. Despite that these sources have been observed at accretion rates where burst oscillations might be expected, we find their behaviour to be not anomalous compared to oscillation sources.

The neutron star low-mass X-ray binary and intermittent millisecond X-ray pulsar HETE J1900.1-2455 returned to quiescence in late 2015, after a prolonged accretion outburst of ~10 yr. Using a Chandra observation taken ~180 d into quiescence we detect the source at a luminosity of ~4.5E31 (D/4.7 kpc)^2 erg/s (0.5-10 keV). The X-ray spectrum can be described by a neutron star atmosphere model with a temperature of ~54 eV for an observer at infinity. We perform thermal evolution calculations based on the 2016 quiescent data and a <98 eV temperature upper limit inferred from a Swift observation taken during an unusually brief (<2 weeks) quiescent episode in 2007. We find no evidence in the present data that the thermal properties of the crust, such as the heating rate and thermal conductivity, are different than those of non-pulsating neutron stars. Finding this neutron star so cold after its long outburst imposes interesting constraints on the heat capacity of the stellar core; these become even stronger if further cooling were to occur.

Using a theoretical model, we track the thermal evolution of a cooling neutron star crust after an accretion induced heating period with the goal of constraining the crustal parameters. We present for the first time a crust cooling model $-\text{ } NSCool\text{ } -$ that takes into account detailed variability during the full outburst based on the observed light curve. We apply our model to KS 1731-260. The source was in outburst for $\sim$12 years during which it was observed to undergo variations on both long (years) and short (days-weeks) timescales. Our results show that KS 1731-260 does not reach a steady state profile during the outburst due to fluctuations in the derived accretion rate. Additionally, long time-scale outburst variability mildly affects the complete crust cooling phase, while variations in the final months of the outburst strongly influence the first $\sim$40 days of the calculated cooling curve. We discuss the consequences for estimates of the neutron star crust parameters, and argue that detailed modelling of the final phase of the outburst is key to constraining the origin of the shallow heat source.