1 - National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, PR China
2 - Research Center for Theoretical Astrophysics, School of Physics, The University of Sydney, NSW 2006, Australia
Received 4 December 2003 / Accepted 19 April 2004
We ascribe the interpretation of the twin kilohertz Quasi Periodic Oscillations (kHz QPOs) of X-ray spectra of Low Mass X-Ray Binaries (LMXBs) to MHD Alfvén wave oscillations in the different mass density regions of the accreted matter at the preferred radius, and the upper kHz QPO frequency coincides with the Keplerian frequency. The proposed model concludes that the kHz QPO frequencies depend inversely on the preferred radius, and that theoretical relation between the upper frequency ( ) and the lower frequency ( ) is , which is similar to the measured empirical relation. The separation between the twin frequencies decreases (increases) with increasing kHz QPO frequency if the lower kHz QPO frequency is more (less) than 400 Hz.
Key words: stars: neutron - X-rays: stars
The launch of the X-ray timing satellite, Rossi X-ray Timing Explorer (RXTE), led to the discovery of Quasi Periodic Oscillations (QPOs) of LMXBs in their X-ray brightness, with frequencies Hz (see van der Klis 2000, for a recent review). Thereafter, much attention has been paid to the kHz QPO mechanism of LMXB; however the proposed models are still far from explaining all detected data. The Z sources (Atoll sources), which are high (less) luminous neutron star (NS) LMXBs (Hasinger & van der Klis 1989), typically show four distinct types of QPOs (van der Klis 2000). At present, these are the normal branch oscillation (NBO) Hz, the horizontal branch oscillation (HBO) Hz, and the kHz QPOs Hz that typically occur in pairs in more than 20 sources, with upper frequency and lower frequency . In 11 sources, nearly coherent burst oscillations Hz have also been detected during thermonuclear type I X-ray bursts; these are considered to be the NS spin frequencies or their first overtone (see, e.g., Strohmayer & Bildsten 2003). Moreover, the existence of a third kHz QPO has also been reported in three low-luminosity sources (Jonker et al. 2000). All of these QPOs but the burst oscillations have centroid frequencies that increase with the inferred mass accretion rate . Furthermore, the frequencies and , as well as the frequencies and , follow very similar relations in five Z sources, and the QPO frequencies of LMXBs and black hole candidates (BHC) have a tight and systematical correlation over three orders of magnitude in frequency (Psaltis et al. 1998, 1999; Belloni et al. 2002).
Various theoretical models have been proposed to account for the QPO phenomenon in X-ray binaries (for a review see, e.g., Psaltis 2000). In the early detection of RXTE, the upper kHz QPO ( ) was simply considered to originate from the Keplerian orbital frequency at the preferred radius, and the lower kHz QPO ( ) is attributed to the beat of this frequency with the stellar spin frequency (Strohmayer et al. 1996; Miller et al. 1998). However, this beat model is inadequate, for the detected frequency separation ( ) decreased systematically with instantaneous (see, e.g., van der Klis 2000). Later on, general relativistic effects were invoked to account for kHz QPOs (Stella & Vietri 1999; Stella et al. 1999; Psaltis & Norman 2000), which can satisfactorily explain the variation in kHz QPO separation . Moreover, the theory of epicyclic parametric resonance in relativistic accretion disks was proposed (Abramowicz et al. 2003), where the twin kHz QPOs occur at the frequency of meridional oscillation and the radial epicyclic frequency in the same orbit, which can explain the frequency ratio 3:2 detected in black hole candidates. Although many other feasible ideas have also proposed, such as the disk seismic model (Wagoner 1999), a two-oscillator model (Osherovich & Titarchuk) and the photon bubble model (Klein et al. 1996), no model has yet explained satisfactorily all observed QPO phenomena of LMXBs until now.
In this paper, the MHD Alfvén wave oscillation model is proposed, and its predictions and comparisons with the well detected sample sources are shown in the figures.
As a phenomenological prescription, we associate the twin kHz QPO frequencies with the Alfvén wave oscillation frequencies (AWOFs) at a preferred radius described in the appendix, where the AWOF with the spherical accretion mass density coincides with the Keplerian orbital frequency and the AWOF with the polar accretion mass density corresponds to a lower frequency. As a preliminary investigation, we are not concerned with the actual mechanisms of producing the kHz QPOs in the X-ray fluxes of LMXBs (see, e.g., Miller et al. 1998; Psaltis 2000). Rather, our main purpose is to stress the consistence between the model and the measured kHz QPO data, and leave the arguments concerning the possibility of the mechanism to a later paper.
At the preferred radius r, the Alfvén
dipole magnetic field strengths at radius r and at the surface
of star with radius R, respectively. Therefore,
|Figure 1: versus plot. The horizontal axis is the lower kHz QPO frequency and the vertical axis is the upper kHz QPO frequency . The kHz QPO data of four detected sample sources are plotted. The model presents a good consistence with the measured data for the NS averaged mass density parameters A=0.6, 0.7 and 0.8, which are shown in the three theoretical curves from bottom to top.|
|Open with DEXTER|
|Figure 2: versus plot. The horizontal axis is the lower kHz QPO frequency and the vertical axis is the twin kHz QPO frequency separation . The kHz QPO data of four detected sample sources are plotted. The model presents a good consistence with the measured data for the NS averaged mass density parameters A=0.6, 0.7 and 0.8, which are shown in the three theoretical curves from bottom to top.|
|Open with DEXTER|
|Figure 3: The ratio between the twin kHz QPO frequencies versus the lower kHz QPO frequency. The three theoretical curves represent the NS mass density parameter conditions with A=0.6, 0.7 and 0.8 from bottom to top. The data of four sample sources are plotted.|
|Open with DEXTER|
|Figure 4: The ratio between the twin kHz QPO frequencies versus the parameter . The horizontal axis is Xand the vertical axis is the ratio . It is found that the theoretical curve is independent of the NS mass density parameter A. The two horizontal lines represent y=5:4 and y=5:3, from bottom to top, respectively, which cover the ratios of the detected kHz QPO data of the four samples.|
|Open with DEXTER|
In conclusion, the consistence between the model and the detections is robust, and the main results are summaried as follows: (1) the twin QPOs are inversely (proportionally) related to the radii (the accretion rate); (2) the theoretical relation between the twin frequencies is ; (3) the separation between the twin frequencies decreases (increases) with increasing the kHz QPO frequency if the lower kHz frequency is more (less) than 400 Hz; (4) the ratio between the twin frequencies is only related to the scaled radius parameter X, and the homogenous kHz QPO frequency distributions of four detected sources indicate that these frequencies are from regions close to the surface of the NS, approximately, from 1.08R to 1.25R. With regard to conclusion (3), for the theoretical - relation with Hz future detection is needed to confirm this prediction.
Thanks are due to T. Belloni, M. Méndez, and D. Psaltis for providing the data files, and discussion with T. P. Li is highly appreciated.
The Alfvén radius
is determined by the condition that
the magnetic ram pressure
should match the plasma momentum energy density (Shapiro & Teukolsky 1983),
the surface field in unit of 108 Gauss,
the accretion rate in unit of 1017 g/s,