J.-M. Hameury1 - D. Barret2 - J.-P. Lasota3 -
J. E. McClintock4 - K. Menou5, - C. Motch1 - J.-F. Olive2 - N.
Webb2
1 - UMR 7550 du CNRS, Observatoire de Strasbourg, 11 rue de
l'Université, 67000 Strasbourg, France
2 - Centre d'Etude Spatiale des Rayonnements, 9 Av. du Colonel
Roche, BP 4346, 31028 Toulouse Cedex 4, France
3 - UMR 7095 du
CNRS, Institut d'Astrophysique de Paris, 98bis Boulevard Arago,
75014 Paris, France
4 - Harvard-Smithsonian Center for
Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
5 -
Department of Astronomy, PO Box 3818, University of Virginia,
Charlottesville, VA 22937, USA
Received 30 september 2002 / Accepted 22 November 2002
Abstract
We report on XMM-Newton observations of GRO J1655-40 and
GRS 1009-45, which are two black hole X-ray transients currently
in their quiescent phase. GRO J1655-40 was detected with a 0.5-10 keV luminosity of
erg s-1 with 20%
accuracy, taking into account both statistical errors and
uncertainty on the distance. This luminosity is comparable to a
previous Chandra measurement, but ten times lower than the 1996
ASCA value, most likely obtained when the source was not yet in a
true quiescent state. Unfortunately, XMM-Newton failed to detect
GRS 1009-45. A stringent upper limit of
erg s-1 was derived by combining data from the EPIC-MOS and PN
cameras.
The X-ray spectrum of GRO J1655-40 is very hard as it can be
fitted with a power law model of photon index
.
Similarly hard spectra have been observed from other systems;
these rule out coronal emission from the secondary or disk flares
as the origin of the observed X-rays.
On the other hand, our
observations are consistent with the predictions of the disc
instability model in the case that the accretion flow forms an
advection dominated accretion flow (ADAF) at distances less than a
fraction (
0.1-0.3) of the circularization radius. This
distance corresponds to the greatest extent of the ADAF that is
thought to be possible.
Key words: accretion, accretion discs - instabilities - black hole physics - X-ray: binaries - stars: individual: GRO J1655-40 - stars: individual: GRS 1009-45
Soft X-ray transients - SXTs (sometimes called X-ray novae) are
semi-detached binaries in which the accreting (primary) star is a
black hole (BH) or a neutron star, and the mass-losing secondary
is usually a late type star. These systems typically brighten in
X-rays by as much as 107 in a week and then decay back into
quiescence over the course of a year. The maximum outburst
luminosities
seen in SXTs are typically
(0.2-1) of the Eddington luminosity
.
Successive outbursts are usually separated by years to decades of
quiescence (see e.g. Tanaka & Shibazaki 1996; Chen et al.
1997, for reviews).
![]() |
D | Ref. | M1 | n21 |
![]() |
![]() |
Reference (X-rays) | |
(hr) | (kpc) | (![]() |
||||||
XTE J1118+480 | 4.1 |
![]() |
1 |
![]() |
0.1 | |||
GRO J0422+32 | 5.1 |
![]() |
2 |
![]() |
1.6 | 30.9 | - | Chandra (13) |
GRS 1009-45 | 6.9 |
![]() |
3 |
![]() |
1.1 | <30.9 | XMM-Newton - this paper | |
A0620-00 | 7.8 |
![]() |
4 |
![]() |
1.9 | 30.6 | 2.07 +0.28-0.19 | Chandra (14) |
0.16 | 30.7 | ROSAT (15) | ||||||
GS 2000+25 | 8.3 |
![]() |
5 |
![]() |
11 | 30.4 | - | Chandra (13) |
XTE J1859+226 | 9.2 |
![]() |
6 |
![]() |
||||
Nova Mus 1991 | 10.4 |
![]() |
7 |
![]() |
2.6 | 31.6 |
![]() |
XMM-Newton (16) |
Nova Oph 1977 | 12.5 |
![]() |
5 |
![]() |
2.8 | <33.0 | ROSAT (17) | |
4U1543-47 | 27.0 |
![]() |
8 |
![]() |
3.5 | <31.5 | Chandra (13) | |
XTE J1550-564 | 37.0 |
![]() |
9 |
![]() |
3.9 | 33.0 |
![]() |
Chandra (14) |
GRO J1655-40 | 62.9 |
![]() |
10 |
![]() |
6.7 | 31.5 |
![]() |
Chandra (14) |
31.8 |
![]() |
XMM-Newton - this paper | ||||||
<3 | 33.1 | 0.7 +0.21-0.4 | ASCA (18) | |||||
V4641 Sgr | 67.6 |
![]() |
11 |
![]() |
0.5 | |||
V404 Cyg | 155.3 | ![]() |
12 |
![]() |
5.4 | 33.6 |
![]() |
Chandra (14) |
11 +3-4 | 33.1 | 2.1 +0.3-0.3 | ASCA (19) | |||||
10 | 33.0 | 1.9 +0.6-0.3 | Beppo-SAX (20) |
Quiescent states of SXTs are very interesting for at least two
reasons. First, they provide the strongest evidence now available
for the existence of stellar-mass black holes. The detection of
the secondary in quiescence allows a determination of the mass
function of the binary system, which is an absolute lower limit on
the primary mass M1. The mass function exceeds 3
in
eight systems; in five others, constraints on the inclination
angle of the system and the secondary mass result in values for
M1 that are in the range 5-10
(Narayan et al.
2002). Thus there are now thirteen SXTs with primary
masses that exceed the maximum stable mass of a neutron star
(
3
).
The second, more mundane reason why quiescent states are interesting is that quiescence provides a strong test of outburst-cycle models. SXT outburst cycles are well described by the disc instability model (Dubus et al. 2001; see Lasota 2001 for a complete review of the disc instability model). In quiescence, the disc is non-steady (a property too often forgotten by too many authors). Its temperature and viscosity are low, and the disc is unable to transport all of the mass supplied by the secondary to the primary; the mass of the disc slowly builds up and the temperature rises finally to the hydrogen ionization temperature. At this moment, the disc becomes thermally and viscously unstable due to strong opacity variations. Propagating heat fronts bring the entire disc into a hot state; in this outburst state, the mass transfer rate is large, and the disc empties until it cannot sustain this regime any longer; it then returns to quiescence.
Even though the disc instability paradigm is widely accepted and
its particular realizations are often successful (requiring
sometimes additional assumptions, see e.g. Esin et al.
2000), the physics of the quiescent state is still rather
controversial. As a prime example, the origin of "viscosity" in
this state is not really known (see e.g. Gammie & Menou
1998; Menou 2002; Lasota 2002), and according
to the simplest version of the disc-instability model (the version
originally used to explain dwarf-nova outbursts; Lasota
2001) the very long SXT recurrence times require unusually
low values of the viscosity parameter .
However, Dubus et
al. (2001) showed that the combined effects of disc
irradiation during outbursts (King & Ritter 1998) and disc
truncation during quiescence (Menou et al. 2000) result
quite effortlessly in long recurrence times for standard values of
the viscosity parameter, making the disc instability model for
SXTs a working and testable hypothesis.
Narayan et al. (1996; see also Lasota et al. 1996
and Narayan et al. 1997a) noticed that spectra of quiescent
SXTs cannot be produced by geometrically thin, optically thick
accretion discs and suggested that the inner flow in such systems
forms an advection dominated accretion flow ADAF. In any case,
according to the disc instability model a disc extending down to
the neutron star surface or the last stable orbit around a black
hole, can be in a cold, neutral state everywhere only for
vanishingly low accretion rates close to the central object
(Lasota 1996): the maximum accretion rate onto the compact
object would be
g s-1, where
is the
Schwarzschild radius (Hameury et al. 1998), much too low
to produce the X-ray luminosity observed in quiescent SXTs.
On the
other hand the ADAF model reproduces well the observed
luminosities, spectra (Narayan et al. 1997a; Quataert &
Narayan 1999) and observed delays between the rises to
outburst in optical and X-rays (Hameury et al. 1997).
As first pointed out by Narayan et al. (1997b), the presence of the radiatively inefficient ADAFs in quiescent SXTs allows one to compare black holes with compact bodies endowed with material surfaces. Because bodies such as neutron stars must re-emit the heat left over in the accretion flow accumulating at their surface they should be brighter than black holes accreting at the same rate since in this case the residual energy is lost forever past the event horizon. X-ray observations of quiescent SXTs showed that this is indeed the case (Narayan et al. 1997b; Menou et al. 1999b; Garcia et al. 2001). However, the ADAF origin of X-rays in quiescent SXTs has been questioned by several authors (see e.g. Bildsten & Rutledge 2001; Nayakshin & Svensson 2001). High quality X-ray observations can solve this controversy.
In this paper, we report on XMM-Newton observations of GRO J1655-50 (AO-1 GTO) and GRS 1009-45 (AO-1 GO).We briefly describe previous observations of these transients in Sect. 2; we present our observations in Sect. 3, and we discuss their implications in Sect. 4.
X-ray observations of quiescent SXTs prior to XMM-Newton and
Chandra have yielded only marginal detections or upper limits,
with very poor constraints on the spectral shape. Among the 13
BHSXTs with known orbital periods, A 0620-00 was the only
one to be detected by ROSAT (McClintock et al., 1995);
ASCA detected GRO J1655-40 (Asai et al. 1998;
Ueda et al. 1998) and V404 Cyg (Narayan et al.
1997a). Finally, among the 5 BHSXTs observed by BeppoSAX
(Campana et al. 2001), only V404 Cyg has been
detected with a 1-10 keV unabsorbed luminosity of 1033 erg s-1.
XMM-Newton and Chandra have already qualitatively changed this
situation, as can be clearly seen from Table 1, which
gives spectral data for most of the 13 SXTs with known orbital
period, and with primary masses exceeding 3 .
The X-ray
luminosities have been computed with the most recent and accurate
distance estimates and associated uncertainties. When no estimate
on the error on distance was available, we arbitrarily assumed it
to be 30% for the distances given by Barret et al.
(1996; compare relative measurements in their Table 1), or
50% in all other cases.
It is now evident that many quiescent SXTs have been detected with luminosities in the range 1031-1033 erg s-1, the faintest being those with the shortest orbital periods. It is also clear that SXTs have very hard spectra with values of the power law spectral index that are often less than 2. Or equivalently, their spectra can be fitted by a bremsstrahlung model with temperatures larger than 10 keV (Kong et al. 2002).
GRO J1655-40 was observed by XMM-Newton on August 30, 2001 with a
40 ksec exposure time. Due to a strong solar flare, only 19 ks
could be used; yet, GRO J1655-40 was clearly detected as one of
the brightest field sources, at a position coincident with that of
the optical counterpart (within 0.3
), with a total source
count rate as determined by the automatic pipeline analysis of
.
SAS version 5.3 was used to calculate the source and background
spectra obtained by extracting data from a circle of approximately
20
radius centered on the nominal source position, and
from an annulus with inner and outer radii
20
and
100
respectively. 74, 73 and 260 counts were
detected in the 20
radius circle in the MOS1, MOS2 and PN
cameras respectively; 618 (MOS1), 403 (MOS2) and 1225 (PN) counts
were detected in the larger area used to determine the background
rate. The spectra were analyzed with Xspec v11. Due to the small
number of counts, in particular in the MOS cameras, the
C-statistics (Cash 1979) was used to determine the best
fits and errors; channels were grouped in such a way that there
were at least 5 counts in each bin. The three cameras, MOS1, MOS2,
and PN were used simultaneously to determine the best fit
parameters.
![]() |
Figure 1:
Top: XMM-Newton spectrum of GRO J1655-40, fitted by a
power law. Bottom: residuals after model subtraction from the
data, in units of 1 ![]() |
Open with DEXTER |
![]() |
Figure 2:
Confidence level at 68, 90 and 99% of the fit given in
Fig. 1. ![]() ![]() |
Open with DEXTER |
Model | ![]() |
![]() |
kT | C statistic | d.o.f. | goodness | F0.5-10 |
(1021 cm-2) | (keV) | 10-14 (erg cm-2s-1) | |||||
power law | 10.0 +1.2-0.7 | 1.54 +1.02-0.72 | 63.8 | 72 | 0.12 | 5.3 | |
power law | 6.7 (fixed) | 1.30 +0.34-0.41 | 64.5 | 73 | 0.11 | 4.8 | |
bremsstrahlung | 6.7 (fixed) | >9 | 64.5 | 73 | 0.12 | 4.1 |
Table 2 summarizes our results. A power law is a good
fit to the results, yielding column density
cm-2, a power law index
,
and normalization
,
poorly constrained because of the large error in
.
If one
uses a hydrogen column density fixed to the value of
cm-2 determined by optical observations (Hynes
et al. 1998), one obtains
,
and
.
All errors given here are
90% confidence limits. Figure 1 shows the EPIC
spectrum for all three instruments, as well as the best fit and
residuals, in units of 1 sigma. Note that these residuals,
although giving a rough estimate of the quality of the fit, cannot
be used to determine the goodness-of-fit by computing a
,
since the C-stat has been used. Note also that, for the sake of
clarity, data in Fig. 1 have been rebinned, so
that more data points were used for performing the spectral
fitting. Figure 2 gives the 68.3, 90 and 99%
confidence contours obtained for the power law index and
normalization.
The total absorbed flux is then 3.96
erg cm-2 s-1 in the
0.5-10 keV range (90% confidence
error bars), corresponding to an unabsorbed flux of
erg cm-2 s-1 and
unabsorbed luminosity of
erg s-1 for a distance
of 3.2 kpc.
We also fitted the spectrum with a thermal bremsstrahlung. As
expected for a hard spectrum with such a photon index, the fit is
good, with a formal best fit value of the temperature kT = 110 keV, indeed unphysical, and very poorly constrained; the 90%
confidence level lower limit on kT is 9 keV, and the 99%
confidence level lower limit is 4.1 keV.
![]() |
Figure 3: Net count rate per 2000 s bins for GRO J1655-40 in the PN camera. Large errors are produced by the strong solar flare. |
Open with DEXTER |
A search for variability yielded negative results for significant variations on time scales of hours. Figure 3 shows the source intensity as a function of time, with 2000 s time bins; the data are consistent with a constant source. The data quality is quite poor during the solar flare that occurs after the first 10 ks, and lasts about 25 ks, with shorts lulls.
GRS 1009-45 was observed on May 30, 2002, during a 20 ksec
exposure. The source was not detected by the automatic pipeline
processing system. Assuming a power law spectrum with photon index
2, and a hydrogen column density
cm-2, as given by optical observations (della Valle et al.
1997), and analyzing simultaneously all three cameras
together, we can set an upper limit for the 0.5-10 keV
unabsorbed flux of
erg cm-2 s-1 (95%
confidence), corresponding to an upper limit on the luminosity of
erg s-1 for a distance of 5 kpc. To determine
this upper limit, we fixed the photon index and column density in
Xspec and calculated the goodness-of-fit for various normalization
factors, with the
statistics (using the C-statistics
yields an almost identical upper limit, as expected). The data are
compatible with a zero flux (we obtain a goodness-of-fit of 63%);
we can set an upper limit on the flux of
erg cm-2 s-1, for which the goodness-of-fit is 0.95, i.e. in
only 5% of the simulated cases did we find a C statistic value
worse than the value actually obtained. This upper limit is almost
independent of the assumed photon power law index.
It has been suggested (Bildsten & Rutledge 2001) that
X-rays detected from quiescent black-hole SXTs are not produced by
accretion, but are instead produced in the secondaries' hot
stellar coronae. There is now growing evidence that this is not
the case because the X-ray flux is too large (especially the relative to the bolometric luminosity of the secondary) and/or
because the spectrum is too hard to be emitted by a stellar corona
(Kong et al. 2002). As we shall see, the flaring disc
model also faces similar difficulties, whereas the ADAF
explanations can still account for the observations.
The spectra of all black hole transients appear to be quite hard, GRO J1655-40 being the hardest, and A0620-00 having the steepest power law index; fits using bremsstrahlung models also lead to very high temperatures. This confirms previous findings (Kong et al. 2002) that coronal emission from the secondary is ruled out. Einstein surveys of late type stars have shown that the coronal temperatures of isolated stars do not significantly exceed 1 keV, and reach at most a few keV in the case of giant RS CVn stars (Schmitt 1990). (Coronal emission is additionally ruled out by the high X-ray luminosities; see below.)
It is also most unlikely that such spectra could be produced by
emission from magnetic loops above the disc (see e.g. Nayakshin &
Svensson 2001) that would very probably result in spectra
similar to those of stellar coronae. Moreover, this model implies
that the accretion disc of quiescent dwarf novae should emit a
comparable soft X-ray flux; these systems are indeed X-ray bright,
with luminosities in the range
1030 - 1032 erg s-1,
but they have deep eclipses indicating that the X-ray emitting
region is very close to the white dwarf; the residual X-ray
luminosity during eclipses is very low (e.g. of the order of
erg s-1 for OY Car, Ramsay et al.
2001).
![]() |
Figure 4:
Quiescent luminosities of SXTs containing black holes
(filled circles: Chandra and XMM-Newton data from
Table 1), or a neutron star (open circles: data from Narayan
et al. 2002, except for SAX J1808.4-3658 for which we
used XMM-Newton data of Campana et al. 2002, who found
the source in deeper quiescence than Wijnands 2002a). Also
shown are the maximum coronal X-ray luminositiy (the "FGM-limit",
Flemming et al. 1989) for a 1 ![]() |
Open with DEXTER |
Figure 4 shows the luminosities of SXTs containing
black holes compared to those in which the accreting object is a
neutron star. We plot the quiescent luminosity (in erg s-1)
as a function of the orbital period of the binary. Here, the data
for BH systems have been taken from Table 1; when several
measurements were available, we selected Chandra or XMM-Newton
values (for GRO J1655-40, the XMM-Newton value given in this paper
is used). The NS data are from Narayan et al. (2002). The
error bars reflect uncertainties on distance only, which are
rather large, and are usually the most important contributor to
the error on luminosity. For SAX J1808.4-3658, we assumed d =
2.5
0.5 kpc (in't Zand et al. 2001); for EXO
0748-676, d = 10 kpc with an arbitrary error of 5 kpc (Gottwald
et al. 1986); for 4U 2129+47, d = 6
2 kpc (Nowak
et al. 2002); for MXB 1659-298, d = 12
2.5 kpc
(Wijnands 2002b); for H 1608-52, d = 3.6 kpc with a
typical error of order 1 kpc (Wachter et al. 2002); for
Cen X-4, d = 1.2
0.3 kpc (Chevalier et al. 1989;
for Aql X-1, d = 5
1 kpc (Rutledge et al. 2001).
Our result confirms the well established fact that systems
containing black holes are, for a given orbital period, much
dimmer than those containing a neutron star (i.e. in which X-ray
bursts have been detected), a property predicted by models in
which the inner disc is truncated and replaced by an ADAF close to
the black hole (Menou et al. 1999b).
The
plot was proposed by Lasota &
Hameury (1998) to assure, when one compares black-hole
systems with those containing neutron stars, that both types of
systems have comparable accretion rates at the outer limit of
the ADAF. Sutaria et al. (2002) pointed out recently that
for various reasons one should expect a "greater spread" of
(the mass-transfer rate) at a given
.
This
is obviously true. However, if SXT outbursts are due to a
dwarf-nova-type disc instability, then in quiescence the accretion
rate at the truncation radius depends only weakly on the
mass-transfer rate (e.g. Menou et al. 1999a; Lasota
2000; Dubus et al. 2001; but see also Menou
2002 who discusses uncertainties linked to our poor
understanding of viscosity in the quiescent state and a possible
alternative solution). Indeed, the accretion rate at the
transition radius can be written as (Lasota 2000)
To obtain Eq. (1) one assumes that according to the disc
instability model the quiescent disc is non-steady and
,
where
is the outer radius of the
geometrically thin accretion disc. The pre-Chandra and
pre-XMM-Newton data (for three systems) satisfied a
relation (dashed line in Fig. 4),
the slope of which suggested that it might result from Eq.
(4) (Lasota 2000).
The new Chandra and XMM-Newton observations show a more
complicated picture. First, the quiescent X-ray luminosity of GRO
J1655-40 is now an order of magnitude lower than it was when
observed by ASCA. In fact the quiescent luminosities of all three
systems (A0620-00, GRO J1655-40 and V404 Cyg) used by Lasota
(2000) to determine the 1.77 slope of the
relation have been revised significantly (see Table
1 and discussion below). In the case of A0620-00 and
V404 Cyg there is little doubt that they are in true quiescence.
Among the new sources, XTE J1550-564 was observed between two
outbursts (see Kong et al. 2002 and references therein) so
that this system was very unlikely to be in its true quiescent
state during the Chandra observations.
The quiescent luminosity might still be described by Eq. (4) but with ft depending on the system. Such a
conclusion is not unexpected since Menou et al. (1999a)
have already noticed that although
for most
systems, V404 Cyg required
.
The value of
should be predicted by the model describing the "evaporation" of
an accretion disc into an ADAF, but we lack such a model.
The flux we obtained for GRO J1655-40 is consistent with that found by Chandra (Kong et al. 2002), the difference being of order of 50% for the unabsorbed flux in the 0.2-10 keV range, and can be essentially accounted for by statistical fluctuations; we do not find significant short time scale variability, which means that these sources are relatively stable. However, both XMM-Newton and Chandra flux measurements are one order of magnitude below that of ASCA (Asai et al. 1998; Ueda et al. 1998). This observation, however, was performed a month before the 1996 outburst and roughly 9 months after the previous (August 1995) eruption. Clearly, the system was not in its quiescent state as it takes more than 30 years (the recurrence time) to refill the disc emptied by the outburst. Matter was then transferred from the secondary at a rate much higher than during the true quiescence (Esin et al. 2000).
Short time scale variability has been detected transients during quiescence containing black holes, such as V404 Cyg (Wagner et al. 1994; Kong et al. 2002), as well as those containing neutron stars, such as Aql X-1 (Rutledge et al. 2002). Such variability is not a surprise, but cannot be significantly constrained in our data due to the low count rate.
We have observed and detected GRO J1655-40 during a 40 ksec
XMM-Newton observation, with a luminosity of
erg
s-1; GRS 1009-45 was undetected, leading to a relatively
small upper limit on luminosity of
erg s-1 for
a distance of 5 kpc. These observations are consistent with the
disc instability model if the accretion disc is truncated at a
radius of 0.1-0.3 times the circularization radius, where the
accretion flow forms an ADAF. The spectrum appears to be quite
hard, much harder than one would expect from coronal emission from
the secondary star or from the accretion disc itself. The quality
of the data is however not sufficient to constrain and distinguish
between various types of optically thin flows in the vicinity of
the black hole.
There is no sign of variability on a time scale of hours, and the flux measured by XMM-Newton from GRO J1655-40 is consistent with the Chandra value (Kong et al. 2002), and about 10 times lower than that which was measured with ASCA by Asai et al. (1998) and Ueda et al. (1998) between two outbursts, indicating that the system was not fully in quiescence at that time, and that the mass transfer rate from the secondary was high.
Acknowledgements
We are grateful to R. Narayan for discussions and comments on the manuscript. This work was supported in part by a grant from the Centre National d'Etudes Spatiales.