Issue |
A&A
Volume 504, Number 1, September II 2009
|
|
---|---|---|
Page(s) | 181 - 184 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/200911944 | |
Published online | 09 July 2009 |
A low luminosity state in the massive X-ray binary SAX J0635+0533
S. Mereghetti - N. La Palombara
INAF, Istituto di Astrofisica Spaziale e Fisica
Cosmica Milano, via E. Bassini 15, 20133 Milano, Italy
Received 25 February 2009 / Accepted 2 June 2009
Abstract
The X-ray pulsar SAX J0635+0533 was repeatedly observed with
the XMM-Newton satellite in 2003-2004. The precise localization
provided by these observations confirms the association of SAX J0635+0533 with a Be star. The source was found, for the first time, in a
low-intensity state, a factor 30 lower than that seen in
all previous observations. The spectrum, well fitted by an
absorbed power law with photon index
1.7 and
cm-2, was compatible with that of the
high state. The low flux did not allow the detection of the
pulsations at 33.8 ms seen in BeppoSAX and RXTE data. In view of the
low luminosity observed in 2003-2004, we reconsider the
peculiarities of this source in both the accretion and
rotation-powered scenarios.
Key words: X-rays: individuals: SAX J0635+0533 - X-rays: binaries
1 Introduction
The X-ray source SAX J0635+0533 was discovered with BeppoSAX in October
1997 (Kaaret et al. 1999) during a search for counterparts of
unidentified gamma-ray sources (Thompson et al. 1995). Its 2-10 keV flux of
erg cm-2 s-1, hard power-law
spectrum (photon index
1.5) extending to 40 keV, and
positional coincidence with a V = 12.8 star of Be spectral type,
immediately suggested classifying SAX J0635+0533 as an accreting high-mass
X-ray binary. The subsequent discovery of X-ray pulsations at 33.8 ms
(Cusumano et al. 2000) in the BeppoSAX data makes this object quite
peculiar and raises some problems for the accretion scenario. In
fact, if the X-ray emission is powered by accretion on the neutron
star surface, the magnetic field in SAX J0635+0533 must be about three
orders of magnitude less than expected in a typical high-mass
X-ray binary to avoid the propeller effect (see, e.g.,
Campana et al. 1995).
The pulse frequencies measured with RXTE in 1999, about two years
after the source discovery, indicated an orbital modulation with a
period of about 11 days and set a lower limit on the long-term
spin-down
s s-1(Kaaret et al. 2000). Such a high
in a rapidly spinning neutron
star implies a large rotational energy loss,
4
erg s-1,
capable of powering the observed X-ray luminosity
without the need of invoking mass accretion. In this
interpretation, SAX J0635+0533 would resemble other binary systems
composed of a fast pulsar orbiting a massive star, such as PSR
B1259-63 (Johnston et al. 1992), in which the X-ray emission is thought to
originate in the shock between the pulsar's relativistic wind and
that of the companion star. The failure to detect radio pulsations
from SAX J0635+0533 (Nicastro et al. 2000) does not rule out this scenario, since
beaming and/or absorption effects might render the radio emission
unobservable.
Here we report the results of XMM-Newton observations, carried out in 2003-2004, during which SAX J0635+0533 was detected at a very low flux level, the lowest ever seen from this source.
2 Observations and data analysis
SAX J0635+0533 was observed by XMM-Newton with ten different pointings between
2003 September 11 and 2004 April 14. The three EPIC focal plane
cameras (Strüder et al. 2001; Turner et al. 2001) were active during these pointings.
The two MOS cameras were operated in the standard full
frame mode (time resolution 2.6 s), in order to cover the whole
30
field-of-view. The pn camera was operated in timing
mode, providing a time resolution of 30
s, but
without imaging information. The thin filter was used for each
observation and focal plane camera.
Table 1: Flux and luminosity values of SAX J0635+0533 in the individual observations.
For each pointing, we retrieved the pps files from the
XMM-Newton archive produced by pipeline processing system,
which is operated by the Survey Science Center. Since
only timing-mode data were obtained with the pn camera, we
concentrated on the MOS data to study the source properties. The
faintness of SAX J0635+0533 during these observations made it impossible
to see the 33.8 ms pulsations in the pn camera data. For each
observation, we looked for possible periods of high instrumental
background caused by flares of soft protons with energies less
than a few hundred keV. To this aim, we selected only single and
double events (PATTERN
4) with energies greater than 10 keV
and recorded in the peripheral CCDs. Then, we set a
countrate threshold for good time intervals (GTI) at 0.5 cts s-1.
By selecting only events within GTIs, we finally
obtained a ``clean'' event list for each MOS data set. The
dates and effective exposure times (after soft-proton rejection)
of the observations are listed in Table 1.
After merging the event lists of the two MOS cameras, we
accumulated an image of the field-of-view comprising all the
observations. This clearly showed a source, with a count rate of
counts s-1 in each camera,
at the coordinates
,
(J2000). This position was obtained
after the astrometric correction of the X-ray coordinates of the
detected sources, based on the positions of five optical
counterparts found in the Guide Star Catalog (Lasker et al. 2008).
The final uncertainty on the source position is 1'' at 90% c.l.
(including statistical and systematic errors). The XMM-Newton localization confirms the association with the star proposed by
Kaaret et al. (1999) on the basis of an estimated 4% probability of
finding by chance a Be star in the 30'' radius error circle
determined with BeppoSAX (see Fig. 1).
We extracted the source spectra by selecting events in a circular
region with a small radius (10'') in order to minimize the
background contribution. The background spectra for the two
MOS cameras were accumulated from large circular areas with no
sources and radii of 100''. We generated ad hoc
redistribution matrices and ancillary files using the SAS tasks
rmfgen and arfgen, respectively. To ensure the
applicability of the
statistics, all spectra were
rebinned with a minimum of 30 counts per bin and fitted in the
energy range 0.3-10 keV using XSPEC 11.3.2. We checked that
separate fits of the two spectra gave consistent results,
therefore we analyzed them simultaneously, in order to increase
the count statistics, imposing common spectral parameters for the
two spectra and introducing a free relative normalization factor
between them, to account for possible differences in the
instrument cross-calibration. We found a 1.02
cross-normalization factor between the MOS2 and the
MOS1 camera.
A good fit was obtained with an absorbed powerlaw
(Fig. 2), yielding a hydrogen column density
cm-2 and a photon
index
.
The absorbed flux in the
energy range 0.2-12 keV is
erg cm-2 s-1,
while the corresponding unabsorbed flux is
erg cm-2 s-1. Acceptable fits
were also obtained with a thermal bremsstrahlung (kT =
9.2 keV) and with a blackbody (
keV).
![]() |
Figure 1:
Optical image of a 2' |
Open with DEXTER |
![]() |
Figure 2:
Top panel: average spectrum of SAX J0635+0533 with the
best-fit power-law model. The spectra of the MOS1 and
MOS2 cameras are shown in black and red, respectively.
Bottom panel: data-model residuals, in units of |
Open with DEXTER |
To study the source variability, we analyzed the ten datasets
separately. For each of them we performed a detailed source
detection in five energy bands (0.2-0.5, 0.5-1, 1-2, 2-4.5,
and 4.5-12 keV), applying the same procedure and parameters as
used by the XMM-Newton SSC to produce the XMM-Newton serendipitous
source catalog (Watson et al. 2009). The source detection was performed
simultaneously on both MOS data sets and using the
corresponding exposure maps, which account for spatial quantum
efficiency variations, mirror vignetting, and effective field of
view. The threshold value for the detection likelihood was set to
5 for the likemin parameter of the SAS task eboxdetect and to 6 for the mlmin parameter of the
SAS tasks emldetect and esensmap. We emphasize
that these values imply a rather loose constraint on the source
significance in order to be detected; they allow detection of
even very weak sources, with a low (i.e. 3)
signal-to-noise ratio. In this way we aim to check that SAX J0635+0533 is
even only marginally detected in any observation. On the other
hand, in our detection procedure we used ``ad hoc'' energy
conversion factors (ECF) to convert the measured count
rates of the detected sources (both in each of the five energy
bands and in the total 0.2-12 keV band) into the corresponding
energy flux. They were derived with the best-fit power-law
model of the source average spectra.
Based on this analysis, we find that SAX J0635+0533 was detected only in six observations, as reported in Table 1. The upper limits on the count rate (at the confidence level corresponding to a detection likelihood L = 6) are obtained from the sensitivity maps at the source position.
The long-term light curve of SAX J0635+0533 is plotted in
Fig. 3. In September-October 2003 the source flux
varied by at least a factor 10. Although the observations are not
continuous, they suggest an outburst lasting about three weeks,
with a rise time of only a few days to a maximum flux of
erg cm-2 s-1, followed by a similarly
rapid decay. A similar flux level was observed again six months
later. Comparing the hardness ratios measured in the different
observations, we found some evidence of a slight spectral
hardening correlated with the source intensity
(Table 1, Fig. 4).
![]() |
Figure 3:
Light curve of SAX J0635+0533. The count rates refer to the
0.2-12 keV energy range and to the sum of 2 MOS. The data of the first two
observations have been merged. The upper limits (obtained with a
threshold in detection likelihood L = 6) correspond to a
|
Open with DEXTER |
![]() |
Figure 4: Hardness ratio versus source flux. The hardness ratio is defined as (H-S)/(H+S), where H and S are the source count rates in the hard ( H = 2-12 keV) and soft ( S = 0.2-2 keV) energy ranges. |
Open with DEXTER |
3 Discussion
The maximum flux we observed for SAX J0635+0533,
erg cm-2 s-1
in the 2-10 keV range, is a factor >30
smaller than that measured at the time of the BeppoSAX discovery in
1997 (
erg cm-2 s-1, Kaaret et al. 1999).
To our knowledge, this is the lowest flux ever reported
for SAX J0635+0533. The upper limits of the XMM-Newton pointings of September
2003 imply an even smaller flux.
The distance of SAX J0635+0533 is not well constrained. The range 2.5-5 kpc was estimated by Kaaret et al. (1999) from the properties of the proposed optical counterpart. We also note that a distance far in excess 5 kpc is unlikely, considering the location of SAX J0635+0533 in the Galactic anti-center direction. In the following discussion, where we consider the two alternative possibilities for the origin of the observed X-ray emission, we conservatively normalize the relevant quantities to a distance d5 of 5 kpc.
3.1 Accretion-powered X-ray emission
Neutron stars accreting from Be star companions constitute the large majority of the high-mass X-ray binary systems present in the Galaxy. Also in view of our improved localization of SAX J0635+0533, it is natural to first discuss this possibility.
Already at the time of its discovery, it was noticed that the
luminosity of SAX J0635+0533 was relatively low compared to classical
Be/neutron star systems. The persistent sources of this class
have X-ray luminosity of 1036-1037 erg s-1. This
is also the luminosity typically reached during the outbursts of
transient Be/neutron stars systems, which comprise the majority of
this population. With the advent of more sensitive observations, a
number of persistent Be binaries with lower luminosity, 1034 erg s-1,
have also been discovered (see, e.g. Reig
& Roche 1999; La Palombara & Mereghetti 2006, 2007).
Our observations indicate for SAX J0635+0533 an average luminosity, a few
1033 d52 erg s-1, much lower than these values,
and provide an upper limit as low as 3
1032d52 erg s-1 in mid September 2003. If we further
consider that the source could well be closer than 5 kpc, we are
faced with an even lower luminosity.
The new data clearly indicate that SAX J0635+0533 is a transient source,
but it differs from the other Be systems for its low luminosity
both during the ``high state'' and in ``quiescence''. The non
detection in September 2003 allows us to set an upper limit on the
mass accretion rate of
g s-1. This
limit applies assuming that the accretion flow proceeds down to
the neutron star surface, which is very unlikely if the neutron
star is indeed rotating at 33.8 ms. In the presence of the
neutron star magnetic field, different scenarios preventing or
reducing the accretion rate onto the neutron star surface can
occur (see, for example, Campana et al. 1998). For such a short
spin period and low luminosity, the direct accretion regime, in
which the magnetospheric radius is smaller than the corotation
radius, hence the magnetic centrifugal barrier is open, requires a
magnetic field less than
108 G. This field is at least
three orders of magnitude lower than expected in a young neutron
star with a Be companion. In fact, all other accreting pulsars of
this class have much longer spin periods. The only exception is
the recurrent transient A 0538-66, which rotates at 69 ms
(Skinner et al. 1982). The pulsations in this systems were only detected
during a bright outburst reaching a luminosity of
1039 erg s-1 (Skinner et al. 1982), implying a magnetic field of
1011 G. Although below average, such a field is not
implausible. Furthermore it is consistent with interpretation of
the unpulsed quiescent luminosity of A 0538-66 (several 1033 erg s-1)
because accretion halted at the centrifugal barrier
(Campana et al. 1995).
If SAX J0635+0533 has a typical magnetic field, the low luminosity
observed with XMM-Newton could come from mass accretion stopped
at the magnetospheric radius. Assuming for simplicity spherically
symmetric accretion, along with a neutron star with mass 1.4
and radius 106 cm, gives in this case an X-ray
luminosity of
B12-4/7
erg s-1 (B12 is the magnetic field
in units of 1012 G and
the accretion rate in
units of 1015 g s-1; see, e.g., Campana et al. 1998).
However, the higher luminosity state observed in the past with
BeppoSAX and RossiXTE cannot be explained in the same way due to
the pulsations with a relatively high pulsed fraction, which are
not expected when the magnetic centrifugal barrier operates.
In conclusion, the difficulties already pointed out in interpreting SAX J0635+0533 as a typical accretion powered Be binary (Nicastro et al. 2000; Kaaret et al. 2000) are reinforced by the low luminosity reported here. Of course, the properties of this system would fit in this scenario better if the fast periodicity were disproved by further observations.
3.2 Rotation-powered X-ray emission
The alternative interpretation of SAX J0635+0533 is that of a
rotation-powered neutron star, whose X-ray emission derives from
the shock between the relativistic pulsar wind and the companion's
wind. The large luminosity difference between our data and the
previous observations could stem from varying shock conditions in
a very eccentric orbit. However, also in this scenario this source
would present some peculiar properties, compared to the
(admittedly few) other systems of this kind. The large variations
seen in September 2003 are difficult to explain if the source was
far from periastron, where no big changes in the shock properties
are expected. Rotation-powered pulsars in interacting binaries,
such as the already-mentioned PSR B1259-63 (Chernyakova et al. 2009,2006) or
the ``black widow'' pulsar PSR B1957+20 (Huang & Becker 2007), have X-ray
efficiencies in the range 10-4-10-2. On the other hand,
the
reported by Kaaret et al. (2000), corresponding to a
rotational energy loss
greater than a few
1038 erg s-1, implies a much lower efficiency for SAX J0635+0533
.
4 Conclusions
Despite SAX J0635+0533 being found by XMM-Newton in a very low-luminosity
state, thanks to the good sensitivity of the EPIC instrument, we
could derive a precise localization that confirms the proposed Be
optical counterpart of this source. We observed large flux
variability when the source was detected in September/October 2003
and derived a stringent upper limit of
d52 erg s-1 for its luminosity when the source was
not detected.
The spectral and flux properties of SAX J0635+0533 are consistent with a variable neutron-star binary powered by accretion from the Be companion or by the loss of rotational energy. Both interpretations imply some peculiarities with respect to other known sources. These result from the very short spin period and possibly high period derivative reported for SAX J0635+0533, which unfortunately we could not confirm owing to the source faintness during our observations.
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 NASA. The XMM-Newton data analysis is supported by the Italian Space Agency (ASI). The Guide Star Catalog II is a joint project of the Space Telescope Science Institute and the Osservatorio Astronomico di Torino. The Digitized Sky Surveys were produced at the Space Telescope Science Institute under US Government grant NAG W-2166.
References
- Campana, S., Stella, L., Mereghetti, S., & Colpi, M. 1995, A&A, 297, 385 [NASA ADS] (In the text)
- Chernyakova, M., Neronov, A., Lutovinov, A., Rodriguez, J., & Johnston, S. 2006, MNRAS, 367, 1201 [NASA ADS] [CrossRef]
- Chernyakova, M., Neronov, A., Aharonian, F., Uchiyama, Y., & Takahashi, T. 2009, [arXiv:0905.3341]
- Cusumano, G., Maccarone, M. C., Nicastro, L., Sacco, B., & Kaaret, P. 2000, ApJ, 528, L25 [NASA ADS] [CrossRef] (In the text)
- Huang, H. H., & Becker, W. 2007, A&A, 463, L5 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Johnston, S., Manchester, R. N., Lyne, A. G., et al. 1992, ApJ, 387, L37 [NASA ADS] [CrossRef] (In the text)
- Kaaret, P., Piraino, S., Halpern, J., & Eracleous, M. 1999, ApJ, 523, 197 [NASA ADS] [CrossRef] (In the text)
- Kaaret, P., Cusumano, G., & Sacco, B. 2000, ApJ, 542, L41 [NASA ADS] [CrossRef] (In the text)
- Lasker, B. M., Lattanzi, M. G., McLean, B. J., et al. 2008, AJ, 136, 735 [NASA ADS] [CrossRef] (In the text)
- Nicastro, L., Gaensler, B. M., & McLaughlin, M. A. 2000, A&A, 362, L5 [NASA ADS] (In the text)
- Skinner, G. K., Bedford, D. K., Elsner, R. F., et al. 1982, Nature, 297, 568 [NASA ADS] [CrossRef] (In the text)
- Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 [NASA ADS] [CrossRef] [EDP Sciences]
- Thompson, D. J., Bertsch, D. L., Dingus, B. L., et al. 1995, ApJS, 101, 259 [NASA ADS] [CrossRef] (In the text)
- Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 [NASA ADS] [CrossRef] [EDP Sciences]
- Watson, M. G., Schröder, A. C., Fyfe, D., et al. 2009, A&A, 493, 339 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
All Tables
Table 1: Flux and luminosity values of SAX J0635+0533 in the individual observations.
All Figures
![]() |
Figure 1:
Optical image of a 2' |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Top panel: average spectrum of SAX J0635+0533 with the
best-fit power-law model. The spectra of the MOS1 and
MOS2 cameras are shown in black and red, respectively.
Bottom panel: data-model residuals, in units of |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Light curve of SAX J0635+0533. The count rates refer to the
0.2-12 keV energy range and to the sum of 2 MOS. The data of the first two
observations have been merged. The upper limits (obtained with a
threshold in detection likelihood L = 6) correspond to a
|
Open with DEXTER | |
In the text |
![]() |
Figure 4: Hardness ratio versus source flux. The hardness ratio is defined as (H-S)/(H+S), where H and S are the source count rates in the hard ( H = 2-12 keV) and soft ( S = 0.2-2 keV) energy ranges. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.