Possible signature of the magnetic fields related to quasi-periodic oscillations 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 the strong gravity regime. 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 for three microquasars, GRS 1915+105, XTE 1550-564 and GRO 1655-40, containing black holes, for magnetized versions of both epicyclic resonance and relativistic precession models and discuss the corresponding constraints of parameters of the model, which are the mass and spin of the black hole and the parameter related to the external magnetic field. The estimated magnetic field intensity strongly depends on the type of objects giving the observed HF QPOs. It can be as small as 10-5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-5}$$\end{document} G if electron oscillatory motion is relevant, but it can be by many orders higher for protons or ions (0.02–1 G), or even higher for charged dust or such exotic objects as lighting balls, etc. On the other hand, if we know by any means the magnetic field intensity, our model implies strong limit on the character of the oscillating matter, namely its specific charge.


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's rotation axis. Due to friction, matter of the accretion disk becomes hot and emits electromagnetic a e-mail: martin.kolos@fpf.slu.cz b e-mail: arman.tursunov@fpf.slu.cz c e-mail: zdenek.stuchlik@fpf.slu.cz radiation, also in X-rays in the 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 useful information regarding the range of parameters of the system [34]. In this connection, the binary systems containing black holes, as compared to neutron star systems, seem to be promising for 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 QPOs with high precision allow us to get useful knowledge about the central object and its background. According to the observed frequencies of QPOs, which cover the range from few mHz up to 0.5 kHz, different types of QPOs were distinguished. Mainly, these are the high frequency (HF) and low frequency (LF) QPOs with frequencies up to 500 Hz and up to 30 Hz, respectively. The HF QPOs in black hole microquasars are usually detected with the twin peaks which have a frequency ratio close to 3:2. This is the case of the Galactic black hole microquasars GRS 1915+105, XTE 1550-564 and GRO J1655-40, which we consider here.
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, disko-seismic 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 the 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 for several reasons, as described below.
Magnetic fields have been detected and measured in nearly all celestial objects. The existence of an accretion disk around a 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 stellarmass 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 though the total charge of an accretion disk is believed to be zero, the presence of a 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 the 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]. Recent studies of the structure of the magnetic fields in galaxies [7] clearly demonstrate the formation of the spiral patterns in almost all considered galaxies, including the Milky Way. Therefore, black holes can also be immersed in the external, large-scale magnetic fields that can have a 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 the 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 in 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 by 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 the specific charge.
The influence of the magnetic field on the charged test particle motion around the black holes has been already studied in the 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 par-ticle 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 X-ray 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 the observed HF QPOs have been studied 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.

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 with the nonzero components of the metric tensor g μν taking in the standard Boyer-Lindquist coordinates the form 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

Uniform external magnetic field
The external axially symmetric and asymptotically homogeneous magnetic field implies the electromagnetic field potential which is a 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 ≡ 2a M 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 the form A α = 1 2 Bξ α (φ) , where B denotes the intensity of the asymptotically uniform magnetic field; ξ (φ) , ξ (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]: 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 solarmass 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 four 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 brevity M = 1, the Wald charge corresponds to Q W = 2aB. The parameter B will be called the magnetic parameter [53]. The estimation of the magnitude of the magnetic parameter shows that the effect of even a 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 under the same conditions, the strength of magnetic the field B has to be 1836 times stronger. A 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 Sect. 4.

Circular obits of charged particles
The motion of a charged particle with charge q and mass m in curved spacetime in the presence of an 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 it 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 two 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 which is obviously independent of the electromagnetic parameters B and Q. Equations (12) and (13) allow us to find two expressions for two unknown quantities u t and u φ . The explicit form of the analytical expressions for u t and u φ cannot be represented in a simple form, however, we can easily solve it numerically. The 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 the 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 equations of motion can then be 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 a magnetized rotating black hole background was studied in detail in [53], where we used 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. The stability of the circular orbits is also discussed in [53] -the harmonic quasi-periodic oscillations are allowed only around the stable circular orbits.

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 Substituting the deviation vector into Eq. (11) and using the first order of an 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 Eqs. (12) and (13). The frequencies for the case with the zero induced charge will be similar to Eqs. (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 Ref. [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 the socalled Larmor angular frequency, Ω L , associated with the pure contribution of an external uniform magnetic field, 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 at large distances from the black hole where the influence of the uniform magnetic field becomes crucial.
The expressions for the fundamental frequencies (20), (21) and (25) are given in dimensionless form. In physical units, one needs to extend the corresponding formulas by the factor c 3 /G M. 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.

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 with a positive charge one can distinguish the following orbits: 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.
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 the case of non-rotating Schwarzschild black holes, the 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 the orbits as given above, we plot the radial profiles of the frequencies ν θ , ν r and ν φ of small harmonic oscillations of charged particle measured by a distant static observer in Fig. 1, for characteristic values of the black hole spin a and magnetic parameter B. The frequencies have Fig. 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 non-magnetic case (last two figures). The second row represents the Kerr black hole with a = 0.7 with magnetic parameter B = ± 0.1. Since the differences between the cases Q = 0 and Q = Q W are negligible, see the discussion of Eqs. (23 and 24), the case Q = Q W is considered only. The last row represents the effect of the nearly extremal Kerr black hole with a = 0.99 for the values of magnetic parameter B = ± 0.1 been calculated using Eq. (26) in the case of small displacements from the equatorial plane. The magnetic field 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 the 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 the ISCO for magnetic parameter B = ± 0.05, in dependence on the black hole spin a, is plotted in Fig. 2. A 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 Eq. (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 magnetic parameter B. For PLO/RLO configurations, with attractive Lorentz force, the frequency ν φ increases with increasing magnetic parameter B. The radial profile of the vertical frequency ν θ does not have a strong dependence on the 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 a Kerr black hole with M = 10 M mass. The change of the mass of the black hole leads to the change only in the 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 parameter B and non-extremal spin a, the difference between the cases Q = 0 and Q = Q W is negligible; see Eqs. (23 and 24). Therefore one can use only the case, e.g., Q = Q W , which is plotted in Fig. 1.
In the Newtonian theory of gravity the fundamental frequencies ν r , ν θ and ν φ of a particle moving around a central object are equal, and the elliptical trajectory is the only possible trajectory for bound orbits. For an uncharged particle orbiting a Schwarzschild black hole, we have ν r = ν θ = ν φ and hence a perihelion shift of the bound quasi-elliptical trajectory can be observed. The fundamental frequency difference in full GR theory allows for the construction of different HF QPOs geodesic models, where the ratio between fundamental frequencies becomes important. In the absence of the magnetic field, the point where any two of fundamental frequencies could cross does not appear in the vicinity of a black hole. The frequencies coincide only at far distances from a black hole, where the spacetime is asymptotically flat. But in the presence of a magnetic field, the radial profiles of frequencies ν r , ν θ and Table 1 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 measurement 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 on the pair of HF QPOs demonstrating the frequency ratio 3:2, common with those observed in other two microquasars 1, the radial profiles of fundamental frequencies in RLO and PLO motion practically cease the change of their profiles with increasing B. Another interesting fact is that the charged particle motion in magnetized black hole background has many similarities with the "string-loop" motion, which was described in [22,42,43].

QPO models and influence of the 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 Table 1. 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 the magnetic field. 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 the 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 the 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): The formula presented 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 Fig. 4 Fitting the observed microquasar parameters with epicyclic resonance and relativistic precession models of HF QPOs (solid lines) in the absence of the magnetic field. As can be seen, none of the most common HF QPOs models is able to explain the observed HF QPOs for all the three microquasars. The plus sign is used for co-rotating, while the minus sign is used for counter-rotating orbits. On the horizontal axis the black hole spin a is given, while on the vertical axis we give the black hole mass M divided by observed upper HF QPOs frequency ν U . Such a mass scaling is using the ν U ∼ 1/M relation (see Eq. (29)) and allows us to compare all three microquasars in one figure. Black boxes are predicted values of black hole mass M and spin a, see Table 1; the thick dashed line is given by the heuristic equation (29) [11, 28,51]. The existence of a simple formula, Eq. (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., the three microquasars considered have completely different spin estimations (see Table  1), while they can be fitted using Eq. (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 the predicted mass M and spin a parameters [11]. As will be demonstrated in the next subsections, the presence of an 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 the magnetic fields potentially affects the predictions of the black hole spin and the mass, given in Table 1. As was mentioned above, the estimates of the parameters M and a given in Table  1 are based on the fitting of the continuum spectrum, which does not take into account the possible interaction of charged particles with the 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 to the position of the inner edge of the accretion disk, which is believed to be located close to the innermost stable circular orbit (ISCO) of the 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 the black hole significantly depends on the magnetic field, as one can see in Fig. 2, or in Figs. 6 and 7 of Ref. [53]. For instance, the presence of such a small magnetic parameter as B = 0.05 in the case of the Schwarzschild black holes shifts the position of the ISCO radius considerably, decreasing it from 6M (for a 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 large values of the magnetic parameter, B 1, the ISCO position can be shifted extremely close to the horizon radius of the black hole. A similar shift of the ISCO position occurs with increasing 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 the ISCO to the black hole horizon for co-rotating 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 1).
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 the 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].

Resonant radii and the fitting technique
The HF QPOs come in pairs 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 Table 1, 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 the 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 a magnetic field the upper and lower frequencies of HF QPOs are functions of the magnetic parameter B, the black hole mass M, the 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 depen-dence of the frequencies on the spin a and magnetic field B is more complicated and hidden inside Ω r , Ω θ , Ω φ functions, as given by Eq. (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 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 Fig. 4. As can be seen in Fig. 4, in the absence of the magnetic field, the calculated frequencies are quite sensitive on the black hole spin a, while the simple heuristic formula of Eq. (29) given by a 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.

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 the spin parameter a. For the magnetic case B = 0, the behavior of the radial profiles of the ν U and ν L frequencies become more complicated, there exist two separate regions: for small values of the radial coordinate = 10.7 while r 2:3 → ∞. Comparing the first row (plotted for B = 0.01) with the second row (plotted for B = 0.1) we can see that the upper ν U and the lower ν L resonance frequencies (at resonance radii) are increasing with increasing the magnetic parameter B. The lower frequency ν L is more influenced by magnetic field presence than the upper one ν L , but this appears due to the influence of the Lorentz force in radial direction r ISCO < r < r 1:1 we have ν U (r ) > ν L (r ); for larger values of the radial coordinate r 1:1 < r < ∞, we have ν U (r ) < ν L (r ). In the absence of the 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 Table 1. 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 the 3:2 and the 2:3 resonances are given. In the 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 the estimations of the magnetic field parameter B in Table 2. The co-rotating particle model (PALO, PLO) can fit the microquasar sources quite well, especially for the 2:3 resonances, where we can fit all three microquasars with one value of the magnetic field parameter B. The counter-rotating particle model (RLO, RALO) is successful in fitting for the RALO case only, the RLO case is completely falsified as regards explaining 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 -the charged particle dynamics in a 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 -the frequencies of the 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 presence of the magnetic field. More detailed consequences of the fitting and discussion of results are given is Sect. 5.

Relativistic precession (RP) model for HF QPOs
A relativistic precession (RP) model was introduced in order to explain the behavior 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  The radial profiles of the upper ν U and lower ν L frequencies for the RP model and the position of the related 3:2 resonance radii are given in Fig. 7.
In the absence of the magnetic field (B = 0), the RP model gives ν φ > ν r . When the magnetic field is taken into account, the limit of the RP model for the 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 the RP model will become negative and has to be redefined. For the RLO/PLO cases, we have ν φ > ν r for any value of the B parameter and the RP modes work properly (Fig. 7).
The RP HF QPOs model modified by the presence of the magnetic field is able to match all observed HF frequencies with estimated black hole mass M and spin a, see Table 1, only in the 3:2 RLO case; see Fig. 8. In the PALO/RALO cases, the RP model will not work properly for large values of the magnetic parameter with B > 0.1, while for the PLO case the magnetic field in the RP model increases the frequencies to values larger than observed. In the RLO scenario, the counter-rotating disk around the central black hole is influenced by the magnetic field with strength B ∼ 0.003. Values of B for every individual microquasars are given in Table 3.

RP model assuming simultaneous observations of HF and LF QPOs
The simultaneous observation of the twin HF QPOs and the LF QPOs (see Table 1) in all three considered microquasars  Table 1. The lines are given by Eqs. (29) and (33) for different values of the magnetic field parameter B, for the non-magnetic case B = 0 the lines are thicker. Only the counter-rotating particle model RLO is successful in fitting all three sources 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]. The 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 the three equations for three independent parameters -the black hole mass M, spin a and magnetic parameter B. For different values of the magnetic parameter B, different combinations of the mass M and spin a can be obtained; see the results in Fig. 9. Only for co-rotating particles (PALO B > 0 and PLO B < 0) the real solution of Eq. (35) can be found. However, the restrictions on the black hole mass M and spin a given by the particle oscillations in the absence of the 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 Table 1. 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 the magnetic field B into account, we are able to fit the mass M of the individual central object given in Table  1. Then the spin a and magnetic field B can be estimated using (35). The results are summarised in Table 4, 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.

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 a charged particle in its relation to the explanation of QPOs observed in microquasars. We have shown that the effect of even a weak magnetic field on the motion of charged particles cannot be neglected and, moreover, in the oscillatory phenomena of charged particles the influence of the 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 the 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 puts limits on the values 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 a 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|BG M/(2mc 4 ), we get the magnetic field strength in Gauss, where the quantities are given in CGS units; see Table 5.
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 the magnetic field in the regions around the three microquasars as given in Table 6. Surprisingly, the range of the magnetic field magnitude estimates coincides with the measurements of the strength of large-scale Galactic magnetic fields. The values of magnetic parameter B sufficient to fit the observational frequencies of different sources are given in Table 2. In all cases, one can predict the magnitude of the magnetic field parameter to be of the order of B ∼ 10 −2 . Assuming that the QPOs appear due to the oscillations of electrons inside the accretion disk one can estimate the magnetic field strength for a particular microquasar with the help of Eq. (36) and Tables 1 and 2. The estimates of the strength of a magnetic field in the vicinity of the three microquasars are given in Table 6. The different class of orbits can correspond to the different configurations of an accretion disk (co-rotating or counter-rotating) and different alignments of the magnetic field lines. As one can see from Table 6, 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 G, 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 the 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 a few of them.
The magnitude of the equipartition magnetic field around the microquasar GRS 1915+105 has been estimated to be ∼ 8 G in [15], based on the assumption that the repeated radio events are of synchrotron 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 gamma-ray burst remnants [56], the authors suggested a 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 synchrotronself-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 particles than electrons, the estimates of the 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 magnetic 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 to 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 for being 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 at least the B parameter) must also be constant through time. The stability of HF QPOs frequencies observed for black hole binaries becomes 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 the magnetic field around neutron star/black hole objects. For neutron stars the magnetic field is assumed to be much stronger and complicated, having a dipole character. On the other hand, our results indicate that for the black hole binaries the Galactic magnetic field could be relevant, which can be expected to be nearly constant on the timescales of years, relevant for QPO observations.

Conclusions
The oscillatory frequencies of charged particles in black hole surroundings filled with an 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 the main results for magnetic ER and RP models as follows: -The influence of the magnetic field B on the charged particle oscillatory motion can mimic the influence of the black hole rotation a for co-rotating particles -the ISCO can be shifted towards the black hole horizon by the presence of the magnetic field, and charged particle oscillatory frequencies can be increased due to the magnetic field, as shown in Fig. 1. -The frequencies given by the magnetized versions of the QPOs geodesic models are only slightly dependent on the black hole spin a -the effect of the 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 Eq. (29), and hence we can fit data from all three considered microquasars with a single QPOs model. -Contrary to the neutral particle motion, the charged particle motion around the black hole immersed in the magnetic field is driven by non-linear equations of motion. Without non-linearity there is no resonance between the ν U and ν L modes of oscillation. Magnetized HF QPOs models will give not only the possibility to fit the values of the black hole spin a and M, but they provide potentially also an explanation of the assumed resonances. -It was demonstrated that the magnetized versions of both the standard ER and the 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) magnitudes 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 co-rotating and counter-rotating disks, while the RP model admits counterrotating 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 the 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 the assumption of a 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 also be solved, assuming magnetic corrections to the spectral continuum model, which can decrease the limits on the black hole spin. -In the magnetized RE or RP models related to the twin HF QPOs, the uniform magnetic field is a 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 the existence of a 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 vicinity of the horizon 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 the 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 Xrays to be 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 a 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 the character of oscillating matter and intensity of the magnetic field.