A&A 380, 516-519 (2001)
DOI: 10.1051/0004-6361:20011390
D. A. Leahy
Department of Physics, University of Calgary, Calgary, Canada T2N 1N4, Canada
Received 19 June 2001 / Accepted 1 October 2001
Abstract
An analysis of RXTE/ASM observations of the unusual binary system LSI +61303
is presented. The X-ray light curve gives a binary period of
days,
and an eccentricity of
(1
).
The eccentricity is determined by
fitting an inverse Compton model to the X-ray light curve.
These are the best determinations to date, as
the optical period and eccentricity have larger uncertainties and
published values of the radio period are variable.
The radio outburst peak is not phase locked to the binary period,
so the radio period is not the same as the binary period.
Key words: stars: neutron - stars: individual: LSI +61303 -
binaries: radio - stars: emission line, Be -
X-rays: stars
Optical data show that the primary is a rapidly-rotating star of
spectral type B0 V (Hutchings & Crampton 1981) or B0-B0.5 III (Paredes & Figueras 1986).
Spectroscopic radial velocity measurements in the optical and UV
confirm the binary nature of the system, and find an orbital period in
agreement with that of the radio outbursts (Hutchings & Crampton 1981). The optical
spectroscopic measurements of the system are however complicated by
the composite and variable nature of the line emission (Hutchings & Crampton 1981),
and have not yielded a good orbital solution. The best fit to these
data imply a high-eccentricity of 0.6-0.8 and a periastron passage
at radio phase of
,
although the data is also consistent
with zero eccentricity. Orbital solutions derived from infrared
photometric observations also imply high eccentricities (0.7-0.8),
with a best-fit periastron passage at
0.5 (Marti & Paredes 1995). For
this orbit and a Be star mass of 10
the semimajor axis of the
orbit is 4-6
cm. Infrared measurements imply the
presence of a dense equatorial disk, and a mass-loss rate in the wind
of 1-4
yr-1 has been estimated from
these observations, assuming an opening angle of 15
for the disk
and an initial wind velocity of 5 km s-1 (Waters et al. 1988).
Originally detected by Einstein, LSI
is a weak X-ray source,
with
ergss-1 (Bignami et al. 1981). ROSAT
observations determined that the X-ray source is variable by a factor
3 on timescales of a day, and has a harder spectrum than that
characteristic of O and B stars (Goldoni & Mereghetti 1995). Soft X-ray monitoring of
LSI +61
303 with ROSAT at nine intervals during a single orbit showed the X-ray
emission to be variable by a factor
10 over the timescale of an
orbit (Taylor et al. 1996). Simultaneous radio monitoring indicated that
the X-ray and radio fluxes may be anti-correlated, which was confirmed
by Harrison et al. 2000.
observations (Leahy et al. 1997), show that the spectrum is
described by a relatively hard power-law, and therefore the emission
mechanism is non-thermal in nature.
LSI +61303 is particularly interesting due to the fact that
it is located within the error box of the strong 100 MeV COS B
source 2CG 135+01 (Bignami & Hermsen 1983). Recent EGRET observations have
narrowed the error region, and currently LSI +61
303 is the most probable
X-ray counterpart (von Montigny et al. 1993).
The luminosity in the EGRET band is
1035 (d/2 kpc)2 erg s-1.
COMPTEL has also detected
emission in the 1-10 MeV band from the same region (van Dijk et al. 1994)
but the error circle also contains the quasar 4U 0241+61.
RXTE/PCA and RXTE/HEXTE observations are presented by Harrison et al. (2000). These
resolve LSI +61
303 from the QSO and show that COMPTEL flux is dominated by
the QSO and not LSI +61
303.
No significant variability has been detected
in the high-energy data, however the observation periods are typically
two weeks and thus sample a large fraction of the binary orbit.
In this paper, RXTE/ASM observations are analysed to give the orbital period, eccentricity and phase of periastron passage.
The RXTE/ASM dwell data and daily-average data were obtained from the
ASM web site (http://xte. mit.edu/ASM_lc.html). The data reduction to
obtain the count rates and errors from the satellite observations was
carried out by ASM/RXTE team, and the procedures are described at
the web site.
The ASM count rates are for the 1.5-12 keV keV band. The data covered
the time period MJD 50088.1 to MJD 51969.2, or about 1881 days.
Most of the RXTE/ASM data points are only upper limits to the source
flux, with data values a factor of several less than the 1
errors.
However, when averaging over many points there is a flux
significantly above zero. For example, averaging the data into 26.4 day bins
yields 71 points: the mean ratio of count rate to 1
error is 3.07,
so that the source is well detected. In the next two subsections, it is
found that the statistical significance of the periodicity found in
power spectrum and epoch folding analysis is high. This is due to the large
amount of data, even though each individual data point is of very low
significance.
The power spectrum was created by averaging the dwell data into
0.05 day bins, subtracting the mean value of all
inhabited bins from each inhabited bin, and taking an FFT.
The dwell data were used in preference to the daily average data
to avoid the effect of averaging over 1day intervals.
Since the interest was only in periods greater than 1-2 days,
and the sampling of the ASM data is poor on short timescales,
the power spectrum was made with 4096 frequencies,
/day
to 1.25/day.
The resulting power at each frequency was divided by the average
power over the whole frequency range.
The
distribution was used to assess the significance of
any peaks in the power spectrum.
The most significant period was at 26.214 days, with a ratio of power/(mean power) of 12.09.
This is statistically significant:
the chance probability is 0.022 in 4096 trials.
Also, the period is nearly the same as published values of the radio period.
The period with next largest power was the adjacent period at 26.426 days.
It had a ratio of power/(mean power) of 8.7, giving a chance probability of
0.68 in 4096 trials.
This is not highly significant by itself, but with the above result confirms
that there is a real periodicity in the data with a period between
26.2 and 26.4 days
An epoch folding analysis was carried out, with the statistic used to assess the reality of variability at any given
trial period. See Leahy et al. (1983) for a description of epoch folding and
the assessment of period uncertainty.
The daily-average data were folded into 24 phase bins in
the period range of 20 days to 30 days, giving a maximum
of 66.6 and a
chance probability of
at a period of 26.448 days.
The dwell data were also epoch folded into 24 bins giving a maximum
of 103.6 and a
chance probability of
at a period of 26.415 days.
The epoch folding results from the dwell data
are more significant than from the daily-average data.
This is not surprising, the lower significance for the daily average data
is probably due to the effect of averaging over 1 day intervals.
Other numbers of phase bins were also tried, since too many phase bins means too
weak a signal in each bin, and too few phase bins results in smoothing the
light curve.
The optimum number of bins was determined by the maximum significance.
Maximum
is not the correct discriminator, as the number of degrees
of freedom (equal to number of bins less one) has a strong effect on the
probability distribution and is different for different numbers of
bins.
10 bins gave a maximum
of 73.9 at 26.410 days and probability of
;
16 bins gave a maximum
of 93.4 at 26.415 days and probability of
;
32 bins gave a maximum
of 101.2 at 26.448 days and probability of
,
Thus the maximum significance is obtained for folding into 16 phase bins.
The
vs. period plot is given in Fig. 1 for 16 bins and two cases.
A wide period range was
folded at periods separated by 1/4 of the Nyquist frequency, and
a narrow period range was folded
at periods separated by 1/32 of the Nyquist frequency
(the latter is offset by +80 in
to avoid confusion of the two lines).
The folded light curve at 26.415 days, the period of maximum
,
is shown in Fig. 2. The solid line is a model fit described below.
The period uncertainty for the dwell data is
days.
The dwell data were split into 6 adjacent sections and epoch folded, to
test for period variability.
There was 1 section of 269 days (MJD 50088.1 to 50357.1,
chosen to be coincident with the time period used by
Paredes et al. 1997), 4 sections
of 300 days (MJD 50357.1 to 51557.1), and 1 section of 312 days
(51557.1 to 51969.2).
The period uncertainty was 1.3 days: this is much larger than for the
full data set due to both the weaker signal and shorter data span.
All of the sections had maximum
from 35.2 to 52.6
and periods consistent with 26.4 days.
Thus there is no evidence for period variability over the whole time period.
![]() |
Figure 1:
Chisquared vs. period for the RXTE/ASM dwell data set folded
into 16 phase bins. A wide period range folded at 1/4 Nyquist period spacing
and a narrow period range folded at 1/32 Nyquist period spacing are
shown (the latter has 80 added to ![]() |
Open with DEXTER |
![]() |
Figure 2: Light curve at 26.415 days for the RXTE/ASM dwell data set folded into 16 phase bins. The solid line is the model fit described in the text. |
Open with DEXTER |
The RXTE/ASM data consist the best long-term X-ray observations of
LSI+61303. The period determined here of 26.42 days is the best
determination of the orbital period of this system. The uncertainty is
dominated by the limit of the time span of observations and is
day. The radio period
shows evidence of being variable (e.g. see Ray et al. 1997), which is likely
due to the radio outbursts not being tightly coupled to the orbital
period. The optical period is the orbital period, but
the only determination (Hutchings & Crampton 1981) has an uncertainty of
days.
The X-ray light curve is produced by inverse Compton emission of stellar photons on the relativistic electrons responsible for the radio synchrotron emission (e.g. see Leahy et al. 1997). The peak of the X-ray emission should be consistently at periastron, unlike the radio emission peak which varies in phase. Thus the peak of X-ray emission and periastron occurs near phase 0.9. For the epoch folding analysis above, the reference for phase zero (T0x) was taken as MJD 50088.1.
To determine the phase of X-ray peak more accurately and also to
estimate the eccentricity
the inverse Compton model was fit to the light curve of Fig. 2. The best fit
model is shown by the solid line. The
resulting best fit eccentricity is 0.30 with a 1
range of
0.26-0.34 and a 2
range of 0.22-0.38.
The eccentricity is consistent with the range
0.2-0.8 allowed by the optical observations
(Hutchings & Crampton 1981).
The best fit phase of periastron passage is 0.916
with a 1
range of 0.90-0.93, and a 2
range of 0.89-0.95.
This yields the epoch of periastron passage of MJD
(2
).
The average radio lightcurve was used to define phase 0 as peak radio flux in
Ray et al. 1997 to yield a reference epoch of JD 2449393.5(0.4).
Taylor & Gregory (1984) defined phase 0.6 as peak radio flux to yield a reference
epoch of JD 2443366.775.
The resulting phase of radio outburst for the two different radio epochs
are at phase, with reference to T0x,
and
.
The uncertainty is primarily due to the uncertainty in the X-ray period,
but the difference is due to the different periods and epochs published
for radio outburst.
So rather than using mean period and epoch for radio outburst, it is preferable to consider cases where radio and X-ray have been observed together.
Two orbital
cycles which have been sampled by observations with ROSAT and RXTE simultaneous
with radio monitoring (Harrison et al. 2000).
These both have the radio outburst occuring
0.4 in orbital phase after
X-ray peak.
The basic nature of the system LSI +61303 is now thought to be as follows:
electrons are accelerated
to relativistic energies probably by a pulsar in orbit around the Be star
primary star.
The acceleration mechanism could be electromagnetic fields in
the vicinity of the pulsar or shock acceleration at the boundary of the pulsar
outflow and the stellar wind from the primary.
The flux of stellar photons
interacts with the power law distribution of relativistic electrons to produce
a broad power law distribution of photons from X-ray through
-ray
energies.
Due to the magnetic field in the vicinity of the pulsar, the same
electrons radiate synchrotron photons in the 1 GHz to >100 GHz range.
The eccentric orbit of the pulsar and its associated
relativistic electron population result
in the moderate modulation of the X-ray through
-ray flux.
The radio outburst is much more highly variable than the X-ray and is not phase locked with the orbit.
The radio outburst follows the X-ray peak by
0.4 in orbital phase,
so that the radio outburst occurs approximately near apastron.
The radio outburst emission has been successfully described as due to a rapidly
expanding cloud of relativistic electrons (Peracaula 1997).
Further work including a detailed model of the radio through -ray emission,
with several new features beyond that presented in Leahy et al. (1997), is presented
in a following paper.
Acknowledgements
This paper made uses of the quick-look results provided by the ASM/RXTE team. DAL acknowledges support from the Natural Sciences and Engineering Research Council.