A&A 488, 249-255 (2008)
DOI: 10.1051/0004-6361:200809506
L. Sidoli1 - N. La Palombara1 - T. Oosterbroek2 - A. N. Parmar3
1 - INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica, via E. Bassini 15, 20133 Milano, Italy
2 - Science Payload and Advanced Concepts Office, ESA, ESTEC, Keplerlaan 1, 2200 AG Noordwijk, The Netherlands
3 - Research and Scientific Support Department of ESA, ESTEC,
Keplerlaan 1, 2200 AG Noordwijk, The Netherlands
Received 4 February 2008 / Accepted 25 May 2008
Abstract
Context. We report on two XMM-Newton observations performed in 2006 of the luminous low-mass X-ray binary XB1832-330, which is located in the galactic globular cluster NGC6652 and is probably an ultracompact binary (orbital period,
,
of <1 h).
Aims. The aim of these observations is to investigate the low-energy absorption towards XB1832-330 and in particular to search for the Ne-rich material local to the binary, which has been suggested as a possible spectral signature of Neon-rich, degenerate companions.
Methods. XMM-Newton observed the source twice, in 2006 September and October. High resolution spectroscopy with the XMM-Newton RGS was used to estimate the ratio of the neutral neon-to-oxygen column densities to search for an anomalous Ne abundance in this X-ray binary.
Results. We find no evidence of anomalous Ne/O ratios, finding Ne/O =
and Ne/O =
(1
uncertainties) in the two observations. These values are consistent with the one in the interstellar medium. Timing analysis of EPIC data suggests possible periodicities at
s and
s, which need confirmation. A by-product of these observations consists of the spatial analysis of the source field, which resulted in the detection of 46 faint sources within the EPIC field of view, which are not present in the Second XMM-Newton Serendipitous Source Catalogue. All these faint sources are probably foreground objects.
Conclusions. We performed the first high spectral resolution observations of XB1832-330, a probable ultracompact binary, without finding any evidence of an anomalous Ne abundance.
Key words: X-rays: individuals: XB 1832-330 - stars: neutron - X-rays: binaries - accretion, accretion disks
XB1832-330 was discovered using HEAO-1 (Hertz & Wood 1985) and is one of the bright low-mass X-ray binaries (LMXBs) located in galactic globular clusters. The source was observed with ROSAT (Predehl et al. 1991; Johnston et al. 1996) and associated with the galactic globular cluster NGC6652. The ROSAT uncertainty region was observed in 2000 May with the High Resolution Camera-Imager (HRC-I) camera on the Chandra X-Ray Observatory (Heinke et al. 2001). Four X-ray sources were detected within 30'' of the cluster centre, including the bright LMXB XB1832-330 (named ``source A'' or CXOGLB J183543.6-325926 in Heinke et al. 2001). The detection of a type I X-ray burst confirmed the neutron star nature of the compact object (in 't Zand et al. 1998), although one of the other fainter X-ray sources detected with Chandra could also be responsible for the burst (Heinke et al. 2001).
The significantly smaller Chandra uncertainty region allowed a secure identification of
the optical counterpart (Heinke et al. 2001)
with an
M V = +3.7 object (HST observations).
A clear 0.08 mag variability in the V and I lightcurves was observed (Heinke et al. 2001), with possible periods of 3300 s,
8000 s (or
16 000 s for ellipsoidal variations), as well as non-periodic flickering. The short temporal coverature did not allow
the three periodicities to be clearly separated.
Heinke et al. also notes that the optical magnitude is compatible with
the donor star being in an ultracompact binary (orbital period,
,
less than about 1 h). If one of the three periodicities is associated with the orbital period,
the secondary could be a white dwarf or a low-mass M dwarf (Heinke et al. 2001).
The X-ray spectroscopy of XB1832-330 has been performed with ASCA (Mukai & Smale 2000) and BeppoSAX (Parmar et al. 2001) and indicates that no simple model provides acceptable fits to the spectrum, suggesting that partial covering absorption, or other spectral complexity plays an important role in modifying the X-ray emission at low energies. Another more recent observation at high energies (above 20 keV) has been performed with IBIS on-board INTEGRAL (Tarana et al. 2006).
![]() |
Figure 1: PN XB1832-330 lightcurves of the two observations in the energy ranges 0.3-2 keV and 2-12 keV. The bottom panels show the hardness ratios (counts between 2-12 keV divided by those between 0.3-2 keV) in the two observations. The binning is 1024 s. |
Open with DEXTER |
The galactic globular cluster NGC6652 is located at a distance of 9.3 kpc.
It has a metallicity of [Fe/H] = -0.9 and a low reddening (
)
(Ortolani et al. 1994), which translates into an absorption,
,
of
cm-2(Predehl & Schmitt 1995).
We report here on the results of two XMM-Newton observations performed in order to investigate the low-energy absorption towards XB1832-330, in particular to search for Ne-rich material local to the binary, which has been suggested as a possible spectral signature of Neon-rich, degenerate companions in binaries (Juett et al. 2001).
XB1832-330 was observed with XMM-Newton on 2006 September 22 (hereafter, obs I)
and on 2006 October 20 (obs II), for a net exposure of
28 ks and
34 ks, respectively (obs-IDs 0400790101 and 0400790301).
Data were reprocessed using version 7.1.0 of the Science Analysis
Software (SAS). Known hot, or flickering, pixels and electronic
noise were rejected. Response and ancillary matrix files
were generated using the SAS tasks rmfgen and arfgen.
The EPIC pn was operated in its Small Window mode, the MOS1 camera was in its Prime Partial Window mode, with the central CCD operated in Timing Mode and with the other CCDs in Imaging Mode. The
MOS2 camera was operated in its Prime Partial Window with all the CCDs in Imaging Mode.
The EPIC MOS and pn observations used the medium thickness filter.
Spectra were extracted using PATTERN 0 for the MOS cameras
(only single events) and PATTERNS from 0 to 4 (single and double events) in the pn.
The RGS source and background events were extracted by making spatial and
energy selections on the event files.
Source counts were extracted from circular regions
of 40
radius for the pn and from a circular corona with an internal radius
of 8
and an external radius of 1' for MOS2. Source counts in
MOS1 timing mode were extracted from
to
,
excluding the central
core with an internal radius of 5 RAW pixels.
The MOS1 and MOS2 data clearly suffered from the effects of pile-up.
We checked iteratively with the SAS task epatplot that the events selected from
annular regions
in the MOS cameras were not significantly affected by pile-up after this spatial selection.
Background counts were obtained from similar sized regions offset from the source positions.
For the MOS1 timing mode, the background region has been taken from an outer CCD
that collected data in imaging mode.
The background did not show evidence of any significant flaring activity.
To ensure applicability of the
statistics, the
net spectra were rebinned such that at least 20 counts per
bin were present and such that the energy resolution was not
over sampled by more than a factor 3.
The XB1832-330 lightcurves and hardness ratios do not show evidence for any large spectral variability (Fig. 1). We therefore extracted averaged spectra from each of the two observations, which were then analysed separately, withouth any further filtering or selection.
Since the three EPIC instruments show significant discrepancies below 0.8 keV when they are compared, with opposite structured residuals if fitted with the same spectral model (Fig. 2), we conservatively decided to limit the spectral analysis to the range 0.8-12 keV, for both the MOS and the pn. Similar problems and discrepancies at low energies have often been observed in other X-ray binaries (e.g. Sidoli et al. 2005; Boirin et al. 2005). The softest energy range was covered with RGS1 and RGS2, which both extend down to 0.4 keV. For each observation separately, the two RGS (0.4-2 keV) spectra were fitted with the EPIC pn and the two MOS (0.8-12 keV) spectra. Free relative normalizations between the instruments were included. The net source count rates are reported in Table 1.
Since the main aim of these XMM-Newton observations was to investigate the low energy absorption towards XB1832-330 and, in particular, the estimate of the Ne/O ratio, we adopted a variable abundance absorption model ( VPHABS in XSPEC), with the elemental abundances set to the ISM values of Wilms et al. (2000), except for those of O, Ne, and Fe, which were set to zero. Their effect on the spectral shape has been replaced with three edges (O-K, Fe-L, Ne-K) with energies fixed at 0.543, 0.706, and 0.869 keV (e.g. Paerels et al. 2001) and edge depths that were allowed to vary. In this way we could also account for a local Fe abundance that could be different from the cosmic value.
We tried different models for the continua of the two observations. A single-component model was never able to fit the data adequately, as already demonstrated by previous observations with ASCA (Mukai & Smale 2000) and BeppoSAX (Parmar et al. 2001). Different combinations of simple models were tried: for the soft component we adopted a blackbody and a disk blackbody model ( DISKBB in XSPEC), while for the high energy component we used a power law, a cutoff power law and a Comptonization emission model ( COMPTT in XSPEC). We tried all possible combinations of these models, obtaining equally good fits and similar spectral parameters. Since there is no evidence in the XMM-Newton energy band of a cutoff (resulting in an unconstrained electron temperature in the COMPTT model or in a high-energy cutoff well above 10 keV), we adopted the power law plus a blackbody model for the continuum.
From the depth of the edges and the elemental cross sections, the element column densities
were measured, hence the Ne/O ratio, which has been compared to that of the
interstellar medium (ISM) value of 0.18, assuming the ISM abundances of Wilms et al. (2000).
The spectral results are reported in Table 2 (for the continuum) and in
Table 3 (for the edge depths,
,
columns,
,
and equivalent hydrogen columns,
,
together with the final Ne/O ratio, for the two observations).
![]() |
Figure 2: Ratio between the data and the model when fitting the MOS1, MOS2, and pn spectra of obs I with a double-component continuum (power law, together with a blackbody) absorbed by a partial covering fraction component, which accounts for the global spectral shape of the X-ray emission. The different (and ``opposite'') residuals in the pn (open circles), compared with the MOS1 (stars) and MOS2 (squares), are evident. A similar structured excess is present in the obs II spectrum. |
Open with DEXTER |
Table 1: Net count rates (in units of s-1) for XB1832-330 during the two observations.
Table 2: Spectroscopy of XB1832-330 with the same model adopted for the two observations (RGS1+RGS2+pn+MOS1+MOS2 data) consisting of a power-law model together with a blackbody, both equally absorbed by a variable absorption model, VPHABS in XSPEC.
Table 3:
Results on the photoelectric absorption towards XB1832-330 (see
Table 2 for the continuum emission parameters) when fitting RGS+pn+MOS spectra.
Quoted uncertainties are at 1
confidence.
We then analyzed the RGS spectra alone, to investigate the fine structure
around the absorption edges more closely. We adopted the same shape for the continua
as found in the broader band spectroscopy, with temperature of the blackbody
component and the power-law photon index fixed to the best-fits reported in Table 2,
but with all the normalizations free. The same method of measuring the edge depths of Ne, O, and Fe was applied to the RGS1 and RGS2 spectra (variable absorption model, together with three edges).
The best fit in both the observations (reduced
(1735 d.o.f.) and
(2070 d.o.f.)) is obtained adding a further absorption feature, a Gaussian line
at energy of
0.53 keV, which is compatible with absorption from interstellar
O I and O II (see e.g. Juett et al. 2004). In Figs. 3 and 4, we show both the net RGS spectra and the structure around the O edge. The absorption line at 0.53 keV has a normalization of
photons cm-2 s-1 (Obs I,
eV) and -
photons cm-2 s-1 (Obs II,
eV). The resulting depths for the absorption edges are the following:
,
,
in Obs I, and
,
,
in Obs II (all errors on the depth of the edges are at 1
).
The resulting ratios Ne/O are
0.39+0.07 -0.08 (Obs I, 1
uncertainty) and
0.27
+0.06 -0.05 (Obs II), which are within 2.6
and 1.8
,
respectively,
so consistent with the ISM value of 0.18. Thus we confirm with high-resolution RGS data alone,
which was found during the broader band analysis, with no evidence of an unusual ratio Ne/O.
![]() |
Figure 3: RGS1 and RGS2 XB1832-330 spectra in the two observations, fitted with the continuum model used to fit the broad-band spectra (see text). The bottom panels show the residuals in units of standard deviations. The missing regions in the RGS plots are due to the failed CCDs. |
Open with DEXTER |
![]() |
Figure 4:
RGS1 and RGS2 XB1832-330 spectra in the two observations, limited to the wavelenght region
around the Oxygen edge. In the left panel, the excesses just above 23 Å could be due to the O III emission line, but it is less than 2![]() |
Open with DEXTER |
From the Lomb-Scargle diagram of the obs I
there is evidence of a
periodicity near
s, with a significance of
about 0.6% (see Fig. 5, left panel).
Folding on this period reveals a smooth sinusodial modulation with an amplitude
of
% (Fig. 6).
This periodicity is not the
same as previously claimed values (Heinke et al. 2001), but periods in the 2-4 h range have been claimed before.
The period is rather long compared to the observation duration,
thus we cannot say much about the long-term stability.
However, it is quite likely (99.3%)
that during this observation an approximately sinusodial modulation with this period was present.
In obs II there is no evidence of the P=9170 s periodicity. However, there is evidence of modulation with a period of
s. From the Lomb-Scargle periodogram this period is significant at the 0.05% level (see Fig. 5, right panel). It is striking that this period is twice the period found in the first observation performed in 2006 September, within the uncertainties.
![]() |
Figure 5: Lomb-Scargle periodograms of the two XMM-Newton observations, 2006 September on the left and 2006 October on the right. |
Open with DEXTER |
![]() |
Figure 6: Lightcurves of the two EPIC observations folded on the periods P = 9170 s (Obs I, upper panel) and P = 18 616 s (Obs II, lower panel). Phase zero is arbitrary. |
Open with DEXTER |
The EPIC images of the XB1832-330 field show the presence of several faint X-ray sources. We therefore performed a detection analysis in order to obtain a source catalogue of this region. We considered only the two MOS cameras, since the pn was operated in Small Window mode and could image a region close to the central target.
The cleaned events files were used to produce MOS1 and MOS2 images in 5 energy ranges: 0.5-1 keV, 1-2 keV, 2-5 keV, 5-10 keV and 0.5-10 keV (total band). For each energy band, a corresponding set of exposure maps (i.e. one for each detector) was generated with the SAS task expmap, to account for spatial quantum efficiency (QE), mirror-vignetting, and field of view variations. To maximize the signal-to-noise ratio (S/N) of our serendipitous sources and to reach lower flux limits, we merged the data (images and maps) from the two MOS cameras for the two pointings separately. After the production of a detector mask (with the task emask), both images and exposure maps were used as a reference for the source detection, which was performed in four steps:
Finally, we used the srcmatch task to crossmatch the three sets of detected sources to check for coincident ones. This list contains all 46 sources whose distribution over the observed sky area is shown in Fig. 7. There we can check that sources #15, #16, #18, #31, and #46 are close to the edge of the EPIC FOV, so we cannot exclude their being spurious detections.
![]() |
Figure 7: Distribution of the 46 detected X-ray sources over the EPIC focal plane. The central CCD is empty since it was operated in Small Window mode and no source detection was performed on its data. Numbers mark the 46 sources listed in Table 4. |
Open with DEXTER |
The main characteristics of these sources are summarised in Table 4.
We list the source position (using the recommended XMM-Newton designations
for serendipitous sources), the position uncertainty, the count rate in the
four fine energy ranges and the two hardness ratios between
the two soft and the two hard ranges.
In the reported position error we not only consider the statistical errors,
associated with the centroiding of the point source,
but also the systematic error of 1'' (Kirsch 2007)
due to the absolute positional accuracy of XMM-Newton.
We also report the total source flux between 0.5 and 10 keV, which
has been calculated from the total count rate in the full energy band.
The applicable count rate-to-flux conversion factor has been calculated
assuming a power-law spectrum with photon index, ,
of 2 and the hydrogen column density in the direction of the globular cluster.
For each source we report whether it was detected in the data of the
single and/or the combined observations.
In this way we find that 11 sources of the first
observation are not detected in the second one,
while 8 sources of the second observation are not detected in the first one.
To verify that this result is due to a real source variability,
we analysed both the position and count rates of these sources.
We find that, in most cases, a source is detected in only one
observation because in the other it either falls out of the imaged
field-of-view (since the two observations are slightly rotated)
or its count rate is just below the detection limit, even if the count rates
measured in both the observations are consistent at a 1
confidence level.
Only in a few cases (highlighted with a ``V'' in Table 4)
we checked that the missing detection is due to a real source variability.
We find that only 14 sources detected in the merged observations
(over a total of 41) are also detected in both the individual observations, while 8 and 6 are only detected in obs I and II, respectively; and the remaining 13 sources are detected neither in the first nor the second observation. Finally, we also note that 3 sources of obs I and 2 sources of
the obs II are not detected in the merged observation.
We searched for possible known counterparts in the SIMBAD database, assuming an uncertainty radius equal to three times the source positional uncertainty. In this way we found no positionally coincident sources. The object known as NGC 6652 10, discovered by ROSAT (Johnston et al. 1996), is at a distance of 10.8'' from our source #4. In Fig. 7 we note that this source is rather bright and, moreover, that there are no other X-ray sources near its position. Therefore, we suggest that it is indeed NGC 6652 10, whose position is now better constrained by the XMM-Newton measurements.
With the same search criteria, we searched for possible counterparts in the
GSC2.3 (McLean et al. 2000), USNO-NOMAD (Monet et al. 2003), and 2MASS (Cutri et al. 2003) catalogues.
As shown in Table 4, with these constraints we found candidate
counterparts for only 10 sources. To evaluate whether they are due to
possible foreground contamination, we have to estimate the probability of a chance
coincidence between an X-ray and an optical source.
This parameter is given by
,
where r is the X-ray
uncertainty-circle radius and
is the surface density of the
optical sources (Severgnini et al. 2005).
In our case, the surface density of the GSC2.3, USNO-NOMAD and 2MASS sources in the
EPIC field-of-view is
,
,
and
arcsec-2,
respectively. Based on the estimated position errors,
we derived P > 10% for all the candidate counterparts;
therefore, we cannot exclude that they are due to random coincidences.
We report here on two XMM-Newton observations of XB1832-330 (located in NGC6652) performed in 2006, about 1 month apart, which are the first high-resolution spectra of this globular cluster source. XB1832-330 is one of the luminous LMXBs ( LX>1036 erg s-1) located in a galactic globular cluster, and it is supposed to be an ultracompact binary, because of the faintness of its optical counterpart (Heinke et al. 2001). Ultracompactness is also suggested by one of the possible periodicities found in optical data (Heinke et al. 2001), if associated with the orbital period.
The 0.3-12 keV spectrum is complex and cannot be fit by a single model, requiring a soft
component (described here by a blackbody with a temperature of 0.6 keV),
together with a hard power law (photon index of 1.5).
The spectral parameters of the soft component are similar to those
obtained by Parmar et al. (2001) during a BeppoSAX observation.
The source also displays a similar 1-10 keV luminosity (with
XMM-Newton we get
erg s-1 in the same energy range).
The XMM-Newton spectroscopy is compatible with what was found by Sidoli et al. (2001)
in ultracompact binaries located in galactic globular clusters.
Indeed, these sources have quite different properties
from all the other bright globular cluster LMXBs with longer orbital periods.
This spectral analogy led Parmar et al. (2001) to
suggest that XB1832-330 could be an ultracompact binary.
The total low-energy absorption resulting from the fits is compatible
with the
optically derived value in the direction of the host globular
cluster of
cm-2 in the first observation, while
it is slightly higher in the second, suggesting an extra-absorption in the line of sight
of
cm-2.
The likely presence of intrinsic neutral
extra-absorption seems to be confirmed by a good
fit when using the partial covering fractional absorption model,
which indicates that XB1832-330 is absorbed by an additional neutral
medium with a covering factor of
30% and an intrinsic
hydrogen column density in the range 3-
cm-2 for the two observations.
From the analysis of ASCA spectra of a few ultracompact binaries, Juett et al. (2001) suggested the presence of an excess absorption of neutral Ne-rich material local to the sources, leading some authors to propose that the donor stars in some ultracompact binary systems are Ne-rich white dwarfs (e.g., Yungelson et al. 2002). If this is the case, the high Ne/O ratio could be a spectral signature of the presence of some types of degenerate companions (e.g. neon-rich CO white dwarf or O-Ne-Mg WD donors). High-resolution spectra taken with XMM-Newton and Chandra confirmed a neutral Ne overabundance in a few cases (e.g., 4U 0614+091, Paerels et al. 2001; 4U 1543-624 and 2S 0918-549, Juett & Chakrabarty 2003). Instead, in other confirmed ultracompact sources, an overabundance of neutral Ne has not been found (e.g. 4U 1850-087, Sidoli et al. 2005; and Juett & Chakrabarty 2005) with high-resolution spectra, contrary to the earlier ASCA measurement (Juett et al. 2001). This indicates a variability in the Ne/O ratio probably due to luminosity variations (see e.g. Juett & Chakrabarty 2003; or Sidoli et al. 2005). On the other hand, several known ultracompact binaries have never shown evidence of unusual Ne/O ratios (e.g. 4U 1820-30, where an He white dwarf donor is present, or 4U 0513-40, Juett & Chakrabarty 2005).
We adopted a variable absorption model, together with three edges corresponding to O-K, Ne-K and Fe-L features in order to measure the column density of the neutral Ne, O, and Fe in the line of sight. The high-resolution XB1832-330 spectra presented here do not show any prominent emission features. Our study of the strength of the edges of neutral Ne and O reveals an Ne/O abundance ratio (see Table 3), which is consistent with the ISM value of 0.18 (Wilms et al. 2000). This implies that there is no evidence of the donor star in XB1832-330 being an Ne-enriched white dwarf.
If we compare the measured elemental column densities
(converted to equivalent H column
densities using the Wilms et al. (2000) ISM abundances) with the
resulting from the overall shape of
the XMM-Newton spectra, there seems to be a
discrepancy, which
could indicate an over-abundance of O and a sub-solar
abundance of Fe. This can be simply explained by the low
Fe abundance of the globular cluster.
We note that we did not assume any uncertainty in the
photoelectric cross-sections, which could be affected by
an uncertainty as large as 30% (Paerels et al. 2001). This could account for the observed
discrepancy. On the other hand, the resulting column density
is also dependent on the model adopted for the continuum.
A Comptonization model, for example,
leads to a lower absorption, because of the roll-over at soft energies.
The timing analysis resulted in a possible periodicity of 2.5 h in the EPIC lightcurve of the first observation and in a second periodicity (about twice what is found in Obs I) during Obs II. The smallest period, if linked to the orbital period of the binary system, would exclude an ultracompact binary (where orbital periods are less than about 1 h). We note, however, that, although the period is significant due to the limited duration of the observations, it is not possible to clearly say if this variability is a stable feature in the X-ray lightcurve or if it is a transient effect not linked to the orbital motion.
![]() |
Figure 8:
Colour-colour diagram of the detected sources. Open circles mark the simulated colours (obtained with an absorbed power-law spectrum, assuming different photon indexes and absorption columns). From top to bottom the photon index is 1, 2, and 3; from left to right different column densities have been assumed: 0.1, 0.5, 1, 3, 6,
![]() |
Open with DEXTER |
Our spatial analysis has led to the compilation
of a catalogue of 46 faint X-ray sources in the XB1832-330 field, most of which probably do not belong
to the globular cluster, as the NGC6652 tidal radius is 4.48
(the 2003 version of the catalogue is described in Harris 1996),
well inside the central CCD of the MOS cameras,
which was not used for source detection.
These sources are not
listed in the second XMM-Newton (2XMM) Serendipitous Source Catalogue (available online at
http://xmmssc-www.star.le.ac.uk/Catalogue/2XMM).
In Fig. 8 we report the colour-colour diagram of
the detected sources. For comparison, we overplot the simulated
colours obtained with an absorbed power-law spectrum, assuming
different photon indices and absorption columns. Most of the sources display a low absorption (
1020 cm-2), lower than the absorption towards NGC6652,
making them likely foreground objects.
Acknowledgements
Based on observations with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA member states and the USA (NASA). This work was supported by contract ASI/INAF I/023/05/0 and ASI/INAF I/008/07/0. We thank the referee (A. Juett) for constructive comments and suggestions which significantly improved our paper.
Table 4: The XMM-Newton Catalogue of sources in the region of X1832-330.