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 International Gamma-Ray Astrophysics Laboratory 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 light curve (normalized 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/Wide Field Camera 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 118,848 observations of the burst sources within the sample time frame. 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 colors, and, for selected sources, the position on the color–color 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.

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 -NSCool -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 ∼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 ∼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.

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 MAXI, Swift, Chandra, and XMM-Newton data of MAXI J0556−332 obtained during and after three accretion outbursts of different durations and brightnesses. We report on new data obtained after outburst III. The source has been tracked up to ∼1800 days after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from ∼333 eV to ∼146 eV was observed during the first ∼1200 days. Outburst III reheated the crust up to ∼167 eV, after which the crust cooled again to ∼131 eV in ∼350 days. 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 ∼17 MeV nucleon −1 and outburst III required ∼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 ∼0 to 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.

With our neutron star crust cooling code 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 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.

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.

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.5 × 10 31 (D/4.7 kpc) 2 erg s −1 (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.

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 rms ≤ 0.10) while oscillations detected in bursts observed at high accretion rates cover a broad spread in amplitudes (0.05 ≤ A rms ≤ 0.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.

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×10 32 to ∼4×10 32 (D/5.8 kpc) 2 erg s −1 between ∼8 to ∼379 days 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. Modeling 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.