A&A 401, 1111-1117 (2003)
DOI: 10.1051/0004-6361:20030201

Hard X-ray behaviour of Her X-1 during second main-on after the anomalous low state

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 $\sim$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


1 Introduction

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$^{\rm s}$ pulsation, 1.7$^{\rm d}$ 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$^{\rm d}$ 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 $\alpha \sim 1$ between 2-20 keV and a sharp cut-off beyond showing a power law fit with $\alpha \sim 4.5$ or a thermal temperature of $kT\sim 5{-}7\:$ keV (Mihara et al. 1991). However, in the broad band observations up to 100 keV, an additional power law component with $\alpha \sim 1.5$ 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 ($\sim$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 $M_\odot$ 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 $35^{\rm d}$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 $.8^{\rm d}/1.6^{\rm d}$ 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 $\sim$10 days during the main-on state and for $\sim$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 $n_{\rm H}$ value which increased to $\sim$ $5\times \rm 10^{23}~\rm cm^2$, 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 ${\dot{P}\over P}\sim 3\times 10^{-6} ~\rm yr^{-1}$. 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.

2 Observations

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 $4.5^\circ \times 4.5^\circ$ 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$\sigma$ sensitivity of the LASE telescope in the entire energy range up to 200 keV is $\sim$ $1 \times 10^{-6}~\rm
cm^{-2}~s^{-1}~keV^{-1}$ for a source observation of 104 s.

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{aah3843f1.eps}\end{figure} 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 $\phi_{1.7}=0.04$ and $\phi_{35}=0.43$ (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 $\phi_{35}\sim 0.23$. The present observations correspond to MJD 51868.28.

3 Results

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 $\sim$$35\sigma$.

In the LASE payload, a 25 $\mu$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 $\chi^2$ deviation of the phase-bins from the phase averaged value. In the present data the best fit period corresponds to 1237.788 $\pm$ 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, $\alpha$) 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. $\sim$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.

 

 
Table 1: Best fit parameters for LASE data.
Model bknpower bremss+power
continuum    
index 1 3.85 $\pm$ 2.4 2.8 $\pm$ 6.3 keV
index 2 1.33 $\pm$ 0.14 1.33 $\pm$ 0.2
cyclotron line    
E1 (abs) 34.24 $\pm$ 2.42 34.82 $\pm$ 3.33
Ew1 4.7 $\pm$ 2.3 4.96 $\pm$ 2.5
Ed1 1.44 $\pm$ 0.91 1.29 $\pm$ 0.79
second harmonics    
$E_2=2\times$E1    
Ew2 12.3 $\pm$ 3.2 13.1 $\pm$ 3.6
Ed2 0.1 $\pm$ 0.11 0.11 $\pm$ 0.12
$\chi^2 (\rm d.o.f.) $ 2.48 (16) 2.68 (16)



  \begin{figure}
\par\vskip 115 mm
\special{psfile=h3843f2.ps vscale=70 hscale=53 voffset=400 hoffset=-70 angle=-90}
\end{figure} 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.

  \begin{figure}
\par\vskip 115 mm
\special{psfile=h3843f3.ps vscale=70 hscale=55 voffset=400 hoffset=-80 angle=-90}
\end{figure} Figure 3: Composite spectral fits to the LASE and PCA data.
Open with DEXTER


 

 
Table 2: Best fit spectral components for PCA+LASE data.
Warm absorber $n_{\rm H} \sim 4.6\pm 2.4 \times 10^{22} $
Iron line 6.57 keV , $\sigma \sim 0.56$
Power law Index $\alpha_1$ $0.94\pm 0.36$
K1 $7.4\pm 2.8 \times 10^{-2}$
Cut off $ \rm 16.6 \pm 1.1 $ keV
Folding energy $89.5 \pm 8.8$ keV
Power law Index $\alpha_2$ $0.75\pm 0.35$
K2 $6.8\pm 2.1 \times 10^{-2}$
cyclotron line  
E1 (abs) 33.4 $\pm$ 1.8
Ew1 10.5 $\pm$ 0.5
Ed1 1.84 $\pm$ 0.06
second harmonics  
$E_2=2\times$E1  
Ew2 30.2 $\pm$ 4.0
Ed2 0.62 $\pm$ 0.06
$\chi^2 (\rm d.o.f.) $ 1.68 (64)



  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{aah3843f4.eps}\end{figure} Figure 4: Comparison of the pre and post ALS hard X-ray observations of Her X-1.
Open with DEXTER

It is seen from the Table that best fit value of $n_{\rm H}\sim 5 \times 10^{22}$ derived for the combined fit is mid way between the zero value for the normal state and the $n_{\rm H}\sim 4\times 10^{23}$ inferred for the low and anomalous low states (Coburn et al. 2001). To illustrate the hard X-ray behaviour in Fig. 4 we plotted a comparison of the present observations with the hard X-ray measurements made during 1997 by the same instrument. PDS data from BeppoSAX in the 2-30 keV band corresponding to the first main-on is also plotted for comparison. It is clearly seen that not only the low energy data are consistent but the hard X-ray spectra above 50 keV are also similar. The only main difference being the energy of the cyclotron absorption feature which occurred at 47.8 $\pm$ 2 keV in the 1997 data (Manchanda 2001).

4 Discussion

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. $ \omega_{\rm M} > \omega_{\rm Kep}$ 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 $\sim$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 $\sim$ $3\times 10^{12}$ 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 $n_{\rm H}$ 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.


  \begin{figure}
\par\includegraphics[width=16.5cm,clip]{aah3843f5.eps}\end{figure} 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  $\rho\: u_{\rm p}\: \Omega\: r^2$, the ordered magnetic energy density in the disk  ${{{r\: B_\phi\: B_{\rm p}}\over 4\pi}}$and the viscous torque $\tau_{{\rm visc}}$. For a dipole field, the $B_\phi$ 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 $60^\circ$ 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 $\sim$ $\Delta\alpha \times 2\pi \sin\chi$, where $\Delta\alpha$ and $\chi$ 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.

References

 
Copyright ESO 2003