A&A 486, 293-302 (2008)
DOI: 10.1051/0004-6361:20078707
J. Barnstedt1 - R. Staubert1 - A. Santangelo1 - C. Ferrigno2 - D. Horns1 - D. Klochkov1 - P. Kretschmar3 - I. Kreykenbohm1,5 - A. Segreto2 - J. Wilms4
1 - Kepler Center for Astro and Particle Physics,
Institut für Astronomie und Astrophysik, Universität Tübingen,
Sand 1, 72076 Tübingen, Germany
2 -
INAF IFC-Pa, via U. La Malfa 153, 90146 Palermo, Italy
3 -
ESA, ESAC, PO Box 78, 28691 Villanueva de la Cañada, Spain
4 -
Dr. Karl Remeis-Sternwarte, Astronomisches Institut, Universität Erlangen-Nürnberg,
Sternwartstr. 7, 96049 Bamberg, Germany
5 -
INTEGRAL Science Data Centre, Chemin d'Écogia 16, 1290 Versoix, Switzerland
Received 19 September 2007 / Accepted 6 May 2008
Abstract
Context. The Galactic Plane Scan (GPS) was one of the core observation programmes of the INTEGRAL satellite. The highly variable accreting pulsar OAO 1657-415 was frequently observed within the GPS.
Aims. We investigate the spectral and timing properties of OAO 1657-415 and their variability on short and long time scales in the energy range 6-160 keV.
Methods. Using standard extraction tools and custom software for extracting INTEGRAL data we analysed energy-resolved light curves with a time resolution of one second - mainly data of the ISGRI instrument. We also analysed phase-averaged broad band spectra - including JEM-X spectra - and pulse-phase resolved spectra of ISGRI.
Results. During the time covered by the INTEGRAL observations, the pulse period evolution shows an initial spin-down, which is followed by an equally strong spin-up. In combining our results with historical pulse period measurements (correcting them for orbital variation) and with stretches of continuous observations by BATSE, we find that the long-term period evolution is characterised by a long-term spin-up overlayed by sets of relative spin-down/spin-up episodes, which appear to repeat quasi-periodically on a 4.8 yr time scale. We measure an updated local ephemeris and confirm the previously determined orbital period with an improved accuracy. The spectra clearly change with pulse phase. The spectrum measured during the main peak of the pulse profile is particularly hard. We do not find any evidence of a cyclotron line, wether in the phase-averaged spectrum or in phase-resolved spectra.
Key words: X-rays: binaries - pulsars: individual: OAO 1657-415
![]() |
Figure 1:
Orbital profile from 2687 individual science
windows in the ISGRI energy range 20-40 keV for source
offset angles <
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The accreting pulsar OAO 1657-415 was discovered by
Polidan et al. (1978) with the Copernicus satellite.
White & Pravdo (1979) detected a 38.22 s pulsation,
whereas the binary orbit and the eclipsing nature of the high mass
binary system were discovered by Chakrabarty et al. (1993). The orbital
period was later improved by Bildsten et al. (1997) from BATSE
measurements to
.
As
OAO 1657-415 lies in a heavily absorbed region, an optical counterpart
could not be identified up to a limit of V>23, but an infrared
counterpart was found by Chakrabarty et al. (2002). The infrared
properties of OAO 1657-415 were found to be consistent with a highly
reddened B supergiant. Analysing ASCA observations
Audley et al. (2006) found a dust-scattered X-ray halo whose
evidence decays through the eclipse. Using this halo they
estimated a distance of
kpc. This is consistent with
a source distance of
kpc derived by
Chakrabarty et al. (2002) from photometric infrared
measurements. The pulsar was found to show a long-term spin-up
from 38.218 s (White & Pravdo 1979) to 37.329 s
(Chakrabarty et al. 2002) with intermediate spin-down periods
(Bildsten et al. 1997) (see Table 2). This
corresponds to an averaged long-time period derivative of
s s-1. For a time span of about
3300 days there are nearly continuous BATSE observations available
with daily pulse period and flux values. Analysing the spin-up and
spin-down episodes of BATSE data, Baykal (1997)
found that the most natural explanation of X-ray flux and angular
acceleration fluctuations is the formation of episodic accretion
disks from stellar wind accretion. Lovelace et al. (1999) developed a
model for magnetic, propeller driven outflows that cause a rapidly
rotating magnetised star accreting from a disk to spin-down.
Another model was introduced by Dai & Li (2006) in
which they assume, in contrast to the model of Lovelace et al. (1999),
that accretion continues during the (soft) propeller stage. This
avoids otherwise occurring singularity problems and also leads to
the fact that no significant changes in the X-ray luminosity occur
during the spin transitions.
The spectral models commonly used for fitting the spectra of OAO 1657-415 usually include an absorbed power law model for the low energy part of the spectrum. For the high energy part an exponential cutoff power law model or a power law with high energy cutoff is used (see Sect. 4.1, for a definition of these models).
Orlandini et al. (1999) searched for cyclotron lines in BeppoSAX spectra of OAO 1657-415 and were able to fit a broad Lorentzian line at 30 keV with a width of 28 keV, using the exponential cutoff model for the continuum. As the width was quite large and not well constrained, they fixed the width at 10 keV, which changed the line position to 36 keV. Using the high energy cutoff continuum model, it was found that adding an absorption line did not significantly improve the fit. So Orlandini et al. (1999) concluded that the spectrum of OAO 1657-415 above 10 keV can be modelled by a series of three power laws with increasing steepness, such that no cyclotron line is needed.
Also, a fluorescence iron line emission was observed in this source with several instruments, e.g. ASCA (Audley et al. 2006), Chandra/HETGS (Chakrabarty et al. 2002), RXTE (Baykal 2000), BeppoSAX (Orlandini et al. 1999) and Ginga (Kamata et al. 1990).
The satellite INTEGRAL (Winkler et al. 2003a) is
equipped with three instruments for the X- and
-ray range: the spectrometer SPI (20 keV-8 MeV)
(Vedrenne et al. 2003), the imager IBIS (15 keV-10 MeV)
(Ubertini et al. 2003), and the X-ray monitor JEM-X (3-60 keV)
(Lund et al. 2003). The imager IBIS consists of the low energy CdTe
detector layer ISGRI (15-1000 keV) (Lebrun et al. 2003) and the
high energy CsI detector layer PICsIT (
0.175-10.0 MeV)
(Labanti et al. 2003). IBIS has a fully coded field of view (FOV) of
and a partially coded FOV of
(50%). In this
paper we have used ISGRI observations of OAO 1657-415, which were mainly
part of the INTEGRAL Galactic Plane Scan
(GPS)
.
Being part of the core program of INTEGRAL, the GPS
observations were performed as a sawtooth pattern along the
accessible part of the Galactic plane with an extension in
latitude of
(Winkler et al. 2003b). A total of 2687 science windows were available. One INTEGRAL science window
(SCW) is a single observation of about 2000 s, so the total time
of data acquisition was roughly
5 Ms. The data used range
from revolution 36 (MJD 52 668) up to revolution 463 (MJD 53 946).
Spectra were obtained also from JEM-X data. Due to a smaller FOV
(fully coded
,
partially (50%) coded
), only
41 science windows (SCWs) from JEM-X1 and 46 SCWs from
JEM-X2 were available (selected from the SCWs used for ISGRI
data, but with the angular distance of the centre of the FOV to
the source restricted to <
). These were JEM-X2
observations of February and March 2003 and JEM-X1 observations
from August 2004 to April 2005. As these 87 SCWs are just a small
fraction of the total SCWs available from ISGRI
observations, only phase-averaged spectra were produced
(OSA 6.0), but no timing analysis was done for JEM-X data.
Generally, the observational data have been reduced using the
Offline Scientific Analysis (OSA) software v. 6.0 provided by the
INTEGRAL Science Data Centre, ISDC (Courvoisier et al. 2003). For
generating spectra and pulse profiles from ISGRI data
(spectral timing), an alternative software provided by
the IASF Palermo (Ferrigno et al. 2007) was used.
For the timing analysis (including spectral timing
analysis) we selected all observations for which OAO 1657-415 was within
from the centre of the FOV (this is recommended by ISDC
since the reliability of the flux determination is reduced when
the source is located in the outer parts of the partially coded
FOV of ISGRI). For all existing science windows we determined the
orbital phase and selected only data in the orbital phase interval
from 0.1 to 0.9, thus excluding data during the eclipse.
Furthermore, we included only those science windows for which the count rate is above 15 s-1 in the energy interval 20-40 keV, corresponding to the ``high state'' (see below). The reason for this is that the timing analysis (the determination of the pulse period and in particular the determination of the pulse phase), needs data with a sufficient signal to noise ratio to yield significant results. All above mentioned restrictions left 842 science windows out of 2687 available (see also Fig. 1, which shows data from all available SCWs). Our results are therefore applicable to this ``high state'' only.
![]() |
Figure 2:
Averaged orbital profile in the ISGRI energy band
20-40 keV. This profile was calculated from 1797 science
windows for which the angular distance to the telescope axis was
less than
![]() |
Open with DEXTER |
Figure 1 shows the mean count rates from
2687 science windows in the ISGRI energy band 20-40 keV as a
function of orbital phase, which was calculated using the
ephemeris of Bildsten et al. (1997) (
and
d). The large scatter in this
figure reflects the strong intrinsic variability of OAO 1657-415.
Figure 2 shows the averaged orbital
profile, which was calculated from 1852 science windows (MJD
52 668-53 946) for which the angular distance of OAO 1657-415 from the
centre of the FOV was less than
.
The dashed line shows
the number of science windows for each data point. In general, we
confirm the asymmetric shape of the profile as it was
measured by Wen et al. (2006) with ASM on board
RXTE. In contrast to these measurements, the
INTEGRAL observations show a remarkable minimum at
phase 0.55, which is visible also in the ASM data averaged for the
time range MJD
52 668-53 946 (Fig. 3).
The data presented by Wen et al. (2006) cover a time
range of roughly 8.5 years (
MJD
50 160-53 240), which is
much wider than the time range we analysed. We conclude therefore
that the dip seen in our data might be a temporary phenomenon.
![]() |
Figure 3: Averaged orbital profile (ASM, 2-10 keV) plotted from quick-look results provided by the ASM/RXTE team. This curve is calculated from 12 830 single dwell values (sum band intensity) covering the used INTEGRAL observation period from MJD 52 668-53 946. The error bars are quadrature averages of the flux errors. The dashed line shows the number of dwells per data point (right axis). |
Open with DEXTER |
The lack of data points around 10-15 counts s-1 (cps) at orbital phases above 0.4 apparent in Fig. 1 may suggest the existence of a bimodal intensity distribution, defining a ``high state'' and a ``low state'' of this source.
The formal uncertainties of the count rates in Fig. 2 are negligible, but an unpredictable scatter is introduced by the fact that only in rare cases the continuous measurements span a complete orbit. Due to the intrinsic variability of the source the available data points do not represent a homogeneously averaged profile. The number of data points (Science Windows) per phase bin used to calculate the average flux varies between 40 and nearly 170. To evaluate the influence of the non-homogeneity of this sampling, we calculated profiles with a constant number of 30 SCWs per phase bin, randomly chosen from the available SCWs. We calculated 10 000 of these profiles and plotted the range between minimum and maximum values for each phase bin as ``error bars'' to Fig. 2. A discussion is given in Sect. 5.
![]() |
Figure 4:
Determination of the centre of eclipse by fitting
straight lines to ingress and egress. The dashed lines denote the
width of the strips (0.03 phase units), which were used to
determine the number of data points falling into the strips left
and right to the fitted straight lines (see text). Only data
points from SCWs with offset angles <
![]() |
Open with DEXTER |
Using the scatter plot shown in Fig. 1 we
have determined the properties of the eclipse. The data points
left and right of the apparent eclipse define a rather sharp
boundary in phase space (orbital phase vs. count rate), leaving
only a few individual outlying points. We fitted these boundaries
by straight lines by maximising the differences in data points
that fall into strips of predefined widths left and right to these
lines (see Fig. 4). For ingress we fitted a line
from 37 cps down to the maximum eclipse intensity of 3 cps
and the egress line from 3 cps to 20 cps. The strip widths were
varied from 0.01 to 0.05 in phase units. No significant change of
the fit parameters was observed for strip widths in the range
0.02-0.04. We therefore used the average values obtained from
the fits with 3 widths of 0.02, 0.03 and 0.04. The intersections
of the ingress and egress lines with the bottom count rate are at
orbital phases 0.899(5) and 1.052(5) (corresponding to
MJD 52 662.84(5) and MJD 52 664.44(5), respectively), such that
we set the centre of eclipse to phase 0.976(5)
(MJD 52 663.64(5)), and the half width of the full eclipse to 0
80(3) (or
). The MJD values are
determined in such a way that they fall into the last orbit before
our observations. The straight lines for ingress and egress are
described by slopes (change of count rates) of 34 cps d-1 and 28 cps d-1, respectively, making the ingress faster than
the egress by a factor of 1.2. We note that the above description
gives upper limits to the speed of ingress or
egress, individual events may be slower.
Table 1: Pulse periods of OAO 1657-415 as measured during the INTEGRAL GPS in 11 observation groups. The date is the centre of the observation time interval.
The pulse period analysis was performed in steps,
successively improving the accuracy of the analysis. This
was necessary because of the strong variation in the pulse period
over the time covered by the INTEGRAL GPS observations. We
started out with an ISGRI light curve in the 20-200 keV energy
range with a time resolution of 1 s (produced by the IASF Palermo
software
). The bin times were first converted
to the solar system bary centre and then corrected for binary
motion using the ephemeris of Bildsten et al. (1997) (T90 =
MJD 48 515.99(5),
,
lt-s,
e=0.104(5) and
). For
each SCW a period search with epoch folding was performed. This
led to first estimates of the pulse periods and already clearly
showed a strong variation with time. As a next step we produced
pulse profiles combining four to five SCWs. A phase connection
analysis of these profiles - by using the mean profile as
a template and fitting the individual profiles to this template -
(Deeter et al. 1981; Nagase 1989; Muno et al. 2002) yielded significantly more
accurate periods. A further refinement was reached using a
smoothing technique: pulse profiles were generated (with the so
far best pulse period) for data sets of at least eight SCWs, where
the start of each data set was shifted by one SCW as compared to
the start of the previous one. This produced ``running mean pulse
profiles'' with good photon statistics, from which the final pulse
periods were again determined by the phase connection technique
(the loss of statistical independence between consecutive
pulse profiles does not pose a problem here). The above described
procedure was performed for each of the eleven groups of
observations (the groups contained 12 to 273 SCWs). The final
periods are listed in Table 1 and plotted in
Fig. 5. Due to the short integration
times it was generally not possible to measure a period derivative
within one group of observations - except for group no. 2, for
which we find a local
of
s s-1 (see Table 1). Phase
connection between the groups was not possible, as the large
gaps in time did not allow the unambiguous determination of the
number of pulsar periods.
![]() |
Figure 5:
Period evolution of OAO 1657-415 as determined in this work
from INTEGRAL ISGRI data (Table 1).
The typical error of 0.004 s shown here was assumed for periods
estimated by the epoch folding method. The actual errors for the
plotted periods estimated by the phase connection method are in
general about a factor of 10 lower. The ![]() |
Open with DEXTER |
The period evolution over the time covered by the
INTEGRAL GPS (Fig. 5) is basically a
triangular function showing an initial spin-down with a mean
s s-1, followed by a spin-up
with a mean
s s-1. As will be
shown below, the long-term mean spin-up in OAO 1657-415 is
s s-1, meaning that the values
measured by INTEGRAL during the spin-down and spin-up
episodes are symmetric with respect to the long-term mean
value of the period derivative.
![]() |
Figure 6: Long-term period evolution of OAO 1657-415. The historic data from the literature (see Table 2) were corrected for binary motion. The BATSE data are from public archive (1, 3) and private communication by Mark Finger (2). The sinusoid with a period of 4.8 yr, superposed on a linear long-term spin-up, is to guide the eye to emphasise a quasi-periodic behaviour. |
Open with DEXTER |
Table 2: Historic pulsar periods, corrected for binary doppler shift. The errors given for the corrected periods are taken from the errors given for the observed periods (except for the value markeda).
Figure 6 shows the long-term evolution of the
pulse periods for all data available to us. There are essentially
three groups of data. First group: the ``historical periods'' as
taken from the literature and summarised in
Table 2. All except one of these reported
periods were without binary correction, since most of them were
determined before the detection of the binary nature of OAO 1657-415.
Therefore we have applied the binary correction using the orbital
parameters of Bildsten et al. (1997). Second group: BATSE data as taken
from public archives,
and a complementing (and partially
overlapping) data set provided through private communication by
Mark Finger. The third group is represented by the data we
processed from INTEGRAL observations presented in this
paper.
The complete set of data in Fig. 6 is described
by a mean spin-up rate of
s s-1. On top of this long-term mean spin-up
deviations are observed with episodes of strong relative
(
)
spin-up and spin-down. In some
of the spin-down episodes even the absolute spin-down
is
quite strong (as in the first half of our measurements
around MJD 53 000). The appearance in time of the deviations is
suggestive of a quasi-periodic behaviour. We do not claim a formal
value for the period, but show in Fig. 6 a model
modulation by a sine curve with an amplitude of 0.13 s and a
period of 1750 d (4.8 yr), which reproduces the main features of
the long-term evolution. One could consider this period
as the characteristical time for a typical spin-down/spin-up
cycle.
After having used the binary ephemeris of Chakrabarty et al. (1993) and Bildsten et al. (1997) for the pulse period analysis (as discussed in Sect. 3.2), we have analysed the non-binary-corrected light curve data from the INTEGRAL GPS in order to establish our own local ephemeris. We again made use of ``running mean pulse profiles'' (from sets of in general eight SCWs, the start of each set shifted by one SCW as compared to the previous one) and the phase connection techniques. For each of the 11 observation groups (see Table 1) a local ephemeris was established, which - except for observation 2 - are generally not very accurate because the observations are short and the intrinsic pulse period shows a strong variation with time (see above). Observation 2, however, samples about 2.5 binary cycles reasonably densely, such that a fit of the pulse arrival time delays (Fig. 7) can be performed and an ephemeris established. We find (keeping the orbital elements as before): T90= MJD 52 663.893(10). Adding all other observations to the analysis yields the same result (with increased uncertainty).
There is no indication for a secular change of the orbital period.
Assuming
,
we can combine our new
ephemeris with values from the past to arrive at a new sidereal
period. Using, e.g., the T90 of Chakrabarty et al. (1993) of
MJD 48 515.99(5), and dividing the time difference to our
ephemeris by 397 orbital cycles, we arrive at
.
Our
is consistent with the value of Bildsten et al. (1997)
within our uncertainty, which is reduced by a factor of 2.5 in
comparison to that of Bildsten et al. (1997).
Table 3 shows the current best set of
orbital parameters.
We have done the same for the observed times of mid-eclipse: if
the elapsed time between our mid-eclipse
(MJD 52 663.64(5), see above) and that of
Chakrabarty et al. (1993) (their Table 1: MJD 48 515.897(72)) is divided
by 397 cycles we find
.
The period found from eclipse timing is consistent
(within uncertainties) with that found using the orbital ephemeris
determined from pulse timing analysis. We point out that this need
not be the case: with the finite eccentricity of e=0.104 apsidal
motion is expected, and
is only equal to
if
is found from averaging over
one (or more) complete apsidal period(s) (see e.g.
Deeter et al. 1987). But so far, the apsidal period of OAO 1657-415 is
not known.
For the time covered by the INTEGRAL GPS observations we
find that the centre of eclipse is earlier than our
T90 by
.
Using the values in
Table 1 of Chakrabarty et al. (1993) we estimate that also in 1991/92 the
eclipse was earlier than T90 by
.
In principle, such observed time differences
can be
used to find the longitude of periastron. To first order
(Deeter et al. 1987), where
is the longitude of periapsis
and e the eccentricity. From the change of
the rate of
apsidal advance (and the apsidal period) can be found.
Unfortunately, neither the data of Chakrabarty et al. (1993) nor our own
data are accurate enough to yield a meaningful constraint.
![]() |
Figure 7: Observed pulse delay times for the second observation period, which covers about 2.5 orbital cycles. The solid line shows the best fit for the orbital parameters (see Sect. 3.4). The residuals (observed - calculated, lower panel) show some variation that may well be due to real pulse period variations on time scales shorter than one orbital period. |
Open with DEXTER |
Table 3:
Current best set of orbital parameters:
T90 and period as estimated in this work,
,
eccentricity and
taken from Bildsten et al. (1997).
![]() |
Figure 8: Energy-resolved pulse profiles. The relative intensity scale is normalised to the average count rate level (horizontal dashed line) in the corresponding energy interval. The percentage values indicate the pulsed fraction. The pulsed fraction of the highest energy interval was set to 100%, as the negative minimum flux values were set to zero. The top panel shows the hardness ratio (see text). The vertical dashed lines denote the four phase intervals for which spectra were extracted. The pulse phase limits of the four intervals are: a=0.025-0.250, b=0.25-0.50, c=0.50-0.75 and d=0.750-1.025. |
Open with DEXTER |
![]() |
Figure 9: Variation of the pulsed fraction with energy. Horizontal bars denote the energy range for which the pulsed fraction was calculated. |
Open with DEXTER |
With the above determined pulse periods, pulse profiles in seven different energy intervals (covering the range 20-160 keV) of 720 SCWs were produced and superposed (using the sharp rise of the main pulse as phase reference - we call this pulse phase 0). The final energy-resolved pulse profiles are displayed in Fig. 8. The pulsation is clearly detected up to the highest energy range (120-160 keV).
The profiles consist of at least three components, the amplitudes
of which decrease during the course of the pulse. For easier
comparison the profiles for the different energy intervals are
shown in the same relative scale (except for the energy range
120-160 keV, where the plot needs a different scale due
to the higher pulsed fraction and the larger errors). The pulsed
fraction - indicated as percentage in Fig. 8
- increases linearly with energy in the range
30-70 keV (Fig. 9). We calculated the
pulsed fraction from the minimum and maximum flux values in the
corresponding energy interval as
.
This is - at least in the case of sinusoidal pulse profiles -
equivalent to the standard definition
,
calculated from the values of the pulsed and the constant flux
parts. It should be noted that negative
values
were set to zero to avoid calculated pulsed fractions larger than 100%, which is the case for the highest energy interval shown in
Fig. 8.
Also obvious from the plot is that the first maximum becomes asymmetric at higher energies, which indicates a harder spectrum in phase interval a (defined in Fig. 8). The top panel of Fig. 8 shows the hardness ratio calculated as ratio of the photon fluxes in the 40-160 keV band to those in the 20-40 keV band.
![]() |
Figure 10: Phase-averaged broadband spectrum for ISGRI, JEM-X1, JEM-X2. The fitted model was HIGHECUT. To avoid artefacts in the residuals of the fit, the transition to the high energy cutoff was smoothed with an additional Gaussian absorption. The fit parameters are given in Table 4. |
Open with DEXTER |
Table 4:
Table of fit parameters for the phase-averaged broad band
spectrum shown in Fig. 10 for two
different spectral models (see text). A systematic
error of 2% was assumed. The uncertainties given are
at 90% confidence level for single parameter errors
(
).
Due to calibration uncertainties below 20 keV and low flux above 160 keV the useable range of the ISGRI spectra was 20-160 keV. The JEM-X spectra were used in the 6-20 keV range, which we consider as the most reliable range. All spectra are background subtracted.
We used the following components for our spectral models:
an absorbed power law model (
,
XSPEC model phabs*powerlaw - in the
following ABSPL), which influences only the JEM-X part of
the spectra. The second component was an exponential
cutoff power law model (
,
XSPEC
model cutoffpl - in the following CUTOFFPL) or a power
law with high energy cutoff
(
if
,
else
;
XSPEC model
powerlaw*highecut - in the following short HIGHECUT).
These components are mainly determined by the ISGRI part
of the spectra.
Being not able to fit the spectra with CUTOFFPL plus
cyclotron absorption as described by Orlandini et al. (1999), we used
in our analysis the HIGHECUT model (Model A in
Table 4). In order to avoid artefacts in the
residuals usually produced by the XSPEC highecut model at
(Kretschmar et al. 1997), the function was smoothed with
a Gaussian absorption gabs as successfully exercised by
Coburn et al. (2002). The centre of this Gaussian was set to
plus 1 keV (from experience with other spectra and
confirmed by the fitting procedure). The width
and the
depth
of the Gaussian were found to be 2.2 keV and 0.9,
respectively, and were subsequently fixed at these values.
The photon index
for the ISGRI part of the phase-averaged
spectra (Fig. 10) was fixed to unity, which is
justified by the JEM-X2 part of the spectrum and is consistent
with other spectral fits reported in literature (see
Table 4). The JEM-X1 spectra show a
significantly different photon index of 1.3 (the
errors for all
values are 0.05). There are several
possible explanations for this behaviour: first, the JEM-X1
spectra cover the orbital phases from 0.25 to 0.65, while the
JEM-X2 spectra are mainly located in the phase intervals
0.1-0.25 and
0.65-0.8. Second, they originate from different
time intervals: JEM-X2 from MJD
52 671-52 869 and JEM-X1 from
MJD
53 232-53 795. A third possibility might be the existence
of uncertainties in the spectral calibration of the two
instruments.
This difference in the JEM-X spectra could also be well
fitted by an additional photon absorption in the JEM-X2 spectra
(ABSPL, see Sect. 4.1). Model B
in Table 4 shows the parameters in this model,
which is formally Model A multiplied by the XSPEC model
phabs. The photon index
was kept identical for
all 3 spectral parts (ISGRI, JEM-X1, JEM-X2). No absorption was
assumed for the ISGRI and JEM-X1 part, while for JEM-X2 cold
matter absorption was fitted. The fit led to a
of
,
which is within the
range of other literature values shown in
Table 4. This finding could imply that the
spectra near the eclipse are more absorbed than in the centre of
the orbital profile. On the other hand it could also mean that at
least part of the strong variability of OAO 1657-415 is due to variable
absorption, which was higher during the JEM-X2 part of the
observation.
![]() |
Figure 11: Phase-resolved ISGRI spectra for the four phase intervals as defined in Fig. 8. The solid line in the upper panel is the model fit of the phase-averaged spectrum as shown in Fig. 10 (model A of Table 4). To enhance the differences between the spectra of the four phase intervals the lower panel gives the ratios between the four spectra and the model of the phase-averaged spectrum in a logarithmic scale. The dotted line in the lower panel corresponds to a fit of spectrum b with a line feature at 49.3 keV. |
Open with DEXTER |
We have investigated ISGRI spectra from the four pulse phase
intervals indicated in Fig. 8. The spectra
are shown in Fig. 11. We fitted all data
using the HIGHECUT model, keeping the power law index
fixed to
.
The results are given in
Table 5.
Table 5:
Fit parameters for the spectral model HIGHECUT,
Model A in Table 4.
In the last line the fit parameters for the second
phase interval are given for the same model, but
multiplied by a cyclotron absorptiona.
The uncertainties given are at 90% confidence level
for single parameter errors (
). No systematic errors were considered.
As cyclotron resonance scattering features are observed in several
pulsars (see discussion in Sect. 5.2) we
tried to find possible evidence of a cyclotron line feature. We
therefore calculated the ratio of the spectra from the four phase
intervals to the phase-averaged spectrum (model A of
Table 4 with ISGRI parameters). In this way,
differences between the individual spectra and features within
these spectra are enhanced. The four individual spectra and the
ratios are shown in Fig. 11. While spectrum c
(pulse phase
0.5-0.75) is nearly identical to the pulse-phase
averaged spectrum (the ratio in the lower panel of
Fig. 11 is around 1.0), spectrum a (phase
0.025-0.25) is clearly harder and spectrum d (phase
0.75-1.025) is clearly softer. Spectrum b (phase
0.25-0.5)
shows a slight ``peculiarity'': there seem to be two components to
the spectrum with a transition between 45 keV and 55 keV, or
alternatively, a depression around these energies. Formally, a
cyclotron line can be fitted with a centroid energy of 49.3 keV.
The
is reduced to 0.9 (d.o.f. = 41) from 1.2
(d.o.f. = 44), resulting in an F-Test value of 1%,
which is usually not considered as sufficient to state
the detection of a cyclotron absorption line - also one
should keep in mind that it is in general not correct to use an
F-Test to test for the presence of a line (Protassov et al. 2002).
In the averaged light curve (Fig. 2) we found a dip at phase interval 0.55, which is confirmed by ASM measurements (Fig. 3) taken in the same time interval as the INTEGRAL observations, but homogeneously sampled. The dip is more pronounced in our INTEGRAL data, but we cannot exclude that this is an effect from the inhomogeneous sampling in the INTEGRAL data (see dashed line in both figures) with a minimum of measurements at the minimum of the profile. In order to evaluate the effects of the sampling we ran a statistical simulation, used to produce the ``error bars'' in Fig. 2. These bars show the range between the minimum and maximum values evaluated from 10 000 profiles calculated from 30 randomly chosen SCWs per phase bin. The dip around orbital phase 0.5 is always present, in line with the ASM results.
We note that the average profile should not be seen as a template for individual light curves, it rather seems to be the result of a random process. Some individual light curves have a maximum in the phase interval 0.2-0.5, others (and more) have a maximum in the interval 0.6-0.85. So the strong variability of the source may be weakly correlated with the orbital phase in such a way that the chance for a ``low state'' is quite high at phase 0.55.
The physical origin of this orbital modulation is still unclear. Topics to be discussed in this context should comprise absorption, scattering, a (partial) blockage of the line of sight to the neutron star, or variations of the accreted wind. All of these mechanisms will be subject to further investigations.
As the pulsar can be considered as a point source compared to its B supergiant companion, the slope of the ingress and egress indicates a soft transition from the stellar atmosphere to the stellar wind region with a decreasing optical thickness for the radiation emitted by the pulsar. The density profile of the outer regions of the supergiant was modelled by Bulik et al. (2005) with an exponential density distribution. They fitted the modelled absorption to the observed folded light curve, which was composed from 480 INTEGRAL SCWs. This model did not account for the asymmetric shape of the orbital profile and should be refined in a further analysis.
The INTEGRAL observations extend the long-term monitoring
of OAO 1657-415 and cover a period with a torque reversal from strong
spin-down to strong spin-up. We have compiled all available pulse
period measurements from the literature and, where necessary,
converted the values to the intrinsic neutron star spin period. We
find that the long-term spin-up is continued and determine a mean
rate of
s s-1. The
spin-down/spin-up behaviour observed by INTEGRAL, together
with the historical data imply a quasi-periodic variation
of the pulse period with a period of roughly 5 years, superimposed
on the long-term mean spin-up.
The orbital elements were determined from pulse timing
analysis. Combining our results with previous data we determine an
orbital period with improved accuracy:
.
We find a difference
between T90 and the centre of eclipse of
,
which is consistent with the previously determined value, and does
not allow a conclusion about a possible apsidal advance.
The long-term evolution of the spin period is similar to that of Her X-1. In this object a quasi-periodic spin-up/down variation is found with a period of about 5 yr (as in OAO 1657-415), superimposed on a long-term spin-up (Staubert et al. 2006). Since in Her X-1 the period variations are correlated with the 35 d turn-on behaviour and the quasi-periodically repeating anomalous low states, this phenomenon is generally connected to the precessing accretion disk and an associated variation in the mass accretion onto the neutron star.
In contrast to Her X-1, however, OAO 1657-415 is a high mass, wind accreting system and we have no direct evidence of an accretion disk. Also no correlation between flux and period change is observed (Inam & Baykal 2000). So the conditions are probably different compared to the low mass X-ray binary Her X-1. OAO 1657-415 is probably better compared to other high mass X-ray binaries (HMXB), e.g. Cen X-3 and Vela X-1. The long-term frequency histories presented by Bildsten et al. (1997) for several pulsars, show a secular spin-up for Cen X-3 with a weak quasi periodic modulation of about 5 to 8 years. For Vela X-1, on the other hand, a secular spin-down is shown, but with a torque reversal around MJD 44 000 (spin-up to spin-down) and possibly also around MJD 49 000 (spin-down to spin-up). Between these dates also a quasi periodic characteristics of the period evolution with a period around 5 years might be inferred. We note here that these quasi-periodic pulse period variations, seen in several accreting X-ray pulsars, have not found much attention in the literature.
In addition to the long-term quasi-periodicity we have found
variations of the pulsar period on time scales shorter than the
orbital period (residuals in Fig. 7).
Bildsten et al. (1997) have discussed power spectra of torque
fluctuations for several objects, including OAO 1657-415. While Her X-1
and Vela X-1 showed a white noise power density spectrum, Cen X-3
shows red noise. For OAO 1657-415 Bildsten et al. (1997) also noted a red noise
spectrum. But a closer look at the power density spectrum shows a
slight increase in power density for frequencies above the orbital
frequency of about
Hz. OAO 1657-415 has the highest
power density at high frequencies of all pulsars presented in
Bildsten et al. (1997), which is consistent with our finding of
pulse period changes on time scales much shorter the orbital
period. Also Blondin et al. (1990) expect fluctuations in the
spin-up/spin-down on short time scales as a result of changes in
the sign and magnitude of the accretion torque.
The energy-resolved pulse profiles show a complex structure with one main peak followed by at least two somewhat weaker peaks. The spectrum of the rising part of the main peak is clearly harder than in the other parts of the pulse.
The change of spectral hardness with pulse phase is
commonly observed. Below 10 keV this behaviour might be modelled
by a variation of an absorbing column density
(e.g. Cen X-3, Santangelo et al. 1999a) but for energies above
10-20 keV cold absorption plays no role anymore. A mixture of
different emission regions, as proposed by Kraus et al. (1996) to
explain energy-resolved pulse profiles of Cen X-3, might explain
the observed profiles, but such detailed theoretical modelling is
beyond the scope of this paper.
As in many other accreting pulsars, e.g. 4U 0115+63 (Tsygankov et al. 2007) or GX 1+4 (Ferrigno et al. 2007), the pulsed fraction increases with energy. Tsygankov et al. (2007) point out that in the accretion column the region emitting harder photons is closer to the neutron star surface and therefore can be more obscured during certain pulse phases than the region higher up emitting photons of lower energy.
The broad band spectrum, composed from ISGRI and JEM-X spectra could be modelled by the HIGHECUT model. We could not find an indication for a cyclotron absorption in the phase-averaged spectrum.
Blondin et al. (1990) have done hydrodynamic simulations of the accretion by a neutron star from stellar wind in massive X-ray binaries. Their work points to episodic changes of the mass accretion rate, which results in corresponding changes of the X-ray luminosity on a time scale of hours. They also find a variation of integrated column density, which varies with orbital phase in such a way that it is highest near eclipse and lowest at orbital phase 0.5. This is in good agreement with our findings for the JEM-X1 and JEM-X2 spectra, which show a higher absorption near eclipse.
Phase-resolved spectra were produced from ISGRI observations for 4 phase intervals. These spectra show significant differences. Phase 0.025-0.25, covering the rising part of the main pulse peak, shows the hardest spectrum, while phase 0.75-1.025 has the softest spectrum. Phase 0.5-0.75 fits nearly exactly the phase-averaged spectrum, while phase 0.25-0.5 shows a feature at about 50 keV, which is reminiscent of a cyclotron absorption line, but with a significance too weak to prove its existence.
Apart from not being able to detect a cyclotron line in our data, it would not be unlikely to find such a line in the spectrum of phase interval 0.25-0.5, as it is common to several pulsars that cyclotron lines are visible only - or at least most prominent - in the falling part of the main pulse, e.g. 4U 0115+63 (Santangelo et al. 1999b), GX 301-2 (Kreykenbohm et al. 2004), Cen X-3 (Burderi et al. 2001), Vela X-1 (Kreykenbohm et al. 2002) and possibly also GX 1+4, which shows a likewise weak absorption feature at 34 keV (Ferrigno et al. 2007).
We have presented a spectral and timing analysis of INTEGRAL observations of OAO 1657-415, which span a time range of about 3.5 years. Though we found our results to be compatible with similar observations of other high mass X-ray binaries, we state that a consistent model for OAO 1657-415 is still missing. Such a model should explain the nature of the strong variability and its possibly weak correlation with the orbital phase. Also the origin of the fluctuations of the pulse period with its observed long-term quasi periodic behaviour should be addressed by this model. Detailed analysis is also required to explain the rather complex structure of the pulse profiles. Further observations with continuous coverage of OAO 1657-415 would berequired to settle the question of a potential spectral feature in the falling part of the main peak of the pulse profile.
Acknowledgements
We thank Mark Finger for complementing the publicly available data on pulse frequency. Part of this work is supported by the German Space Agency (DLR) under contracts 50 OG 9601 and 50 OG 0501. We also thank the anonymous referee for valuable comments and suggestions.