Possible signature of magnetic fields related to quasi-periodic oscillation observed in microquasars

The study of quasi-periodic oscillations (QPOs) of X-ray flux observed in the stellar-mass black hole binaries can provide a powerful tool for testing of the phenomena occurring in strong gravity regime. The high frequency (HF) QPOs usually come in pairs of double peaks which have a frequency ratio close to $3:2$. In addition to HF QPOs, some sources display simultaneous existence of the low frequency (LF) QPOs in Fourier power spectra. We demonstrate that the explanation of these phenomena can be well related to the epicyclic oscillations of charged particles in accretion disks orbiting Kerr black holes immersed in external large-scale magnetic fields. Magnetized versions of the standard geodesic models of QPOs can explain the observationally fixed data from the three microquasars. We perform a successful fitting of the HF QPOs observed in three microquasars, GRS 1915+105, XTE 1550-564 and GRO 1655-40, containing black holes, for magnetised versions of both epicyclic rezonance and relativistic precession models and discuss the corresponding constraints of parameters of the model, which are the mass and spin of black hole and the parameter related to the external magnetic field. Assuming the main source of synchrotron radiation producing X-rays are the relativistic electrons, we estimate the magnetic field in the vicinity of the black hole in the three sources to be of order $10^{-5}$ Gs which can serve as possible signature of the Galactic magnetic field magnitude. For heavier particles (protons, ions) larger magnetic fields are necessary for fitting the data.


I. INTRODUCTION
Microquasars are binary systems composed of a black hole and a companion (donor) star; matter floating from the companion star onto the black hole forms an accretion disk and relativistic jets -bipolar outflow of matter along the black hole -accretion disk rotation axis. Due to friction, matter of the accretion disk becomes to be hot and emits electromagnetic radiation, also in X-rays in vicinity of the black hole horizon.
Applying the methods of spectroscopy (frequency distribution of photons) and timing (photon number time dependence) for particular microquasars, one can extract a useful information regarding the range of parameters of the system [34]. In this connection, the binary systems containing black holes, being compared to neutron star systems, seem to be promising due to the reason that any astrophysical black hole is thought to be a Kerr black hole (corresponding to the unique solution of general relativity in 4D for uncharged black holes which does not violate the no hair theorem and the weak cosmic censorship conjecture) that is determined by only two parameters: the black hole mass M and the dimensionless spin |a| ≤ 1.
One of the promising tools to probe the phenomena occurring in the field of the black hole candidates is the study of the quasi-periodic oscillations (QPOs) of the X-ray power density observed in microquasars. The current technical possibilities to measure the frequencies of After the first detection of QPOs, there were many attempts to fit the observed QPOs, and different models have been proposed, such as the hot-spot models, diskoseismic models, warped disk model and many versions of resonance models, developed in the framework of general relativity or alternative theories of gravity. It is particularly interesting that the characteristic frequencies of HF QPOs are close to the values of the frequencies of test particle, geodesic epicyclic oscillations in the regions near the innermost stable circular orbit (ISCO) which makes it reasonable to construct the model involving the frequencies of oscillations associated with the orbital motion around Kerr black holes. However, until now, the exact physical mechanism of the generation of HF QPOs is not known, since none of the models can fit the observational data from different sources. Surprisingly, the situation changes considerably if one assumes the oscillating matter to be electrically charged, and takes into account the influence of external magnetic fields in the curved black hole background. This is the main aim of the present paper, such assumptions are relevant and astrophysically motivated due to several reasons described below.
Magnetic fields have been detected and measured in nearly all celestial objects. The existence of accretion disk around black hole which usually contains highly conducting plasma can lead to the appearance of the regular electromagnetic fields in the vicinity of a black hole. Jets from stellar-mass compact objects such as the black hole microquasars are believed to be generated due to the interaction of charged matter with the gravitational field of a black hole combined with surrounding electromagnetic fields, rather than the radiation pressure. Even thought the total charge of an accretion disk is believed to be zero, the presence of plasma implies the existence of local charges in small limited regions along an accretion disk which may be influenced by the external magnetic fields [12]. Moreover, the collimation of jets at large distances supports the idea of the existence of large-scale magnetic fields in the interstellar medium.
This implies that the role of magnetic fields in the vicinity of astrophysical black holes cannot be neglected. Typically, the Galactic large-scale magnetic fields are weak, as can be estimated from the intensity of total synchrotron radiation producing the X-rays, or by methods including the Faraday rotation of the plane of polarization of a starlight in optical range propagating in a magnetized plasma, or by the Zeeman effect in the vicinity of stars. The average equipartition strength of the magnetic field in spiral galaxies is estimated to be of order ∼ 10 −5 G. The strength of the magnetic field increases closer to the center of the Galaxy, achieving from tens up to hundred Gauss around supermassive black holes [13]. The recent studies of the structure of magnetic fields in galaxies [7], clearly demonstrate the formation of the spiral patterns in almost all considered galaxies, including Milky Way. Therefore, black holes can be also immersed into the external, large scale magnetic fields that can have globally complicated structure, but at large distances from its source, in a finite element of space, its character can be considered as locally uniform.
We concentrate our attention on the particular and simplified case of weak magnetic fields that are considered to be asymptotically uniform at spatial infinity, described by the well known Wald solution for weakly magnetized black holes [54]. The condition of the weakness of magnetic field implies the fact that its stress-energy tensor does not violate the geometry of the black hole spacetime. In this sense the test field approximation can be applied when the strength of the magnetic field in the vicinity of a black hole with mass M fulfills the following condition (see, e.g., [16]) Obviously, for astrophysical black holes the weak field condition (1) is perfectly satisfied. However, as we will see in the following sections, for the motion of charged test particles the magnetic field plays a crucial role due to the large factor of the specific charge q/m entering the equations.
In the present paper, we consider the motion of charged test particles around a Kerr black hole immersed into an external asymptotically homogeneous magnetic field having the field lines aligned with the black hole rotation axis. We look especially for the existence and properties of the harmonic or quasi-harmonic oscillations of charged particles in the magnetized black hole backgrounds. The quasi-harmonic oscillations around a stable equilibrium location and the frequencies of these oscillations are then compared with the frequencies of the HF and LF QPOs observed in the three particular microquasars GRS 1915+105, XTE 1550-564, and GRO 1655-40 [27,51]. We demonstrate that the effect of the external magnetic field, which is usually neglected in most of the QPO models, allows us to fit the observed QPO frequencies for different sources in a single model. Assuming the detected X-rays in QPOs are generated due to the synchrotron radiation of electrons, we estimate the magnetic field strength necessary and sufficient for the fitting of the QPOs frequencies. We show that the magnetic field strength in this case is of the order which is comparable to the magnitude of the Galactic magnetic field. In principle, the oscillating matter could be constituted from heavier particles, protons, ions, charged dust, or some exotic objects similar to ball lighting, proposed in [49]; for such possible sources, the magnetic field should be significantly stronger than in the electron case, in dependence on specific charge.
Influence of the magnetic field on the charged test particle motion around black holes has been already studied in literature [16,23,24,26]. Charged particle motion in an uniform magnetic field, and related high frequency oscillations, have been studied for Schwarzschild black holes in [22], locally measured angular frequencies for the magnetized Kerr black hole with the dynamics of charged particles were studied in [53]. A complementary simple model of creation of charged jets has been introduced in [45]. For charged particle motion in the Reissner-Nordström and Kerr-Newman backgrounds see [5,10,31,32,41]. The effect of the dipole magnetic field on the oscillatory charged particle motion has been studied in [3,4]. The dipole magnetic field configuration is assumed to be much more relevant for neutron star Xray binaries, while for black hole X-ray binaries we can use the uniform configuration of the magnetic field. The influence of the magnetic field on the modes of the diskoseismic oscillations of an accretion disk orbiting a black hole, and its relation to observed HF QPOs, has been studies in [19].
Throughout the present paper, we use the spacelike signature (−, +, +, +), and the system of geometric units in which G = 1 = c. However, for expressions having astrophysical relevance we use the physical constants explicitly. Greek indices are taken to run from 0 to 3.

II. CHARGED PARTICLE MOTION IN THE FIELD OF MAGNETIZED KERR BLACK HOLE
We start from the general description of the motion of a charged test particle in the field of a rotating black hole in the presence of the external magnetic field. The geometry of the rotating black hole is given by the Kerr metric ds 2 = g tt dt 2 + 2g tφ dtdφ + g φφ dφ 2 + g rr dr 2 + g θθ dθ 2 , (2) with the nonzero components of the metric tensor g µν taking in the standard Boyer-Lindquist coordinates the form where Here, M is the gravitational mass of the black hole and a is its spin parameter. The physical singularity is located at the ring with r = 0, θ = π/2. The outer horizon is located at A. Uniform external magnetic field The external axially symmetric and asymptotically homogeneous magnetic field implies the electromagnetic field potential which is solution of the vacuum Maxwell equations taking in general the form [54] Here Q refers to the induced electric potential which exists due to the rotation of the black hole giving the contribution to the Faraday induction. This causes the process of a selective accretion of charged particles into the rotating black hole. However the stage of the selective accretion ends shortly in the most of the astrophysical scenarios when the black hole gets the induced charge Q = Q W ≡ 2aM B. The charge Q W is called the Wald charge. Finally, the four-vector potential of an electromagnetic field after the process of selective accretion completed takes form A α = 1 2 Bξ α (φ) , where B denotes the intensity of the asymptotically uniform magnetic fiels; ξ (φ) , ξ (t) denote the axial and time Killing vector fields. Hereafter, we will consider two limiting cases with the induced charge Q of a black hole: Q = 0 and Q = Q W . So, the four vector potential of electromagnetic field for these two cases takes the form One can neglect the gravitational effect of the induced charge, if the following condition is satisfied [53] Q << This condition holds well in our approach since the maximal induced charge of a magnetized rotating black hole, namely Q W , is 12 orders of magnitude less than Q G for the solar mass black holes. The electromagnetic field tensor reads F µν = A ν,µ − A µ,ν . For the potential (6), and in the equatorial plane θ = π/2, one can write the nonvanishing independent components of F µν as where the upper and lower signs correspond to the cases with Q = 0 and Q = Q W , respectively. The expressions for the Maxwell tensor F µν in the equatorial plane will be useful in the next section for the calculation of the frequencies of the charged particle epicyclic oscillations. In general, the antisymmetric tensor F µν has 4 nonzero components in the magnetized rotating black hole case. In order to describe the dynamics of the charged particle motion, we introduce the following dimensionless parameters Putting hereafter for the shortness M = 1, the Wald charge corresponds to Q W = 2aB. The parameter B we will be called the magnetic parameter [53]. The estimation of the magnitude of the magnetic parameter shows that the effect of even weak magnetic field cannot be neglected for the charged particle motion due to the large value of the specific charge q/m of the test particle. For example, for electrons in the vicinity of a stellar mass black hole M ≈ 10M ⊙ with the magnetic parameter B ≈ 1, the magnetic field strength corresponds to the value of B ≈ 0.002 G; for protons in the same conditions, the strength of magnetic the field B has to be 1836 times stronger. More detailed analysis of the magnetic parameter in astrophysical scenarios related to the generation of QPOs observed in the black hole microquasars will be given in Sec. IV.

B. Circular obits of charged particles
The motion of a charged particle with charge q and mass m in curved spacetime in the presence of electromagnetic field can be treated by the Lorentz equation where u µ = dx µ /dτ is the four-velocity of the particle. The regular motion of charged particles in magnetized black hole backgrounds is always bounded in radial direction which makes important to study the cyclic motion of the particles. In particular, the circular motion of a particle is possible in the equatorial plane θ = π/2, which follows from the properties of the symmetry of the geometry and background electromagnetic field. The four velocity of the circular motion has only 2 nonvanishing components, u µ = {u t , 0, 0, u φ }. This implies that for the equatorial circular motion, the radial component of Eq. (11) can be written in the form where the upper sign corresponds to the case with Q = Q = 0, and the lower sign to the case with Q = Q W . The normalization condition g µν u µ u ν = −1 gives the second equation for the nonzero components of the four-velocity in the form (13) which is obviously independent of the electromagnetic parameters B and Q. Eq. (12) and (13) allow us to find two expressions for two unknown quantities u t and u φ . The explicit form of analytical expressions for u t and u φ cannot be represented in a simple form, however we can easily solve it numerically. Results of the numerical calculations will be demonstrated in the following section. This form of the equations of motion will be useful for the calculation of frequencies of the quasi-harmonic oscillations of charged particles along circular orbits in the field of magnetized black holes.
One can find the equations of motion in a different way using the Hamiltonian formalism. Due to the symmetries of the Kerr geometry and the axial symmetry of an external magnetic field, the components of four velocity are related to the conservation of the following quantities where E = E/m and L = L/m are identified as the specific energy and specific angular momentum of the charged particle, respectively. This gives us a right to write the Hamiltonian of the charged particle motion in the form [55] where P µ is a generalized (canonical) four-momentum related to the four-velocity by the relation Equations of motion can be then written, in the form of the Hamilton equations where the affine parameter ζ = τ /m is related to the proper time of the particle τ . The dynamics of charged particles in magnetized rotating black hole background was studied in detail in [53], where we use the Hamiltonian formalism combined with the so called formalism of forces which allows us to find the analytical solutions for the motion of charged particles. Stability of the circular orbits is also discussed in [53] -the harmonic quasiperiodic oscillations are allowed only around the stable circular orbits.

III. HARMONIC OSCILLATIONS AS PERTURBATION OF CIRCULAR ORBITS
In order to describe the oscillatory motion of charged particle we use the perturbation of the equations of motion around the stable circular orbits. If a charged test particle is slightly displaced from the equilibrium position x µ 0 corresponding to a stable circular orbit located at the equatorial plane, the particle will start to oscillate around the stable orbit realizing thus epicyclic motion governed by linear harmonic oscillations. Following [2], one can introduce the deviation vector ξ µ (τ ) = x µ (τ ) − x µ 0 (τ ). Substituting the deviation vector into the equation (11) and using the first order of expansion, we get the equation for ξ µ in the form where Ω r , Ω θ are the epicyclic frequencies of the particle oscillations near the black hole, measured by a distant static observer. The radial and vertical oscillatory frequencies for the case Q = Q W = 2aB take the explicit form where the introduced coefficients read while u φ and u t are the components of the four velocity of the oscillating particle given by the solution of Eq. (12) and (13). The frequencies for the case with the zero induced charge will be similar to the Eq. (20) and (21). The difference between the frequencies of the oscillations in the presence, Ω (W) , and the absence, Ω (0) , of the Wald charge are given by the following expressions Note that the epicyclic frequencies (20) and (21) are given for the oscillations measured from the rest at infinity. Locally measured frequencies of the radial and latitudinal harmonic oscillations were presented in our previous paper [53] in terms of the energy and angular momentum of a charged particle at the circular orbit. In addition to the epicyclic frequencies Ω r and Ω θ , there exist also the Keplerian frequency, Ω K , and so called Larmor angular frequency, Ω L , associated with the pure contribution of an external uniform magnetic field, that are given by the relations The Larmor angular frequency depends only on the strength of the magnetic parameter B and the redshift factor u t . It is fully relevant in large distances from the black hole where influence of the uniform magnetic field becomes to be crucial. The expressions for the fundamental frequencies (20), (21) and (25) are given in the dimensionless form. In the physical units, one needs to extend the corresponding formulas by the factor c 3 /GM . Then the frequencies of the charged particle radial and latitudinal harmonic oscillations, and the Keplerian frequency, measured by the distant observers in Hz, are given by where i ∈ {r, θ, φ}. The radial profiles of the frequencies ν i are shown in Fig. 1, the analysis of the properties of these frequencies and their relations is given in the following subsection.

A. Frequencies of charged particle oscillations
The properties of the motion of charged particles around a black hole immersed in an uniform magnetic field can be separated into four qualitatively different classes [53]. For the particles a with positive charge one can distinguish the following orbits i. Prograde anti-Larmor orbits (PALO) with L > 0, B > 0. Particle is co-rotating. Magnetic field lines co-oriented with the rotation axis of the black hole. The Lorentz force is repulsive, i.e., directed outwards the black hole.
ii. Retrograde Larmor orbits (RLO) with L < 0, B > 0. Particle is counter-rotating. Magnetic field lines cooriented with the rotation axis of the black hole. The Lorentz force is attractive, i.e., directed towards the black hole.
iii. Prograde Larmor orbits (PLO) corresponding to L > 0, B < 0. Particle is co-rotating. Magnetic field lines counter-oriented to the rotation axis of the black hole. The Lorentz force is attractive, i.e., directed towards the black hole.
iv. Retrograde anti-Larmor orbits (RALO) corresponding to L < 0, B < 0. Particle is counter-rotating. Magnetic field lines counter-oriented with the rotation axis of the black hole. The Lorentz force is repulsive, i.e., directed outwards the black hole.
In case of the non-rotating Schwarzschild black holes, PALO/RALO and RLO/PLO classes coincide. In the absence of a magnetic field, for rotating black holes, the classes PALO/PLO and RALO/RLO coincide.
Classifying orbits as given above, we plot the radial profiles of the frequencies ν θ , ν r and ν φ of small harmonic oscillations of charged particle measured by distant static observer in Fig. 1, for characteristic values of black hole spin a and magnetic parameter B. The frequencies have been calculated using formula (26) in case of small displacements from equatorial plane. The magnetic field presence increases the value of the radial (horizontal) frequency ν r and the largest increase of the frequency ν r can be observed in the PALO/RALO cases. It is also notable that ν r does not tend to zero for large values of radial coordinate r; this is caused by the role of the uniform magnetic field at large distances from the black hole. The radial frequency ν r vanishes at the innermost stable circular orbit (ISCO), The position of ISCO for magnetic parameter B = ±0.05, in dependence on black hole spin a, is plotted in Fig.  2. The detailed analysis of the properties of the ISCO radii for different values of the magnetic field parameter has been presented in [53]. Below the ISCO radius, the particle starts to spiral down into the black hole. The Keplerian frequency ν φ related to the orbital coordinate velocity by relation (25) behaves differently for different classes of the circular orbits. For PALO/RALO configurations, with repulsive Lorentz force, ν φ frequency (as well as orbital u φ velocity [53], decreases with increasing the magnetic parameter B. For PLO/RLO configurations, with attractive Lorentz force, the frequency ν φ increases with increasing the magnetic parameter B. The radial profile of the vertical frequency ν θ does not have strong dependence on magnetic parameter B. Since the charged particle motion in vertical direction is not affected by the magnetic field, the Lorentz force acting on a charged particle is oriented in horizontal direction. The frequencies in Fig. 1. have been calculated for Kerr black hole with M = 10 M ⊙ mass. The change of the mass of the black hole leads to the change only in scale of the plots, but not their shapes due to (26). The frequencies of oscillations depend also on the induced charge parameter Q; however, for small values of param- Figure 1: Radial profiles of the frequencies of small harmonic oscillations ν θ , νr and ν φ of charged particle around Kerr black hole with mass M = 10M⊙ and spin a in external magnetic field measured by static distant observer. The first row represents the influence of magnetic parameter B in the absence of rotation (first two figures) and rotation in the nonmagnetic case (last two figures). Second row represents the Kerr black hole with a = 0.7 with magnetic parameter B = ±0.1. Since the differences between Q = 0 and Q = QW cases is negligible, see discussion of eq. (23-24), the case Q = QW is considered only. Last row represents the effect of nearly extremal Kerr black hole with a = 0.99 for the values of magnetic parameter B = ±0.1.
eter B and non-extremal spin a, the difference between Q = 0 and Q = Q W cases is negligible, see equations (23)(24). Therefore one can use the only, e.g., Q = Q W case, which is plotted in Fig. 1.
In the Newtonian theory of gravity the fundamental frequencies ν r , ν θ and ν φ of a particle moving around central object are equal to each other and the elliptical trajectory is the only possible trajectory for bound orbits. For an uncharged particle orbiting Schwarzschild black hole, we have ν r = ν θ = ν φ and hence the perihelion shift of bound quasi-elliptical trajectory can be observed.
Fundamental frequency difference in full GR theory allows construction of different HF QPOs geodesic models, where the ratio between fundamental frequencies become important. In the absence of magnetic field, the point where any two of fundamental frequencies could cross each other does not appear in the vicinity of a black hole. The frequencies coincide only in far distances from a black hole where the spacetime is asymptotically flat. But in the presence of magnetic field, the radial profiles of frequencies ν r , ν θ and ν φ cross each other in the vicinity of a black hole as demonstrated in Fig. 1. Concentrating our attention on the stable circular perturbations, we will further focus on small values of the magnetic parameter B < 1. But it is worth to note, that for large values of B >> 1, the radial profiles of fundamental frequencies in RLO and PLO motion practically ceases the change of their profiles with increasing B. Another interesting fact is that the charged particle motion in magnetized black hole background have many similarities with the "string-loop" motion which was described in [22,42,43].

IV. QPO MODELS AND INFLUENCE OF MAGNETIC FIELDS
The charged particle oscillations around circular orbits, studied in the previous section, suggest interesting astrophysical application, related to LF QPOs and HF QPOs observed in many Galactic Low Mass X-Ray Binaries (LMXB) containing neutron stars [6] or black holes [9,33]. We have selected three microquasars, GRS 1915+105, XTE 1550-564 and GRO 1655-40, where the central object is a black hole candidate and for which the masses M and spins a are estimated, see Tab. I. The limits on the values of the mass and the spin of the black hole using the spectral continuum fitting, which is independent of the QPO methods and does not include the effect of external large-scale magnetic fields, can be found in [34,36,37]. In the following subsection, we discuss the general constraint methods of black hole parameters from QPOs and the possible influence of external magnetic fields on the predictions of the rotation and mass parameters of black holes. In the later subsections we describe the general technique useful for the QPO fittings and particularly consider two most common QPO models, known as epicyclic resonance and relativistic precession models, modifying them accordingly by taking into account the effects of magnetic field.

A. General properties of QPOs and constraints of black hole parameters
The peaks of LF QPOs with frequencies f low and the twin peaks of the HF QPOs with upper f U and lower f L frequencies are sometimes observed in the Fourier power spectra. In the microquasars, i.e., LMXB systems containing a black hole, the twin HF QPOs appear at the fixed frequencies that usually have nearly exact 3 : 2 ratio [27]. The observed high frequencies are close to the orbital frequency of the marginally stable circular orbit representing the inner edge of the accretion disks orbiting black holes; therefore, the strong gravity effects are believed to be relevant for the explanation of HF QPOs [51].
The models of twin HF QPOs involving the orbital motion of matter around black hole can be generally separated into four classes: the hot spot models (the relativistic precession model and its variations [39,49], the tidal precession model [25]), resonance models [48,50,51] and disk oscillation (diskoseismic) models [29,35]. These models were applied to match the twin HF QPOs and the LF QPO for the microquasar GRO J1655-40 in [47]. Of course, the models can be applied also for intermediate massive black holes [44].
Unfortunately, none of the models recently discussed in literature, based on the frequencies of the harmonic geodesic epicyclic motion, is able to explain the HF QPOs in all three microquasars simultaneously, assuming that their central attractor is a black hole [52]. There is no generally accepted QPOs model for micorquasars, the observed HF QPO frequencies can be nearly fitted by the following heuristic formula (see Fig. 3) Presented formula is useful not only for the black holes with few solar masses (including three microquasars presented in our paper), but also for objects with completely different mass magnitudes such as active galactic nuclei [11,28,51]. Existence of a simple formula (29) nearly fitting the observational data implies the presence of the general mechanism in black hole physics related to the generation of QPOs. Moreover, if such a mechanism exists, it has to be almost independent (or weakly dependent) on the black hole spin, since, e.g., considered three microquasars have completely different spin estimations (see Tab. I), while can be fitted using the formula (29), independent of the spin parameter a. As one can see from Fig. 4, the straight dashed line corresponding to Eq. (29) fits all the three microquasars with the different values of the spin parameter a. It is quite challenging to find a single model within which one can fit the frequencies with predicted mass M and spin a parameters [11].
As it will be demonstrated in the next subsections, the presence of external uniform magnetic field can change the situation substantially allowing us to perform the fitting of all the three microquasar sources of HF QPOs. However, one has to point out that taking into account the effects of magnetic fields can potentially affect the predictions of the black hole spin and the mass, given in Tab. I. As it was mentioned above, the estimates of the parameters M and a given in the Table I are based on the fitting of the continuum spectrum which does not take into account the possible interaction of charged particles with external magnetic field. While the mass of the central object can be estimated using several different methods, e.g. observing the motion of the companion star, the prediction of the black hole spin parameter is quite sensitive on the position of the inner edge of an accretion disk, which is believed to be located close to the innermost stable circular orbit (ISCO) of test particle. Particularly, in the spectral continuum models the ISCO location is one of the most important characteristics influencing the spin parameter constraints. However, the ISCO position for the charged particle moving around black hole significantly depends on the magnetic field, as one can see in Fig. 2, or in Fig. 6 and 7 of our recent paper [53]. For instance, the presence of such a small magnetic parameter B = 0.05 in case of the Schwarzschild black holes, shifts the position of the ISCO radius considerably, decreasing it from 6M (for neutral particle) to 5M . This occurs due to the Lorentz force acting on a charged particle, which shifts the ISCO position towards the black hole horizon. For the large values of magnetic parameter B >> 1 the ISCO position can be shifted extremely close to the horizon radius of the black hole. Similar shift of the ISCO position occurs with increasing the spin parameter of the black hole. The ISCO radius coincides with the horizon radius in the extremely rotating black hole case, a = 1. This implies that the uniform magnetic field can in some sense mimic the black hole spin a in the sense that both increasing a and B shift ISCO to the black hole horizon for corotating particles, see Fig. 2. Thus, the implication of a magnetic field on the spectral continuum model could lead to the new predictions of black hole spin parameter a, shifting it to lower values, a mag < a nomag , than the one we have in the 'non-magnetic' spectral continuum fitting ( Table I).
The influence of strong magnetic fields (with B >> 1) on the profiled Fe spectral lines has been studied in [57] where the authors calculated the splitting of the profiled Fe spectral lines through the Zeeman effect. In case of a weak magnetic field, the modified profiled Fe spectral line originate from different angular velocity distribution inside the accretion disk, and the shift of the inner edge of the accretion disk towards the black hole horizon [17,18].  Table I: Observed twin HF QPO data for three microquasars and the restrictions on mass and spin of the black holes located in them, based on measurements independent of the HF QPO measurements given by the optical measerument for mass estimates and by the spectral continuum fitting for spin estimates [34,36]. Note that in the GRS 1915+105 microquasar more HF QPOs are observed [8]. Here we concentrate attention to the pair of HF QPOs demonstrating the frequency ratio 3:2, common with those observed in other two microquasars.

B. Resonant radii and the fitting technique
The HF QPOs come in pair of two peaks with upper f U and lower f L frequencies in the timing spectra. For the QPOs from black hole microquasars given in Tab. I, the frequency ratios f U : f L are very close to the fraction 3 : 2. Observation of this effect in different non-linear systems indicates the existence of the resonances between two modes of oscillations. In case of geodesic QPO models, the observed frequencies are associated with the different linear combinations of the particle fundamental frequencies ν r , ν θ and ν φ . In the presence of magnetic field the upper and lower frequencies of HF QPOs are the functions of magnetic parameter B, black hole mass M , spin a and the resonance position r, M, a, B), ν L = ν L (r, M, a, B).
It is worth to note that the frequencies ν U and ν L are inversely proportional to the mass M of a black hole, while the dependence of frequencies on the spin a and magnetic field B is more complicated and hidden inside Ω r , Ω θ , Ω φ functions, as given by equation (26). The resonant models of the twin HF QPOs assume a particular parametric or non-linear forced resonance of the oscillatory modes of the accretion disk [2]. In order to fit the frequencies observed in HF QPOs with the black hole parameters, one needs first to calculate the so called resonant radii r 3:2 (r 2:3 ) ν U (r 3:2 ) : ν L (r 3:2 ) = 3 : 2, ν U (r 2:3 ) : ν L (r 2:3 ) = 2 : 3.
(31) Resonant radii r 3:2 (or r 2:3 ) in general case are given as the numerical solution of higher order polynomial in r, for given values of spin a and magnetic field B parameters. Since the equation (31) is independent of the black hole mass explicitly, the resonant radius solution also has no explicit dependence on the black hole mass and techniques introduced in [49] can be used. Substituting the resonance radius into the equation (30), we get the frequencies ν U and ν L in terms of the black hole mass, spin and magnetic field.
The inverse dependence of the frequencies on the black hole mass allows us to write the mass M in terms of frequency, spin and magnetic field. The example of this dependence in the absence of magnetic field is shown in the Fig. 4. As it can be seen in the Fig. 4, in the absence of magnetic field, the calculated frequencies are quite sensitive on the black hole spin a, while the simple heuristic formula (29) given by dashed line is not dependent on black hole spin a.
In the following, we will examine magnetic field contribution to the two most relevant geodesic QPO models: epicyclic resonance and relativistic precession models.

C. Epicyclic resonance (ER) model
The epicyclic resonance (ER) models [1,51] consider resonance of axisymmetric oscillation modes of accretion discs. The frequency commensurability is thus crucial ingredient of the resonant models, and a particular case of this commensurability occurs for the parametric (internal) resonant phenomena that become strongest in the case of the 3 : 2 frequency ratio [49,51]. The simplest variant of the resonant model is the resonant epicyclic model where the two resonant modes are identified with the radial and vertical epicyclic frequencies [51].
Radial profiles for upper ν U and lower ν L frequencies for ER model are given in Fig. 5. In the non-magnetic case B = 0 (geodesic motion in Kerr spacetime), we have ν θ (r) > ν r (r) for any value of spin parameter a. For the magnetic case B = 0, the behaviour of the radial profiles of the ν U and ν L frequencies become more complicated, there exist two separate regions: for small values of radial coordinate r ISCO < r < r 1:1 we have ν U (r) > ν L (r); for larger values of radial coordinate r 1:1 < r < ∞, we have ν U (r) < ν L (r). In the absence of magnetic field B = 0 the radii r 1:1 where two frequencies coincide (ν U (r) = ν L (r)) is located at infinity r → ∞. For magnetic case B = 0, there exist two resonant radii r 3:2 and r 2:3 where one can produce the observed frequency ratios according to Tab. I. Thus the radii having importance for our study are lined up in the following way: r ISCO < r 3:2 < r 1:1 < r 2:3 . Fits of the three microquasar sources, HF QPOs modified ER model of charged particle oscillations in the field of magnetized black holes, are given in Fig. 6. Four different configurations of orbits (PALO, RLO, PLO and RALO) for both 3:2 and 2:3 resonances are given. In absence of the magnetic field, the ER model is not able to explain the observed HF QPOs in all the three microquasars, however, in the magnetic case it is possible -see magnetic field parameter B estimations in Tab. II. The co-rotating particle model (PALO,PLO) can fit the microquasar sources quite well, especially for 2:3 resonances, where we can fit all three microquasars with one value of magnetic field parameter B. Counter-rotating particle model (RLO,RALO) is successful in fitting for the RALO case only, the RLO case is completely falsified to explain the GRO 1655-40 source for any value of B, while the GRS 1915+105 source can be fitted by the RLO case only partially.
Increasing the magnetic field parameter B, the charged particle frequency will change. For the RLO and PLO cases, such changes become smaller and smaller with magnetic field B increasing -charged particle dynamics in magnetic field approaches the "string-loop" limit [22]. Moreover, when B is increased, the predicted frequencies become more independent on the black hole spin a -frequencies of charged particle oscillation become more influenced by the magnetic field than gravity for large B. The existence of new 2:3 resonance radii, at large radii, is also caused by the magnetic field presence. More detailed consequences of the fitting and discussion of results is given is Section V.

D. Relativistic precession (RP) model for HF QPOs
Relativistic precession (RP) model was introduced in order to explain the behaviour of the QPOs observed in X-ray neutron star and black hole binaries [30,38,40]. The RP model belongs to the class of geodesic QPO models, with the frequencies of ν U , ν L for HF QPOs defined as Radial profiles of the upper ν U and lower ν L frequencies for the RP model, and position of the related 3 : 2 resonance radii, are given in Fig. 7.
In the absence of magnetic field (B = 0), the RP model gives ν φ > ν r . When magnetic field is taken into account, the limit of the RP model for PALO/RALO configurations corresponds to the emergence of ν r = ν φ radii. So for large B there exist radii at which ν φ < ν r and the lower frequency ν L given by (33) in RP model will become negative and has to be redefined. For RLO/PLO cases, we have ν φ > ν r for any value of B parameter and the RP modes work properly.
The RP HF QPOs model modified by the presence of the magnetic field is able to match all observed HF    (29) and (33) for different values of magnetic field parameter B, for non-magnetic case B = 0 the lines are thicker. Only counter-rotating particle model RLO is successful in fitting all three sources.  The simultaneous observation of the twin HF QPOs and the LF QPOs (see Tab.I) in all three considered microquasars enables us to obtain the stringent restrictions on the mass and dimensionless spin of the central black hole, if we assume that all of these QPOs arise at a given radius of the accretion disc [30]. RP model for simultaneous HF and LF QPOs has frequencies ν U , ν L defined as Now, one can use the simultaneously observed values of HF and LF QPOs f U , f L , f low frequencies and identify them with formula for RP model (34) through three equations However, the restrictions on the black hole mass M and spin a given by the particle oscillations in the absence of magnetic field, B = 0, are in contradiction with the values obtained from the spectral fits for GRO J1655-40 and GRS 1915+105 source [46]. For B = 0, only the XTE 1550-564 source fits mass M and spin a restrictions from Tab. I. Note that in accord with the idea of [46], we could remove the discrepancy with fitting to the limits on the BH spin a given by the spectral measurements by breaking the assumption that the twin HF QPOs and the LF QPO arise at the same radius.
Taking magnetic field B into account, we are able to fit mass M of individual central object given in Tab. I. Then the spin a and magnetic field B can be estimated using (35). The results are summarised in Tab. IV, where we can observe that the RP model for HF and LF QPOs with magnetic field is predicting lower values of the spins a than the predictions due to the spectral fits.

V. DISCUSSION
We have studied the influence of an external weak magnetic field, approximated as asymptotically uniform with field lines oriented along the black hole rotation axis, on the motion of charged particle in its relation to the explanation of QPOs observed in microquasars. We have shown that the effect of even weak magnetic field on the motion of charged particles cannot be neglected and moreover, in the oscillatory phenomena of charged particles the influence of magnetic field is dominant over the influence of the black hole rotation. We have demonstrated that in order to fit the observational frequencies of HF QPOs with the epicyclic resonance model, the effect of external magnetic field can be sufficient.
In this article we assume the magnetic field to be uniform. Real magnetic fields around microquasar black holes and their accretion disks are far away from being completely regular and uniform, the Wald uniform magnetic field solution is used as a useful approximation. In the case of QPOs models, it is enough to assume uniformity of the magnetic field in the region where the oscillatory motion occurs. Fitting of the observed QPOs with the charged particle oscillations put limits on the values electron proton Fe+ charged dust B = 0.004 10 −5 G 0.02 G 1 G 10 8 G of the magnetic parameter B, which however contains, together with the field strength, also the specific charge of the oscillating test particles. This implies that in order to make proper estimation of the magnetic field values in the vicinity of particular microquasars, one needs to identify first the type of oscillating matter. In our approach, the "charged particle" can represent matter ranging from electron to some charged inhomogeneity orbiting in the innermost region of the accretion disk. The specific charges q/m for any of such structures will then range from the electron maximum to zero. Recalling the physical constants in the dimensionless magnetic parameter as B = |q|BGM/(2mc 4 ), we get the magnetic field strength in Gauss where the quantities are given in CGS units, see Tab. V. Assuming the oscillating particles to be electrons (as the main candidate source of synchrotron radiation producing X-rays), we get the estimations of the strength of magnetic field in the regions around the three microquasars as given in Table VI. Surprisingly, the range of the magnetic field magnitude estimates coincides with the measurements of the strength of large-scale Galactic magnetic fields. GRO   The values of magnetic parameter B sufficient to fit the observational frequencies of different sources are given in Tab. II. In all cases, one can predict the magnitude of magnetic field parameter to be of the order of B ∼ 10 −2 . Assuming that the QPOs appear due to the oscillations of electrons inside accretion disk one can estimate the magnetic field strength for a particular microquasar with the help of Eq. (36) and Tables I and II. The estimates of the strength of a magnetic field in the vicinity of the three microquasars are given in Table VI. The different class of orbits can correspond to the different configurations of an accretion disk (corotating or counterrotating) and different alignments of the magnetic field lines. As one can see from Table VI, the estimated values of the magnetic field in all cases for all microquasars with different configurations of accretion disks are of the same order, 10 −5 ÷ 10 −4 Gauss which suggests the interpretation of the origin of these fields as the Galactic magnetic field. There are many indications that the magnetic field must be present in the interstellar medium, and our estimates are in accord with these observations. There were also several attempts to estimate the strength of magnetic field in the vicinity of the three microquasars considered in the present paper. The estimates were using different methods for each source, independently on the QPO models. Here are few of them.
Magnitude of the equipartition magnetic field around the microquasar GRS 1915+105, has been estimated to be ∼ 8 G in [15], being based on the assumption that the repeated radio events are synchrotron in origin.
Fitting the flux data of the X-ray jets from the microquasar XTE J1550-564, with the trans-relativistic external shock model based on the analogy with the gammaray burst remnants [56], the authors suggested the magnetic field strength of the value 5×10 −4 G. Similar values have been confirmed in [21] from the observation of both jets of XTE J1550-564.
For the microquasar GRO J1655-40, the magnetic field strength has been estimated in [20], using the synchrotron-self-Compton model where the strength is obtained to be in the range between 0.05 G and 49 G, depending on the parameters of the model.
If one assumes the HF QPOs to be produced due to the oscillations of heavier than electrons particles, the estimates of magnetic field strength in our model will increase accordingly. For protons surrounding a black hole with 10 M ⊙ , the magnetic field strength with average value of the magnetic parameter B = 0.05 gives a value of 0.2 G. Some QPOs were detected in the modulations of an iron line of spectra. Thus, if one assumes that the QPOs come from the oscillations of iron atoms, one can get the magnetic field of 12 G, for the black hole with 10 M ⊙ and equipartition magnetic parameter B = 0.05. The physical interpretation of the X-ray emission from heavy particles, such as ions, can be realized by considering, e.g., the nuclear-Compton scattering processes. As one can see, the estimates of the strength of the magnetic field in our models do not disprove the independent mag-netic field measurements. In the case of LF QPO, which was simultaneously observed with HF QPOs, the fitting by the RP model gives the average magnetic parameter B = 0.005 for all sources. Using similar estimates way, as described above we get the strength of a magnetic field of an order of (1 − 10) × 10 −6 G assuming the synchrotron emission to be produced by electrons.
There are still some open questions and unresolved issues in the presented magnetic HF QPOs model. QPOs frequencies f U , f L observed in microquasars are notorious to be constant through time. They remain constant for different accretion disk regimes, where the accretion ratesṀ are assumed to be different. For magnetic HF QPOs models this means that the external uniform magnetic field (or et last B parameter) must also be constant through time. The stability of HF QPOs frequencies observed for black hole binaries, become apparent with comparison of HF QPOs observed in neutron star binaries -the neutron star HF QPO frequencies are substantially changing in time. Maybe, this tells something about the structure of magnetic field around neutron star/black hole objects. For neutron stars the magnetic field is assumed to be much stronger and complicated, having dipole character. On the other hand, our results indicate that for the black hole binaries the Galactic magnetic field could be relevant that can be expected nearly constant on year timescales relevant for QPO observations.

VI. CONCLUSIONS
The oscillatory frequencies of charged particles in black hole surroundings filled with asymptotically uniform magnetic field can be well related to the frequencies of QPOs observed in the microquasars GRS 1915+105, XTE 1550-564, GRO 1655-40. We can summarise main results for magnetic ER and RP models as follows: • Influence of magnetic field B on charged particle oscillatory motion can mimic the influence of the black hole rotation a for co-rotating particles -the ISCO can be shifted towards to the black hole horizon by magnetic field presence, and charged particle oscillatory frequencies can be increased due to the magnetic field, as shown in Fig. 1.
• Frequencies given by the magnetized versions of the QPOs geodesic models are dependent only slightly on the black hole spin a -the effect of weak magnetic field on the charged particle oscillations is stronger than the corresponding contribution of the black hole rotation. In the magnetized QPOs model, we obtained the formula ν ∼ 1/M , similar to the Eq. (29), and hence we can fit data from all the three considered microquasars with a single QPOs model.
• Contrary to the neutral particle motion, the charged particle motion around black hole immersed into magnetic field is driven by non-linear equations of motion. Without non-linearity there is no resonance between ν U and ν L modes of oscillation. Magnetized HF QPOs models will give not only the possibility to fit the values black hole spin a and M , but provide potentially also explanation of the asumed resonances.
• It was demonstrated that both the magnetized versions of the standard ER and RP models of twin HF QPOs can explain the observationally fixed data from the three microquasars, GRS 1915+105, XTE 1550-564 and GRO 1655-40, if we assume nearly equal (of the same order) magnitude of the external (near-uniform) magnetic field, and limits on the black hole mass M (and dimensional spin a) given by independent optical (spectral continuum) measurements. The assumed magnetic field magnitude is in agreement with estimates of the Galactic magnetic field for the electron case.
• The magnetic ER model admits both the corotating and counter-rotating disks, while the RP model admits counter-rotating disks only. Moreover, in the RP model containing naturally also the explanation of the simultaneously observed LF QPO, we have found a contradiction with the limits on spin in the microquasars GRS 1915+105 and GRO 1655-40, if the creation of both HF QPOs and the related LF QPOs occurs at the same radius. As proposed in [46], this controversy could be solved by breaking assumption of common location of the twin HF QPOs and LF QPOs. Problem of lower predicted black hole spins a, given by the RP model with the LF QPOs, could be also solved, assuming magnetic corrections to the spectral continuum model, which can decrease limits the on the black hole spin.
• In the magnetized RE or RP models related to the twin HF QPOs, the uniform magnetic field is necessary ingredient only in the vicinity of the black hole horizon where the epicyclic oscillations have to occur. Therefore, the state of the accretion disk should be relatively regular, corresponding to quiescent disks with sufficiently low accretion flow, allowing for existence of nearly uniform magnetic field near the horizon. This really corresponds to the so called hard states where the HF QPOs are observed [14]. In the so called soft states, related to high accretion rates and relatively irregular disks, no HF QPOs are observed, thus we can expect that the uniformity of the magnetic field in the horizon vicinity becomes strongly violated by the strong irregular accretion flow.
The charged test particle oscillatory frequencies are sensitive to any magnetic field presence. The applicability of the magnetized HF QPO ER and RP models for QPOs observed in completely different sources allows us to conclude that the proposed models with inclusion of the effect of magnetic field can be considered as one of the possible explanations of the HF QPOs occurring in the field of microquasars. Assuming the main source of synchrotron radiation producing X-rays are the relativistic electrons, we estimate the magnetic field in the vicinity of the black hole in the three sources to be of order 10 −5 G which can serve as possible signature of the Galactic magnetic field magnitude. For heavier particles (ions, charged hotspots) much larger magnetic fields are necessary for fitting the data. The model also gives limits on character of oscillating matter and intenzity of the magnetic field.