arxiv:astro-ph/ v1 24 May 2003

Transcription

1 Mon. Not. R. Astron. Soc., () Printed 6 September 8 (MN LATEX style file v.) An analytic model of rotating hot spot and khz QPOs in X-ray binaries arxiv:astro-ph/35478v 4 May 3 Ding-Xiong Wang,3, Ren-Yi Ma, Wei-Hua Lei and Guo-Zheng Yao Department of Physics, Huazhong University of Science and Technology, Wuhan,4374,China Department of Physics, Beijing Normal University, Beijing 875, China 3 Send offprint requests to: D.-X. Wang (dxwang@hust.edu.cn) Accepted INTRODUCTION ABSTRACT With the existence of a magnetic field connecting a Kerr black hole (BH) to its surrounding disc, energy and angular momentum can be transferred from the BH to the disc (Blandford 999; Li ; Li a, hereafter Lia; Wang, Xiao & Lei, hereafter WXL), which can be regarded as one of the variants of the Blandford-Znajek (BZ) process proposed two decades ago (Blandford & Znajek 977). Henceforth this energy mechanism is referred to as the magnetic coupling (MC) process. Recently Wilms et al. () found that XMM-Newton observation of the nearby bright Seyfert galaxy MCG reveals an extremely broad andredshifted Fe Kα line indicatingits origin from thevery most central region of the accretion disc, and they suggested that the rotating energy of a BH is extracted by the magnetic field connecting the BH or plunging region to the disc. Later Li pointed out that the magnetic coupling of a spinning BH with a disc can produce a very steep emissivity with index α = , which is consistent with the above XMM-Newton observation (Li b, hereafter Lib). This result can be regarded as the observation signatures of the existence of the MC process. Very recently we obtained the same result in a more detailed MC model, which is consistent with the XMM-Newton observation with reference to a variety of parameters of the BH-disc system (Wang, Lei & Ma 3, hereafter WLM). Another relevance of MC model to the observation is An analytic model of rotating hot spot is proposed to explain kilohertz quasiperiodic oscillations (khz QPO) in X-ray binaries, which is based on the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc. The hot spot in the inner region of the disc is produced by energy transferred from a spinning BH with non-axisymmetric magnetic field. The frequency and energy band of the hot spot turns out to be related to six parameters, by which the strength and the azimuthal angular region of the bulging magnetic field on the BH horizon, the mass and spin of the BH, and the power-law index of the magnetic field varying with the radial coordinate on the disc are described. In addition, the correlation of the fluctuation of the bulging magnetic field with the widths of QPO frequency is discussed. Key words: accretion, accretion discs black hole physics kilohertz quasi-periodic oscillations (khz QPO) in X-ray binaries. As argued by van der Klis (), khz QPOs in X- ray binaries probably originate from the inner edge of an accretion disc with a BH of stellar-mass, since millisecond is the natural timescale for accretion process in these regions. Recently Li suggested that the non-axisymmetric MC of a spinning BH of stellar-mass with a disc might be used to explain the khz QPOs (Li, hereafter Li). Motivated by Li s suggestion we propose an analytic model of a rotating hot spot to explain the khz QPOs, where the hot spot arises from the energy transferred to the inner region of the relativistic disc by the closed field lines of a nonaxisymmetric magnetic field on the BH horizon. Six parameters are involved in our model. Among them three are used to describe the strength and azimuthal angular region of the bulging magnetic field on the BH horizon, another two are used to describe mass and spin of the Kerr BH, and one is the power-law index of the magnetic field varying with the radial coordinate on the disc. We give a detailed discussion on a rotating hot spot for fittingkhz QPOs in terms of these parameters. This paper is organized as follows. In Section we show that the position of a hot spot is determined by nonaxisymmetric magnetic field with a mapping relation between the poloidal angular coordinate on the horizon and the radial coordinate on the disc. It turns out that the hot spot is limited very close to the inner edge of the disc for the

2 RADIAL POSITION OF ROTATING HOT SPOT θ θ Black Hole r in r out Disk Figure. The poloidal magnetic field connecting a spinning BH with its surrounding disc f (φ) +δ f (φ) +δ φ φ (a) (b) Figure. Azimuthal profile of the non-axisymmetric magnetic field on the BH horizon appropriate values of the BH spin and the power-law index of the magnetic field on the disc. In Section 3 we prove that the hot spot can reach the energy level of emitting X-ray, provided that the bulging magnetic field on the horizon is strong enough to reach 8 Gauss. Furthermore the rotatingfrequencyofahot spotis expressedin terms ofkeplerian angular velocity at the place where r F is maximized, where F is the sum of radiation flux due to disc accretion and that arising from the MC between the spinning BH and the disc. In Section 4 we fit khz QPO widths by assuming the fluctuation of the bulging magnetic field on the horizon. Finally, in Section 5, we give a summary for our model. Throughout this paper the geometric units G = c = are used. π π φ φ In order to facilitate the discussion of the rotating hot spot we make the following assumptions. (i) The magnetic field is assumed to connect the BH with the surrounding disc as shown in Fig.. The radii r in and r out are the radii of inner and outer boundaries of the MCregion, respectively,and θ andθ are thecorresponding poloidal angular coordinates on the horizon. (ii) The toroidal profile of the magnetic field is nonaxisymmetric on the BH horizon, being expressed by a function f (φ) of the azimuthal coordinate φ as shown in Fig.. The magnetic field is assumed to vary in a power-law as the radial coordinate of the disc as suggested by Blandford (976). (iii) The disc is both stable and perfectly conducting, and the closed magnetic field lines are frozen in the disc. (iv)the disc is thinandkeplerian, locating in theequatorial plane of the BH with the inner boundary at the last stable circular orbit. Generally speaking, the magnetic field on the horizon is non-axisymmetric. Unfortunately, a reliable profile of the magnetic field on the horizon has not been given due to lack of knowledge about it. As a simple analysis we assume the profile function f (φ) as shown in Fig.(a), which is expressed by B H (φ) = BH f (φ), { +δ, < φ < φ; f (φ), φ φ π. where BH is root-mean-square of themagnetic fieldover the poloidal angular coordinate from θ to θ. As shown in Fig.(a), the parameter δ is used to describe the strength of the bulging magnetic field in the region < φ < φ, which is the azimuthal angular range of the hot spot on the disc. The drawback of the pattern given in Fig.(a) lies in the uncontinuity of the magnetic field at φ = and φ. In fact, function f (φ) can be modified by a smooth joint-function as shown in Fig.(b), and we shall modify the model of rotating hot spot by using a more realistic profile function in future. The radial range of the hot spot is determined by a mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the disc as follows: cosθ = where G(a ;ξ,n) = () G(a ;ξ,n)dξ, () ξ n χ ms +a χ 4 ms ξ +a χ 6 ms ξ 3 (+a χ 4 ms +a χ 6 ms)( χ ms ξ +a χ 4 ms ξ ). The angle θ = π/ is assumed in derivation. Equation () associated with equation (3) is the revised version of the mapping relation derived in WXL. Parameter ξ r/r ms is a dimensionless radius of the disc in terms of the radius of the last stable circular orbit, r ms = Mχ ms, and a J / M is dimensionless parameter of the BH spin defined by the BH mass M and angular momentum J. Following Blandford (3)

3 An analytic model of rotating hot spot and khz QPOs in X-ray binaries ξout, ξc.5.5 ( out, ξ ' out ) (.3594, ) ξout, ξc.5.5 ( out, ξ ' out ) (.3594, ) ξout, ξc ( out, ξ ' out ) (.3594, ) (a) (b) (c) Figure 3. The curves of ξ c(dot-dashed line) and ξ out versus a for < a < with n =.,.5 and 3 in solid, dashed and dotted lines, respectively. The values of θ are given as (a)θ = π/, (b)θ = π/6 and (c)θ = π/4. (976) the normal components of magnetic field, B z, on the disc is assumed to vary with ξ as follows. B z ξ n. (4) In deriving equations () and (3) the boundary condition is used as given in WLM, i.e. πr HB = π D(r ms)b z (r ms), (5) where B is the normal component of magnetic field at horizon, and D(r ms) is the cylindrical radius at r ms, D(r ms) = Mχ ms +χ 4 msa +χ 6 msa. (6) Equation (5) implies the conservation of magnetic flux corresponding to the two loops of the same width, and the ratio of B to B z (r ms) varies from.8 to 3 for < a <. Considering that the strength of the magnetic field at the horizon is likely greater than that in the disc by numerical simulation (Ghosh & Abramowicz 997 and the references therein), we think that the relation (5) is more reasonable than the boundary condition B = B z(r ms) given in WXL. Substituting θ = θ and ξ = ξ out into equation (), we obtain an integral equation for determining the outer boundary ξ out r out/r ms of the MC region as follows: cosθ = out G(a ;ξ,n)dξ. (7) From equation (7) we find that ξ out(a ;n,θ ) behaves as a non-monotonic function of a for different values of n and θ as shown in Fig.3. In WLM the transfer direction of energy and angular momentum between a rotating BH and its surrounding disc is discussed in detail by using the co-rotation radius r c, which is defined as the radius where the angular velocity Ω D of the disc is equal to the BH angular velocity Ω H. The parameter ξ c(a ) r c/r ms is a function of a determined by the following equation, Ω D/Ω H = (+q) a where [ ( ) ] 3 ξχms +a =, (8) Ω D = M(χ 3 +a ), ΩH = q = a. a M(+q), In equation (9) χ r/m = ξχ ms is a dimensionless radial parameter on the disc, and χ ms is defined as χ ms rms/m. The MC region is divided by r c into two parts: the inner MC region (henceforth IMCR) for < ξ < ξ c and the outer MC region (henceforth OMCR) for ξ c < ξ < ξ out. The parameter ξ c decreasing monotonically with a is shown by the dot-dashed line in Fig.3, and the transfer direction of energy and angular momentum is described as follows. (i) For.3594 < a < the transfer direction is from the BH to the disc with ξ c <. (ii) For a < a out the transfer direction is from the disc to the BH with ξ c > ξ out, where a out is the BH spin corresponding to the intersection of the curve ξ c(a ) with the curve ξ out(a,n,θ ) as shown in Fig.3. (iii) For a out < a <.3594 the transfer direction is from the BH to OMCR with ξ c < ξ < ξ out, while it is from IMCR to the BH with < ξ < ξ c. Since a hot spot is produced by the energy transferred from the BH to the disc, its radial range is expressed by { ξc < ξ < ξ out, for a out < a <.3594 (casei) () < ξ < ξ out, for.3594 < a < (caseii). From Fig.3 we find that the outer boundary ξ out is very close to ξ in = in both case I and case II, provided that the power-law index n is small, such as n =.,.5. This result seems insensitive to the values of θ. Therefore the hot spot is confined to a small radial region very near to the inner edge of the disc. 3 RADIATION ENERGY AND ROTATING FREQUENCY OF ROTATING HOT SPOT Based on conservation of energy and angular momentum the basic equations for BH evolution in the coexistence of disc accretion and the MC process are written as (9)

4 4 dm/dt = E ms Ṁ D P MC, () dj/dt = L ms Ṁ D T MC. () where P MC and T MC are MC power and torque in the MC process, respectively. In WXL we derived P MC and T MC for axisymmetric magnetic field, by using an improved equivalent circuit as follows: P axis MC = B H M a T axis MC = 4 B H M 3 a (+q) β( β)g(a ;ξ,n) csc θ ( q) dξ, (3) ( β)g(a ;ξ,n) csc θ ( q) dξ. (4) Equations (3) and (4) can be extended to the case for non-axisymmetric magnetic field as follows: P naxis MC T naxis MC = λp axis MC, (5) = λt axis MC, (6) where the superscripts axis and naxis indicate the quantities for axisymmetric and non-axisymmetric magnetic field, respectively. The parameter λ is expressed by λ = [(+δ)ε+( ε)] = (+δε), (7) where ε φ/π is a parameter indicating the azimuthal region of the bulging magnetic field on the horizon expressed by equation (). Equations (5) and (6) can be derived in the same way as in WXL, provided that the magnetic flux between two adjacent magnetic surfaces Ψ = B Hπ l is replaced by Ψ = BH [(+δ) φ+(π φ)] l = B H π (+δε) l, (8) where equation (8) is given for the bulging magnetic field expressed by equation (). Obviously, equations (5) and (6) reduce to axisymmetric ones for either δ or ε with λ =, and they also become axisymmetric ones for ε with λ = (+δ). Based on equations (3) (6) we can express the rate of extracting energy and angular momentum from the BH to the hot spot as follows: P HS MC = (λ )P axis MC +εp axis MC = (λ+ε )P axis MC, (9) T HS MC = (λ )T axis MC +εt axis MC = (λ+ε )T axis MC. () From equations (5), (6), (9) and () we find that the ratio of PMC HS topmc naxis is the same as that of TMC HS totmc naxis, i.e. R(δ,ε) = PHS MC P naxis MC = ε (+δε), = THS MC = λ+ε T naxis λ MC () where the ratio R(δ,ε) is used to express the fractional energy and angular momentum attributable to the hot spot in the MC process. The curves of R(δ,ε) varying with the parameters δ and ε are shown in Fig.4. From Fig.4 we find that the ratio R(δ,ε) is independent of the parameters M, a and n, and it increases almost linearly with the parameters δ and ε. For an appropriate angular region of the bulging magnetic field on the horizon, such as ε =., we have the ratio varying from. to.8 R R δ (a) ε Figure 4. The curves of R(δ,ε) (a) versus δ for < δ <.5 with ε =.,.5 and (b) versus ε for < ε <. with δ =.5,.3 in solid and dashed lines, respectively. for a small value range of the strength of the bulging magnetic field, < δ <.5. These results show that the effects of the bulging magnetic field on the energy and angular momentum transferred to the hot spot are significant. Since the hot spot is located very near to the inner edge of the disc, it should be optically thick. So we can define its effective radiation temperature (T HS) eff as a blackbody spectrum: (T HS) eff = [( F DA +F HS MC)/ σ ] /4 () where σ is the Stefan-Boltzmann constant. The fluxes F DA and FMC HS are radiated from the hot spot due to disc accretion and the MC process, respectively, and are expressed as follows (Page & Thorne 974; Lia): F DA = Ṁ D 4πM χ 4 ms dω D ξdξ ( E + Ω DL +) (b) ( FMC HS = dω D ξdξ E + Ω DL +) ( E + Ω DL +) dl + dξ dξ, (3) ( E + Ω DL +) H HSξdξ. (4)

5 An analytic model of rotating hot spot and khz QPOs in X-ray binaries HS F MC / F, FDA / F HS F MC / F, FDA / F.5..5 HS F MC / F, FDA / F ξ ξ ξ (a) (b) (c) Figure 5. The curves of F DA /F (dot-dashed lines) and F HS MC /F for < ξ < ξ out and ε =., δ =.5 with n =.,.5 and 3. in solid, dashed and dotted lines, respectively. (a) a =.3594; (b) a =.8; (c) a =.998. In this paper we only discuss the hot spot produced by the bulging magnetic field on the BH with spin greater than a =.3594, which corresponds to case II in equation (). In equations (3) and (4) E + and L + are specific energy and angular momentum of accreting particles, respectively, and are expressed respectively by (Novikov & Thorne 973) E + = and χ +a χ 3 ( 3χ +a χ 3 ) / (5) L + = Mχ( a χ 3 +a χ 4 ). (6) ( 3χ +a χ 3 ) / In equation (4) H HS is the flux of angular momentum transferred from the BH to the hot spot by the bulging magnetic field, being related to TMC HS by 4πεrH HS = r φh HS = T HS MC/ r (7) Combining equations (), (3), (4), () and (7), we have { A(a,ξ,n,δ,ε)ξ n, < ξ < ξ out H HS/H =, (8), ξ > ξ out where A(a,ξ,n,δ,ε) = (+δ+δ ε)a ( β)(+q) F πχ ms[csc θ ( q)] A(a,ξ) and F A(a,ξ) = +a χ 4 ms ξ +a χ 6 ms ξ 3 (+a χ 4 ms +a χ 6 ms)( χ ms ξ +a χ 4 ms ξ ) (9) H = B H M =.48 3 B 4m BH g s. (3) In equation (3) B 4 and m BH are BH and M in the units of 4 Gauss and one solar mass, respectively. Since the magnetic field on the horizon is brought and held by its surrounding magnetized disc (Macdonald & Thorne 98; Thorne, Press & Macdonald 986), there must exist some relations between B H and ṀD. As a matter of fact these relations might be rather complicated, and would be very different in different situations. One of them is given by considering the balance between the pressure of the magnetic field on the horizon and the ram pressure of the innermost parts of an accretion flow (Moderski, Sikora & Lasota 997), i.e. B H/ (8π) = Pram ρc ṀD /( 4πr H ), (3) From equation (3) we assume the relation as Ṁ D = B H M (+q) /. (3) The accretion rate Ṁ D is still axisymmetric in spite of the existence of the bulging magnetic field on the horizon, provided that the extra angular momentum transferred into the disc is totally radiated away from the disc. Incorporating equations (3), (3), (4) and (8), we have F DA/F = f DA(a,ξ) (+q) f DA(a,ξ) = 8πχ 4 ms(e + Ω D L + ) (33) and with F HS MC/ F = f MC, f MC = dω D ξdξ MdΩ D (E + Ω D L + ) ξdξ ( E + Ω DL +) dl + dξ dξ ( E + Ω DL +) A(a,ξ,n,δ,ε)ξ n dξ, (34) F = B H c = B 4 erg cm s. (35) Inspecting equations (33) (35), we obtain the following results: (i) Both F DA and FMC HS are proportional to BH, and are independent of the BH mass due to equation (3). (ii) The bulging magnetic field on the horizon only affects FMC HS rather than F DA, since Ṁ D is independent of the bulging magnetic field. By using equations (33) and (34) we have the curves of F DA and FMC HS varying with ξ as shown in Fig. 5. From Fig.5 we obtain the following results:

6 6 Table. Radiation energy and frequency of a rotating hot spot produced by the MC process and disc accretion with B 4 =, δ =.5 and ε =.. a.6 n ( ) T max HS (K) eff Emax HS (kev) ξmax m BHν HS (Hz) (i) The relation FMC HS > F DA always holds for.3594 < a <, implying that the contribution from the MC process always dominates over that from disc accretion for the above value range of the BH spin. (ii) Both FMC HS and F DA approach zero as ξ ξ in =, and both vary non-monotonically with ξ. The flux FMC HS varies much more sharply than F DA does, attaining its peak value between the inner and outer boundaries of the hot spot. (iii) The peak of FMC HS is closer to the inner edge of the disc than that of F DA. The more is the value of a, the greater is the peak value, and the closer is the peak to the inner edge. Substituting equations (33) and (34) into equation (), we have the effective temperature and the corresponding radiation energy of the hot spot as follows: (T HS) eff = T [f DA +f MC] /4, (36) E HS k B (T HS) eff = E [f DA +f MC] /4, (37) where T = (F /σ) / B / 4 K, (38) and E = k B (F /σ) /4 4.4 B / 4 kev. (39) The maxima of (T HS) eff and E HS corresponding to the different values of the BH spin and the power-law index are shown in Table. It is easy to prove that the maxima of either (T HS) eff or E HS increase monotonically with the increasing BH spin. For the fast spinning BH with spin.6 < a <.998 we conclude that the radiation energy of the hot spot reaches.6.7kev, provided that the bulging magnetic field reaches B 4 or 4 Gauss, while this value range becomes.6.7kev for the bulging magnetic field reaching 4 B 4 or 8 Gauss. So the hot spot can have energy reaching soft to hard X-ray in the MC process, provided that a fast-spinning BH has the bulging magnetic field of 4 8 Gauss on the horizon. Since the magnetic field lines are frozen in the disc plasma, the rotating frequency of the hot spot is regarded as QPO frequency. As suggested by Nowak and Lehr (998), the QPO frequency can be worked out by calculating Keplerian angular velocity at the place where r F attains its B Figure 6. The contours of E HS of constant values (kev) for.3594 < a <.998 and 3 < B 4 < 4 with δ =.5, ε =. and n =.,.5, 3. in solid, dashed and dotted lines, respectively. maximum,andf is thesum off MC andf DA.Wecan define a dimensionless function as follows: F QPO r F / rmsf = ξ ( / ) F DA/F +FMC HS F (4) Incorporating equations (33), (34) and (4), we obtain ξ max corresponding to the maximum of function F QPO, and the QPO frequency ν HS can be calculated by ν HS ν (ξ 3/ maxχ 3 ms +a ), (4) where we have ν HS Ω D/π ξ=ξmax and ν (m BH) Hz. The values of ξ max and ν HS corresponding to different values of BH spin and power-law index are shown in Table, where ξ max turns out to be very close to the inner edge of the disc. For the fast-spinning BH with spin.6 < a <.998 we have ν HS varying from.4khz/m BH to khz/m BH with m BH > 3. In order todescribe the variation of E HS and ν HS of the rotating hot spot more clearly we plot the contours of E HS and ν HS in terms of the concerning parameters as shown in Figs 6 and 7, respectively. From Fig. 6 we find that less magnetic field is required as the BH spin increases for the given value of E HS. Furthermore the more is the BH spin, the more steeply the contours

7 An analytic model of rotating hot spot and khz QPOs in X-ray binaries m BH 5 dν HS / ν HS Figure 7. The contours of ν HS of constant values (Hz) for.3594 < a <.998 and 3 < m BH < 3 with δ =.5, ε =. and n =.,.5, 3. in solid, dashed and dotted lines, respectively (a) of E HS decline. This result implies that the MC process becomes very efficient in transferring energy from a very fast spinning BH to the hot spot. From Fig. 7 we find that more BH mass is required as the BH spin increases for the given values of ν HS. Furthermore the more is the BH spin, the more contours of ν HS will be crossed over for the given BH mass. This result implies that ν HS increases very rapidly as the BH spin a approaches unity for the given BH mass. The BH spin acts as a bridge to link the contours of E HS and ν HS in Figs 6 and 7. This linkage is helpful to fit observations of khz QPOs in X-ray binaries in some scope. 4 FLUCTUATION OF PARAMETERS AND WIDTHS OF QPO FREQUENCIES It is shown that the widths of observed frequencies of kilohertz QPOs in X-ray binaries are rather narrow, usually varying from several to some dozens of hertz, e.g. as shown in Tables -3 given by van der Klis (). However, it is difficult for us to explain the narrow range of QPOs by the rotating hot spot produced by the stationary bulging magnetic field on the horizon. As a matter of fact we can fit the narrow range of QPO widths by the fluctuation of δ and ε in our model. The argument is given as follows. From equation (4) we find that QPO frequency is determined by ξ max, and the latter depends on the maximum of function F QPO in equation (4), which consists of the sum of F DA/F and FMC/ HS F. As mentioned above that the peak of FMC/ HS F is located closer to the inner edge of the disc than that of F DA/F, so we infer that the narrow range of QPO widths might be produced by the variation of FMC/ HS F, provided that the parameters δ and ε are fluctuating rather than constant. The curves of dν HS/ν HS versus a for the given fluctuation of δ and ε are shown in Fig.8. From Fig.8 we have the following results. (i) The value of dν HS/ν HS varies between.4 to. with δ =.3 and ε = for.54 < a <.998. This dν HS / ν HS Figure 8. The curves of dν HS /ν HS versus a for.54 < a <.998, (a) n = 3, ε =., ε =, δ = 5, δ =.3 (solid line), δ =.3 (dashed line), (b) n = 3, δ =.5, δ =, ε =., ε =. (solid line), ε =. (dashed line). result implies that the range of QPO widths is of dozens of hertz, and could be produced by the fluctuation of the strength of the bulging magnetic field of the BH as shown in Fig.8(a). (ii) The effect of the variation of ε on the range of khz QPO width is about one order of magnitude less than that of δ as shown in Fig.8(b). The values of dν HS/ν HS corresponding to different values of the concerning parameters are listed in Table. As shown in Fig.8(a) and Table we can fit the range of khz QPO width by a variety of parameters in our model, such as a, n, ε, δ, ε and δ. 5 SUMMARY In this paper an analytic model of rotating hot spot is proposed to explain khz QPOs in X-ray binaries by using the MC of a rotating BH with the disc. Six parameters are used in our model. (i) The parameters δ, ε and B 4 are given for describing the bulging magnetic field on the horizon by equation (). (b)

8 8 Table. The values of dν HS /ν HS corresponding to different values of the concerning parameters with ε =., ε = and δ =.5. a n m BH ν HS dν HS /ν HS δ =. δ =. δ =. δ = Specifically speaking, the parameter B 4 is used to indicate root-mean-square of the magnetic field over the poloidal angular coordinate from θ to θ. The parameters δ and ε are in charge of the fractional energy and angular momentum attributable to the hot spot by equation (). The fluctuation of δ and ε can affect the range of khz QPO widths. (ii)theparameters m BH anda aregivenfordescribing the mass and spin of the BH, which determine the frequency of QPOs by equations (4). (iii) The parameter n is used to indicate the concentration of magnetic field on the central region of the disc, and it is involved in equation (7) to confine the radial width of the hot spot. In addition, it is involved in equation (4) to determine ξ max. We can find that the index n plays a role of fine-tuning energy and frequency of the hot spot and the range of khz QPOs widths as shown in Table, Fig.6, Fig.7 and Table, respectively. Nowak, M. A., Lehr, D. E., 998, in Theory of Black Hole Accretion Disks, ed. M. A. Abramowicz, G. Bjornsson, & J. E. Pringle, (Cambridge: Cambridge Univ. Press), p.33, preprint (astro-ph/984) Page D. N., Thorne K. S., 974, ApJ, 9, 499 Thorne K. S., Price R. H., Macdonald D. A., 986, Black Holes: The Membrane Paradigm, Yale Univ. Press, New Haven van der Klis, M.,, ARA&A, 38, 77 Wang D. X., Xiao K., Lei W. H.,, MNRAS, 335, 655 (WXL) Wang D. X., Lei W. H., Ma R. Y., 3, MNRAS (in press), preprint (astro-ph/335) (WLM) Wilms J. et al., MNRAS, 38, L7 Acknowledgments. This work is supported by the National Natural Science Foundation of China under Grant No. 734 and No. 53. The anonymous referee is thanked for his suggestion about the effects of the parameters δ and ε on the rotating hot spot. REFERENCES Blandford R. D., 976, MNRAS, 76, 465 Blandford R. D., Znajek R. L., 977, MNRAS, 79, 433 Blandford R. D., 999, in Scllwood J. A., Goodman J., eds, ASP Conf. Ser. Vol. 6, Astrophysical Discs: An EC Summer School, Astron. Soc. Pac., San Francisco, p.65 Ghosh P., Abramowicz M. A., 997, MNRAS 9, 887 Li L. -X., ApJ, 533, L5 Li L.-X. in X-rayEmissionfromAccretion onto Black Holes, ed. T. Yaqoob and J. H. Krolik (published electronically at JHU/LHEA workshop, June -3, (Li) Li L. -X. a, ApJ, 567, 463 (Lia) Li L. -X. b, A&A, 39, 469 (Lib) Macdonald D., Thorne K. S., 98, MNRAS, 98, 345 Moderski R., Sikora M., Lasota J.P., 997, in Relativistic Jets in AGNs eds.m. Ostrowski, M. Sikora, G. Madejski & M. Belgelman, Krakow, p. Novikov, I. D., & Thorne, K. S., 973, in Black Holes, ed. Dewitt C (Gordon and Breach, New York) p.345