A&A 365, L267-L272 (2001)
K. Wu1,2 - R. Soria1 - M. J. Page1 - I. Sakelliou1 - S. M. Kahn3 - C. P. de Vries4
Send offprint request: K. Wu
1 -
Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, RH5 6NT, UK
2 -
Research Centre for Theoretical Astrophysics, School of Physics, University of Sydney, NSW 2006, Australia
3 -
Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027, USA
4 -
Space Research Organization of The Netherlands, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands
Received 2 October 2000 / Accepted 10 November 2000
Abstract
We report the results of preliminary analysis
of the XMM-Newton EPIC and RGS observations
of the candidate black-hole binary 3
between February and June 2000.
The observations covered
both the soft and the hard X-ray spectral states.
The hard-state spectra were dominated by a power-law component
with a photon index
.
The soft-state spectra consisted of a thermal component
with a multi-colour disk temperature
of
keV
and a power-law tail with
.
The model in which
the X-rays from 3 in the high-soft state
are powered by a strong stellar wind
from a massive companion
is not supported by
the small line-of-sight absorption
(
cm-2)
deduced from the RGS data.
The transition from the soft to the hard state
appears to be a continuous process
associated with the changes in the mass-transfer rate.
Key words: accretion, accretion disks - black hole physics - binaries: general - stars: early type - stars: winds, outflows - X-rays: general
3 (Leong et al. 1971) is a persistent X-ray source in the
Large Magellanic Cloud (LMC).
It has an orbital period of 1.70 d,
and a mass function
(Cowley et al. 1983).
The optical brightness of
indicates that the system has a massive companion
(van Paradijs et al. 1987).
The companion is often classified as a B3V star
(Cowley et al. 1983), although a B5IV spectral type has also been
suggested (Soria et al. 2001).
The non-detection of eclipses in the X-ray curve
implies that the orbital inclination of the system
is
(Cowley et al. 1983).
The inferred mass of the compact star is
(Paczynski 1983),
thus establishing that the system is a black-hole candidate (BHC).
Most BHCs show soft and hard X-ray spectral states.
Their X-ray spectrum in the soft state
generally consists of both a thermal and a power-law component.
The thermal component can be fitted by a blackbody spectrum
with a temperature
1 keV, and
it is interpreted as thermal emission from the inner accretion disk.
The power-law component is believed to be comptonised emission
from a disk corona (Sunyaev & Titarchuk 1980)
or from the high-speed infalling plasma
near the black-hole event horizon (Titarchuk & Zannias 1998).
The photon index
of the power law is
2.5-4.
In the hard state, the thermal component is insignificant.
The spectrum is dominated
by a flat (
), extended power-law component.
3 is found in the soft state most of the time. LMC X-1, another high-mass BHC in the LMC, has been seen in the soft state only. In contrast, Cyg X-1, the high-mass BHC in our Galaxy, tends to be in the hard state for the majority of the time. Such preference for either the hard or the soft state is not evident in the low-mass transient BHCs (e.g. GRO J1655-40). It is still unclear how mass transfer occurs in 3: whether via Roche-lobe overflow or capture of the stellar wind from the companion. Given the fact that a B3V star is unable to drive a strong stellar wind or to fill its Roche lobe if it is in a binary system with a 1.7-d orbital period, it is difficult to explain the observed X-ray luminosity of the system.
In this letter, we report on the XMM-Newton observations of 3 between 2000 February and June. The system was in an unusual hard state in April, and our observations span the transition from the soft state to the hard state and then back to a softer state. Making use of these new X-ray spectral data, we attempt to shed light on the question of mass transfer in the system (see also Soria et al. 2001).
3 was observed with
the European Photon Imaging Camera (EPIC; Turner et al. 2001),
the Reflection Grating Spectrograph (RGS; Brinkman et al. 2001) and
the Optical Monitor (OM; Mason et al. 2001) on board XMM-Newton,
during the Performance Verification Phase
between 2000 February and June.
Here we present part of the EPIC and RGS observations.
The OM observations are presented in Soria et al. (2001).
The logs of the EPIC-MOS1/PN and RGS observations used for this study
are shown in Tables 1 and 2 respectively.
| Revolution/ID | Instrument | Exposurea | Modeb |
| Rev0028/201 | PN | 08.1 ks | small window |
| Rev0028/301 | PN | 0.3 ks | full PN |
| MOS1 | 08.9 ks | partial RFS | |
| Rev0030/501 | MOS1 | 05.9 ks | partial W5 |
| Rev0041/101 | PN | 0.7 ks | full PN |
| Rev0041/401 | MOS1 | 03.5 ks | partial W2 |
| MOS1 | 03.6 ks | partial W5 | |
| MOS1 | 03.4 ks | partial W4 | |
| Rev0045/101 | PN | 01.6 ks | full PN |
| Rev0045/301 | MOS1 | 06.6 ks | partial W3 |
| Rev0066/101 | PN | 26.6 ksc | small window |
| Rev0092/201 | PN | 24.2 ks | small window |
| MOS1 | 20.3 ks | partial W4 |
| Revolution/ID | Instrumenta | Exposure | count rateb |
| Rev0028/201 | RGS2 (s1) | 17.6 ks |
|
| Rev0030/501 | RGS2 (tot) | 12.1 ks |
|
| Rev0045/101 | RGS1 (s1) | 06.4 ks |
|
| RGS2 (s1) | 06.4 ks |
|
|
| Rev0045/201 | RGS1 (s1) | 06.7 ks |
|
| RGS2 (s1) | 06.8 ks |
|
|
| Rev0045/301 | RGS1 (s1) | 03.0 ks |
|
| RGS1 (s3) | 03.0 ks |
|
|
| RGS1 (s5) | 03.0 ks |
|
|
| RGS2 (s2) | 00.3 ks |
|
|
| RGS2 (s4) | 03.0 ks |
|
|
| RGS2 (s6) | 00.3 ks |
|
|
| Rev0066/101 | RGS1 (s7) | 44.3 ks |
|
| RGS2 (s8) | 44.3 ks |
|
|
| Rev0092/101 | RGS1 (s2) | 51.9 ks |
|
| RGS2 (s2) | 51.9 ks |
|
|
| Rev0092/201 | RGS1 (s4) | 25.6 ks |
|
| RGS2 (s5) | 25.6 ks |
|
| a The labels "s#'' identifies multiple exposures during the |
| same revolution; for more details see the Preferred |
| Observation Sequence files of each revolution. | b In photons s-1. |
All EPIC exposures were taken with the "medium'' filter;
Some of the data suffer from pile-up
when the source was observed in the bright state (before Rev0066).
MOS exposures taken in the "full window'' mode are worst affected,
and so are excluded from the analysis.
"Partial window'' MOS exposures and "full window'' PN exposures
are, however, less affected.
In our analysis, we have removed the central pixels
(which are most affected by pile-up)
from the extraction regions,
for the partial MOS and full PN exposures that are affected.
The corresponding normalisations of the fitted spectral models
for these observations are therefore smaller than the true values,
and they are put in brackets in Table 3.
We are aware that
the photon index of the power-law component
is affected by pile-up.
For the Rev0028, 0030, 0041 and 0045 observations,
we estimate that the fitted values
are
0.15-0.20 smaller (harder) than the real value.
The "small window'' PN exposures
are not significantly affected by pile up,
and they are more reliable
for the determination of the power-law photon indices.
The data were processed using the September 2000 release of the SAS, except for the PN small window mode data, which were processed with the October 2000 SAS release. All spectral fits presented here were performed with the most up-to-date response matrices available at the end of October 2000 (mos1_medium_all_qe17_rmf3_tel5_15.rsp for EPIC-MOS1 and epn_fs20_sY9_medium.rmf for EPIC-PN). We are aware that the PN response matrix assumes a nominal event threshold of 20 adu, whereas a threshold of 23 adu was used in the Rev0028 and Rev0041 observations. This may have a small effect on the low-energy end of the spectrum and make the determination of the column density from the early PN data less reliable. We used the RGS spectra to obtain a better determination of the column density.
![]() |
Figure 1:
(First panel from the top): the XMM-Newton EPIC observations of 3, marked by solid vertical lines,
are shown and compared with the RXTE/ASM light curve. The three RXTE/PCA observations that confirmed the hard
state of the system (Boyd et al. 2000) are marked by dotted vertical lines. (Second panel): the
0.35-2.5 keV light curve from the RGS observations. (Third panel): the evolution of the fit temperature of the
thermal component in the EPIC-PN and MOS spectra. The open circles are the means of the fit disk-blackbody
temperatures
|
| Open with DEXTER | |
We consider a conventional model
consisting of an absorbed multi-colour disk blackbody
plus a power law,
wabs
zvfeabs
(diskbb+powerlaw) in XSPEC,
to fit the EPIC data.
As the fitting process is limited
by the reliability of the response matrices currently available,
more sophisticated models or the inclusion of additional features
(such as lines) are not warranted.
The generic models used here are therefore sufficient
for this preliminary analysis.
For the MOS observations, we consider only the channels
in the 0.3-8.0 keV energy range;
for the PN observations we extend the energy range
to 0.3-12.0 keV.
The data are binned with the FTOOLS task grppha,
such that the number of counts in each group of channels
is larger than 40.
This improves the signal-to-noise ratio
for energies
keV,
while it does not affect low-energy channels.
We have also checked that different
binning criteria give consistent results.
No systematic error has been added to the data.
We consider the line-of-sight Galactic absorption toward
LMC X-3 to be
cm-2
following Wilms et al. (2000), and
use it as the value for the (fixed) first absorption component.
The second absorption component is then
the intrinsic photoelectric absorption
within the binary system and the LMC.
We also fix both the iron abundance and the total metal abundance
to be 0.4 times the solar values, to account for the
lower metallicity of the LMC (Caputo et al. 1999).
The best fit parameters are shown in Table 3.
| Observation |
|
|
|
|
|
Rev0028/201 | PN |
|
0.93+0.01-0.01 | 4.32+0.24-0.25 | 2.71+0.09-0.08 |
|
1.31 (784) | Rev0028/301 | PN |
|
0.93+0.04-0.05 | (6.26+1.38-1.15) | 2.45+0.53-0.44 |
|
1.02 (192) | MOS1 (s7) |
|
0.88+0.06-0.06 | (0.65+0.24-0.15) | 2.53+0.41-0.35 |
|
0.75 (170) | Rev0030/501 | MOS1 (s7) |
|
0.86+0.02-0.02 | (6.10+0.74-0.65 | 2.45+0.14-0.14 |
|
1.28 (282) | Rev0041/101 | PN |
|
0.59+0.07-0.06 | (10.91+4.13-5.17) | 2.57+0.22-0.27 |
|
1.01 (259) | Rev0041/401 | MOS1 (s3) |
|
0.62+0.04-0.03 | 16.80+5.25-4.36 | 2.47+0.08-0.08 |
|
1.12 (260) | MOS1 (s5) |
|
0.64+0.03-0.02 | 16.65+3.85-3.37 | 2.33+0.08-0.09 |
|
1.22 (258) | MOS1 (s7) |
|
0.58+0.02-0.02 | 30.38+6.62-6.12 | 2.41+0.09-0.09 |
|
1.25 (258) | Rev0045/101 | PN |
|
0.76+0.03-0.02 | (9.00+1.46-1.33) | 2.32+0.14-0.17 |
|
1.09 (454) | Rev0045/301 | MOS1 (s7) |
|
0.74+0.02-0.02 | (8.93+1.22-1.13) | 2.28+0.12-0.14 |
|
1.02 (271) | Rev0066/101 | PN |
|
1.89+0.10-0.06 |
|
1.08 (54) | Rev0092/201 | PN |
|
0.22+0.01-0.01 | 766.2+81.6-62.0 | 2.05+0.02-0.02 |
|
1.15 (1000) | MOS1 (s7) |
|
0.26+0.01-0.01 | 380.1+58.3-50.1 | 1.97+0.03-0.04 |
|
1.49 (376) |
a The Galactic line-of-sight absorption (
cm-2) has been subtracted.
We assumed a metallicity of 0.4 times the solar value, i.e., Z = 0.008.
| Observationa |
|
|
|
|
|
|
Rev0028/201 | RGS2 (s1) |
|
0.55+0.02-0.02 | 183.8+17.1-15.6 | 2.7 |
|
1.38(485) | Rev0045/101 | RGS1 (s1) |
|
0.63+0.07-0.06 | 60.73+25.47-22.18 | 2.5 |
|
1.04(545) | RGS2 (s1) |
|
0.49+0.03-0.03 | 195.0+5.8-38.4 | 2.5 |
|
1.14(488) | Rev0045/201 | RGS1 (s1) |
|
0.67+0.08-0.06 | 51.16+22.14-19.01 | 2.5 |
|
1.08(543) | RGS2 (s1) |
|
0.50+0.03-0.03 | 178.5+23.7-40.9 | 2.5 |
|
1.05(484) | Rev0045/301 | RGS1 (s1) |
|
0.53+0.05-0.05 | 127.4+54.8-41.7 | 2.5 |
|
1.12(543) | RGS1 (s3) |
|
0.70+0.18-0.10 | 42.99+14.98-24.12 | 2.5 |
|
0.91(541) | RGS1 (s5) |
|
0.53+0.06-0.04 | 119.0+53.3-38.5 | 2.5 |
|
1.08(540) | RGS2 (s4) |
|
0.45+0.03-0.03 | 252.7+91.0-74.7 | 2.5 |
|
1.03(487) | Rev0066/101 | RGS1 (s7) |
|
0.07+0.03-0.03 | 3311+4171-3056 | 1.9 |
|
1.11(44) | Rev0092/101 | RGS1 (s2) |
|
0.33+0.01-0.02 | 211.2+29.6-25.3 | 2.05 |
|
1.27(538) | RGS2 (s2) |
|
0.27+0.02-0.02 | 478.1+103.8-81.9 | 2.05 |
|
1.25(485) | Rev0092/201 | RGS1 (s4) |
|
0.29+0.03-0.03 | 252.0+83.9-52.9 | 2.05 |
|
1.08(537) | RGS2 (s5) |
|
0.27+0.03-0.03 | 500.0+65.1-52.7 | 2.05 |
|
1.15(485) |
| a We do not show fit parameters for observations without enough counts or without a reliable |
| energy calibration. We assumed a metallicity of 0.4 times the solar value, i.e., Z = 0.008. |
|
b The Galactic line-of-sight absorption (
|
| c Fixed. |
The reduced
are about one in most fits.
The large
in some fits
(e.g. the MOS1 observations in Rev0092/201)
are probably due
to the uncertain response and effective-area calibrations.
In particular, the feature seen at
0.5 keV
in the PN spectra (Figs. 2-4) is due to uncertainties
in the charge transfer efficiency correction for the small window mode.
Values of
are obtained
when the 0.4-0.6 keV energy range is excluded from the fit.
We consider the RGS data
to search for possible emission lines,
and to constrain the value of the column density.
The RGS spectra are fitted
with the same model used in the analysis of the PN data.
As the RGS spectra cover only the 0.35-2.5 keV energy range,
the thermal and the power-law components cannot be constrained
simultaneously.
We therefore fix
to the values determined from the EPIC fits
(with a correction to the value of
from Rev0045
to take the pile-up into account).
Because of the low count rate,
the data in the Rev0066 observations
are binned in groups of 100 channels each.
The error-weighted means of the fit values of
are 0.52 keV and 0.29 keV
for observations before and after Rev0066 respectively.
We do not find strong evidence of emission lines.
The column density
was generally below 1021 cm-2,
except for Rev0066,
where
is not well constrained.
After the line-of-sight column density to the LMC is subtracted,
the value of the column density for the observations
before Rev0066 is
cm-2,
and the value after Rev0066 is
cm-2.
(These values are smaller than those
obtained from the fits to the EPIC data.)
Given that the error may be dominated by systematic effects,
the error obtained in the fits may not truly represent
the statistical weights of the fit parameters.
Apart from the conclusion that
the spectra do not show strong column absorption,
we are unable to determine whether the variations in
are correlated with the X-ray luminosity.
Previous RXTE/PCA observations (Wilms et al. 2000)
have shown that
when the source was bright
(with RXTE/ASM count rates
ct s-1),
the disk blackbody component was prominent, and
its fit temperature
keV.
There was also a power-law tail with
in the spectra.
The temperature
appeared to decrease
when the X-ray luminosity decreased,
while the normalisation parameter
remained approximately constant.
When the RXTE/ASM count rate dropped below
0.6 ct s-1,
was reduced to
0.7 keV.
The transition was accompanied by the hardening of the power-law component.
The photon index
became
2.0-3.0.
When the RXTE/ASM count rate decreased
below the 0.3-ct s-1 level,
the disk blackbody component was not detected,
and the spectrum was a power law with
.
![]() |
Figure 2: The top panel shows the data and fit EPIC-PN spectra of the Rev0028/201 observation (2000 February 2). An absorbed diskbb+powerlaw model is used in the fit (see Table 3). The residuals are shown in the bottom panel |
| Open with DEXTER | |
![]() |
Figure 3: Same as Fig. 2 for the Rev0066/101 EPIC-PN observation (2000 April 19) |
| Open with DEXTER | |
![]() |
Figure 4: Same as Fig. 2 for the Rev0092/201 EPIC-PN observation (2000 June 10) |
| Open with DEXTER | |
The soft-to-hard transition was also seen in our data.
The spectra obtained in the Rev0028-Rev0045 observations
show a disk blackbody component with
-1.0,
and a power-law component with
.
The RXTE/ASM count rate was below 2 ct s-1
during the observations.
The spectral properties are similar to
those observed previously
when the system had similar RXTE/ASM count rates.
The spectrum obtained in the Rev0066/101 observations
is dominated by a power law with
.
The photon index is consistent with that observed
in the previous hard states
(Wilms et al. 2000), while a lower value of
was obtained by Boyd et al. (2000) on 2000 May 7.
The
of the fit which
includes a thermal disk black-body component
(best fit
keV)
is 58.486 for 52 degrees of freedom.
If the thermal component is not included,
then we obtain
for 54 degrees of freedom.
The thermal component is therefore insignificant.
The system was at the transition
from the hard state back to the soft state
during the Rev0092 observations (see Fig. 1).
The power-law component had steepened,
with
(for the PN data).
The disk blackbody component reappeared, with a temperature
keV, significantly lower than
that before Rev0066.
In summary, between 2000 February and June LMC
X-3 underwent a transition
from a soft to a hard state,
and then in the process of returning to the soft state.
During the soft-to-hard transition
the fit temperature
of the disk blackbody component decreased,
and the power law component became harder.
As the system started to return to the soft state,
the disk blackbody component became more prominent.
The power law appeared to be steeper than that
obtained from the RXTE/PCA observations near the middle of the
faint state, yet the disk blackbody temperature was still well below
the values of the previous soft state.
A similar conclusion is also obtained independently
from an analysis of the optical/UV properties of the system
(Soria et al. 2001).
The OM data obtained in the Rev0066
suggest that the companion is a B5 subgiant
instead of a B3 main-sequence star.
Roughly 3% of the X-rays would be intercepted
by the companion, so that
the rate at which energy is deposited into its atmosphere
can be
erg s-1 in the soft state.
This rate is larger than
the intrinsic luminosity of the companion.
If the companion in 3 is indeed a subgiant star
its tenuous envelope is susceptible to irradiation heating.
The soft-to-hard transitions
seen in the RXTE observations in 1997/1998
and the XMM-Newton observations in 2000 may be caused by
variations in the rate of mass overflow
from the Roche lobe of the subgiant companion.
Wilms et al. (2000)'s interpretation of the decreases in the RXTE/ASM count rate as evidence for transitions from soft to hard state is consistent with our data. We further propose that the decrease in the X-ray luminosity is caused by the decrease in the fraction of the Roche lobe filled by the companion star. When the companion is detached from its critical Roche surface, mass transfer will be dominated by a focused wind.
It is worth noting that the three known high-mass BHCs
are all persistent X-ray sources
which show preferential X-ray spectral states.
While Cyg X-1 tends to be in the hard state,
LMC X-1 and 3 are more often found in the soft state.
Recent studies
(e.g. Igumenshchev et al. 1999;
Beloborodov & Illarionov 2000)
have shown that
accretion of matter with low angular momentum
will give rise to hard X-rays instead of soft X-rays.
The relative angular momentum of the accreting matter
is smaller for wind accretion than for Roche-lobe overflow.
Cyg X-1 has a 33-
O-type companion
(Giles & Bolton 1986), which has a strong stellar wind.
The companion stars in LMC X-1 and 3
are less massive (
)
B stars,
whose stellar wind is much weaker.
Therefore, we suggest a unified scenario
which relates the mode of mass transfer (wind or Roche-lobe overflow)
to the spectral state preferentially observed
in these three high-mass BHC binaries (hard or soft respectively).
The subgiant companion of 3
may occasionally underfill its Roche lobe
because of feedback irradiative processes
or instabilities in its envelope.
This leads to the residual accretion of the (focused) wind matter
which has relatively low angular momentum.
The BHC 3 was observed in 2000 February-June,
with the XMM-Newton EPIC and RGS,
throughout a soft-hard transition.
The system was apparently in the process
of returning to the soft state in 2000 June.
The hard-state spectra are dominated by a power-law component
with a photon index
.
The soft-state spectra consist of a thermal component
with an inner-disk temperature
of
0.9 keV
and a power-law tail with
-2.7.
The line-of-sight absorption deduced from the EPIC and RGS data
is
cm-2.
Our observations therefore do not support the wind accretion
model for this system in the soft state.
The transition from the soft to the hard state
appears to be a smooth process
associated with the changes in the mass-transfer rate.
Acknowledgements
This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA member states and the USA (NASA). We thank Keith Mason for his comments. KW acknowledges a PPARC visiting fellowship.