A&A 401, 1111-1117 (2003)
DOI: 10.1051/0004-6361:20030201
R. K. Manchanda
Tata Institute of Fundamental Research, Mumbai, India
Received 16 July 2002 / Accepted 13 February 2003
Abstract
The binary X-ray pulsar Her X-1 shows a variety of time
variations both in luminosity and its pulse features. For the first
time an about 600 days long anomalous low state was seen in the light curve
of ASM on board RXTE during the period MJD 51200-51800.
We observed the source in the 20-200 keV energy band
during its second main-on in the post anomalous state on
November 21, 2000. The source was found to be extremely bright
at high energies and a total of 26 000 photons from the source
were detected during the 35 min snapshot observation. The cyclotron absorption line at
34 keV
and its second harmonics were clearly detected in the data. We present the
spectral details and discuss its implication
in terms of the geometrical properties of the source.
Key words: accretion, accretion disks - X-rays: stars - pulsars: individual - Her X-1
The binary X-ray pulsar Her X-1 is a
template system for the low mass X-ray binaries group.
The source exhibits a variety of temporal variations like
the 1.24
pulsation, 1.7
binary period, 35 day cyclic variation of high
and low fluxes during which source remains on for about 5 binary cycles,
Mid-on emission during the 35 d cycle, and the possible 1.6
periodicity
in the intensity dips seen clearly from the source outside the
eclipse phase (Tananbaum et al. 1972; Deeter et al. 1981).
The X-ray spectrum of Her X-1 in the low energy region below 30 keV
is characterized by a power law with
index
between 2-20 keV and a sharp cut-off beyond
showing a power law fit with
or a thermal temperature
of
keV (Mihara et al. 1991). However, in the broad band observations
up to 100 keV, an additional power law component with
is
required to best fit the data (Manchanda et al. 1973, 1984; Soong et al. 1990).
A prominent cyclotron absorption line feature between 35-55 keV distinguishes
its spectrum from the other known LMXB pulsating X-ray sources.
The high energy line
feature is interpreted as arising due to cyclotron absorption
in the intense magnetic field region (
1012 G) near the magnetic
poles of the rotating neutron star (Truemper et al. 1978). A second
cyclotron absorption feature corresponding to the second harmonics has been
observed on occasions (Manchanda 2001). However,
there is no common consensus as yet on the
line width and the centroid energy. The data suggest
an apparent correlation with the binary phase, 35 day cyclic phase as well as the source
luminosity (Staubert et al. 1983; Voges et al. 1982; Tueller et al. 1984; Mihara et al. 1991, 1998;
Gruber et al. 1998; Soong et al. 1990; Dal Fiume et al. 2000).
A binary system containing a 1.3
neutron star
with its optical companion HZ Her fits quite well the observations both in
optical and in the
X-ray band (Doxsey et al. 1973). Pulse features and the bulk spectral properties of the
source can be explained by the mass accretion on to the neutron star from the
accretion disk, which in turn is fed by the mass flow from the companion star.
However, there is no universally accepted model for explaining the
on-off cycle of the source. A variety of mechanisms have been
proposed in literature (Roberts 1974; Petterson 1975; Gerend & Boynton 1976; Priedhorsky & Holt 1987; van Kerkwijk et al. 1998). The poor quality
of the 35 day clock leads to some of the difficulties in invoking precise
mechanisms (Oegelman 1987). The origin of the fourth
clock
is still unknown (Crosa & Boynton 1980; Truemper et al. 1986; White et al. 1995).
Since its discovery, a regular 35 day on-off cycle in the source intensity has been
the key characteristic of the source, and hence its geometrical models.
The source intensity remains high for 10 days during the main-on state and
for
5 days during the mid-on state half way between the cyclic off state. The
cycle regularity is clearly visible in the light curve of the source in the
BATSE archival data since 1991 in the hard X-ray energy band of 20-100 keV and in the ASM data since 1994 in the soft X-ray band.
A low intensity behaviour of Her X-1 at low energy was first observed in 1983
from the EXOSAT data during June to August (Parmar et al. 1985). No contemporary data is available
in the hard X-ray band. However, the non detection of
the source by ASCA in August 1993 in the soft X-ray band (Mihara & Soong 1994)
corresponded to only a reduced activity of the source in the hard X-ray band as
is seen from the BATSE data.
The source entered an anomalous low intensity state in March 1999
which persisted for about 18 months.
As seen in the BATSE and the RXTE light curves, this recent anomalous phase was
extended over the entire X-ray energy band of 0.1 to 200 keV or possibly beyond. The ASM and BATSE
light curve show a gradual drop in the X-ray luminosity of the sources until the anomalous
low state (ALS) where the normal turn-on did not occur.
Optical spectra of the HZ Her taken during ALS showed prominent emission lines
features of CIII/NIII and He II at comparable strength to the normal main-on observations
(Margon et al. 1999). Since these emission line features are produced due to
Bowen pumping caused by X-ray irradiation of the companion star, it firmly
establishes that even though no X-ray flux was detected from the source,
the X-ray photon production in the emitting regions of the source continued
at the same level as in the normal main-on state. The spectral data taken before the source
ingress to the ALS, show that only significant variation was observed in the
value which increased to
,
all other parameters
remained unchanged within the statistical uncertainties. The pulse period near to the ALS was
measured to be 1237.77485 msec. A comparison of pulse period obtained during the preceding
normal state, indicated that the source had switched to a rapid spin-down mode
compared to its normal spin-up behaviour. (Parmar et al. 1999; Levine & Corbet 1999; Coburn et al. 2000).
Using BeppoSAX observations during the first main-on after the anomalous low state,
Oosterbroek et al. (2001) have
reported the evidence that the rapid spin down of the source continued during
the entire ALS period and the source pulse period changed to the value observed in 1984/5.
The measured pulse period of 1237.7697 corresponds to a
.
The pulse profiles and the spectral behaviour of the source in the 0.1 to 30 keV energy
band was found to be consistent with the earlier observations during high state.
In this paper, we present the hard X-ray observations in the extended energy region of 20-200 keV made in the middle of second main-on phase as seen in ASM light curve after the anomalous low state. This observation was planned for the mid-on phase in the 35 day cycle based on the old ephemeris. The source was found to be quite bright at hard X-ray energies. The observed time averaged spectrum is hard with a clear indication of cyclotron features. We have also used the RXTE archival data form the PCA detectors taken on the same day to constrain the model parameters.
The X-ray observations were carried out using LASE payload at a ceiling altitude of 2.5 g cm-2 during a balloon flight conducted on November 21, 2000 from Hyderabad, India. The newly developed balloon-borne payload code named LASE consists of an X-ray telescope made of three identical modules of specially developed large area scintillation detectors having both passive and active shielding and fitted on a fully steerable alt-azimuth mount. One of the primary goal of the LASE experiment is to look for periodic and quasi-periodic components in the hard X-ray region from LMXB and AGN sources in the hard X-ray energy region up to 200 keV (Manchanda 1998).
Each detector module is a stand alone unit and consists of
a combination of thin and thick sodium iodide
detector arranged in a back-to-back configuration and
has a geometric area of 400 cm2.
The field of view of each module is
and is
defined by a slat collimator made from specially designed
sandwiched material of Lead, Tin and Copper. The response matrix or
each module is generated by using the preflight calibration data at
different energies. Each module is also provided with a Am241 calibration source
for gain calibration during the flight using ground command. The details of the
detector design, associated electronics and the control system are given
elsewhere (D'Silva et al. 1998). A 3
sensitivity of the LASE
telescope in the entire energy range up to 200 keV is
for a source observation of 104 s.
![]() |
Figure 1: X-ray light curve of Her X-1. The arrow in the figure represents the epoch of present observation. |
Open with DEXTER |
Her X-1 was observed for a continuous stretch of 35 min between 0630 and 0705 UT on November 21st and the source background was measured for 20 min each
before and after the observation. The off-source pointing location was
carefully selected blank field from the known X-ray source catalog.
Present observations correspond to binary phase
and
(old ephemeris). Using the X-ray light curve
from the all sky monitor on-board RXTE for the
35-day turn-on corresponding to the present observations, the derived
phase is about
.
The present observations correspond to MJD 51868.28.
The present observations were timed to make a snap-shot spectral measurement
of the source near the short mid-on in the 35 days cycle. In Fig. 1, we have plotted
an expanded view of the X-ray light curve of the source taken from the all sky monitor
on board RXTE. The arrow in the figure indicates the epoch of present observations.
It is clearly seen from the figure that present data correspond to the main-on due to a
delayed onset of the source as noted by Oosterbroek et al. (2001). It is also seen from the
figure that in the post-anomalous state, the second main-on represents the normal behaviour
during which the source intensity remained high for about 10 days,
while the first main-on was of much shorter duration. A total excess of 26 000 photons
due to Her X-1 were recorded in the three detectors during the present observations.
This corresponds to a combined statistical significance of .
In the LASE payload, a 25 sec temporal resolution for each accepted event
is achieved by referencing the arrival time with the telemetry frame start pulse.
To determine the pulse period, we computed a table of arrival times for all photons in
the entire energy range. After making the barycentric conversion we binned the data and
searched for a
coherent pulsations around 1237.7 msec using XRONOS folding routine which maximizes the
deviation of the phase-bins from the phase averaged value.
In the present data the best fit period corresponds to 1237.788
0.01 msec
and is quite consistent with the value of 1237.7697 msec, observed during first mid-on in the
post-anomalous state of the source. The present data
supports the continued rapid spin-down of Her X-1 during the ALS as inferred
from the BeppoSAX data (Oosterbroek et al. 2001).
To derive the time averaged source spectrum, the source contribution was binned in 20 energy channels between 19.5 to 200 keV band. Each bin in the grouped data contained more than 600 source photons. The ASCII data files were then converted to
FITS format and analyzed using XSPEC spectral analysis model fitting routine. The FITS files of
the detector response matrix, atmospheric transmission and the detector resolution function
constitute the other key input parameters during the spectral fitting.
A broken power law as well as a composite function
with thermal and non thermal terms represent the continuum data equally well. Therefore,
the model fits were made using both the broken power law and the composite model (kT, )
for the source continuum along with the cyclotron absorption features.
The data folded with best fit broken power law model is shown
in Fig. 2. The cyclotron absorption features were fitted with the Mihara's cyclotron
function as provided in the XSPEC analysis tools. The del-chi in figure indicates the quality
of the fit to the data and the improvement in the chi-square value of the fit was 55% with the
cyclotron absorption lines.
The best fit values of various
parameters for both models are summarized in Table 1.
The data presented in Fig. 1, clearly demonstrates that present
observations corresponded to the main-on phase rather than the near mid-on
expected on the basis of mean 35 day ephemeris. From a detailed analysis of
the on-state data since 1996, Oosterbroek (2001) have derived that
main-on in the post anomalous state are delayed by 7 days i.e.
0.2 in phase.
We can therefore, safely assume that the delay in the on-set was caused by the anomalous
low state of the source. However,
such large delay is very difficult to account in any reasonable geometrical
model of the source. The measured pulse period
in our data is consistent with the BeppoSAX measurements of rapid spin down
of the source during ALS. Detection of the source with high
significance during our observations clearly suggests that Her X-1 had regained
all of its characteristics by the second main-on. The spectral data plotted in Fig. 2 show not only that the source was quite bright at energies
above 30 keV, but also that the cyclotron absorption features are clearly
seen in the data.
Model | bknpower | bremss+power |
continuum | ||
index 1 | 3.85 ![]() |
2.8 ![]() |
index 2 | 1.33 ![]() |
1.33 ![]() |
cyclotron line | ||
E1 (abs) | 34.24 ![]() |
34.82 ![]() |
Ew1 | 4.7 ![]() |
4.96 ![]() |
Ed1 | 1.44 ![]() |
1.29 ![]() |
second harmonics | ||
![]() |
||
Ew2 | 12.3 ![]() |
13.1 ![]() |
Ed2 | 0.1 ![]() |
0.11 ![]() |
![]() |
2.48 (16) | 2.68 (16) |
![]() |
Figure 2: Spectral fit for the LASE data. See Table 1 for details of various spectral components. |
Open with DEXTER |
In order to further constrain the spectral parameters in the low energy region we have used
archival data form PCA observation
on board RXTE made on Nov. 21, 2000 at 0341-0553 UT. Fortunately, this observation was made only 4 hour before our observation in hard X-ray band and therefore, must represent a
very high degree of correlation to the high energy data. Figure 3 shows the
PCA and LASE data and the best fit value of various components. Since the PCA data above
40 keV have large uncertainties even for bright sources, we have not utilized the data from higher energies. The model
component used for the combined fit are a leaky warm absorber similar to the one proposed from the ASCA
data (Endo et al. 2000) and consists of a Gaussian iron line, power law with a high energy cut-off, a second power law component
and the cyclotron absorption line features as defined by the Mihara function. The values of different best-fit
parameters are given in Table 2. The addition of the second harmonic improved the value chi-square from 3.98 to 1.68 per degree of freedom. The most important fact to note is that no arbitrary scaling
or additive terms of any kind were necessary to make a combined fit to the data from
the two instruments. This clearly
points to the very well understood response matrix of the LASE payload.
![]() |
Figure 3: Composite spectral fits to the LASE and PCA data. |
Open with DEXTER |
Warm absorber |
![]() |
Iron line | 6.57 keV ,
![]() |
Power law Index ![]() |
![]() |
K1 |
![]() |
Cut off |
![]() |
Folding energy |
![]() |
Power law Index ![]() |
![]() |
K2 |
![]() |
cyclotron line | |
E1 (abs) | 33.4 ![]() |
Ew1 | 10.5 ![]() |
Ed1 | 1.84 ![]() |
second harmonics | |
![]() |
|
Ew2 | 30.2 ![]() |
Ed2 | 0.62 ![]() |
![]() |
1.68 (64) |
![]() |
Figure 4: Comparison of the pre and post ALS hard X-ray observations of Her X-1. |
Open with DEXTER |
In the standard picture of accretion in close X-ray
binary systems, the all important transfer of mass from the
primary to the compact companion takes place by Roche lobe overflow
or by a stellar wind and secular changes in the spin period
are associated with the magnitude of the accretion torque (Ghosh & Lamb 1979).
For the 35 day on-off cycle, the proposed geometries include, precession
of a slave-disk, warp disk etc. (Roberts 1974; Petterson 1975; Pringle 1996; Shakura et al. 1999).
A reduced accretion and a reduction in the inclination of the accretion disk has been proposed
to be the cause for the ALS by Coburn et al. (2000). The reduced accretion leads to the
decrease in the spin-up rate or even spin down in extreme circumstance, however,
such a model does
not account for the delay in the 35 day cycle which simply amounts to resetting the
clock at an arbitrary epoch. It is to be noted that the effective change in the pulsar period
consists of the sum of its normal spin up rate and the observed magnitude of the spin down rate.
The rapid
spin down as observed Her X-1, therefore, not only requires reduction in the on-going mass accretion
but also an efficient pathway to loose angular momentum. The propeller effect is
commonly assumed as the likely cause for reduced accretion
(Illarionov & Sunyaev 1975; Davies et al. 1979). However, this condition
requires that the angular velocity of the star be larger than the Keplerian value of the
angular frequency
i.e.
such that spinning neutron star makes a centrifugal barrier and thereby
mechanically expels the material. For such a condition to occur, one requires either a very fast spinning
pulsar or extremely high magnetic field of
1015 G.
On the other hand the heated wind model for
mass loss requires higher temperature of the disk, which in turn implies a larger X-ray luminosity of the source (Illarionov & Kompaneets 1990).
The observed X-ray and optical data clearly limits three vital parameters in the case of Her X-1.
First, the magnetic field of
G as inferred from the cyclotron absorption line,
second, the X-ray production remained at the normal level during the ALS
phase and third, the spin down of the source was rapid. A spin down
model based on torque reversal due to
warping of the accretion disk has been explored by van Kirkwijk et al. (1998). Apart from a large
number of implausible assumptions necessary to explain the observed data,
authors note that once the instability which leads to warping
of the inner part of the disk is set, it is difficult to re-normalize the system.
As mentioned earlier, about the large value
of
observed during the ALS, it is tempting to relate the anomalous
low state with an extended normal off-mode behaviour during the 35 day cycle. However,
the observation of the phase difference of the main-on in the post ALS state as shown in the Fig. 5, based on the old ephemeris remains unresolved. In addition, though the non-detection of
the X-ray flux from the source can be explained by obscuration in
this model, the observation of rapid spin down of the source during the ALS would require an
independent set of assumptions.
In Fig. 5 we have plotted the long term light curve of the source in the pre and post ALS
using ASM data. The dotted lines in the figure represent the expected on-set of the 35 day cycle.
We have initialized the data visually at MJD 50807 corresponding to start of main-on. It is to be
noted from the figure that unlike earlier observation of fluctuations on the mean ephemeris for 35 day period (Oegelman 1987), the present phase shift appears to be of permanent nature to date.
Thus in the absence of any consistent geometrical model to explain both the spectral and temporal behaviour of the source, we propose that the emission of a strong jet from the Her X-1 led to the observed anomalous state of the source by causing an increased period of obscuration in the line of sight for about 18 months. A magneto-centrifugal wind has been discussed as the underlying mechanism for the spurting jets in which the ionized particles within the accretion disk are flung out along magnetic field lines that ultimately twist into a helical shape along the rotational axes of the disk. Orbital energy from the accretion disk is dissipated in propelling the polar jets (Königl 1989). Detailed models for the origin and self collimation of such jets from the surface of Keplerian accretion disks surrounding the central objects like YSO and black-holes have been discussed in literature (Pelletier & Pudritz 1992). Similarly, the ejection of moderately relativistic jets powered by rotational energy has also been proposed based upon the observations of superluminal jets from micro-quasars (Zhang et al. 1997; King & Kolb 1999). Such escaping jets not only lead to a large loss of gravitating mass but also efficiently extract specific angular momentum from the system i.e. the disk as well as the central object. Observationally, episodes of jet ejection have been observed for GRS 1915+105 (Mirabel & Rodriguez 1994), SN 1987a (Manchanda et al. 1989) and Cyg X-3 (Mioduszewski et al. 1998). In addition, pre-flare quenching in hard X-ray band of 20-100 keV has been noted for the large flare (Hjellming & Han 1995). Therefore, an episode of jet ejection, will affect the geometry of the precessing accretion disk and also force a rapid spin down in the source. The escaping material will disrupt the inner accretion disk structure, which is dynamically stabilized by the precessing and the motion of the accreted material, thus leading to not only obscuration of the source but also it will cause a shift in the 35 day ephemeris, since the amplitude of the perturbation will die down at an arbitrary time scale independent of the established temporal behaviour of the source. It is seen from Fig. 5 that the phase shift of about 7 days seen during the first main-on (Oosterbroek et al. 2001) is a persistent feature of the source and is consistent with the conclusion of re-establishing the stability of the system randomly.
![]() |
Figure 5: Long term light curve of Her X-1. The dotted line represent the expected turn-on dates based on the old ephemeris. |
Open with DEXTER |
In the present model we need to address the question, as to what triggers such transient
spurting of the jet. The
steady state angular momentum equation of an axisymmetric accretion disk threaded by the magnetic field
and undergoing viscous stress depends upon the angular momentum of the accreting matter
,
the ordered magnetic energy density in the disk
and the viscous torque
.
For a dipole field, the
is taken as Bz.
We propose that poloidal magnetic field in the accretion disk grows
on arbitrary time scale depending upon the accretion rate and the development of instabilities in the
disk. When sufficiently strong, this component can then drive the ejection of a jet.
It has been shown that poloidal magnetic field may be dragged radially inward by the accreting gas until its dissipative escape, driven by the gradient in magnetic pressure, limits its growth. Provided that the magnetic field makes an angle of less than
with the radius vector at the disk in Keplerian motion, it will be energetically favorable for the gas to leave the disk in a centrifugally driven wind (Blandford 1989; Königl & Ruden 1992).
A large inclination of the magnetic axis with respect to the spin axis in the case of Her X-1
further helps the obscuration scenario for the anomalous low state,
set off by the disruption of the precessing disk due to
the escape of the material in the form of a jet.
The high inclination of the magnetic axis with respect to the spin axis in Her X-1 can be inferred from the following.
The low energy pulse profile of Her X-1 shows a double peaked main pulse and a
extremely weak inter-pulse but at higher energies it is a the sharp peaked pulse profile.
This clearly suggests that orientation of the spin axis is such that
pulsed X-rays are emitted from the funnel of one of the magnetic pole having
an hollow cone structure. Since the duty cycle of the
pulse is approximately given by
,
where
and
are the beaming angle and the inclination of the magnetic pole
to the rotation axis, a large duty cycle of the pulse clearly require large
angle of inclination. In fact even in the normal main-on mode, the possibility that the
X-ray emission cone almost grazes the disk has been inferred from
the possible scattering of the beam necessary to account for the slight
asymmetry in the pulse profile as well as
explaining the 6.7 keV iron-line seen in the Her X-1 spectrum (Bai 1980; Mihara et al. 1991).
Therefore, a perturbation of the precessing disk can easily lead to
the observed flux anomaly from the source.
The present model can also account for the fluctuations in the mean ephemeris for the 35 day period observed
during the pre ALS data, if we assume that
a low intensity variable magneto-centrifugal wind from the accretion disk is a normal feature of the
source at all times.
Acknowledgements
It is pleasure to thank sarvashri J. A. R. D'Silva, P. P. Madhwani, N. V. Bhagat, B. G. Baggage and Ms. N. Kamble for mechanical and electronic fabrication of the payload. I also thank Mr. V. Rana and Dr. B. Paul for their help in reducing the RXTE data. RXTE teams are gratefully acknowledged for the archival data.