Testing Ghasemi-Nodehi-Bambi metric parameters with quasi-periodic oscillations

QPOs are seen as peak features in the X-ray power spectral density of stellar mass black holes and neutron stars, of which frequencies can be measured with high precision. These QPO frequencies are believed to be related to fundamental frequencies of test particles, which are mainly determined by the background metric. We consider the metric introduced in Ghasemi-Nodehi and Bambi (Eur Phys J C 76:290, 2016). The fundamental frequencies in this metric are sensitive to some of the metric parameters but insensitive to other parameters, which means, the differences in fundamental frequencies in this metric and those in the Kerr ones can be significant for small changes of some but not all parameters around the Kerr value. By comparing with the QPO observations of GRO J1655-40, we find that only one parameter of the Ghasemi-Nodehi-Bambi metric can be strongly constrained, but other parameters cannot. We also use nested algorithm to investigate whether better constraints on the Ghasemi-Nodehi-Bambi metric parameters can be obtained from QPO observations of multiple objects by simulations. We find that four parameters can be strongly constrained while other parameters cannot . Our results suggest that QPOs may be important tools for testing the Kerr metric.


Introduction
Einstein theory of General Relativity (GR) was born more than 100 years ago and successfully passed test in the weak gravity field regime [1,2,3,4,5]. According to GR, astrophysical black holes (BHs) are described by the Kerr metric with only two parameters, i.e., mass and spin [6,7]. However, the nature of the Kerr BH is still to be verified. Important test of GR in the strong gravity field regime can be from direct observational confirmation of the Kerr metric, via either electromagnetic wave [10,11] or gravitational wave [8,9] observations. The electromagnetic tests in the literature mainly make use of the iron line method [30] and/or continuum fitting method [31] by utilizing X-ray data [28,29]. Black hole shadow also offers a unique tool to probe the event horizon and nature of BHs, which has been demonstrated by recent observations of M87 through the Event Horizon Telescope (EHT) [32,33].
There have also been a number of attempts to generalize the Kerr solution by phenomenological parameterization [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. The idea is mostly to write a more general metric that includes the Kerr one as its special case, similar to the PPN of the Schwarzschild metric [2]. Generally, parameters of a email: mghasemin@ipm.ir the non-Kerr metric may be degenerated in testing strong gravity fields. Therefore, the Kerr metric cannot be verified only by the detection of its characteristic properties, and it would be proofed only when all deformation parameters of the non-Kerr metric are required to vanish by observations.
Quasi-periodic oscillations (QPOs), peak features in the X-ray power spectral density of stellar mass BHs and neutron stars, are recently recognized as promising tools to test GR in the strong gravity field regime [38,39,40,41,42,43,44]. QPOs have been detected in four systems GRO J1655-40, XTE J1550-564, GRS 1915+105, and H1743-322, they are not sensitive to the properties of accretion flows. The exact physical mechanism responsible for the production of the QPOs is not clear, though a number of different scenarios have been proposed, including relativistic precession model [34], diskoseismology models [35], resonance models [36], and p-mode oscillations of an accretion torus [37]. In most scenarios, QPO frequencies are related to the characteristic orbital frequencies of test particles, which are determined only by the BH metric. The main advantage of using QPOs with respect to other methods is that the QPO frequencies can be measured with high accuracy. Thus, QPOs can provide a unique tool to probe the spacetime and geometry around the central compact objects.
One of us proposed a new parameterization to Kerr metric in Reference [26] (hereafter the GB metric). The Kerr case would be recovered by setting all the parameters involved in the GB metric equal to 1. Any parameter different from 1 indicates a deviation from the Kerr case. Assuming the GB metric, we have studied the BH shadow in the same paper and X-ray reflection spectroscopy in [45]. BH shadow studies show that only 4 parameters involved in the GB metric have weak impact on the shadow shape while other parameters have almost no impact, and thus it will be difficult to constrain the GB metric by using BH shadow measurements. In the paper [45], one of us studied the X-ray reflection spectroscopy resulting from accretion systems to constrain the parameters in the GB metric. Our results show that all parameters except one can be potentially well constrained with near future X-ray missions. In the present paper, we investigate the impact of GB parameters of Ref. [26] on QPOs. Our results show significant parameter degeneracy in constraining the GB parameters.
The content of this paper is as follows. In Section 2, we provide the theoretical framework of our calculations. Results and discussions are written in Section 3. Summary and conclusions are presented in section 5. Throughout the paper, we use units in which G N = c = 1.

Theoretical Framework
Three fundamental frequencies characterize the equatorial orbit of a test particle rotating around a BH, i.e., the Keplerian frequency ν φ , radial epicyclic frequency ν r , and vertical epicyclic frequency ν θ . These frequencies only depend on the BH metric and radius of orbit.
The general case can be studied by the line element that is for a generic stationary and axisymmetric spacetime and the metric components are independent of t and φ.
There are two conserved quantities associated to t and φ, specific energy at infinity, E, and z-component of specific angular momentum at infinity, L z . The t and φ components of 4-velocity of a test particle are [46] One may introduce V eff from the conservation of the rest mass g µνẋ µẋν = −1, where Geodesic equation can also be written as Circular orbit requires thatṙ =θ =r = 0, so the radial component is This equation gives us the orbital angular velocity where the sign is (+/−) for corotating/counterrotating orbits and orbital frequency can be calculated as We find E and L z since −E = g ttṫ + g tφφ and L z = g tφṫ + g φφφ , Substituting E and L z in equation (4) together with metric component gives us the V eff . From equation (4) we consider small perturbation around circular orbit along radial and vertical directions at linear order. If δ r and δ θ are the small displacements around the mean orbit, r = r 0 + δ r and θ = π/2 + δ θ , and neglecting terms O(δ 2 r ) and O(δ 2 θ ), we have following differential equations where The radial epicyclic frequency and the vertical epicyclic frequency are, respectively, ν r = Ω r /2π and ν θ = Ω θ /2π. In this paper, we consider relativistic precession model that recent studies seem to support it [47,48]. This model first proposed to explain the frequency of the QPOs in power spectral density of neutron star and then is extended to the QPOs of the stellar mass BH [34]. This model simply relates the observed frequencies with the orbital motion of a test particle. One can find the periastron precession frequency, ν p , and nodal precession frequency from fundamental frequencies, These frequencies can be measured observationally and compared with the fundamental frequencies computed from the background metric. This provides a unique opportunity with high accuracy to probe the strong gravity regime and test possible deviations from the Kerr solution of GR.
In Ref [26] one of us introduced 11 parameters, b i as follows, Kerr case recovers for all b i = 1. Any differences from b i = 1 show deviations from general relativity. b 1 can be set as one because it is coefficient of mass, b 1 M and also b 3 = 1 in the same way asymptotic specific angular momentum is b 3 a. We also have b 6 close to one from solar system experiment. So we set b 1 = b 3 = b 6 = 1 and we do not consider these three parameters in our QPO calculations.

Results
QPOs showing in the X-ray power spectral density of stellar mass BHs and neutron stars provide precise information about the background metric. QPOs are frequencies measured for a BH system are believed to be determined by the BH metric but not the accretion flow, though it may model dependent. In this paper, we adopt the relativistic precession model and assume that QPO frequencies are related to the characteristic orbital frequencies of test particles. These characteristic frequencies are solely determined by the BH metric. We also assume that the GB metric describes the BH metric, for which the Kerr metric is a special case. We first calculate frequencies ν φ , ν r , ν θ , ν p , and ν n for the GB metric with different settings of the metric parameters. Figures 1 and 2 show ν φ , ν p , and ν n as a function of radius r for GB BHs with mass 5.4M , dimensionless spin parameter a * = a/M = 0.9, but different b i (not too much different from the Kerr case b i = 1). As seen from these two Figures, different settings of b 2 or b 4 (e.g., b 2 or b 4 = 1.2 and 5) may lead to substantial difference in the resulting ν n but negligible differences in ν φ and ν p (top left and right panels of Fig. 1); different settings of b 7 , b 9 , or b 10 (e.g., 1.2 and 5) may lead to significant difference in the resulting ν p but not ν φ and ν n (bottom-right panel of In the Kerr metric, the three fundamental frequencies can be derived analytically as, The Schwarzschild case would recover by imposing a = 0, i.e., In the Kerr metric ν θ ≥ ν r , and in the Schwarzschild metric ν φ = ν θ > ν r . In our calculation for b i = 1.2, similar to the Kerr case and Schwarzschild case all ν θ s are greater than ν r , but for the case b 2 = 5, ν r ≥ ν θ . This is possible, for example, if the ISCO is marginally vertically stable, ν θ = 0 at the ISCO. Thus, we have ν r > ν θ . In order to show differences with respect to the Kerr case, we plot Figure 3. This figure shows the impact of different GB parameters on the shape of ν φ (top left panel), ν p (top right panel), ν n (bottom panel). For all cases, the mass and dimensionless spin parameter are set to 5.4 M and 0.9, respectively. It is obvious that ν φ curves for GB metric with b i = 5 (i = 2, 4, 5, 7, 8, 9, 10, and 11, respectively) are almost the same as that for the Kerr case with all b i = 1 (top-left panel), while ν p curves for the cases with b 7 = 5, b 9 = 5, or b 10 = 5 are significantly different from those for other cases (almost identical; top-right panel), and ν n curves for the cases with b 2 = 5, b 4 = 5, b 9 = 5, or b 10 = 5 are significantly different from those for other cases (almost identical; bottom panel).
According to Figures 1, 2 and 3, we conclude that the GB metric can be distinguished from the Kerr one if any of the GB parameters b 2 , b 9 , b 10 is substantially different from 1. The reason is that some of the resulting fundamental frequencies from such a GB metric are substantial different from those of the Kerr one. The GB parameters b 4 and b 7 may also be constrained but with less significance since the differences in the fundamental frequencies are smaller.
We note here that b 11 may also be constrained since the QPO frequencies may be only significantly affected by b 11 when BH spin is extremely high, which is discussed later in Section 3.2 and can be seen from Figure 5.
Other GB parameters like b 5 and b 8 cannot be constrained because of the negligible differences in fundamental frequencies for different settings of those b i .

QPO observations for GRO J1655-40
Rossi X-ray Timing Explorer (RXTE) mission observations of GRO J1655-40 have revealed that it has two high frequency QPOs and one low frequency QPO [47,49]. These frequencies are measured with high accuracies, i.e., ν 1 = 441 +2 −2 Hz, ν 2 = 298 +4 −4 Hz, and ν 3 = 17.3 +0.1 −0.1 Hz. If we assume that these frequencies correspond to ν φ , ν p , and ν n , respectively, and the central BH is described by the Kerr metric, these frequencies should be determined by three unknown parameters M , a, and r as seen from Equations (19)- (21). Therefore, strong constraints should be able to obtain on the mass M and dimensionless spin a of the BH.
If the BH metric in GRO J1655-40 is described by the GB metric rather than the Kerr metric, the determined QPO frequencies may be used to put constraint on the GB parameters b i s. For simplicity, we assume that only one of the b i s is different from that of the Kerr case, i.e., b i = 1 and b j =1 with j = i, in order to get a meaningful constraint. For each given i (2, 4, 5, 7, 8, 9, 10, or 11), we calculate χ 2 as follows to obtain a constraint on the b i where ν φ , ν n , and ν p are, respectively, the orbital, nodal, and periastron frequency. The numbers 441, 298, and 17.3 are the QPO frequencies measured from the X-ray observations of GRO J1655-40, and 2, 4, and 0.1 are the corresponding 1-σ errors σ i here. We set M = 5.4 M here, which is obtained from independent optical observations [50]. Here we adopt the χ 2 statistics to obtain such constraints. We first minimize χ 2 here over parameter r and then fix r to obtain constraint on b i . 1 Figure 4 shows the contour levels for parameters b 2 (top-left panel), b 9 (topright panel), and b 10 (bottom panel), respectively. As seen from this Figure, b 2 can be well constrained to a value of 2.2 +0.395 −0.523 . The minimum of χ 2 is 0.02 and ∆χ 2 for 1σ, 2σ, and 3σ confidence levels are 1, 4, and 9, respectively. The fitting results suggest that the BH in GRO J1655-40 may be better described by the GB metric with b 2 = 1, but the 1 Set r = rISCO, we may get a much worse fitting. significance is less than 2σ and a Kerr metric is also compatible with the data. This figure also shows degeneracies between b 9 and a and between b 10 and a.

Mock QPO observations
The QPO observations of a single object have already enabled some constraints on the GB metric. If there were QPO observations for many more objects, it may be possible to get better constraints. In this subsection, we investigate whether better constraints can be obtained from QPO measurements of multiple objects and whether the GB metric can be distinguished from the Kerr metric. To do this, we first generate 10 Kerr BHs, each with a randomly assigned mass and spin. We also generate another mock GB BH sample, each with a randomly assigned mass and spin but a fixed b 2 = 5. For both mock samples, we calculate ν φ , ν p , and ν n for each sample object. We assume that these frequencies are correspondingly the QPO frequencies that can be measured with typical errors of 1%, similar to that of GRO J1655-40, and thus we get mock QPO observations. We adopt the Bayesian method to investigate whether whether b 2 can be strongly constrained as follows.
For a given set of n observations d, from Bayes' theorem: For QPO observations of each object j, χ 2 j can be obtained similarly as before: where ν j φ , ν j p , and ν j n are the frequencies from the GB metric with given b i , (ν j φo , ν j po , ν j no ) and (σ j ν φo , σ j νp o , σ j νn o ) are the mock frequencies and its corresponding uncertainties, respectively.
We first mock the observation data by assuming the Kerr metric, and fit the data with the generalized GB metric. We only vary one of the GB parameters b i each time. For each of those cases that only one b i is considered to be free, there are totally 21 parameters θ = {b i , M 1 , ...M 10 , a 1 , ..., a 10 } in the model.
In Mock data production, we also assume r = r ISCO , for simplicity. We sample the posterior with nested sampling algorithm using dynesty [51]. The prior for all the parameters are set to be flat with b i ∈ [0, 7] 2 , M j ∈ [4,11]M and a j ∈ [0.01, 0.99] with j = 1, ..., 10.    The result shows that b 2 , b 9 , and b 10 can be well constrained as expected from previous analysis shown by Figure 1, 2, and 3. Note that b 7 cannot be well constrained due to little frequency difference at ISCO induced by the change of b 7 .
However, b 4 also cannot be well constrained though significant differences in ν can be found for cases with different b 4 (see Fig. 1), which may be due to strong degeneracies with spin a and mass M .
The constraint on b 11 is quite good, despite the small frequencies difference induced by the change of b 11 shown in Figure 1. The reason may be seen from Figure 5, which shows the impact of b 10 (left panel) and b 11 (right panel) on QPO frequencies at ISCO as a function of the BH spin. Apparently, the impact of b 10 gradually increases with increasing spin, while the impact of b 11 is significant only when the BH spin is high (> 0.9). and become comparable with b 10 at high spin. In the mock QPO data, there is only one mock BH with a = 0.93, which is main contributor to the constraint on b 11 . As we see for other parameters like range of bi can be adopted, which may give a looser constraint on bi compared to the present one. b 5 in Fig. 1 and b 8 in Fig. 2 are overlap with each other and it is hard to constrain.
To summarize, we have shown that the GB parameters b 2 , b 9 , b 10 , and b 11 can be well constrained if using QPO frequency measurements with a precision of 1% for even only ten mock BH systems. For such a case, the reconstructed values of these b i parameters, and the mass and spin of each mock BH are listed in Table 1. If the number of measurements is larger, then better constraints may be obtained. One should note that the QPO frequencies are assumed to be originated at r ISCO when generating the mock data and performing the fitting.
In realistic cases, the QPO phenomenon may not occur exactly at r ISCO , and this uncertainty may reduce the constraint power. As a comparison, in our study of GRO J1655-40 above, we minimized χ 2 over r to first obtain the radius for the QPO frequencies to occur. We defer detailed study on that the QPO phenomenon occurs at r different from ISCO to a future work.
We also consider mock data by assuming the GB metric but not the Kerr ones. For each one b i , we generate the mock data by setting this b i to 5, and fit the data        with the GB metric. we only vary one b i parameter each time. Similar to the Kerr case, only b 2 , b 9 , b 10 , and b 11 can be relatively well constrained, though in most cases the constraints are not as good as those in the Kerr case. The reconstructed parameters can be found in Table 2.
We may compare these results with those using BH shadow and iron line methods as follows. BH shadow [26] shows very weak dependent on the parameters b 2 , b 8 , b 9 and b 10 . Parameters b 4 , b 5 , b 7 and b 11 do not produce any impact on shadow boundary and thus cannot be constrained. According to the iron line studies in [45], all GB parameters except b 11 can be constrained with near fu-ture X-ray mission, though b 5 and b 8 are relatively more difficult to constrain. Different methods may be complementary to each other and put constraints on different deformation GB parameters, therefore, they can be combined together to help break the degeneracies and finally pin down the metric of black holes.

Discussions
In this paper, we simply assume that QPOs are related to the fundamental frequencies at the same radius in the Table 1. Inferred values of the GB parameter from the mock observations generating by assuming the Kerr metric.
We also note that there are various works in literature on constraining alternative metrics for black holes. Below we give a short summary for them and compare some of them with ours presented in this paper.
Reference [38] checked braneworld Kerr BH when bulkspace influence is described by a single, brany tidal charge parameter, b. Similar to Kerr-Newman solution in GR in which the square of the electric charge Q 2 is replaced by a tidal charge b. The behavior of radial and vertical epicyclic frequencies is qualitatively similar to Kerr and brany Kerr BHs. But there are strong differences in the case of naked singularities. The vertical epicyclic frequency could be even lower than radial one. In the structure of radial profile the number of local extrema could be higher in comparison with standard Kerr naked singularity. Also, the radial epicyclic frequency has no zero point for some special family of brany naked singularity. Reference [39] introduced the Jahannsen and Psaltis (JP) metric and studied their QPOs. In their quasi Kerr space time, quadrupole moment is a free parameter in addition to Mass and spin. They show that for moderate spin, the Keplerian frequency is independent of small deviations of the quadrupole moment from Kerr value. They also showed that the epicyclic frequencies shows significant variations. Reference [41] study QPOs in the space time of rotating braneworld gravity. Such a BH carries a tidal charge as imprint of the extra fifth dimension. The metric contains β parameter as tidal charge. They found that large enough value of the positive tidal charge is not supported by observations of high frequency QPOs. But for large enough negative tidal charge, the braneworld BH are similar to the Kerr case. They also showed that over-rotating braneworld and extreme Kerr BH can not be distinguished in the high frequency QPOs observations. Reference [42] considered Bardeen BH metric where parameter g in the metric can be introduced as the magnetic charge of a nonlinear electromagnetic field. He also considered JP metric where metric has an infinite number of deformation parameter k . He fixes the mass and found the spin parameter a/M = 0.279 and g/M < 0.56 at the 68% C.L.. For second metric he found a/M = 0.27 and 3 = 0.5 at the 68% C.L.. In the reference [43], they studied Einstein-Dilaton-Gauss-Bonnet (EDGB) theory with QPOs. α/M 2 characterizes the theory. They considered low spinning BH in this work and showed the LOFT observation can put constrain on the parameter α/M 2 of EDGB. Reference [44] considered Konoplya and Zhikendao metric. The metric is obtained by deforming the Kerr metric by adding a static deformation to mass. They found constraint on spin parameter and deformation parameter δr/r Kerr . If we compare our work with these above mentioned work, we measure QPO frequencies and use MCMC method in addition to directly check data with observations.

Conclusions
In this paper, we investigate whether QPOs can be used to constrain the parameters of the GB metric introduced in Ref. [26] and possible deviations from the Kerr one.
We adopt the general relativistic precession model to relate the QPO frequencies to those fundamental ones, including the periastron and nodal precession frequencies, of test particle at the inner edge of accretion disks around black holes. We find that at least some of the QPO frequencies resulting from the GB metric with b 2 , b 4 , b 9 , b 10 , or b 11 significantly deviating from 1 can be substantially different from those from the Kerr case, though there are no such differences for other b i s.
Considering the case of GRO J1655-40 with three accurately measured QPO frequencies, we find that the GB parameter b 2 is required to be 2.2 +0.395 −0.523 but other GB parameters cannot be well constrained. Although such a constraint on b 2 seems to indicate that a GB metric fits to the QPO data better, but the significance is smaller than 2-σ. By generating mock samples of 10 black hole systems and each black hole with three QPO frequency measurements, we use the nested algorithm to reconstruct the input b parameter(s) by using these mock samples. We find that b 2 , b 9 , b 10 , and b 11 can be well reconstructed, while b 5 , b 7 , and b 8 cannot, as expected. However, b 4 cannot also be well constrained due to parameter degeneracies among b 4 , M , and a. The results presented in this paper suggest that strong constraints on a few GB parameters can be obtained by using the QPO frequency measurements of only about ten or more black hole systems.