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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_{11}$$\end{document}. The parameter b4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_4$$\end{document} 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 observaa e-mail: mghasemin@ipm.ir (corresponding author) tions 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 [30], respectively. The study of quasi-periodic 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 [31][32][33][34].
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 the 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 considering that I have a 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 [30]. 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 Sect. 3. Section 4 shows our χ 2 calculations in this paper. Section 5 is for the results and discussions. Summary and conclusions are in Sect. 6.

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 [35]. The thermal photons from the accretion disk can interact with hot electrons in the corona. Because of inverse Compton scattering the corona becomes an Xray 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 disk-corona model [36,37]. However, a different geometry can be considered [38,39]. Here I consider the simple lamppost geometry. The disk-corona model is described by the parameter 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 [40]. 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 [40]: 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 Sect. 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: 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 . (3) 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. [41] 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 the 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 [33] 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 [33] 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 [33], 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 [42,43]. 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.

Results and discussion
The I study contour levels as discussed in Sect. 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: 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: 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, a reverberation mapping can constrain the b i parameters except for b 11 and also the b 4 is harder to constrain. The power of a reverberation measurement for constraining deviations from GR is clear.
In Fig. 6, I plot the response function for different b i . We see that the responses are not similar to the Kerr case, so we expect that these parameters can be constrained. From the figure, the response shifts to earlier time. We will better see 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 [40].

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  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 a 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  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.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All results can be obtained from the equations presented in this manuscript and reference therein.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .