Iron line reverberation mapping in Ghasemi-Nodehi-Bambi background

The reverberation associated with the iron line is the time lag between direct photons from the corona and the photons reflected from the disk. The resulting line spectrum is called the 2D transfer function. The shape of the 2D transfer function is determined by the geometry of spacetime and the properties of BH. In a paper (Ghasemi-Nodehi and Bambi in EPJC 76:290, 2016), the authors have proposed a parametrization. This parametrization is aimed to test the Kerr nature of astrophysical black hole candidates. In this paper, I provide a reverberation mapping of the Ghasemi-Nodehi-Bambi metric in order to constrain the parameter of spacetime. All parameters can be constrained with the exception of b11. The parameter b4 is harder to constrain too.


Introduction
Einstein proposed the theory of General Relativity (GR) over a century ago [1]. While largely successful in the weak-field tests [2,3,4,5], the validity of GR in the strong field regime is still unexplored. For instance, the motions of pulsars are sensitive to strong-gravitational self-field effects [6] According to the predictions of GR the formation of astrophysical BHs from gravitational collapse is well understood by the Kerr spacetime paradigm. The Kerr spacetime is a stationary, axisymmetric, asymptotically flat solution of vacuum Einstein equations. It is only determined by the mass and spin parameters [7,8]. Because of the presence of a highly ionized environment the charge reaches equilibrium and its value is negligible.
There are different scenarios in which the metric around astrophysical black holes (BHs) may deviate from the Kerr case. The problem is that there is a degeneracy between the parameters of Kerr spacetime and theories deviating from the Kerr case. In order to verify the Kerr solution of GR, it is not only enough to detect the Kerr properties by observations but also any deviations from the Kerr case should be ruled out. Recently, there have been a lot of efforts to study confirmation of GR and constraining the deviation from GR using both electromagnetic radiation and gravitational waves [9,10,11,12,13]. The recent detection of gravitational waves from the coalescence of BH binaries can explore the dynamical strong field regime, but the available data do not constrain significantly the possible deviation from GR [14,15,16]. The constraining power of gravitational waves has been shown in [17].
In Ref. [28] the authors proposed a parametrized metric in which one can recover the Kerr case by setting the deformation parameters equal to one. We already studied BH shadow and X-ray reflection spectroscopy of this metric in [28] and [31], respectively. The study of quasiperiodic oscillation QPO observations of this metric is in progress. In the present paper, we study the reverberation mapping of the metric.
In the framework of the corona-disk model, reverberation is associated with the iron line is the time lag between direct photons from the corona and the photons reflected from the disk. The resulting line spectrum is called the 2D transfer function. The shape of the 2D transfer function is determined by the geometry of spacetime and the properties of a BH. Thus, an accurate measurement of the 2D transfer function can be used as a probe of the spacetime geometry [32,33,34,35].
In this paper, first, I have studied the impact of the deformation parameters b i on the shape of the 2D transfer function. Then, using a minimum χ 2 approach, I study the contour levels to constrain the parameters. Our reference BH is a Kerr BH and total photon count is 10 3 , which is for current observations. Except for b 4 and b 11 , the parameters of the Ghasemi-Nodehi-Bambi (GB) metric can be constrained. b 4 and b 11 introduce a parameter degeneracy. However, higher spin values can remove the degeneracy of parameter b 4 . I also checked the contours for the reference being a non-Kerr GB BH. The results are similar to the cases with the reference Kerr BH. I also checked if it is possible to constrain the height of the corona con-sidering that I have spin parameter from an independent measurement. As is shown, the height of the corona can be constrained except for the case b 11 . I also have studied the BH shadow of the GB metric [28] and X-ray reflection spectroscopy of the metric [31]. In shadow studies, the parameters b 2 , b 8 , b 9 and b 10 leave a small signature in the shadow boundary and b 4 , b 5 , b 7 and b 11 do not produce any specific signature in the shadow shape. In iron line studies, in the presence of the correct astrophysical model, 200 ks observations with future observational facilities such as LAD/eXTP can constrain all the Kerr parameters except for b 11 . The impact of b 11 on the iron line profile is extremely weak. I also try to provide possible constraints on the parameters using QPO, which is still in progress.
The content of the paper as follows. Section 2 is devoted to the iron line reverberation mapping. The geometry of spacetime is discussed in section 3. Section 4 shows our χ 2 calculations in this paper. Section 5 is for the results and discussions. Summary and conclusions are in section 6.
2 Iron line reverberation mapping I use the accretion disk-corona model; the disk is on the equatorial plain orthogonal to the BH spin. The accretion disk emits as blackbody locally and as multicolor blackbody when integrated radially. The so-called corona is a hotter, usually optically thin electron cloud enshrouding the accretion disk. Its exact geometry is unknown, also in this regard some work has been done [36]. The thermal photons from the accretion disk can interact with hot electrons in the corona. Because of inverse Compton scattering the corona becomes an X-ray source with a power-law spectrum. The corona works as a point source located on the axis of the accretion disk just above the BH. This arrangement is known as the lamppost geometry of the diskcorona model [37,38]. However, a different geometry can be considered [39,40]. Here I consider the simple lamppost geometry. The disk-corona model is described by parameter the h as the height of the corona above the disk in addition to the parameters of the BH spacetime. Furthermore, the inner edge is at ISCO. A photon of the corona enters the disk and may produce fluorescence emission line also referred to as the reflection component. The strongest line is the iron Kα line at ∼ 6.4 keV. Here I only consider the iron line. The coronal flux received by the disk obeys a power-law, as r −q , where r is the disk radius; here I consider q = 3, which recovers the Newtonian limit at large distances but it might be different at small radii, r ≈ h.
Here I consider reverberation associated with the iron line. It is the time lag between direct photons from the corona and the photons reflected from the disk. The resulting line spectrum, which is a function of both time and photon energy, is called the 2D transfer function. The shape of the 2D transfer function is different for different geometries and BH spacetimes. In addition to the fundamental properties of BH, it also depends on the height of the corona above the disk, h, and the inclination angle of the disk with respect to the observer's line of sight, i.
The time delay or lag is caused by the difference in light travel time between primary emission and reprocessed emission. I also calculate the frequency dependence and the energy dependence of the lag. First, I calculate the response function for different b i . To plot the response, it is assumed that the rest frame spectrum is simply a δfunction iron line at energy 6.4 keV. Then, in order to calculate the frequency dependence of the lag, I follow the approach considered in [41]. First, I take the Fourier transform of the transfer function. The transfer function in the frequency domain is where ψ(τ ) is the transfer function in the time domain. The phase difference, φ, is as follows [41]: The time lag is φ/2πf . For the energy dependence of the lag I calculate the response of the disk at every energy. To obtain the energy dependence of the lag one should consider a frequency range to plot. The results are discussed in section 5.

Ghasemi-Nodehi-Bambi spacetime
In Ref. [28] one proposed a new parametrization to the Kerr metric. As another parametrization we want to constrain possible deviations from the Kerr solution of GR.
Here we recover the Kerr case when all deformation parameters are equal to 1. Meanwhile, in other metrics the deformation parameters are additive and reduce to the Kerr case for vanishing deformation parameters. We want to see how mass and spin in metric component deform the spacetime. We introduce 11 new parameters in front of any mass and/or spin term. Any deviation from 1 deforms spacetime more or less from that of the prediction of GR. The metric is as follows: ar sin 2 θ r 2 + b 4 a 2 cos 2 θ dtdφ + r 2 + b 5 a 2 cos 2 θ r 2 − 2b 6 M r + b 7 a 2 dr 2 + r 2 + b 8 a 2 cos 2 θ dθ 2 + r 2 + b 9 a 2 + 2b 10 M a 2 r sin 2 θ r 2 + b 11 a 2 cos 2 θ sin 2 θdφ 2 .
The metric reduces to the Kerr metric for the b i = 1 for all i. Here we set b 1 = b 3 = b 6 = 1. b 1 is equal to 1 because it is the coefficient of the mass and b 3 = 1 in the same way; b 3 a is the asymptotic specific angular momentum. b 6 is close to 1 from solar system experiments. We do not consider these three parameters in our reverberation mapping calculations. The primary aim of the parametrization to this metric was to see how each part of the metric provides a signature on observations. Moreover, Ref. [42] provides an example of the usage of this metric. One considered a special class of quintessential Kerr black holes. This class generated a modification to Kerr geometry that apparently is an extension of the modification of the Kerr geometry represented by Ghasemi-Nodehi-Bambi (GB). This modification is because of a special class of quintessential fields. One discusses shadow and spectral line of this modification to the Kerr geometry.

χ 2 calculation for comparison of Kerr and GB background
Here I follow the approach of [34] for comparison of the reverberation transfer functions of the Kerr and non-Kerr backgrounds. I first consider a primary model with spin a * , parameter b i , viewing angle i, emissivity profile q, and height of the corona h. For the 2D transfer function I use the notation [34] for the photon flux number density in the energy bin [E j , E j + ∆E] and in the time bin [t k , t k + ∆t k ] . The secondary model with parameters a * , b i , i , q , h (n jk = n(a * , b i , i , q , h )) would be compared with the primary model by introducing a normalized negative log-likelihood [34], where α = j,k n jk j,k n 2 jk /n jk .
The corresponding χ 2 is N L. N is the number of detected photons. Here I consider N = 10 3 , which is for high quality observation today. Furthermore, the ∆E here is 50 eV and ∆t = M . If I consider M = 10 6 M , ∆t would be about 5s.
I simulate the 2D transfer function with N photons. The simulation is an extension of the code described in [43,44]. I also added Poisson noise to my data. I treat my simulations as mocked data and apply the χ 2 ∼ N L approach to comparing the data. The contours show 1 − σ, 2 − σ and 3 − σ levels. I study contour levels as discussed in section 4. In the first study, the reference model is a Kerr black hole with a * = 0.6, h = 10, q = 3 and i = 55 • . In all cases the total photon count is N = 10 3 . According to contour studies the reverberation mapping can constrain b i except for b 4 and b 11 . The parameters b 4 and b 11 are degenerate. This degeneracy means that Kerr results can be reproduced by the parameters b 4 and b 11 of the GB metric. But in the case of b 4 the degeneracy would be removed by considering a higher spin value, 0.9, but the case for b 11 is still degenerate even for a high spin value and also for a non-Kerr reference model. The values for the b i are as follows:

Results and discussion
Second, I also study the contours with the GB reference model for b i = 1. The reference BH is a GB BH with spin 0.6, the viewing angle 55 • and b i = 5 for each case. Similar to the case with the Kerr BH as a reference, except for b 4 and b 11 , all other parameters can be constrained. The value for b i = 5 is as follows: Next, I consider a spin value of 0.9 from independent observation. The reference height of the corona is set to 10. I try to see if I can constrain the height of the corona. The height of the corona can be constrained in all cases, with the exception of b 11 . The case for b 11 seems more challenging. The value for the height of the corona is as follows: h = 10 +0.001 −0.002 (9) As an example, the plot for the contour level from the analysis of the 2D transfer function with GB BH as the reference BH with b 2 = 5, a * = 0.6 and i = 55 • is drawn in Fig. 4. Also, b 2 vs h contour levels from the analysis of the 2D transfer is shown in Fig. 5. The reference model is a Kerr BH with a * = 0.9, h = 10 and i = 55 • As we see reverberation mapping can constrain the b i parameters except for b 11 and also the b 4 is harder to constrain. The power of reverberation measurement for constraining deviations from GR is clear. In Fig. 6, I plot response function for different b i . We see the responses are not similar to the Kerr case, so we expect these parameters can be constrained. From the figure, the response shifts to earlier time. We will better see this difference in the frequency and energy dependence of the lag.
I have also drawn the lag vs frequency in Fig. 7. As we see the lag starts from a value and then oscillates around zero. There is significant difference from the Kerr case but the cases for b 2 , b 4 , b 5 , b 7 , b 8 , b 11 are very similar. The cases for b 9 and b 10 have less lag at low frequencies than the other parameters.
The energy dependence of the lag is shown in Fig. 8. There is less lag than in the Kerr case for all parameters. The blue part shifts to lower energies if we compare with the Kerr case. The red wings of the parameters b 2 , b 4 , b 5 , b 7 , b 8 , b 11 are longer than the parameters b 9 and b 10 . Also b 9 and b 10 have a smaller lag than the other ones. For the effect of other parameters such as the BH spin, inclination, height of the corona and so on one can refer to [41].

Summary and conclusions
Iron line reverberation mapping concerns the time lag between coronal photons and photons reflected from the disk. This lag depends on the light travel distances. Thus, this provides us with the opportunity to probe spacetime geometry in strong gravity regimes. Here I consider the lamppost coronal geometry; I have an additional height of the corona just above the BH as my parameters. Furthermore, parametrization of the Kerr BH is one way to study deviations from Kerr BH. In this paper I apply reverberation mapping studies to the GB parametrized metric. In GB background I have eight Kerr parameters in addition to the mass and spin parameters of the BH. By the opportunity to study the reverberation associated with the iron line, our GB parameters can be constrained; there is an exception for the parameter b 11 , which introduces degeneracy with the Kerr case. The parameter b 4 is also harder         to constrain. I already studied this metric by the BH shadow, and the time-integrated iron line. Time-integrated iron line studies can constrain our Kerr parameters with the exception of b 11 . Notice that our iron line study deals with future observational facilities, but in a reverberation mapping the photon count 10 3 is for current high quality data. The boundary of the shadow can be slightly altered by the parameters b 2 , b 8 , b 9 and b 10 but there are no signatures on the shadow shape for the parameters b 4 , b 5 , b 7 and b 11 . The response function, frequency dependence and energy dependence of the lag also discussed in the paper. In conclusion, using available current data, the power of the reverberation mapping in the study of strong gravity regimes is clear.