Probing the Crust of the Neutron Star in EXO 0748–676
X-ray observations of quiescent X-ray binaries have the potential to provide insight into the structure and the composition of neutron stars. EXO 0748–676 had been actively accreting for over 24 yr before its outburst ceased in late 2008. Subsequent X-ray monitoring revealed a gradual decay of the quiescent thermal emission that can be attributed to cooling of the accretion-heated neutron star crust. In this work, we report on new Chandra and Swift observations that extend the quiescent monitoring to 5 yr post-outburst. We find that the neutron star temperature remained at 117 eV between 2009 and 2011, but had decreased to 110 eV in 2013. This suggests that the crust has not fully cooled yet, which is supported by the lower temperature (95 eV) measured 4 yr prior to the accretion phase in 1980. Comparing the data to thermal evolution simulations reveals that the apparent lack of cooling between 2009 and 2011 could possibly be a signature of convection driven by phase separation of light and heavy nuclei in the outer layers of the neutron star.
Subject headings:accretion, accretion disks — binaries: eclipsing — stars: individual (EXO 0748–676) — stars: neutron — X-rays: binaries
Transient neutron star low-mass X-ray binaries (LMXBs) are excellent laboratories for increasing our understanding of the structure and the composition of neutron stars, and how matter behaves under extreme physical conditions. In these binary star systems a neutron star is accompanied by a late-type star that overflows its Roche lobe and transfers matter to an accretion disk. This matter is rapidly accreted onto the neutron star during outburst episodes, whereas little or no matter reaches the compact primary during quiescent intervals.
These accretion cycles have a profound effect on the interior properties of neutron stars. They cool during quiescence as they lose thermal energy via photons emitted from their surface and neutrinos escaping from their crust and core (e.g., Yakovlev et al., 2003; Page et al., 2006; Steiner & Reddy, 2009; Schatz et al., 2014). However, neutron stars can re-gain thermal energy during accretion outbursts.
The accretion of matter compresses the crust of a neutron star, which causes successive electron captures, neutron emissions and pycno-nuclear fusion reactions. Together, these processes deposit an energy of 2 MeV per accreted nucleon in the crustal layers (e.g., Haensel & Zdunik, 1990a, b; Gupta et al., 2007; Steiner, 2012). This energy is thermally conducted both towards the stellar core and surface, and can effectively maintain the interior temperature of the neutron star at K. This temperature is set by the energy injected during its historic accretion activity and the efficiency of the neutrino cooling processes (e.g., Brown et al., 1998; Colpi et al., 2001; Yakovlev et al., 2003; Wijnands et al., 2013).
During quiescent episodes, thermal X-rays from the surface of the neutron star may be detected. This allows a measurement of its temperature, which can encode valuable information about its interior properties. Of particular interest are observations obtained shortly after the cessation of an outburst; heating due to accretion may lift the temperature of the crust well above that of the stellar core and the subsequent cooling may be observable once back in quiescence (Wijnands et al., 2001; Ushomirsky & Rutledge, 2001; Rutledge et al., 2002).
Indeed, dedicated X-ray monitoring of six LMXBs (KS 1731–260, MXB 1659–29, XTE J1701–462, EXO 0748–676, IGR J17480–2446, and MAXI J0556–332), revealed that the temperature of the neutron star decreased for years following the cessation of accretion, consistent with the heating/cooling paradigm (e.g., Wijnands et al., 2002, 2004; Cackett et al., 2006, 2008, 2010a; Degenaar et al., 2009, 2011a, 2011b, 2013b; Díaz Trigo et al., 2011; Fridriksson et al., 2010, 2011, J. Homan et al., in preparation). Comparison with thermal evolution simulations has yielded valuable insight into the thermal and transport properties of neutron star crusts (Shternin et al., 2007; Brown & Cumming, 2009; Degenaar et al., 2011a; Page & Reddy, 2013; Turlione et al., 2013).
Despite these successes, interpretation of the quiescent thermal emission and crustal cooling is complicated by the question whether accretion onto the neutron star fully comes to a halt. Searching for (strong) non-thermal emission in the quiescent X-ray spectrum, irregular quiescent X-ray variability, or optical/UV signatures of the quiescent accretion stream can shed light on whether residual accretion occurs (see e.g., Cackett et al., 2010b, 2011, 2013b; Degenaar & Wijnands, 2012; Bernardini et al., 2013, for recent studies).
1.1. Exo 0748–676
The neutron star LMXB EXO 0748–676 was discovered almost three decades ago (Parmar et al., 1985). The detection of X-ray eclipses indicates that the binary is viewed at high inclination (), and led to a measurement of the orbital period ( hr; Parmar et al., 1986; Wolff et al., 2009). The source displays thermonuclear X-ray bursts, which allows for a distance determination ( kpc; e.g., Galloway et al., 2008).
EXO 0748–676 was first detected in outburst in 1984 with EXOSAT (Reynolds et al., 1999), and was serendipitously detected in quiescence with Einstein in 1980 (Parmar et al., 1986; Garcia & Callanan, 1999). The source remained in outburst for 24 yr and during this time the flux was moderately stable with occasional excisions to higher and lower fluxes. However, its activity suddenly ceased in 2008 September (Wolff et al., 2008a, b; Hynes & Jones, 2008; Torres et al., 2008). Subsequent monitoring with Swift, Chandra and XMM-Newton revealed a relatively hot neutron star that gradually cooled over time (Degenaar et al., 2009, 2011b; Díaz Trigo et al., 2011).
In this work we report on new X-ray observations of EXO 0748–676 to further monitor how the accretion-heated crust cools, and to search for signs of continued low-level accretion. We also re-analyze the Einstein data obtained in 1980 to measure the pre-outburst temperature of the neutron star. We then compare the entire data set to crust cooling simulations.
2. Observations and Data Analysis
Table 1 gives an overview of all new Chandra and Swift data of EXO 0748–676. A list of earlier X-ray observations obtained during the quiescent phase can be found in Degenaar et al. (2009), Degenaar et al. (2011b) and Díaz Trigo et al. (2011). To benefit from the latest calibration updates and to ensure a homogeneous analysis, these quiescent X-ray observations were re-reduced and re-analyzed in this work.
|Swift||90420001||2010 May 1||8.5|
|Swift||90420002||2010 Jun 22||7.7|
|Swift||90420003||2010 Aug 22||6.5|
|Swift||90420004||2010 Aug 24||4.0|
|Chandra||11060||2010 Oct 20||27.2|
|Swift||90420005||2010 Oct 21||10.3|
|Swift||90420006||2010 Dec 23||4.0|
|Swift||90420007||2010 Dec 24||6.3|
|Swift||90420008||2011 Feb 22||10.5|
|Swift||31272051||2011 Jun 28||6.0|
|Swift||31272052||2011 Jun 29||4.0|
|Chandra||12414||2011 Jul 2/3||38.1|
|Swift||31272053||2013 Mar 23||7.8|
|Chandra||14663||2013 Aug 1||42.9|
Note. – The listed exposure times and count rates are not corrected for eclipses. The count rates refer to the 0.3–10 keV energy band. Errors represent 90% confidence levels.
2.1. New Chandra Observations
We obtained three new Chandra observations of EXO 0748–676 between 2010 October and 2013 August (Table 1). The setup was similar to previous quiescent observations of the source, using the S3 chip of the Advanced CCD Imaging Spectrometer (ACIS; Garmire et al., 2003). The ACIS-S3 CCD was operated in a 1/8 sub-array during the 2010 observation and in a 1/4 sub-array during the 2011/2013 observations. We reduced the data using the ciao package (ver. 4.5). All observations were free from background flaring.
We extracted source events using a circular region of radius centered on EXO 0748–676, and a annulus was used to estimate the background (excluding a circular region with a radius centered on a faint point source). Count rates and light curves were extracted using the task dmextract. Source and background spectra, as well as the corresponding response files, were created using the meta task specextract. We used grppha to group the spectra into bins of at least 20 photons. About 3800–6400 net source events were collected for each Chandra observation.
The light curves obtained from the individual observations clearly show the presence of eclipses at times consistent with the ephemeris of Wolff et al. (2009). During the eclipses, the X-ray emission from the neutron star is temporarily blocked by the companion. To calculate the time-averaged X-ray fluxes excluding the eclipsed epochs, we therefore reduced the exposure times of all Chandra observations by 500 s per eclipse (which is the observed duration of the eclipses in quiescence; Bassa et al., 2009).
Note. – Quoted uncertainties represent 90% confidence intervals. represents the total unabsorbed model flux and gives the thermal flux (0.5–10 keV). The ratio of these two fluxes, , represents the fractional contribution of the thermal emission to the total model flux. represents the unabsorbed thermal flux in the 0.01–100 keV energy band. The 0.5–10 keV total model luminosity () and the thermal bolometric luminosity () are for a distance of kpc.
2.2. New Swift Observations
EXO 0748–676 has been observed with the X-Ray Telescope (XRT; Burrows et al., 2005) onboard Swift many times since it transitioned to quiescence (see Wolff et al., 2008a, b; Degenaar et al., 2009, 2011b). Here we report on 11 new observations that were carried out between 2010 May and 2013 March (see Table 1). All XRT data were obtained in the photon-counting mode.
The Swift observations were reduced using the heasoft suite (ver. 6.13). We first processed the raw XRT data using the xrtpipeline. Employing xselect, we then extracted source events using a circular region with a radius of , which optimized the signal to noise ratio for the observed count rates (cf. Evans et al., 2007). A surrounding annular region of radius was used as a background reference. We obtained 50–350 net source counts per observation.
The Chandra and XMM-Newton observations provide superior spectral quality compared to the short Swift/XRT pointings. However, the Swift data offer unique dense sampling that allows us to closely track intensity variations occurring on a time scale of days–weeks. For the Swift data we therefore focused our attention on the long-term quiescent light curve. The count rates extracted for each observation were corrected for losses due to bad pixels and bad columns. Due to the relatively low count rates, it was not possible to identify eclipses in the data. We therefore checked the good time intervals against the ephemeris of Wolff et al. (2009); if (part of) eclipses were expected, we reduced the exposure times accordingly.
2.3. Archival XMM-Newton Observations
Following its return to quiescence, EXO 0748–676 was observed with the European Photon Imaging Camera (EPIC; Strüder et al., 2001; Turner et al., 2001) onboard XMM-Newton on four occasions (Bassa et al., 2009; Díaz Trigo et al., 2011). We reduced and analyzed these data using the sas package (ver. 13.0). After reprocessing with emproc and epproc, and removing background flaring events, we extracted spectra and light curves using the task evselect. A circular region with a radius of was used for the source and a -radius circular region placed on an empty part of the CCD served as the background.
Response files were generated using arfgen and rmfgen. The spectra and response files of the three detectors (MOS1, MOS2, and PN) were combined with the tool epicspeccombine.
2.4. Pre-outburst Einstein Observation
EXO 0748–676 was serendipitously detected with the Imaging Proportional Counter (IPC) onboard the Einstein observatory on 1980 May 22 (ObsID 7708; Parmar et al., 1986; Garcia & Callanan, 1999). The source intensity during the ks exposure was , whereas the local background is estimated at . We created spectra and response files using XSelect. Source events were collected from a region of 20 pixels (), and background events using a source-free region of twice that size. Extrapolating the ephemeris of Wolff et al. (2009) back to the time of the Einstein observation suggests that the source may have been eclipsed for s. We therefore reduced the exposure time by this amount.
There is little reported information in the literature about cross-calibration between the Einstein/IPC and the current generation of X-ray instruments (e.g., the Crab was too bright for the IPC; Kirsch et al., 2005). A study of nova-like variables yielded strong consistency with ROSAT/PSPC results, lending credence to the low-energy response of the IPC (Verbunt et al., 1997). Systematic uncertainties are estimated to be at the 10%-15% level (e.g., Fabricant et al., 1984; David et al., 1990). To account for (cross-)calibration uncertainties, we included a conservative 20% systematic error for the Einstein data in our spectral fits.
2.5. Spectral Fitting Procedures
The Chandra, XMM-Newton and Einstein spectra (each with their own response and background files) were fitted simultaneously in the 0.3–10 keV energy range using xspec (ver. 12.8; Arnaud, 1996). The Einstein data and the last three XMM-Newton observations were all dominated by background noise above an energy of keV (see also Díaz Trigo et al., 2011). For these observations we therefore excluded the data above 3 keV.
Expanding on previous work (Degenaar et al., 2009, 2011b; Díaz Trigo et al., 2011), we concentrated on fitting the data to a combination of a neutron star atmosphere model (nsatmos; Heinke et al., 2006), and a power-law (pegpwrl; to model any possible non-thermal emission). To account for the neutral hydrogen absorption along the line of sight, , we included the tbabs model adopting the vern cross-sections and wilm abundances (Verner et al., 1996; Wilms et al., 2000).
Since EXO 0748–676 is viewed at a relatively high inclination and hence the absorption along the line of sight could be variable (see, e.g., Cackett et al., 2013a), we allowed to change between the different observations. Only for the Einstein spectrum we could not obtain good constraints on . For this data set we therefore fixed to the average value obtained for the Chandra and XMM-Newton spectra (Table 2).
For the nsatmos model we fixed the normalization at unity, i.e., we assumed that the emitting area was the same during all observations and corresponded to the entire neutron star surface. When fixing the neutron star mass and radius to and 10 km, no acceptable fit could be obtained (see also Degenaar et al., 2009, 2011b; Díaz Trigo et al., 2011). Both parameters were therefore left free to find the best fit values (assuming that and did not change over time, they were tied between the different data sets). The source distance was not well constrained when left to vary freely and drove the neutron star mass and radius to unrealistic values. We therefore fixed the distance to kpc in all spectral fits (Galloway et al., 2008). The neutron star temperature was allowed to vary freely.
For the non-thermal emission component it was not possible to constrain the power-law index for each observation individually. We therefore chose to fix this parameter for the entire data set to the value obtained from the first XMM-Newton observation, which provided the best constraints (; Degenaar et al., 2011b). The power-law normalization was free to vary.
From the fit results we calculated the effective neutron star temperature as seen by an observer at infinity, . Here, is the gravitational redshift factor, with being the Schwarzschild radius, the gravitational constant and the speed of light. Using cflux in xspec, we determined the unabsorbed 0.5–10 keV flux and the thermal bolometric flux. The latter was estimated by extrapolating the nsatmos model component to the 0.01–100 keV energy range. All fluxes were converted to luminosities by assuming kpc.
|Fit Parameter (unit)||Value|
|Exponential decay, base level fixed|
|Decay time, (days)|
|Constant offset, (eV)||94.6 fixed|
|Exponential decay, base level free|
|Decay time, (days)|
|Constant offset, (eV)|
|Power-law decay, base level fixed|
|Constant offset, (eV)||94.6 fixed|
|Power-law decay, base level free|
|Constant offset, (eV)|
|Power-law decay, no constant offset|
|Constant offset, (eV)||fixed|
Note. – The data was fit to an exponential decay of the form , and a power-law decay shaped as . The fixed constant offset of eV corresponds to the temperature inferred from the pre-outburst Einstein observation. The power-law fit without constant offset is included to allow for a comparison with the literature. In all fits the start of the cooling curve, was set to 2008 September 5 (MJD 54714; Degenaar et al., 2009).
3.1. X-Ray Spectral Evolution
Simultaneously fitting the Chandra, XMM-Newton and Einstein spectra to a combined neutron star atmosphere and power-law model as described in Section 2.5, resulted in a reduced chi-squared value of for 1950 degrees of freedom (dof) with a p-value of . The best fit yielded a neutron star mass of , and a radius of km. The uncertainty in dominates the errors in all other parameters (see also the discussion in Díaz Trigo et al., 2011). For the error calculation we therefore fixed at the best fit value (whereas was still free). The results are summarized in Table 2.
We find that the values of obtained for the XMM-Newton data are systematically lower than for the Chandra observations. This likely arises due to cross-calibration uncertainties (Kirsch et al., 2005; Degenaar et al., 2011b; Tsujimoto et al., 2011).
The first Chandra and first XMM-Newton observation (obtained within 2 months after the outburst) required the presence of a hard spectral component, although its contribution to the total 0.5–10 keV flux was small (15%; Degenaar et al., 2011b; Díaz Trigo et al., 2011). This non-thermal component is not statistically required in subsequent observations, although it may still account for up to 10% of the 0.5–10 keV flux (Table 2). The quiescent spectra of EXO 0748–676 are thus strongly dominated by soft, thermal emission throughout the quiescent phase (see also Section 4.3).
The neutron star temperature gradually decreased by 11 eV within the first year after the outburst (2008 October till 2009 June), but showed little change in the subsequent two years. Indeed, requiring the temperatures of the 6 observations obtained between 2009 June and 2011 July to be the same results in a good fit ( for 1955 dof, ), with an average temperature of eV. It therefore appeared that the neutron star crust had fully cooled well within two years of entering quiescence (Degenaar et al., 2011b; Díaz Trigo et al., 2011). However, our Chandra observation obtained in 2013 August ( yr after the outburst) shows a decrease in temperature by 7 eV compared to 2009–2011 (Table 2). This is illustrated by Figure 1, which shows the evolution of the temperature over time. Forcing the 2013 temperature to be the same as during the previous 6 observations results in a poor fit ( for 1956 dof, ).
Since we found little variation in the absorption along the line of sight, the drop in temperature observed in 2013 appears to be genuine, and hence indicates continued or accelerated cooling of the neutron star (see Section 3.2). This is supported by the fact that the temperature inferred from the pre-outburst Einstein observation (95 eV) is lower than in our latest Chandra observation (see Table 2 and Figure 1). We note that different values of , , , and may shift the absolute temperatures by a few () eV. However, the observed relative temperature change, i.e., the crust cooling curve, is not affected by these systematic uncertainties.
3.2. Crust Cooling Curve Fits
To characterize the temperature evolution of EXO 0748–676 and compare it with that of other sources, we fitted the crust cooling curve to an exponential decay of the form . Here, is a normalization, the -folding time, and represents the quiescent base level. The reference time is set to the presumed end of the outburst, 2008 September 5 (MJD 54714; Degenaar et al., 2009). Fixing eV to the 1980 level does not provide an acceptable fit (yielding for 9 dof, ). The fit improves when the base level is allowed to vary ( for 10 dof, ), but the obtained value of eV is higher than observed with Chandra in 2013 ( eV). An exponential decay therefore may not be a good representation of the crust cooling curve of EXO 0748–676 (dashed and dashed-dotted lines in Figure 1).
We also fitted the data to a power-law function of the form , which is the theoretically expected shape for a cooling crust (e.g., Brown & Cumming, 2009). Fixing eV to the Einstein level yields a good fit with a decay index of ( for 9 dof, ). With the constant offset allowed to vary we found eV and ( for 10 dof, ). The base level we obtain is consistent within the errors with the Einstein measured temperature. The power-law decay fits are indicated by the solid and dotted lines in Figure 1.
The results of the exponential and power-law decay fits are summarized in Table 3. In the literature the crust cooling curves are often fit to a power-law decay without a constant offset. Since the inclusion of a base level changes the resulting slope significantly, we also reference a fit without a constant level to allow for a direct comparison with other sources.
3.3. Swift Quiescent X-Ray Light Curve
The 2008–2013 Swift/XRT count rate light curve of EXO 0748–676 is shown in Figure 2. The source intensity changed gradually from to within the first year after the outburst, but showed little variation thereafter. Simple decay fits with an exponential ( for 56 dof) or a power-law function ( for 56 dof) leveling off to a constant suggests that the quiescent light curve flattened 1 yr post-outburst (Figure 2). This may indicate an episode of relatively constant intensity, as was also hinted by our analysis of the Chandra/XMM-Newton spectral data (Section 3.1). The Swift spectral data is not of sufficient quality to test whether the last data point (obtained in 2013 March) supports the lower temperature seen during the 2013 August Chandra observation.
There is also more stochastic variability among the data points, although the 90% error bars largely overlap. Since we corrected the count rates for dead zones on the CCD and the possible occurrence of eclipses (Section 2.2), we suspect that the variations are due to photon statistics (this is perhaps supported by the fact that the largest outliers in the light curve concern short observations, i.e., which collected a small number of photons). It is clear that there are no strong intensity flares as are sometimes seen in quiescent neutron star LMXBs (see Section 4.3).
3.4. Thermal Evolution Simulations
We briefly explored thermal evolution simulations to gauge the properties of the neutron star in EXO 0748–676. Using the physical model described in Brown & Cumming (2009) and Medin & Cumming (2014), we calculated source-specific cooling curves assuming an outburst duration of 24 yr, and an outburst mass-accretion rate of (e.g., Degenaar et al., 2011b). The model parameters that are then adjusted to match the data are the core temperature , and the impurity parameter . The latter parametrizes how organized the structure of the ion lattice is and hence determines the thermal conductivity (e.g., Itoh & Kohyama, 1993; Brown & Cumming, 2009).
In recent years, evidence has accumulated that there is more heat generated in the crusts of neutron stars than is currently accounted for by nuclear heating models (e.g., Brown & Cumming, 2009; Degenaar et al., 2011a, 2013a; Schatz et al., 2014). We therefore allowed for the inclusion of an additional heat source placed at a column depth of . Table 4 gives an overview of the model calculations.
We find that reproducing the observed high temperatures of EXO 0748–676 requires a crust impurity parameter of , and an additional heat source of MeV nucleon. This model produces the data reasonably well (dashed curve in Figure 3). It suggests that the source has so far moved along a continuous cooling track and will reach its pre-outburst quiescent level many years from now.
Interestingly, the occurrence of a possible plateau of slow cooling in EXO 0748–676 starting 1 yr post-outburst (Sections 3.1 and 3.3) resembles recent calculations of Medin & Cumming (2014). These authors showed how cooling curves are affected by a convective heat flux that arises when light and heavy nuclei in the crust separate out in liquid and solid phases (Horowitz et al., 2007; Medin & Cumming, 2011). Therefore, we also performed model calculations with compositionally-driven convection taken into account. This introduces one extra fit parameter, , which is the mass fraction of oxygen at the base of the liquid ocean (for details, see Medin & Cumming, 2014).
We find that the inclusion of convection leads to a plateau of slow cooling between 150–750 days post-outburst, broadly consistent with the data (solid curve in Figure 3). This arises because the compositionally-driven convection transports heat inward, temporarily slowing the cooling in the crustal layers where the phase separation occurs. We note that the model including convection is not statistically preferred over the non-convective case.
|Parameter (unit)||No Convection||Convection|
Note. – represents the impurity parameter, the additional heat energy (placed at a column depth of ), the mass fraction of oxygen at the base of the liquid ocean, and the core temperature. We assumed an outburst duration of 24 yr, and mass-accretion rate of . The equation of state used in these simulations results in and km (for details, see Medin & Cumming, 2014).
4.1. Crustal Cooling in EXO 0748–676
Our new Chandra and Swift observations of EXO 0748–676 extend the quiescent monitoring to 4.9 yr after the cessation of its very long (24 yr) active period. We find that the neutron star temperature gradually decreased during this time, consistent with expectations for cooling of the accretion-heated neutron star crust. In the first year of quiescence, between 2008 and 2009, the temperature decreased from to eV. It then hovered around eV for at least yr till 2011, but our most recent observation obtained in 2013 indicates a further decrease in temperature to eV.
Despite the high inclination of the binary, there are no indications that the lower temperature in 2013 is due to a changing absorption column density, such as possibly seen in MXB 1659–29 (Cackett et al., 2013a). Whereas the apparent lack of temperature evolution after 2009 led to the suggestion that the neutron star crust restored equilibrium with the core (Degenaar et al., 2011b; Díaz Trigo et al., 2011), the new data presented in this work suggests that cooling is still ongoing and that a further decrease in temperature may be expected.
4.2. A Signature of Convection?
The possible “plateau” of stalled cooling starting yr post-outburst is reminiscent of the crust cooling curve of XTE J1701–462. That source too appeared to level off within 2 yr of entering quiescence (Fridriksson et al., 2011), but Page & Reddy (2013) predicted that after a temporary plateau an accelerated temperature decay would occur, which seems to be borne out by more recent observations (J. K. Fridriksson et al., in preparation). XTE J1701–462 experienced a relatively short (1.6 yr) but very bright accretion phase (an average flux near the Eddington limit). As a result of this vigorous heating, the temperature in the crust likely did not reach a steady state but rather had double peaked profile, which would naturally give rise a plateau. This is in sharp contrast to EXO 0748–676, which was active for 24 yr at relatively low X-ray flux (5% of Eddington), implying that the crust had ample time to reach a steady state profile (cf. Brown & Cumming, 2009).
Another mechanism that may cause a plateau in the cooling curve is a convective heat flux driven by chemical separation of light and heavy nuclei in the outer layers of the neutron star (Horowitz et al., 2007; Medin & Cumming, 2011, 2014). Inclusion of the inward heat transport by compositionally-driven convection in the model calculations for EXO 0748–676 leads to an episode of slow cooling that is broadly consistent with the observations. The crust cooling curve of EXO 0748–676 may thus bear an imprint of this process, although the data can also be satisfactory modeled without convection. Perhaps another possibility is that a crustal shell of rapid neutrino cooling as recently identified by Schatz et al. (2014) is connected to the period of stalled cooling. This process is highly temperature-sensitive and may be related to the fact that a plateau appears to be seen only in the two hottest crust-cooling neutron stars EXO 0748–676 and XTE J1701–462. However, this could also be an observational bias, since these two sources were more intensely monitored than the others (see Degenaar et al., 2011b, for a comparison). Further theoretical investigation is required to grasp the implications of this newly identified cooling process on neutron star crust cooling curves.
It is of note that the model calculations of EXO 0748–676 require rather high values for the impurity parameter (), and the additional heat ( MeV nucleon), to keep the crust hot as long as observed. In contrast, the crust cooling curves of KS 1731–260, MXB 1659–29, XTE J1701–462, and IGR J17480–2446 suggested an impurity parameter of order unity (Brown & Cumming, 2009; Degenaar et al., 2011a; Page & Reddy, 2013; Medin & Cumming, 2014), consistent with expectations from molecular dynamics simulations (Horowitz et al., 2008). In fact, taking into account allowed ranges in mass, radius, and accretion rate, Brown & Cumming (2009) set an upper limit of for KS 1731–260 and MXB 1659–29. The higher value that we find here could imply that the crust of EXO 0748–676 has a more impure (i.e., less organized) structure than the other neutron stars, although it is unclear why that would be the case. Moreover, the obtained value of is sensitive to other model parameters.
If we allow for a higher mass-accretion rate, e.g., , the crust temperature rises and therefore we require a lower impurity parameter (), and less extra heat ( MeV nucleon). This mass-accretion rate is higher than inferred from X-ray observations (; e.g., Degenaar et al., 2011b; Díaz Trigo et al., 2011), but not implausible. There are large uncertainties in determining the accretion rate from X-ray observations, in particular for high-inclination systems such as EXO 0748–676 when part of the central X-ray source may be blocked from our line of sight, causing to be underestimated. Nevertheless, even for this higher accretion rate remains considerably larger than found for the other sources. Another possible way of keeping the crust in EXO 0748–676 hot for a long time is residual accretion during quiescence.
4.3. On the Possibility of Quiescent Accretion
Our interpretation of the observations of EXO 0748–676 in the crustal heating/cooling framework relies on the assumption that accretion onto the neutron star stopped when the source transitioned to quiescence. It is not straightforward to test this hypothesis with observations. Low-level accretion may generate a thermal emission spectrum like that of a cooling neutron star (Zampieri et al., 1995; Soria et al., 2011). However, the measured temperature would then reflect that of the stellar surface that is continuously heated by residual accretion and masks the interior temperature of the neutron star. We therefore searched for signatures of quiescent accretion in EXO 0748–676.
X-ray monitoring with Swift has revealed X-ray flares in several quiescent neutron stars, e.g., XTE J1701–462, Aql X-1, Cen X-4, KS 1741–293, GRS 1741.9–2853, and SAX J1750.8–2900 (e.g., Bernardini et al., 2013; Degenaar & Wijnands, 2009, 2013; Degenaar et al., 2012b; Fridriksson et al., 2011; Wijnands & Degenaar, 2013; Coti Zelati et al., 2014). During these X-ray flares the XRT count rate increased for several days by a factor of 2 to even 10 for some of these sources. A corresponding hardening of the X-ray spectrum is observed and suggests that these flares are possibly caused by a spurt of low-level accretion. EXO 0748–676 was monitored with Swift roughly once per month for 10 ks between 2008 and 2011. We did not detect any irregular X-ray variability or flaring events such as seen in other sources. Regular Swift monitoring has therefore not revealed any indications of ongoing low-level accretion in EXO 0748–676. However, accretion flares appear to be short-lived events (lasting days), and could therefore easily be missed (e.g., Degenaar & Wijnands, 2009, 2013; Fridriksson et al., 2011; Wijnands & Degenaar, 2013; Coti Zelati et al., 2014).
The first Chandra and XMM-Newton observations of EXO 0748–676 (obtained in 2008, within 2 months after the outburst end) both showed the presence of non-thermal emission, albeit contributing only 15% to the total unabsorbed 0.5–10 keV flux (whereas this is 50% in some other neutron stars such as SAX J1808.4–3658, Swift J1749.4–2807, and EXO 1745–248; e.g., Heinke et al., 2009; Degenaar & Wijnands, 2012; Degenaar et al., 2012a). Optical spectroscopic and photometric observations performed shortly after the transition to quiescence hinted the presence of an accretion disk that could allow for continued accretion onto the neutron star (Bassa et al., 2009; Hynes & Jones, 2009). However, optical spectroscopy and Doppler tomography performed one year later did not show evidence for an accretion disk any more (Ratti et al., 2012). Any contribution from non-thermal X-ray emission also remained low at this time (Section 3.1). Finally, there are no dips or other features in the quiescent X-ray light curves that might evidence the presence of a residual accretion stream or remnant disk (see also Díaz Trigo et al., 2011).
We conclude that there are no obvious signs of ongoing accretion in the quiescent state of EXO 0748–676, particularly not 1 yr after the outburst ended. Given the optical signatures of a quiescent accretion disk and the presence of non-thermal X-ray emission we cannot exclude, however, that matter was falling onto the neutron star shortly after the outburst appeared to have ended.
- affiliation: Hubble Fellow
- slugcomment: Received 2014 March 9; accepted 2014 June 26; published 2014 ???
- For all four XMM-Newton observations we found that the fits results for the combined PN/MOS spectra were consistent with fits to the separate PN and MOS spectra.
- We note that in Degenaar et al. (2011b) it was found that the first XMM-Newton observation had an elevated temperature compared to adjacent Chandra and Swift observations. In those spectral fits was fixed between the different data sets. Here we leave free, which does not yield an elevated temperature for the first XMM-Newton observation (see Table 2 and Figure 1).
- Degenaar et al. (2011b) noted that the Einstein flux reported by Garcia & Callanan (1999) was consistent within the errors with that inferred from the 2010 Chandra data. However, Garcia & Callanan (1999) used a different physical model to fit the spectrum, which may introduce biases. Fitting the data in tandem with the post-outburst Chandra and XMM-Newton observations suggests that the pre-outburst temperature was lower than currently seen, provided the caveats mentioned in Section 2.4.
- Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17
- Bassa, C. G., Jonker, P. G., Steeghs, D., & Torres, M. A. P. 2009, MNRAS, 399, 2055
- Bernardini, F., Cackett, E. M., Brown, E. F., et al. 2013, MNRAS, 436, 2465
- Brown, E. F., Bildsten, L., & Rutledge, R. E. 1998, ApJ, 504, L95
- Brown, E. F., & Cumming, A. 2009, ApJ, 698, 1020
- Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Science Reviews, 120, 165
- Cackett, E. M., Brown, E. F., Cumming, A., et al. 2013a, ApJ, 774, 131
- —. 2010a, ApJ, 722, L137
- Cackett, E. M., Brown, E. F., Degenaar, N., et al. 2013b, MNRAS, 433, 1362
- Cackett, E. M., Brown, E. F., Miller, J. M., & Wijnands, R. 2010b, ApJ, 720, 1325
- Cackett, E. M., Fridriksson, J. K., Homan, J., Miller, J. M., & Wijnands, R. 2011, MNRAS, 414, 3006
- Cackett, E. M., Wijnands, R., Linares, M., et al. 2006, MNRAS, 372, 479
- Cackett, E. M., Wijnands, R., Miller, J. M., Brown, E. F., & Degenaar, N. 2008, ApJ, 687, L87
- Colpi, M., Geppert, U., Page, D., & Possenti, A. 2001, ApJ, 548, L175
- Coti Zelati, F., Campana, S., D’Avanzo, P., & Melandri, A. 2014, MNRAS, 438, 2634
- David, L. P., Arnaud, K. A., Forman, W., & Jones, C. 1990, ApJ, 356, 32
- Degenaar, N., Brown, E. F., & Wijnands, R. 2011a, MNRAS, 418, L152
- Degenaar, N., Patruno, A., & Wijnands, R. 2012a, ApJ, 756, 148
- Degenaar, N., & Wijnands, R. 2009, A&A, 495, 547
- —. 2012, MNRAS, 422, 581
- Degenaar, N., & Wijnands, R. 2013, in IAU Symp. 290, Feeding Compact Objects: Accretion on All Scales, ed. C. M. Zhang, T. Belloni, M. Méndez, & S. N. Zhang (Cambridge: Cambridge Univ. Press), 113
- Degenaar, N., Wijnands, R., Cackett, E. M., et al. 2012b, A&A, 545, A49
- Degenaar, N., Wijnands, R., & Miller, J. M. 2013a, ApJ, 767, L31
- Degenaar, N., Wijnands, R., Wolff, M. T., et al. 2009, MNRAS, 396, L26
- Degenaar, N., Wolff, M. T., Ray, P. S., et al. 2011b, MNRAS, 412, 1409
- Degenaar, N., Wijnands, R., Brown, E. F., et al. 2013b, ApJ, 775, 48
- Díaz Trigo, M., Boirin, L., Costantini, E., Méndez, M., & Parmar, A. 2011, A&A, 528, 150
- Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 379
- Fabricant, D., Rybicki, G., & Gorenstein, P. 1984, ApJ, 286, 186
- Fridriksson, J. K., Homan, J., Wijnands, R., et al. 2010, ApJ, 714, 270
- —. 2011, ApJ, 736, 162
- Galloway, D. K., Özel, F., & Psaltis, D. 2008, MNRAS, 387, 268
- Garcia, M. R., & Callanan, P. J. 1999, AJ, 118, 1390
- Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, Jr., G. R. 2003, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 4851, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. J. E. Truemper & H. D. Tananbaum, 28
- Gupta, S., Brown, E. F., Schatz, H., Möller, P., & Kratz, K.-L. 2007, ApJ, 662, 1188
- Haensel, P., & Zdunik, J. L. 1990a, A&A, 229, 117
- —. 1990b, A&A, 227, 431
- Heinke, C. O., Jonker, P. G., Wijnands, R., Deloye, C. J., & Taam, R. E. 2009, ApJ, 691, 1035
- Heinke, C. O., Rybicki, G. B., Narayan, R., & Grindlay, J. E. 2006, ApJ, 644, 1090
- Horowitz, C. J., Berry, D. K., & Brown, E. F. 2007, Phys. Rev. E, 75, 066101
- Horowitz, C. J., Dussan, H., & Berry, D. K. 2008, Phys. Rev. C, 77, 045807
- Hynes, R., & Jones, E. 2008, ATel, 1816
- Hynes, R. I., & Jones, E. D. 2009, ApJ, 697, L14
- Itoh, N., & Kohyama, Y. 1993, ApJ, 404, 268
- Kirsch, M. G., Briel, U. G., Burrows, D., & et al. 2005, in Proc. of the SPIE, Vol. 5898, UV, X-Ray, and Gamma-Ray Space Instrumentation for Astronomy XIV, ed. O. H. W. Siegmund, 22–33
- Medin, Z., & Cumming, A. 2011, ApJ, 730, 97
- —. 2014, ApJ, 783, L3
- Page, D., Geppert, U., & Weber, F. 2006, Nuclear Physics A, 777, 497
- Page, D., & Reddy, S. 2013, Physical Review Letters, 111, 241102
- Parmar, A. N., White, N. E., Giommi, P., & Gottwald, M. 1986, ApJ, 308, 199
- Parmar, A. N., White, N. E., Giommi, P., et al. 1985, IAU Circ., 4039
- Ratti, E. M., Steeghs, D. T. H., Jonker, P. G., et al. 2012, MNRAS, 420, 75
- Reynolds, A. P., Parmar, A. N., Hakala, P. J., et al. 1999, A&AS, 134, 287
- Rutledge, R. E., Bildsten, L., Brown, E. F., et al. 2002, ApJ, 580, 413
- Schatz, H., Gupta, S., Möller, P., et al. 2014, Nature, 505, 62
- Shternin, P. S., Yakovlev, D. G., Haensel, P., & Potekhin, A. Y. 2007, MNRAS, 382, L43
- Soria, R., Zampieri, L., Zane, S., & Wu, K. 2011, MNRAS, 410, 1886
- Steiner, A. W. 2012, Phys. Rev. C, 85, 055804
- Steiner, A. W., & Reddy, S. 2009, Phys. Rev. C, 79, 015802
- Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18
- Torres, M. A. P., Jonker, P. G., Steeghs, D., & Seth, A. C. 2008, ATel, 1817
- Tsujimoto, M., Guainazzi, M., Plucinsky, P. P., et al. 2011, A&A, 525, A25
- Turlione, A., Aguilera, D. N., & Pons, J. A. 2013, ArXiv:1309.3909
- Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27
- Ushomirsky, G., & Rutledge, R. E. 2001, MNRAS, 325, 1157
- Verbunt, F., Bunk, W. H., Ritter, H., & Pfeffermann, E. 1997, A&A, 327, 602
- Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487
- Wijnands, R., & Degenaar, N. 2013, MNRAS, 434, 1599
- Wijnands, R., Degenaar, N., & Page, D. 2013, MNRAS, 432, 2366
- Wijnands, R., Guainazzi, M., van der Klis, M., & Méndez, M. 2002, ApJ, 573, L45
- Wijnands, R., Homan, J., Miller, J. M., & Lewin, W. H. G. 2004, ApJ, 606, L61
- Wijnands, R., Miller, J. M., Markwardt, C., Lewin, W. H. G., & van der Klis, M. 2001, ApJ, 560, L159
- Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914
- Wolff, M. T., Ray, P. S., & Wood, K. S. 2008a, ATel, 1736
- Wolff, M. T., Ray, P. S., Wood, K. S., & Hertz, P. L. 2009, ApJS, 183, 156
- Wolff, M. T., Ray, P. S., Wood, K. S., & Wijnands, R. 2008b, ATel, 1812
- Yakovlev, D. G., Levenfish, K. P., & Haensel, P. 2003, A&A, 407, 265
- Zampieri, L., Turolla, R., Zane, S., & Treves, A. 1995, ApJ, 439, 849