A&A 386, 910-915 (2002)
DOI: 10.1051/0004-6361:20020281
A. N. Parmar - T. Oosterbroek - L. Boirin - D. Lumb
Astrophysics Division, Research and Scientific Support Department of ESA, ESTEC, Postbus 299, 2200 AG Noordwijk, The Netherlands
Received 12 November 2001 / Accepted 22 February 2002
Abstract
We report the discovery of narrow X-ray absorption
features from the dipping low-mass X-ray binary X1624-490 during an
XMM-Newton observation in 2001 February. The features are identified
with the K absorption lines of Fe XXV
and Fe XXVI and have
energies of
keV and
keV
and equivalent widths (EWs) of
-7.5+1.7 -6.3 eV and
-16.6+1.9 -5.9 eV, respectively.
The EWs show no obvious dependence on orbital phase, except
during a dip,
and correspond to a column of
1017.3 Fe atom cm-2.
In addition,
faint absorption features tentatively identified with
Ni XXVII K
and Fe XXVI K
may be present.
A broad emission feature at
6.58+0.07 -0.04 keV with an
EW of
78+19 -6 eV
is also evident. This is probably
the 6.4 keV feature reported by earlier missions since fitting
a single Gaussian to the entire Fe-K region gives an energy
of
6.39+0.03 -0.04 keV.
A deep absorption feature is present during the dip with an energy
consistent with Fe XXV K
.
This is the second dipping LMXRB source from
which narrow Fe absorption features have been observed. Until recently
the only X-ray binaries known to exhibit narrow X-ray absorption lines
were two superluminal jet sources and it had been suggested that these
features are related to the jet formation mechanism.
It now appears likely that ionized absorption features
may be common characteristics of accreting systems with accretion disks.
Key words: accretion, accretion disks - stars: individual: X1624-490 - stars: neutron - X-rays: general
The low-mass X-ray binary (LMXRB) source X1624-490 exhibits irregular dips in
X-ray intensity that repeat every orbital cycle of
21 hr (Jones & Watson 1985).
This long orbital period means that as well as being the most
luminous dip source, X1624-490 also has the largest stellar
separation and is often referred to as the "Big Dipper''.
It is generally accepted that these dips are due to obscuration in
the thickened outer regions of an accretion disk (White & Swank 1982).
The depth, duration and spectral properties of the dips
vary from source to source and from cycle to cycle (see
Parmar & White 1988 and White et al. 1995 for reviews).
The X1624-490 X-ray lightcurve exhibits prominent flares, but no X-ray bursts
have been observed (e.g., Baucinska-Church et al. 2001;
Smale et al. 2001).
The X-ray continuum of X1624-490 may be modeled using
blackbody and Comptonized components which are
absorbed by a column of
atom cm-2(Ba
ucinska-Church et al. 2000). There is evidence for
the presence of a dust scattering halo (Angelini et al. 1996).
During dips the point-like blackbody is rapidly obscured and the
extended Comptonized component, presumably an
accretion disk corona, is progressively covered by an extended absorber
(Ba
ucinska-Church et al. 2000).
XMM-Newton is proving to be very successful in improving our
understanding of the nature of emitting and absorbing regions in
LMXRB. Reflection Grating Spectrometer
(RGS) observations of the eclipsing and dipping system EXO0748-676,
revealed the presence of a rich variety
of emission lines identified with O VIII Ly,
O VII He-like complex, Ne X Ly
,
Ne IX He-like complex
and N VII Ly
(Cottam et al. 2001).
In addition, photo-electric absorption
edges of both O VIII and O VII were detected
as well as narrow radiative recombination continua of O VIII
and O VII near their respective absorption edges.
The line features are significantly broadened
with velocity widths of 1000-3000 km s-1, depending on the
line, and no obvious velocity shift (<300 km s-1).
There is no obvious dependence of the line
properties on orbital phase (even during X-ray eclipse).
European Photon Imaging Camera
(EPIC) and RGS spectra of the eclipsing and dipping system
MXB1658-298 have revealed the presence of narrow
resonant absorption features identified with
O VIII K
,
K
,
and K
,
Ne X K
,
Fe XXV K
,
and Fe XXVI K
,
together with a broad Fe emission feature at
6.47+0.18 -0.14 keV
(Sidoli et al. 2001).
Again, the properties of these features show no obvious
dependence on orbital phase, even during dipping intervals.
Until recently, the only X-ray binaries known to exhibit such narrow
X-ray absorption
lines were two superluminal jet sources and it had been suggested that
these features are related to the jet formation mechanism. This now
appears unlikely, and instead their presence may be
related to the viewing angle of the system, or they may be a common
feature of LMXRB systems as suggested by the recent detection of a
narrow (
eV) absorption feature with an energy of
keV from
GX13+1 (Ueda et al. 2001). Here we report the discovery
of narrow X-ray absorption features from highly
ionized Fe in the XMM-Newton spectrum of X1624-490.
![]() |
Figure 1: The EPIC PN 2-10 keV lightcurve of X1624-490 with a binning of 64 s. The upper panel shows the hardness ratio (counts between 4-10 keV divided by those between 2-4 keV) with a binning of 256 s. Time is hours of 2001 February 11. |
Open with DEXTER |
The XMM-Newton Observatory (Jansen et al. 2001) includes three 1500 cm2 X-ray telescopes each with an EPIC at the focus. Two of the EPIC imaging spectrometers use MOS CCDs (Turner et al. 2001) and one uses a PN CCD (Strüder et al. 2001). The region of sky containing X1624-490 was observed by XMM-Newton on 2001 February 11 23:23 UT to February 12 23:06 UT as part of the Performance Verification Programme. In order to minimize the effects of pile-up the PN and one MOS were operated in their Small Window modes. The other MOS CCD was operated in its normal Full Frame mode and the data are strongly affected by pile-up and not discussed further. Due to the large amount of absorption there are too few counts in the 0.35-2.5 keV RGS to perform a useful analysis. The OM was not operated during the observation. Since useful data on X1624-490 was only obtained with one MOS, we concentrate on the analysis of PN data, and use the lower count rate MOS data to check for consistency.
Raw data products were extracted from
the public XMM-Newton archive and then reprocessed using
version 5.2 of the Science Analysis Software (SAS),
before being further filtered using XMMSELECT. For the MOS
only X-ray events corresponding to patterns 0-12 were selected.
Known hot, or flickering, pixels and electronic noise were rejected
using the SAS.
Source events were extracted from circular regions of 60
radius centered on X1624-490.
Background spectra were obtained from a total of
60 ks
observations of several fields in the normal Full Frame modes.
Background subtraction, however, is not critical for such a
bright source.
All quoted count rates and fluxes are corrected for
instrumental deadtime.
![]() |
Figure 2: The PN spectrum of the X1624-490 persistent emission and the best-fit blackbody and power-law continuum model together with a broad emission line. The lower panel shows the ratio of model to data counts when the normalizations of the line features are set to zero. |
Open with DEXTER |
![]() |
Figure 3:
The 5.5-8.5 keV residuals from the fit shown in
Fig. 2.
A broad emission feature at 6.58 keV
and 2 narrow features identified with Fe XXV K![]() ![]() ![]() ![]() |
Open with DEXTER |
The 2-10 keV X1624-490
lightcurve obtained from the PN during the 2001 February
observation is shown in Fig. 1 with a binning of
64 s. The upper panel
shows the PN hardness ratio (counts in the energy range
4-10 keV divided by those between 2-4 keV) with a binning of
256 s. A dip started after 2001 February 12 05:47 UT
and had a duration of <3 hours.
During the dip the source exhibited rapid, irregular changes in
count rate and hardness ratio.
Such behavior
is typical of X1624-490 (see e.g., Baucinska-Church et al. 2000;
Smale et al. 2001).
During the persistent emission some flaring is evident, as investigated in detail
by Smale et al. (2001) and
Ba
ucinska-Church et al. (2001) who show that it
is primarily a high-energy effect, not visible in RXTE Proportional Counter
Array 2-6 keV data, but becoming
dominant at energies
10 keV.
The persistent emission spectrum of X1624-490 was first investigated by
excluding source events between 2001 February 12 05:47 and 09:06 UT
during the dip. This selection results in a
2.0-10 keV PN count rate of 89 s-1 and an exposure of
52 ks. This count rate is close to the 100 s-1 level,
above which pile-up effects become important (Strüder et al. 2001).
At high count rates the energies of multiple events are added and
assigned to a single event. This results in harder spectral
continua than in the unaffected case, but should not produce additional
line features.
In the PN Small Window mode only a 63
63 pixel region of the
central CCD is read out with a time resolution of 5.7 ms.
In order to investigate the effects of pile-up a second
spectrum was accumulated with the region contained
within the central 11
radius of the
point spread function (PSF) excluded. This results in the
loss of 45% of the counts.
Both spectra were
rebinned to oversample the full width half-maximum of the energy
resolution by a factor 3 and to have additionally a minimum of 20 counts per bin to allow use of the
statistic.
In order to account for systematic effects a 2%
uncertainty was added quadratically to each spectral bin.
The photo-electric absorption cross sections of Morrison & McCammon
(1983) are used throughout.
All spectral uncertainties are given at 90%
confidence, unless indicated otherwise.
The overall continuum was modeled using absorbed blackbody and
power-law components.
A power-law was chosen rather
that a cutoff power-law or Comptonization model since the cutoff
energy of 12 keV (Ba
ucinska-Church et al. 2000)
is too high to be important in the EPIC energy range.
This combination models the overall shape of the continuum moderately
well and gives
s of 517.5 and 483.6 for 173
degrees of freedom (d.o.f.) for the total and PSF excluded
spectra, respectively.
A comparison of the best-fit parameters
shows that, as expected, the total spectrum is significantly harder than
the PSF excluded one. We therefore use the PSF excluded spectrum
to determine the continuum parameters.
Two broad spectral features remain in both spectra centered on 2.5 and 3.5 keV.
The first feature is almost certainly
due to an incorrect modeling of the Au mirror edges,
while the second feature, with a peak-to-peak amplitude of 5% is likely also
due to remaining PN calibration uncertainties and
is not visible in the MOS spectrum.
Examination of the remaining fit
residuals in both spectra
shows that broad positive residuals are present near
6.6 keV, as well as two narrow negative residuals at slightly higher energies.
The properties of these features are consistent between the
two spectral fits. Since pile-up is not expected to significantly affect
narrow spectral features, we use the results obtained from the fit
to the total spectrum to investigate the properties of any discrete
spectral features.
The addition of a
keV
broad emission line at
6.39+0.03 -0.04 keV with
an equivalent width (EW) of
33+6 -5 eV reduced the
to
299.7 for 170 d.o.f. When modeled in this way, this feature is
consistent with that seen by
previous missions such as ASCA which measured an energy of
6.4+0.1 -0.3 keV and an EW of
13+13 -9 eV
(Asai et al. 2000).
The F-statistic value of 41.2 indicates that
the probability of such a decrease occurring by chance is
.
The narrow features were modeled as Gaussian absorption lines.
Including the first feature at
keV with an EW
of
-16.6+1.9 -5.9 eV reduced the
to 249.0 for 167 d.o.f.
This decrease in
corresponds to an F-statistic of 11.3. The probability of such
a decrease occurring by chance is
.
Including a second narrow feature at
keV with an EW of
-7.5+1.7 -6.3 eV
reduced the
to 218.9 for 164 d.o.f. This decrease in
corresponds to an F-statistic of 7.5. The probability of such
a decrease occurring by chance is 10-4.
When the two absorption lines are included in the fit, the properties of the
emission
feature change to a best-fit energy of
6.58+0.07 -0.04 keV, a width (
)
of
eV and an EW of
78 +19 -6 eV.
Figure 2 shows the best-fit spectral fit
model for the X1624-490 total spectrum. The best-fit
parameters derived from the fit to the PSF core excluded spectrum
consist of a
1.18+0.01 -0.02 keV
blackbody with a radius of
9.60+0.04 -0.10 km
(for a distance of
15 kpc), together with a power-law with a photon index,
,
of
.
Both components are modified by absorption equivalent
to
atom cm-2.
The best-fit continuum and line model parameters are given in
Table 1 and the residuals in the region around
the Fe features are shown in Fig. 3.
Figure 3 also shows evidence for the presence
of two additional narrow absorption features
at
7.83+0.03 -0.05 and
8.28+0.07 -0.04 keV
with EWs of
-1.4+1.4 -1.5 and
-1.7+1.1 -2.2 eV,
respectively (at 68% confidence).
The measured energies are consistent with those
of Ni XXVII K
at 7.806 keV and Fe XXVI K
at 8.26 keV, both features
seen in the ASCA spectrum of GRS 1915+106 by Kotani et al. (2000).
The reductions in
when these features are added to the best-fit
spectral model are however not significant.
We caution also that
this is a spectral region where the EPIC calibration is still relatively
uncertain, and so we regard the presence, and proposed identifications,
of the features as tentative.
We have searched for absorption edges in the PN spectrum
due to neutral Fe, Fe XXV, and Fe XXVI expected
at 7.12, 8.83 and 9.28 keV.
The 95% confidence upper-limit to the optical depth of an edge
at 7.12 keV is 0.009 and there is some evidence for the presence
of edges at 8.83 and 9.28 keV with optical depths of
and
0.011+0.013 -0.011.
Component | Parameter | Value |
![]() |
![]() |
|
L (erg s-1) | 3.8 1037 | |
![]() |
195.6/164 | |
Blackbody | kT (keV) | 1.18+0.01 -0.02 |
Radius (km) | 9.60+0.04 -0.10 | |
Power-law | ![]() |
![]() |
Normalization at 1 keV |
![]() |
|
Fe emission |
![]() |
6.58+0.07 -0.04 |
feature | ![]() |
![]() |
![]() |
11.7+2.9 -0.9 | |
EW (eV) | 78+19 -6 | |
Fe XXV abs |
![]() |
![]() |
feature | ![]() |
<50 |
![]() |
-1.2+0.3 -1.0 | |
EW (eV) | -7.5+1.7 -6.3 | |
Fe XXVI abs |
![]() |
![]() |
feature | ![]() |
<56 |
![]() |
-2.3+0.3 -0.8 | |
EW (eV) | -16.6+1.9 -5.9 |
![]() |
Figure 4: Residuals in the 5.5-8.5 keV energy range when a power-law and blackbody continuum is fit to the shallow (upper panel) and deep (lower panel) dip spectra. The complex nature of the remaining structure is evident. |
Open with DEXTER |
In order to investigate whether the properties of the absorption features
depend on orbital phase, the persistent emission interval was divided
into three approximately equal duration parts, one before the dip, and
two after (see Fig. 1). The approximate phases covered
are 0.57-0.88, 0.03-0.37 and 0.37-0.70 using the ephemeris of
Smale et al. (2001), which has an uncertainty of 0.03 at the
time of the XMM-Newton observation (phase 0.0 corresponds to the dip
centers).
The same continuum model as used in Sect. 3.2
was fit to the 3 PN spectra and the
EWs of the two Fe absorption features determined. For the deeper
7.00 keV feature, EWs of
-23.4+5.2 -7.0,
-24.1+5.3 -6.4, and
-15.7+4.6 -8.5 eV were
obtained. Thus, there does not appear to be any obvious dependence of the
EW of the Fe XXVI feature on orbital phase, at least outside of
dipping intervals. A similar analysis, when applied to the shallower
Fe XXV feature is not so constraining, but gives a similar result.
The properties of the absorption features during the dip were
next investigating by accumulating PN spectra during intervals of
deep and shallow dipping activity. Events were selected by
making a lightcurve with 16 s resolution of the dipping interval. Time
intervals corresponding to PN count rates above and below 37 s-1were then derived.
Deep and shallow dip spectra were accumulated using
these time intervals.
The 2-10 keV count rates are 23 and 61 s-1 for exposures of
1047 and 480 s.
In order to characterize the continua, the same blackbody and
power-law model as in Sect. 3.2 with all
spectral parameters allowed to vary was used to give
values of 198.3 and 172.2 for 167 and 164 d.o.f., for the
shallow and deep dip spectra, respectively.
It is well known that detailed modeling of the spectral
changes during dips requires multiple components
which undergo different amounts of absorption and covering (e.g.,
Church & Ba
ucinska-Church 1995). However, the simple
approach adopted here allows us to investigate the properties of
any narrow absorption features present.
Figure 4 shows the residual in the energy range
5.5-8.5 KeV obtained in this way for the shallow and deep dipping
intervals. There is a great deal of spectral structure
remaining. The shallow dip spectrum
shows evidence for two deep features which may be modeled
as absorption lines with energies of
and
keV
(Table 2).
The lower energy feature may be Fe XXV K
resonant
absorption, while the energy of the second feature does not correspond
to that of any likely line or edge feature.
The same two features are also
present in the deep dip spectrum, together with an
emission feature at
6.47+0.09 -0.15 keV. This may be the
same Fe emission feature as seen in the persistent emission spectrum.
The deep Fe XXVI K
feature at 7.00 keV
seen in the persistent emission spectrum may also be present in both
dip spectra (see Fig. 4),
but is much shallower than the feature
identified with Fe XXV K
.
This suggests that during
dipping the additional absorbing material is less strongly ionized
than that responsible for the absorption features seen in the persistent
emission. A detailed
study of such features may allow the ionization state, velocity and
temperatures of the material responsible for the dips to be probed.
Shallow dipping | |||
E (keV) | ... |
![]() |
![]() |
EW (eV) | ... |
![]() |
-151+44 -35 |
![]() |
... | ![]() |
![]() |
Deep dipping | |||
E (keV) | 6.47+0.09 -0.15 |
![]() |
7.74+0.06 -0.14 |
EW (eV) |
![]() |
-78+43 -70 | -120+40 -150 |
![]() |
370+390 -110 | ![]() |
140+150 -70 |
![]() |
Figure 5: Fe column densities derived from the curve of growth analysis for the Fe XXV and Fe XXVI absorption lines observed in the persistent emission. The solid lines show the values derived using the best-fit values of the EWs and the shaded areas show the 90% confidence regions. The dotted line indicates the Fe column density equivalent to one Thomson depth. |
Open with DEXTER |
Narrow X-ray absorption lines were first detected from the superluminal
jet sources GRO J1655-40 (Ueda et al. 1998;
Yamaoka et al. 2001) and GRS 1915+105 (Kotani et al. 2000;
Lee et al. 2002). ASCA observations of GRO J1655-40 revealed the
presence of absorption features due to Fe XXV and Fe XXVI which
did not show any obvious dependence of their EWs on orbital phase.
GRO J1655-40 has been observed to undergo deep absorption dips
(Kuulkers et al. 1998) consistent with observing the source
at a high inclination angle of -
(e.g., Frank et al. 1987). ASCA observations of GRS 1915+105
revealed, in addition, absorption features due to Ca XX,
Ni XXVII and Ni XXVIII. A recent Chandra HETGS observation
of this source revealed absorption edges of Fe, Si, Mg, and S
as well as resonant absorption features from
Fe XXV and Fe XXVI and possibly Ca XX
(Lee et al. 2002).
Until recently, it was possible that these absorption features were
peculiar to superluminal jet sources and related in some way to the jet
formation mechanism. With the discovery of narrow absorption features from
the LMXRBs GX13+1 (Ueda et al. 2001),
MXB1658-298 (Sidoli et al. 2001) and now
X1624-490, this appears not to be the case, and as proposed
by Kotani et al. (2000) ionized absorption features
may be common characteristics of systems accreting via a disk.
In order to estimate the column densities of Fe that produced
the absorption lines observed in the persistent emission
we performed a curve of growth analysis
following that of Kotani et al. (2000).
This allows the EWs of the lines to be converted into column densities
depending on the kinematic temperature of the absorbing material.
The kinematic temperature includes contributions from thermal
motions as well as any bulk motions.
An Fe abundance of
was assumed
(Anders & Grevesse 1989).
Figure 5 shows the calculated columns for
the Fe XXV and Fe XXVI features
for a range of assumed temperatures.
As expected, a lower kinematic temperature requires a higher ion
column density, and hence a higher hydrogen column density.
With the range of assumed temperatures, the derived Fe column densities
are always less than the
corresponding
to one Thomson optical depth (
atom cm-2). Such a high
column (shown as a horizontal line in Fig. 5)
would result in the absorption lines being
strongly diminished. Examination of Fig. 5
shows that the total Fe column density must be
1017.3 cm-2 for the range of assumed temperatures.
The total line EW observed from X1624-490 is a factor 3
lower than that observed
from the eclipsing, dipping source MXB1658-298 by
Sidoli et al. (2001) also using XMM-Newton. In addition, the
ratio of H-like to He-like line EWs of
2.2 in X1624-490 is higher than in
MXB1658-298 where a ratio of
1.3 was determined.
This suggests that there is relatively more Fe XXVI than Fe XXV
present in the absorber in X1624-490, compared to MXB1658-29.
There could be a number of reasons for these differences. (1)
The overall metal abundance could be higher in MXB1658-298
than in X1624-490, resulting in deeper absorption features. (2) The
geometries of the absorbing regions could be different, resulting
in different fractions of scattered photons entering the line of
sight. (3) As implied by the larger ratio of Fe XXVI to Fe XXV EW in X1624-490 compared to MXB1658-298 the ionization parameter
of the absorbing system may be higher in this system and much of the Fe
in the absorber may be fully ionized. (4) Since MXB1658-298 shows
dips and eclipses, while X1624-490 only shows dips, it is possible
that MXB1658-298 is viewed closer to the orbital plane than
X1624-490. This could mean that the line of sight may pass through more
material. It will be very interesting to obtain column
densities of absorption features in other LMXRB (dip-)sources to
investigate which of the above possibilities is true.
Asai et al. (2000) report on a spectral survey of Fe-K emission
lines from 20 LMXRB with ASCA. Significant lines are detected from about
half the sources. The mean line energy is
keV with a full width half-maximum
of
0.5 keV. Neither the line energy, EW nor width appear
to depend on source luminosity or type. When the X1624-490 persistent
emission PN spectrum is fit with a single Gaussian emission feature, the
best-fit energy of
6.39+0.03 -0.04 keV is significantly below
this mean value. However, when the two absorption features are included,
then the best-fit energy of
6.58+0.07 -0.04 keV is fully
consistent with the mean energy found by ASCA for the LMXRB sample.
This difference well illustrates the dependence of derived spectral
results on the chosen model. It also suggests that other LMXRB with
Fe-line energies well below the average in the ASCA sample may exhibit
absorption features. Examination of Table 5 of Asai et al. (2000)
shows that the line observed from XB 1916-053 has the lowest energy
in the sample of
5.9+0.2 -0.1 keV. Interestingly, XB 1916-053 is
another dip source (White & Swank 1982), supporting the idea that
the presence of absorption features may depend on viewing angle.
Acknowledgements
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). L. Boirin acknowledges an ESA Fellowship. We thank T. Kotani for making his curve of growth software available.