Testing Ghasemi-Nodehi–Bambi spacetime with continuum-fitting method

The continuum-fitting method is the analysis of the thermal spectrum of the geometrically thin and optically thick accretion disk around stellar-mass black holes. A parametrization aiming to test the Kerr nature of astrophysical black holes was proposed in Ghasemi-Nodehi and Bambi in EPJC 76:290, 2016. The metric contains 11 parameters in addition to the mass and spin parameters. One can recover the Kerr case by setting all parameters to one. In this paper, I study the continuum-fitting method in a Ghasemi-Nodehi–Bambi background. I show the impact of each of the parameters on the spectra. I then employ χ2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2$$\end{document} studies and show that using the continuum-fitting method all parameters of Ghasemi-Nodehi–Bambi spacetime are degenerate. However, the parameter b9\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_9$$\end{document} can be constrained in the case of a high spin value and the Ghasemi-Nodehi–Bambi black hole as reference. This degeneracy means that the spectra of the Kerr case cannot be distinguished from spectra produced in Ghasemi-Nodehi–Bambi spacetime. This is a problem as regards measuring the spin of astrophysical black holes and constrain possible deviations from the Kerr case of General Relativity.


Introduction
General Relativity (GR) was born more than a century ago and is still used for describing the gravitational field and geometry of spacetime [1]. The theory has been extensively explored in the weak field regime, such as in solar system experiments and radio pulsars observations [2][3][4][5]. There are many theoretical models that make the same prediction as general relativity in the weak field regime but they deviate in strong field regimes. The strong field regime can be either studied with gravitational waves (GWs) or with electromagnetic radiation. Dynamical systems have a signature on GW but an electromagnetic radiation method can be studa e-mail: mghasemin@ipm.ir (corresponding author) ied when spacetime is static. Effects such as the deviation from geodetical motion can be studied by electromagnetic radiation-related methods.
GR predicts astrophysical black holes (BHs) are Kerr BHs with spin and mass parameters [6,7]. There are different scenarios that deviate from the Kerr case but there is degeneracy in measurement of their parameters and spin parameters of Kerr spacetime. One way to test the Kerr hypothesis is to parametrize the Kerr metric and try to constrain deviations from Kerr spacetime. There are several parametrization in the literature [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22]. In this paper, I consider the Ghasemi-Nodehi-Bambi (GB) metric, proposed in [21]. The Kerr case would be recovered by setting all GB parameters to one. For an electromagnetics-based test of the strong gravity regime, one can use X-ray continuum spectra and iron line emission. One can also study BH shadow observations with an event horizon telescope at mm/sub-mm wavelength. In spite of large systematic errors of quasi-periodic oscillations (QPOs), one can use QPO frequencies to probe strong field regimes as well. Iron line reverberation mapping can be used as well.
Continuum spectra, iron lines and QPOs have already been observed by X-ray missions including RXTE, Chandra, Nustar, ASCA, Suzaku, and XMM-Newton.
The shadow [21], X-ray reflection spectroscopy [23], QPO [33] and iron line reverberation [33] of the GB metric have already been studied. In the current paper, I study the continuum-fitting method to check possible constraints of the GB parametrization.
The continuum-fitting method is used to study the thermal spectrum of geometrically thin and optically thick accretion disks [24]. One can only apply this method to stellar-mass BHs. This is because the spectrum of the supermassive BHs is in the optical/UV band where dust absorption limits the capability of accurate measurement. The temperature of the disk depends on the mass of the compact object. The spectrum of the thin disks around stellar-mass black holes is in the soft X-ray band. Currently, assuming the Kerr BH as astro-physical BH, the spin parameter, a * = J/M 2 , of the BH can be measured by this method. If one considers non-Kerr BHs, this technique can measure deviations from the Kerr solution [25,26].
In this paper, I consider a GB background metric. I simulate thermal spectra of the thin accretion disk and show the impact of the deformation parameters on the spectra. By decreasing (increasing) the deformation parameter b 7 the spectra get softer (harder) and by decreasing (increasing) the parameters b 4 , b 9 , b 10 the spectra get harder (softer). The impact of the parameters b 2 , b 5 , b 8 , b 11 is very close to the Kerr one. The impact of the other parameters was studied before. References [25,27,28] provide the impact of all other parameters. By increasing (decreasing) the mass the spectra get harder (softer). The mass accretion rate changes the temperature of the disk. The distance of the source alters the normalization of the spectrum. For higher (lower) inclination angle, the area of the observed emitting region decreases (increases). If the spin parameter increases the spectra get harder.
Then, in order to compare results with the Kerr case, I use a minimum χ 2 approach. I plot the related contour plots. The results show that using continuum-fitting method the GB parameters are degenerate and cannot be distinguished from the Kerr case. However, the parameter b 9 may be constrained in the case of a high spin value and the GB BH as a reference.
The structure of the paper is as follows. In Sect. 2, I provide details of the continuum-fitting method. Section 3 is devoted to a description of the metric of the spacetime. Results and discussion are presented in Sect. 4. Finally, a summary and conclusions are in Sect. 5.

The continuum-fitting method
The thermal spectrum can be written in terms of the photon flux number density as measured by a distant observer [25,[27][28][29][30]. The photon flux number density is given by where I obs is specific intensity of radiation, E obs is the photon energy, and ν is the photon frequency measured by a distant observer. dΩ obs = dX dY/D 2 , X and Y are coordinate position of the photon on the sky as seen by a distant observer, D is the distance of the source. In order to include all relativistic effects, one needs to compute the photon trajectory from the emission point in the disk to thre image plane of the distant observer (detection point). We have where I e is the local specific intensity of the radiation emitted by disk and w the redshift factor, which are Regarding Eq. (3), since the disk is in thermal equilibrium the emission is blackbody-like. One can define the effective temperature, T eff , F(r ) = σ T 4 eff , where σ is the Stefan-Boltzmann constant and F(r ) is the time averaged energy flux from Novikov-Thorne model [31,32]. The Novikov-Thorne model can be used to describe geometrically thin and optically thick accretion disks around the black hole. The assumption is that the disk is in the equatorial plane. The gas of the disk moves on a nearly geodesic circular orbit. F(r ) can be written as whereṀ is the mass accretion rate, and F(r ) is E, L z , and Ω are, respectively, the conserved specific energy, the conserved z-component of the specific angular momentum, and the angular velocity for equatorial circular geodesics. r in is assumed to be at the innermost stable circular orbit (ISCO). Actually the disk's temperature near the inner edge is high and non-thermal effects are not negligible. For this reason one can introduce a hardening factor or a color factor, f col and the color temperature is T col (r ) = f col T eff . In Eq. (3), ν e is the photon frequency, h is Planck's constant, k B is the Boltzmann constant and Υ is a function of the angle between the wavevector of a photon emitted by the disk and the normal of the disk surface. Here I consider Υ = 1 for isotropic emission.
Regarding the w in Eq. (4), E e = hν e , ν being the photon frequency measured by a distant observer, k α is the 4momentum of the photon, u α obs = (−1, 0, 0, 0) is the 4velocity of the distant observer, and u α e = (u t e , 0, 0, Ωu t e ) is the 4-velocity of the emitter. From Liouville's theorem we have I e (ν e )/ν 3 e = I obs (ν obs )/ν 3 . Finally, the photon flux number density can be written as where A 1 and A 2 are One finds from the normalization condition g μν u μ e u ν e = −1 and thus where λ = k φ /k t is a constant of motion along the photon's path. All relativistic effects are encoded in the redshift factor w.

Ghasemi-Nodehi-Bambi spacetime
Reference [21] proposed a new parametrization to the Kerr metric. The purpose of this parametrization is to constrain possible deviations from a Kerr solution of GR. The Kerr case would be recovered when all deformation parameters are equal to one. While in other metrics the deformation parameters are additive and reduce to the Kerr case for vanishing deformation parameters. In this parametrization, 11 new parameters are introduced in front of any mass and/or spin term. Any deviation from one deforms the spacetime more/less than that of prediction of GR. The metric is as follows: Mar sin 2 θ r 2 + b 4 a 2 cos 2 θ dtdφ + r 2 + b 5 a 2 cos 2 θ r 2 − 2b 6 Mr + b 7 a 2 dr 2 + r 2 + b 8 a 2 cos 2 θ dθ 2 + r 2 + b 9 a 2 + 2b 10 Ma 2 r sin 2 θ r 2 + b 11 a 2 cos 2 θ sin 2 θ dφ 2 . (11) Here we set b 1 = b 3 = b 6 = 1. b 1 is equal to one because it is the coefficient of mass and b 3 = 1 in the same way; b 3 a is the asymptotic specific angular momentum. b 6 is close to one from a solar system experiment. I do not consider these three parameters in my continuum-fitting method calculations. The impact of all other parameters is studied in [25,27,28]. Now, in order to compare my results, I use the χ 2 approach. I compare the spectra produced in Kerr spacetime with one expected from the GB background. The χ 2 is defined as follows: where summation is performed over n sampling energies The brown, orange and lighter orange are for 1σ, 2σ and 3σ (68%, 95% and 99.7%) confidence level, respectively. If the contour plots would be closed, this means we can constrain the parameters and there is no degeneracy, but, here,  Figs. 2 and 3. So we cannot constrain the parameters. In other words, there is degeneracy between measurements of the parameters b i and measurement of the spin parameter, so that the parameters b i cannot be constrained. Degeneracy means one cannot distinguish the properties of these spectra from spectra of a disk around a Kerr BH.
This degeneracy is a problem as regards measuring the spin of astrophysical BHs. So far, spin measurement is based on the Kerr BH hypothesis but as we see from the result of this study with the continuum-fitting method, the same electromagnetic observation (thermal spectra here) can be reproduced by non-Kerr spacetime as an example of GB spacetime with different spin and deformation parameters. I also consider the case for higher values of the spin, such as 0.9, as reference model but still the results are degenerate. Then we consider the reference BH not to be a Kerr BH. I consider a GB BH as reference BH in a contour study with spin 0.9 and b i = 18. In this case, only the parameter b 9 may be constrained because the contour looks closed, as we see in Fig. 4. The parameter b 10 is harder to constrain is this case as well.
As the result of other tests with electromagnetic radiation, considering triplet frequency observations of GRO J1655-40, we find that the GB parameter b 2 is 2.2 +0.395 −0.523 but the other parameters cannot be constrained [33]. If we assume we have 10 observations in the future, we can constrain more parameters, we find that b 2 , b 9 , b 10 and b 11 can be well constrained, while b 4 , b 5 , b 7 and b 8 cannot. With X-ray reflection spectroscopy, known as the iron line method, 200 ks observations with future observational facilities such a LAD-eXTP  11 . An iron line reverberation mapping study of GB spacetime can constrain all parameters except b 11 . b 4 also is harder to constrain [33]. As we see, the BH shadow can constrain some parameters but can not constrain others. Iron line studies either with time-integrated analysis or with reverberation can constrain all with exception of b 11 . But a mock observation of 10 BHs can constrain b 11 in addition to the other parameters. Unfortunately, due to degeneracy of the parameters in the continuum-fitting method, we cannot further restrict GB parameters using this method. This might be a general problem as regards measuring the spin of the BH.
A combination of different observational methods may help to constrain the deformation parameters and break the degeneracy of the parameters.

Summary and conclusions
The continuum-fitting method is a method to study thermal spectra. The disk should be geometrically thin and optically thick. Also this study is valid for stellar-mass black holes.
Furthermore, parametrization of the Kerr BH hypothesis is one way to study deviations from the predictions of GR and provide constraints on the deviation from GR.
In this work I use the continuum-fitting method to check whether it can constrain the GB metric parameters introduced in [21]. The GB metric contains 11 deformation parameters, which the Kerr case recovers when all deformation parameters are set to one. I first compute the thermal spectra for different b i parameters. My results show that by increasing parameter b 7 the spectra get harder. Also by increasing the parameters b 4 , b 9 , b 10 the spectra get softer and impact of the parameters b 2 , b 5 , b 8 , and b 11 is very close to the Kerr one.
I then consider my simulation as observational data and try to study possible constraints of the deformation parameters of GB spacetime using the χ 2 approach. My contour studies show that all deformation parameters are degenerate. However, for the case of a high spin value with GB BH as reference, the parameter b 9 may be constrained. This degeneracy means the observation with specific spin and the deformation parameter can reproduce the spectra of the Kerr case. This is a problem as regards measuring the spin of astrophysical black holes and constrain possible deviations from GR. Also as a result of this degeneracy, in order to verify the Kerr hypothesis, any deviation from the Kerr case should vanish.
Different observational tests such as the BH shadow method, X-ray reflection spectroscopy, iron line reverberation mapping and QPO introduce different impacts of the parameters. Combination of different methods may help to break the degeneracies.