Thermonuclear Bursts Research Articles

The very faint hard state of the persistent neutron star X-ray binary SLX 1737–282 near the Galactic Centre

We report on a detailed study of the spectral and temporal properties of the neutron star low mass X-ray binary SLX 1737-282, which is located only ~1degr away from Sgr A. The system is expected to have a short orbital period, even within the ultra-compact regime, given its persistent nature at low X-ray luminosities and the long duration thermonuclear burst that it has displayed. We have analysed a Suzaku (18 ks) observation and an XMM-Newton (39 ks) observation taken 7 years apart. We infer (0.5-10 keV) X-ray luminosities in the range 3-6 x10^35erg s-1, in agreement with previous findings. The spectra are well described by a relatively cool (kTbb = 0.5 keV) black body component plus a Comptonized emission component with {\Gamma} ~1.5-1.7. These values are consistent with the source being in a faint hard state, as confirmed by the ~ 20 per cent fractional root mean square amplitude of the fast variability (0.1 - 7 Hz) inferred from the XMM-Newton data. The electron temperature of the corona is >7 keV for the Suzaku observation, but it is measured to be as low as ~2 keV in the XMM-Newton data at higher flux. The latter is significantly lower than expected for systems in the hard state. We searched for X-ray pulsations and imposed an upper limit to their semi-amplitude of 2 per cent (0.001 - 7 Hz). Finally, we investigated the origin of the low frequency variability emission present in the XMM-Newton data and ruled out an absorption dip origin. This constraint the orbital inclination of the system to 65 degr unless the orbital period is longer than 11 hr (i.e. the length of the XMM-Newton observation).

Thermonuclear burst oscillations and the dense matter equation of state

AbstractMatter in neutron star cores reaches extremely high densities, forming states of matter that cannot be generated in the laboratory. The Equation of State (EOS) of the matter links to macroscopic observables, such as mass M and radius R, via the stellar structure equations. A promising technique for measuring M and R exploits hotspots (burst oscillations) that form on the stellar surface when material accreted from a companion star undergoes a thermonuclear explosion. As the star rotates, the hotspot gives rise to a pulsation, and relativistic effects encode information about M and R into the pulse profile. However the burst oscillation mechanism remains unknown, introducing uncertainty when inferring the EOS. I review the progress that we are making towards cracking this long-standing problem, and establishing burst oscillations as a robust tool for measuring M and R. This is a major goal for future large area X-ray telescopes.

Thermonuclear Bursts with Short Recurrence Times from Neutron Stars Explained by Opacity-driven Convection

Abstract Thermonuclear flashes of hydrogen and helium accreted onto neutron stars produce the frequently observed Type I X-ray bursts. It is the current paradigm that almost all material burns in a burst, after which it takes hours to accumulate fresh fuel for the next burst. In rare cases, however, bursts are observed with recurrence times as short as minutes. We present the first one-dimensional multi-zone simulations that reproduce this phenomenon. Bursts that ignite in a relatively hot neutron star envelope leave a substantial fraction of the fuel unburned at shallow depths. In the wake of the burst, convective mixing events driven by opacity bring this fuel down to the ignition depth on the observed timescale of minutes. There, unburned hydrogen mixes with the metal-rich ashes, igniting to produce a subsequent burst. We find burst pairs and triplets, similar to the observed instances. Our simulations reproduce the observed fraction of bursts with short waiting times of ∼30%, and demonstrate that short recurrence time bursts are typically less bright and of shorter duration.

The 2015 outburst of the accretion-powered pulsar IGR J00291+5934: INTEGRAL andSwiftobservations

IGR J00291+5934 is the fastest-known accretion-powered X-ray pulsar, discovered during a transient outburst in 2004. In this paper, we report on Integral and Swift observations during the 2015 outburst, which lasts for $\sim25$ d. The source has not been observed in outburst since 2008, suggesting that the long-term accretion rate has decreased by a factor of two since discovery. The averaged broad-band (0.1 - 250 keV) persistent spectrum in 2015 is well described by a thermal Comptonization model with a column density of $N_{\rm H} \approx4\times10^{21}$ cm$^{-2}$, a plasma temperature of $kT_{\rm e} \approx50$ keV, and a Thomson optical depth of $\tau_{\rm T}\approx1$. Pulsations at the known spin period of the source are detected in the Integral data up to the $\sim150$ keV energy band. We also report on the discovery of the first thermonuclear burst observed from IGR J00291+5934, which lasts around 7 min and occurs at a persistent emission level corresponding to roughly $1.6\%$ of the Eddington accretion rate. The properties of the burst suggest it is powered primarily by helium ignited at a depth of $y_{\rm ign}\approx1.5\times10^9$ g cm$^{-2}$ following the exhaustion by steady burning of the accreted hydrogen. The Swift/BAT data from the first $\sim20$ s of the burst provide indications of a photospheric radius expansion phase. Assuming this is the case, we infer a source distance of $d = 4.2 \pm 0.3$ kpc.

THE ACCRETION RATE DEPENDENCE OF BURST OSCILLATION AMPLITUDE

ABSTRACT 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 analyzed previously by Muno et al., 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 ( ) while oscillations detected in bursts observed at high accretion rates cover a broad spread in amplitudes ( ). 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 the fact that these sources have been observed at accretion rates where burst oscillations might be expected, we find their behavior not to be anomalous compared to oscillation sources.

THE FERMI–GBM THREE-YEAR X-RAY BURST CATALOG

ABSTRACT The Fermi Gamma-ray Burst Monitor (GBM) is an all-sky gamma-ray monitor well known in the gamma-ray burst (GRB) community. Although GBM excels in detecting the hard, bright extragalactic GRBs, its sensitivity above 8 keV and its all-sky view make it an excellent instrument for the detection of rare, short-lived Galactic transients. In 2010 March, we initiated a systematic search for transients using GBM data. We conclude this phase of the search by presenting a three-year catalog of 1084 X-ray bursts. Using spectral analysis, location, and spatial distributions we classified the 1084 events into 752 thermonuclear X-ray bursts, 267 transient events from accretion flares and X-ray pulses, and 65 untriggered gamma-ray bursts. All thermonuclear bursts have peak blackbody temperatures broadly consistent with photospheric radius expansion (PRE) bursts. We find an average rate of 1.4 PRE bursts per day, integrated over all Galactic bursters within about 10 kpc. These include 33 and 10 bursts from the ultra-compact X-ray binaries 4U 0614+09 and 2S 0918-549, respectively. We discuss these recurrence times and estimate the total mass ejected by PRE bursts in our Galaxy.

THE DENSE MATTER EQUATION OF STATE FROM NEUTRON STAR RADIUS AND MASS MEASUREMENTS

ABSTRACT We present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of ∼ 2 M ⊙ ?> neutron stars, we place strong and quantitative constraints on the properties of the equation of state between ≈ 2 – 8 ?> times the nuclear saturation density. We find that around M = 1.5 M ⊙ , ?> the preferred equation of state predicts radii between 10.1 and 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential.

Observational constraints on neutron star masses and radii

Precise and reliable measurements of the masses and radii of neutron stars with a variety of masses would provide valuable guidance for improving models of the properties of cold matter with densities above the saturation density of nuclear matter. Several different approaches for measuring the masses and radii of neutron stars have been tried or proposed, including analyzing the X-ray fluxes and spectra of the emission from neutron stars in quiescent low-mass X-ray binary systems and thermonuclear burst sources; fitting the energy-dependent X-ray waveforms of rotation-powered millisecond pulsars, burst oscillations with millisecond periods, and accretion-powered millisecond pulsars; and modeling the gravitational radiation waveforms of coalescing double neutron star and neutron star -- black hole binary systems. We describe the strengths and weaknesses of these approaches, most of which currently have substantial systematic errors, and discuss the prospects for decreasing the systematic errors in each method.

Temperature measurement during thermonuclear X-ray bursts with BeppoSAX

We have carried out a study of temperature evolution during thermonuclear bursts in LMXBs using broad band data from two instruments onboard BeppoSAX, the MECSand the PDS. However, instead of applying the standard technique of time resolved spectroscopy, we have determined the temperature in small time intervals using the ratio of count rates in the two instruments assuming a blackbody nature of burst emission and different interstellar absorption for different sources. Data from a total of twelve observations of six sources were analyzed during which 22 bursts were detected. We have obtained temperatures as high as ∼3.0 keV, even when there is no evidence of photospheric radius expansion. These high temperatures were observed in the sources within different broadband spectral states (soft and hard).

EVIDENCE FOR ENHANCED PERSISTENT EMISSION DURING SUB-EDDINGTON THERMONUCLEAR BURSTS

The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts involves subtraction of the pre-burst emission as background. This approach implicitly assumes that the persistent flux remains constant throughout the burst. We reanalyzed 332 photospheric radius expansion bursts observed from 40 sources by the Rossi X-ray Timing Explorer, introducing a multiplicative factor $f_a$ to the persistent emission contribution in our spectral fits. We found that for the majority of spectra the best-fit value of $f_a$ is significantly greater than 1, suggesting that the persistent emission typically increases during a burst. Elevated $f_a$ values were not found solely during the radius expansion interval of the burst, but were also measured in the cooling tail. The modified model results in a lower average value of the $\chi^2$ fit statistic, indicating superior spectral fits, but not yet to the level of formal statistical consistency for all the spectra. We interpret the elevated $f_a$ values as an increase of the mass accretion rate onto the neutron star during the burst, likely arising from the effects of Poynting-Robertson drag on the disk material. We measured an inverse correlation of $f_a$ with the persistent flux, consistent with theoretical models of the disc response. We suggest that this modified approach may provide more accurate burst spectral parameters, as well as offering a probe of the accretion disk structure.

THE PECULIAR GALACTIC CENTER NEUTRON STAR X-RAY BINARY XMM J174457-2850.3

The recent discovery of a milli-second radio pulsar experiencing an accretion outburst similar to those seen in low mass X-ray binaries, has opened up a new opportunity to investigate the evolutionary link between these two different neutron star manifestations. The remarkable X-ray variability and hard X-ray spectrum of this object can potentially serve as a template to search for other X-ray binary/radio pulsar transitional objects. Here we demonstrate that the transient X-ray source XMM J174457-2850.3 near the Galactic center displays similar X-ray properties. We report on the detection of an energetic thermonuclear burst with an estimated duration of ~2 hr and a radiated energy output of ~5E40 erg, which unambiguously demonstrates that the source harbors an accreting neutron star. It has a quiescent X-ray luminosity of Lx~5E32 erg/s and exhibits occasional accretion outbursts during which it brightens to Lx~1E35-1E36 erg/s for a few weeks (2-10 keV). However, the source often lingers in between outburst and quiescence at Lx~1E33-1E34 erg/s. This unusual X-ray flux behavior and its relatively hard X-ray spectrum, a power law with an index of ~1.4, could possibly be explained in terms of the interaction between the accretion flow and the magnetic field of the neutron star.

The neutron star transient and millisecond pulsar in M28: from sub-luminous accretion to rotation-powered quiescence

The X-ray transient IGR J18245–2452 in the globular cluster M28 contains the first neutron star (NS) seen to switch between rotation-powered and accretion-powered pulsations. We analyse its 2013 March–April 25 d long outburst as observed by Swift, which had a peak bolometric luminosity of ∼6 per cent of the Eddington limit (LEdd), and give detailed properties of the thermonuclear burst observed on 2013 April 7. We also present a detailed analysis of new and archival Chandra data, which we use to study quiescent emission from IGR J18245–2452 between 2002 and 2013. Together, these observations cover almost five orders of magnitude in X-ray luminosity (LX, 0.5–10 keV). The Swift spectrum softens during the outburst decay (photon index Γ from 1.3 above LX/LEdd = 10−2 to ∼2.5 at LX/LEdd = 10−4), similar to other NS and black hole transients. At even lower luminosities, LX/LEdd = [10−4–10−6], deep Chandra observations reveal hard (Γ = 1–1.5), purely non-thermal and highly variable X-ray emission in quiescence. We therefore find evidence for a spectral transition at LX/LEdd ∼ 10−4, where the X-ray spectral softening observed during the outburst decline turns into hardening as the source goes to quiescence. Furthermore, we find a striking variability pattern in the 2008 Chandra light curves: rapid switches between a high-luminosity ‘active’ state (LX ≃ 3.9 × 1033 erg s−1) and a low-luminosity ‘passive’ state (LX ≃ 5.6 × 1032 erg s−1), with no detectable spectral change. We put our results in the context of low-luminosity accretion flows around compact objects and X-ray emission from millisecond radio pulsars. Finally, we discuss possible origins for the observed mode switches in quiescence, and explore a scenario where they are caused by fast transitions between the magnetospheric accretion and pulsar wind shock emission regimes.

Strong neutrino cooling by cycles of electron capture and β− decay in neutron star crusts

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometre, 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, β(-) decay, involving neutron-rich nuclei at a typical depth of about 150 metres, 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. The 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.

Double-peaked thermonuclear bursts at the soft–hard state transition in the Rapid Burster

Long suspected to be due to unstable accretion events, the type II bursts exhibited by the Rapid Burster (RB, or MXB 1730-335) still lack an exhaustive explanation. Apart from type II bursts, the transient RB also shows the better-understood thermonuclear shell flashes known as type I bursts. In search of links between these two phenomena, we carried out a comprehensive analysis of all $\textit{Rossi X-ray Timing Explorer}$ observations of the RB and found six atypical type I bursts, featuring a double-peaked profile that is not due to photospheric radius expansion. The bursts appear in a phase of the outburst decay close to the onset of the type II bursts, when the source also switches from the high/soft to the low/hard state. We also report the discovery of a simultaneous low-frequency quasi-periodic oscillation present in the persistent emission as well as in the burst decaying emission. We discuss several scenarios to understand the nature of the peculiar bursts and of the accompanying oscillation, as well as their possible relation with each other and with the type II burst phenomenon. We favour a model consisting of two accretion channels, one polar and one equatorial, in a source viewed at low inclination. We also, for the first time, clearly identify the atoll nature of the RB.

DETERMINING NEUTRON STAR MASSES AND RADII USING ENERGY-RESOLVED WAVEFORMS OF X-RAY BURST OSCILLATIONS

Simultaneous, precise measurements of the mass $M$ and radius $R$ of neutron stars can yield uniquely valuable information about the still uncertain properties of cold matter at several times the density of nuclear matter. One method that could be used to measure $M$ and $R$ is to analyze the energy-dependent waveforms of the X-ray flux oscillations seen during some thermonuclear bursts from some neutron stars. These oscillations are thought to be produced by X-ray emission from hotter regions on the surface of the star that are rotating at or near the spin frequency of the star. Here we explore how well $M$ and $R$ could be determined by generating, and analyzing using Bayesian techniques, synthetic energy-resolved X-ray data that we produce assuming a future space mission having 2--30 keV energy coverage and an effective area of 10 m$^2$, such as the proposed \textit{LOFT} or \textit{AXTAR} missions. We find that if the hot spot is within 10$^\circ$ of the rotation equator, both $M$ and $R$ can usually be determined with an uncertainty of about 10% if there are $10^6$ total counts from the spot, whereas waveforms from spots within 20$^\circ$ of the rotation pole provide no useful constraints. These constraints can usually be achieved even if the burst oscillations vary with time and data from multiple bursts must be used to obtain 10$^6$ counts from the hot spot. This is therefore a promising method to constrain $M$ and $R$ tightly enough to discriminate strongly between competing models of cold, high-density matter.

EVIDENCE FOR ACCRETION RATE CHANGE DURING TYPE I X-RAY BURSTS

The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts involves subtraction of the pre-burst emission as background. This approach implicitly assumes that the persistent flux remains constant throughout the burst. We reanalyzed 332 photospheric radius expansion bursts observed from 40 sources by the Rossi X-Ray Timing Explorer, introducing a multiplicative factor fa to the persistent emission contribution in our spectral fits. We found that for the majority of spectra the best-fit value of fa is significantly greater than 1, suggesting that the persistent emission typically increases during a burst. Elevated fa values were not found solely during the radius expansion interval of the burst, but were also measured in the cooling tail. The modified model results in a lower average value of the χ2 fit statistic, indicating superior spectral fits, but not yet to the level of formal statistical consistency for all the spectra. We interpret the elevated fa values as an increase of the mass accretion rate onto the neutron star during the burst, likely arising from the effects of Poynting–Robertson drag on the disk material. We measured an inverse correlation of fa with the persistent flux, consistent with theoretical models of the disk response. We suggest that this modified approach may provide more accurate burst spectral parameters, as well as offering a probe of the accretion disk structure.

A bright thermonuclear X-ray burst simultaneously observed withChandraand RXTE

The prototypical accretion-powered millisecond pulsar SAX J1808.4-3658 was observed simultaneously with Chandra-LETGS and RXTE-PCA near the peak of a transient outburst in November 2011. A single thermonuclear (type-I) burst was detected, the brightest yet observed by Chandra from any source, and the second-brightest observed by RXTE. We found no evidence for discrete spectral features during the burst; absorption edges have been predicted to be present in such bursts, but may require a greater degree of photospheric expansion than the rather moderate expansion seen in this event (a factor of a few). These observations provide a unique data set to study an X-ray burst over a broad bandpass and at high spectral resolution (lambda/delta-lambda=200-400). We find a significant excess of photons at high and low energies compared to the standard black body spectrum. This excess is well described by a 20-fold increase of the persistent flux during the burst. We speculate that this results from burst photons being scattered in the accretion disk corona. These and other recent observations of X-ray bursts point out the need for detailed theoretical modeling of the radiative and hydrodynamical interaction between thermonuclear X-ray bursts and accretion disks.

THEFERMI-GBM X-RAY BURST MONITOR: THERMONUCLEAR BURSTS FROM 4U 0614+09

Thermonuclear bursts from slowly accreting neutron stars (NSs) have proven difficult to detect, yet they are potential probes of the thermal properties of the neutron star interior. During the first year of a systematic all-sky search for X-ray bursts using the Gamma-ray Burst Monitor (GBM) aboard the Fermi Gamma-ray Space Telescope we have detected 15 thermonuclear bursts from the NS low-mass X-ray binary 4U 0614+09, when it was accreting at nearly 1% of the Eddington limit. We measured an average burst recurrence time of 12+/-3 d (68% confidence interval) between March 2010 and March 2011, classified all bursts as normal duration bursts and placed a lower limit on the recurrence time of long/intermediate bursts of 62 d (95% confidence level). We discuss how observations of thermonuclear bursts in the hard X-ray band compare to pointed soft X-ray observations, and quantify such bandpass effects on measurements of burst radiated energy and duration. We put our results for 4U 0614+09 in the context of other bursters and briefly discuss the constraints on ignition models. Interestingly, we find that the burst energies in 4U 0614+09 are on average between those of normal duration bursts and those measured in long/intermediate bursts. Such a continuous distribution in burst energy provides a new observational link between normal and long/intermediate bursts. We suggest that the apparent bimodal distribution that defined normal and long/intermediate duration bursts during the last decade could be due to an observational bias towards detecting only the longest and most energetic bursts from slowly accreting NSs.

MILLIHERTZ QUASI-PERIODIC OSCILLATIONS AND THERMONUCLEAR BURSTS FROM TERZAN 5: A SHOWCASE OF BURNING REGIMES

We present a comprehensive study of the thermonuclear bursts and millihertz quasi-periodic oscillations (mHz QPOs) from the neutron star (NS) transient and 11 Hz X-ray pulsar IGR J17480-2446, located in the globular cluster Terzan 5. The increase in burst rate that we found during its 2010 outburst, when persistent luminosity rose from 0.1 to 0.5 times the Eddington limit, is in qualitative agreement with thermonuclear burning theory yet opposite to all previous observations of thermonuclear bursts. Thermonuclear bursts gradually evolved into a mHz QPO when the accretion rate increased, and vice versa. The mHz QPOs from IGR J17480-2446 resemble those previously observed in other accreting NSs, yet they feature lower frequencies (by a factor ~3) and occur when the persistent luminosity is higher (by a factor 4-25). We find four distinct bursting regimes and a steep (close to inverse cubic) decrease of the burst recurrence time with increasing persistent luminosity. We compare these findings to nuclear burning models and find evidence for a transition between the pure helium and mixed hydrogen/helium ignition regimes when the persistent luminosity was about 0.3 times the Eddington limit. We also point out important discrepancies between the observed bursts and theory, which predicts brighter and less frequent bursts, and suggest that an additional source of heat in the NS envelope is required to reconcile the observed and expected burst properties. We discuss the impact of NS magnetic field and spin on the expected nuclear burning regimes, in the context of this particular pulsar.

THE MASS AND RADIUS OF THE NEUTRON STAR IN THE BULGE LOW-MASS X-RAY BINARY KS 1731–260

Measurements of neutron star masses and radii are instrumental for determining the equation of state of their interiors, understanding the dividing line between neutron stars and black holes, and for obtaining accurate statistics of source populations in the Galaxy. We report here on the measurement of the mass and radius of the neutron star in the low-mass X-ray binary KS 1731-260. The analysis of the spectroscopic data on multiple thermonuclear bursts yields well-constrained values for the apparent angular area and the Eddington flux of the source, both of which depend in a distinct way on the mass and radius of the neutron star. The binary KS 1731-260 is in the direction of the Galactic bulge, allowing a distance estimate based on the density of stars in that direction. Making use of the Han & Gould model, we determine the probability distribution over the distance to the source, which is peaked at 8 kpc. Combining these measurements, we place a strong upper bound on the radius of the neutron star, R <= 12 km, while confining its mass to M <= 1.8 M_sun.