Qualifying ringdown and shadow of black holes under general parametrized metrics with photon orbits

The motion of photons around black holes determines the shape of shadow and match the ringdown properties of a perturbed black hole. Observations of shadows and ringdown waveforms will reveal the nature of black holes. In this paper, we study the motion of photons in a general parametrized metric beyond the Kerr hypothesis. We investigated the radius and frequency of the photon circular orbits on the equatorial plane and obtained fitted formula with varied parameters. The Lyapunov exponent which connects to the decay rate of the ringdown amplitude is also calculated. We also analyzed the shape of shadow with full parameters of the generally axisymmetric metric. Our results imply the potential constraint on black hole parameters by combining the Event Horizon Telescope and gravitational wave observations in the future.


Introduction
One of the main problems of gravity theories is to test the theory in strong field regime with high accuracy. However, the LIGO-Virgo experiments on gravitational wave observation [1][2][3] and observation of black hole shadow image at the center of M87 elliptic galaxy [4] give us the opportunity to develop new tests of general relativity and modified/alternative theories of gravity in the strong field regime. Further improvements of the experiments and observations will give more precise results and opportunity to obtain constraints on different theories of gravity.
Up to now the general relativity proposed by Albert Einstein in 1916 is considered as the main theory of gravity which is justified by the different tests in the weak (e.g. Solar system tests [5,6]) and strong field (e.g. gravitational a lisong20@mails.ucas.ac.cn b ahmadjon@astrin.uz c Corresponding author: wbhan@shao.ac.cn waves [7,8], shadow of M87 [9-12, 114, 116] regimes with high accuracy. However, general relativity meets number of problems connected with the spacetime singularity at the origin of the solutions of the field equations, renormalization, compatibility with quantum field theory and etc. Despite the attempts to resolve the singularity problem introducing conformal transformation [13][14][15][16][17], coupling with nonlinear electrodynamics [18][19][20], higher dimensional corrections [21] and etc, no single unique theory has been found to resolve all fundamental problems of general theory of relativity. Thus, we need to deal with the big number of the modified/alternative theories of gravity to resolve the fundamental problems of general relativity and describe the current observational and experimental data. There is a belief that astrophysical black holes are described by the Kerr spacetime (at least by Schwarzschild spacetime when the effect of rotation is negligible). However, the so-called Kerr hypothesis has not been yet well confirmed by current experimental and observational data. The new parameters of solutions within modified or alternative theories of gravity representing the deflections from Kerr spacetime may mimic the effects of the spin parameter of the Kerr black holes [22][23][24]. Due to this fact one needs to develop more independent tests of the gravity theory [25,26] together with further improvements of the astrophysical instrumentation [27,28].
On the other hand big number of alternative and modified theories of gravity and corresponding solutions of field equation describing the compact gravitating object creates the difficulties associated with resolving the cross mimicking of the parameters of theories. One of the attempt to resolve this issue is to use the parametrization of the spacetime metric. The parametrization of the spacetime around rotating black hole may help to cover big number of the solutions of different gravity models. Several ways of parametrization of the spacetime describing the rotating black hole have been proposed by different authors [29][30][31]. So-called KRZ metric was propsed by Konoplya, Rezzolla, and Zhidenko [30] where they have suggested parametrization of the spacetime of rotating black hole. They have proposed a parametric frame to describe the spacetime of the axisymmetric black holes which have the Killing vector η = (0, 0, 0, 1). The most interesting point of the KRZ metric is that it contains many significant parameters. These parameters associated with physical properties of the black hole and we briefly discuss them in the text of the paper. Particularly, when all the parameters of KRZ parametrization are equal to zero, the specetime metric reduces to Kerr one.
The most distinguishable feature of the metric theories of gravity is the light deflection and other effects connected with photon motion in the curved spacetime. Some authors studied the photon surfaces [32] through the photon motion. The gravitational lensing was the one of the first consequences of the general relativity and has been discovered by Einstein. The study of the light deflection can be used to study either gravitational object or distant source. The review of the gravitational lensing effect can be found in Refs [33][34][35][36]. The effects of the plasma on gravitational lensing in different spacetimes have been studied in [37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54][55][56]. The most basic gravitational lensing is the Schwarzschild lensing and it have been studied in [57,58], some papers also studied the cosmic censorship hypothesis (CCH) with the gravitational lensing [59,60].
Additionally, the properties of circular orbits of photons around black holes reflect on the ringdown signal from the merger of binary black holes [61,62]. The latter are quasinormal modes (QNMs) which describe the end state of a black hole-black hole merger. Therefore the gravitationalwave emission at late times can be well described by the properties of null geodesics on unstable circular orbits at the black hole's light ring [63,64]. The direct calculation of QNMs of the KRZ metric is difficult, however, one may use the alternative method to calculate the frequency and Lyapunov exponent of unstable circular orbits leading to the features of ringdown of waveforms [61][62][63][64]. Some authors work on the spherical and static parametrized RZ metric and proposed the higher order WKB method to reduce the difficult of computing the QNMs [65]. On the other hand the standard metric perturbations of the Schwarzschild black hole has been studied in the pioneer work of Regge and Wheeler [66] and Zerilli [67]. Later the QNM of the black holes have been calculated using the perturbative method in different works (see, e.g. [68][69][70][71][72][73] and reference therein).
The recent [74][75][76][77] and future observation of the black hole shadow by Event Horizon Telescope (EHT) using very long baseline interferometry (VLBI) technique can be used to explore the the gravity in the strong field regime around supermassive black hole (SMBH). At the same time one may test the gravity theories using the observational data from the black hole shadow. The shadow of SMBH has been theoretically studied in Refs. [43,52, within the different gravity models. Here we plan to study shadow of the black holes described by the parametric spacetime metric proposed in [30].
The paper is organized as follows: Sect. 2 is devoted to briefly review of the motion of massive and massless particles in the KRZ space-time and construction of the ray tracing algorithm necessary to investigate the shadow. We also study the frequency of photon orbits (Sect. 3) and the Lyapunov Exponent (Sect. 4) of light ring. Sect. 5 describes the ray-tracing code used to construct the shadow of the KRZ metric. In Sect. 6, we consider the shadow cast by the KRZ space-time for observer at infinity. Finally, in Sect. 7 we summarize the obtained results. Throughout the paper we use a space-like signature (−, +, +, +), a system of units in which G = c = 1. Greek indices run from 0 to 3, Latin indices from 1 to 3.

Photon motion
In this section we explore the parametrized KRZ metric proposed in [30] and investigate the photon motion around compact object described by KRZ metric. The lowest-order metric expression of the KRZ parametrization has the following form: wherer = r/M,ã = a/M and the other metric functions are defined as [105]: In this paper we use the following parameters defined as [105]: where r 0 is the radius of the event horizon in the equatorial plane and δ i (i = 1, 2, 3, 4, 5, 6, 7, 8) is the dimensionless parameter describing the corresponding deformation of the parameter in the metric (1). Particularly, δ 1 corresponds to the deformation of g tt , δ 2 and δ 3 correspond to the deformations of spin, δ 4 and δ 5 correspond to the deformations of g rr , δ 6 corresponds to the deformation of the event horizon.
In the case when δ i = 0 the KRZ one (1) reduces to Kerr metric andã = 0 reduces the Kerr metric to Schwarzschild one. The stationary and axisymmetric KRZ metric is independent of t and φ coordinates which leads to existence of timelike and spacelike Killing vectors. Consequently, these two Killing vectors correspond to two conserved quantities: the energy E and the z-component of the angular momentum L z of test particle. The conserved energy and angular momentum of the test particle can be expressed as: where the overhead dot represents the derivative with respect to the affine parameter (proper time for a massive particle). One can thus express the equation of motion of test particles with these two conserved quantities. Substituting Eqs. (17)- (18) into the normalization condition of the four-velocity u α u α = −1 for a massive particle, where u α = (ṫ,ṙ,θ ,φ ) is the 4-velocity, one may obtain the following equation for the motion in the equatorial plane (θ = π/2): Similarly one can consider the orbits of photons around black hole. For the photon orbits the normalization condition of the four-velocity take the form u α u α = 0. Considering the orbit in the equatorial one may get the following expression: 1. We expand the normalization equation u α u α = −1 or u α u α = 0. 2. We substitute the equations for the conserved quantities E and Φ (17)- (18) into the normalization equation. 3. We rewrite the normalization equation by the two conserved quantities in a form similar to the equation in Newtonian mechanics.
The equations of radial motion containing the effective potential for particle and photon have the following form: where V eff and W eff are the effective potential for particle and photon orbits, respectively. The expressions for the V eff and W eff depend on parameters of the chosen spacetime metric. However, application of these steps to the KRZ metric becomes problematic due to its complicate form. Thus here we use a method described in Ref. [106]. For the light ring (LR) we have the following condition where the V eff has the following form One may now easily introduce new potential functions rewriting the Eq. (24)in the following form where the σ = 1/b. b = |Φ/E| is the inverse impact parameter. Newly introduced effective potential H ± of photon orbits on the orthogonal 2-space has the following form where D ≡ g 2 tφ − g tt g φ φ . The circular orbits of photon obey the condition ∂ r H ± =0, where the ± sign is associated to the two direction of the rotation. In this paper for our analysis we consider the case corresponding to the sign "+". Now we can calculate the circular orbits of photon in the equatorial plane i.e. θ = π/2 through the method described in [106]. Since we consider the motion in the equatorial plane and use the condition ∂ r H ± =0, one can find out that the circular orbits of photon depend on δ 1 , δ 2 andã only. We cannot get the exact analytical relationship of these three variables δ 1 , δ 2 andã, because the relationship of the three variables δ 1 , δ 2 andã in the equation ∂ r H ± =0 are complicated. So we try to perform a numerical fit to get the fitting equation F r (δ 1 , δ 2 ,ã). After a series of fittings, we get the equation F r (δ 1 , δ 2 ,ã) (see Appendix A). Fig. 1 shows the numerical results and fitting results with different variables. From Fig. 1 we can see that the fitting results fits well with the numerical results. With the increase of spin a and δ 2 , the radius of photon orbits decrease.
After obtaining the radius of the photon circular orbit, one can explore the photon motion. The geodesic equations for null geodesics has the following form: where λ is the affine parameter, Γ µ ντ are Christoffel Symbols defined as From Eqs. (27)- (28), one can easily get a differential equations including the terms Using the relation for radius of photon orbits and the Eqs. (17), (18) and (20), one can determine the whole set of initial conditions.

Frequency of the photon orbits
As we have mentioned before, the QNM of a perturbative black hole can be determined by the unstable circular orbits of photons around the black hole. The frequency of QNM is related the orbital frequency of the light ring [61]. Following the procedure described in Ref. [64], we calculate the frequency of light ring in order to quantify the frequency of ringdown waveforms from KRZ black holes.
In this subsection, we will explore the frequency of the photon orbits in the equatorial plane around black hole described the KRZ spacetime metric. One can calculate the frequency of photon orbits around Kerr black hole in the following way. Using the normalization condition of the fourvelocity u α u α = 0 one may easily write the radial equation of motion in the form: where σ m ≡ sign(Φ) and photon effective potential W eff (r, b, σ m ) has the following form: Solving the equation ∂W eff / ∂ r| r cir = 0 and using the value for r cir one can easily get the the angular velocity Ω = dφ /dt and consequently get information about the frequency. We may apply the same method to calculate the frequency for the KRZ spacetime.
From the Eqs. (18) and (20) one may get the following equation Comparing the Eq. (29) with Eq. (31), one can see that the Eq. (31) does not contain the term 1/b 2 . This is due to fact that the parameter b is independent of the radial coordinate r. Consequently, one can include the parameter b into the expression for W eff in order to calculate the frequency. Now we solve the equation using the value r cir obtained using the fitting and get the frequency of the circular photon orbits in KRZ spacetime in the equatorial plane. Since we consider the equatorial plane the frequency of photon orbits depend only on δ 1 , δ 2 andã under the condition δ 4 = 0. It is easy to find the frequency of the photon orbits for the different fixed values of δ 1 , δ 2 and a. However, the analytical expression describing the relation of radius with δ 1 , δ 2 andã cannot be found explicitly due to complicated view of the metric functions. Here we tried to perform a numerical fit to obtain the equation F ω (δ 1 , δ 2 ,ã) for the frequency of photon orbits. After a number of fittings, we have obtained the expressions F ω (δ 1 , δ 2 ,ã) (see Appendix B). Fig. 2 shows the numerical and fitting results for the frequency of photon orbits. Although the data shows some differences on the graphs, the relative error between the data is found to be less than one percent after calculation. From Fig. 2 we can find that when δ 1,2 is as large as 0.5, the frequency changes less than 10%. Therefore, only when these non-Kerr parameters are large enough, we can read the derivation from the frequency of ringdown waveform.

The Lyapunov Exponent of Light Ring
The Lyapunov exponent (LE) is the key indicator of chaos in dynamical systems. Interestingly, the LE of the unstable circular orbit of photon corresponds to the decay rate of the ringdown amplitude from a perturbed black hole [62]. The Lyapunov coefficient characterizes the rate of divergence of nearby null geodesics. Based on the orbital frequency and LE of light-ring, analytical black hole binary merger waveforms are constructed in [64], and the waveform amplitude decays as where A p is the amplitude when the congruence converges. In this section we compute the LE of the the photon circular motion around a generally axisymmetric black hole which is described by the KRZ metric. Authors of Ref. [107] have considered two particle's orbital motions with small differences in their initial conditions, then using the metric functions they have calculated LE. However, in this paper we use another approach proposed by McWilliams [64]. We consider the perturbation in the radial direction and use the following expression for LR radius [64]: where t p is the time when the congruence converges, ε is a small dimensionless order-counting parameter and r lr is the light ring radius. ϒ in Eq. (32) is a function defined as: where γ is the Lyapunov exponent. Since we consider the motion in the equatorial plane we consider the effect of the paramaters δ 1 and δ 2 in the KRZ metric Eq. (1). Using the Eqs. (32)-(33) one may explore the dependence of LE from the metric parameters. For the appropriate time we plan to plot three-dimensional graphics and explore the influence of the parameters δ 1 and δ 2 of the metric Eq. (1) on LE.
We have chosen the different orbits having separation in theṙ of the order 10 −5 with the other initial conditions to be the same. In Fig. 3 we present the 3-D dependence of LE from the parameters δ 1 and δ 2 . From the Fig. 3 one can easily see that when δ 1 increases from 0 to 1, LE will have two maximum values and there is no obvious downward trend in the overall dependence. When δ 2 increases from 0 to 1 LE have a single maximum value and the overall trend of the graph is down. From the results of Fig. 3 one may speculate that this is related to the initial conditions: theṫ 0 anḋ φ 0 decrease with increasing δ 1 , theṫ 0 andφ 0 increase with increasing δ 2 . Moreover, the angular momentum l increases with the decrease of δ 1 and δ 2 . However, the angular momentum l decreases faster with the increase of δ 2 . Finally, we conclude that the results show in Fig. 3 is related to the interaction of these initial conditions affected by the parameters of KRZ metric. From Fig. 3, we can find that the magnitude of LE is very sensitive to the KRZ parameters δ 1 and δ 2 . Comparing with the frequency changes due to non-Kerr parameters, the decay rate of ringdown can constrain δ 1 and δ 2 better. However, from the figure, the variation of LE is very complicated and hard to be fitted. This creates a problem to construct a parameterized ringdown waveform.

Ray-tracing code for photons
The information from the the distant source in the KRZ spacetime comes through the study and analyze of the light ray from them to the observer. Through the light ray one can get an image of the source and consequently get some information. In this work we study the the trajectories of photons in the KRZ spacetime using the ray-tracing code described in the Ref. [108]. The code describes the trajectories of photons near the black hole.
The evolution of the photon's position with different components: the tand φ -components can be obtained by the first-order differential equations (17)- (18). Then one can rewrite p t and p φ in terms of two parameters: the normalized affine parameter λ = E/λ For the remaining rand φ -components of the photon's position in the KRZ space-time, we can use the second-order geodesic equations with the normalized affine parameter and the Christoffel symbols Γ σ µν as: In this way we can get the system of equations that the ray-tracing code can be used for KRZ spacetime.
We suppose that the massive source described by the KRZ spacetime is located at the origin of the reference frame and coordinate system. We choose the mass of the object M = 1 since it does not affect the shape of the shadow. We assume that the observer's screen is located at a distance away the source of d = 1000, the azimuthal and polar angles are γ rt and 0, respectively. The celestial coordinates (α, β ) on the observer's sky are related to polar coordinates r scr and φ scr on the screen by α = r scr cos(φ scr ) and β = r scr sin(φ scr ). Since we only know the positions and momenta of the photon in the screen, we should solve the geodesic equations from the screen to the source. The photons depart from the screen with a four-momentum perpendicular to the screen and other initial conditions. The method assumes that the screen at spatial infinity, only the photons which moving perpendicular to the screen at a distance d could influence the infinite screen.
The initial position and four-momentum of each photon in the KRZ spacetime are given as [109] and Using of the equation (20) of the photon four-momentum to be zero one can find the component (dt/dλ ) i . The conserved quantity b, which is involved in Eqs. (34) and (35), is calculated from the initial conditions of E and L z .
The initial conditions of the code on the screen is defined in the following way. The confines of the location of the compact source is found inside 0 ≤ r scr ≤ 20, and the value of the φ scr in the range 0 ≤ φ scr ≤ 2π with step of π/180. The confines is the border between the photons that are captured by the compact source and the photons that are able to escape to infinity. The photons are considered as captured by the compact source if they cross the surface r = r surf + δ r with δ r = 10 −3 , where r surf is the radius of the horizon. Then the confines is amplified in to an accuracy of δ r scr = 10 −3 to accurately determine the shadow boundary with the value of r scr for the corresponding value of φ scr . This method allows one to accurately calculate the shadow produced by light ray in the KRZ parametrized metric with high accuracy with respect to sampling the entire screen.

The shadow of the KRZ metric
In this section, we plan to study the apparent shape of the compact object shadow under the KRZ spacetime. We can use the celestial coordinates α and β [17] to describe the shadow of the compact object described by the KRZ spacetime where r 0 is the distance between the massive source and observer and θ 0 is the inclination angle between the observer lens axis and the normal of observer's sky plane (see Fig. 4). In order to describe the dependence of the shadow's shape with different deformation parameters, we will use the coor-   dinate independent formalism proposed in [110]. The shape of the shadow is described by the horizontal displacement from the center of the image D, the average radius of the sphere R , and the asymmetry parameter A. Since the KRZ spacetime is axially symmetric, the parameter D is always identically equal to zero. In the papers Ref. [111,112] authors proposed different ways to describe the shape of the shadow. However, the results of the different approaches are similar to each other. The average radius R is the average distance of the boundary of the shadow from its center, which is defined by where R(ϑ ) ≡ (α − D) 2 + β (α) 2 1/2 , D = 0 and ϑ ≡ tan −1 [β (α)/α)].
The asymmetry parameter A is the distortion of the shadow from a circle and defined as The shadow of the compact object in KRZ spacetime for the different values of metric parameters is shown in Figs. 5-8. From the Figs. 5-8 one may come to the following conclusions: -From Fig. 5 (a) one can see that with the increase of rotation parameter a one side of the shadow goes away from the center while other one comes closer. -From Fig. 5 (b) one can observe that the presence of the parameters δ 1 and δ 2 force the shape of shadow to be more flatter. -The Fig. 5 (c) shows the same effect similar to one caused by the rotation parameter a: the result is obvious because  the effect of inclination angle between the observer lens axis and the normal of observer's sky can be converted by projection of the spin a.
-From Fig. 5 (d), one can see that the parameters δ 1 and δ 2 don't affect inclination as much as they spin parameter does. -The Fig. 6 (a) shows the shadow is more small when δ 1 increases. Since δ 1 related to g tt one may conjecture as g tt increases the size of shadow decreases just like as the particle or photon moves around the black hole with different energy then captured by the black hole . -The Fig. 6 (b) shows that with the increase of δ 2 the shadow becomes more asymmetric and the center of the image will be shifted. -The Fig. 6 (c) shows the similar effect as one due to parameter δ 2 . Here the interesting point is that with increase of the value of parameter δ 3 the additional protrusions will appear in the image of the shadow. Since the effect of the δ 2 and δ 3 is similar to the one due to spin parameter a, they may be related to rotational deformations of the spacetime metric. -The Fig. 6 (d)-(e) show that the size of shadow becomes bigger with the increase either δ 4 or δ 5 . Due to this reason δ 4 and δ 5 may be related to deformations of g rr . -The Fig. 6 (f) shows that as the value of δ 6 moves away from 0 the shadow spreading inward/outward is more obvious, and the graphs become similar to the egrosphere's outer boundary. This may be due to fact that the δ 6 is related to deformations of the event horizon.
Summarizing we may conclude that the parameters δ 1 and δ 4 can change the size of the shadow but with opposite effect, the parameters δ 2 makes the shadow deviate from the center point, the parameters δ 5 also can change the size of the shadow but not obvious and a little stretched, the parameters δ 3 and δ 6 can change the contour shape.
From Fig. 7 we can see when δ 7 increases, the change of the shadow is small, the small picture in the upper right corner of the picture is the enlarged view of the graph. It shows the minor difference corresponding to the different values of parameter δ 7 . From Fig. 8 one can observe similar results corresponding to the parameter δ 8 . Since the effect of the δ 7 and δ 8 is very small, in some places (see, e.g. [31]) authors neglect the parameters δ 7 and δ 8 and put the Eq. (7) to be 1. Although the simplification has very weak impact on the shadow, it is mathematically misleading, because the metric used in the Ref. [31] does not reduce to the Kerr metric when the δ i = 0.
We also want to know how the KRZ metric can reproduce the exact metric, so we use the Kerr-Sen(KS) metric and Einstein Dilaton Gauss-Bonnet(EDGB) metric [115]. Fig. 9 shows the difference between the KS/EDGB metric and KRZ metric. Through this picture we can know KRZ metric has a good approximation to KS and EDGB metric. Fig. 10 shows the dependence of R and A as a function of the spin parameter (Kerr metric) and δ i (KRZ metric) for the fixed values of inclination angle θ 0 = π/4 and the spin parameter a = 0.5. From Fig. 10 (b) one can easily see that the δ 1 and δ 4 have a greater impact on R , but with different trend: when δ 1 increases the value of R decreases and when δ 4 increases the value of R increases. This is due to the fact that δ 1 corresponds the deformation of g tt , while δ 4 corresponds the deformation of g rr . On the other hand δ 2 and δ 3 have the same trend: when their values increase the R decreases. From Fig. 10 (d), we may also see that the δ 2 have a greater impact on A. Other δ i 's have small influence on A. From Fig. 10 (b)-(d), one may see that the trend of R and A with δ 4 , δ 5 always opposite to effect of δ 2 , δ 3 and δ 6 . Finally, with the increase of δ 7 and δ 8 , the value of R and A change very slowly, which also shows that δ 7 and δ 8 are adjustment parameters and have weak effect on the KRZ metric.

Conclusion
In the present work we have studied the photon motion in the equatorial plane around the generally axisymmetric black hole which is described by the parametrized KRZ metric. From the properties of the circular orbits of photon, we can quantify the frequency and decay rate of ringdown GW signals by the orbital frequency and Lyapunov exponent of the light-ring. At the same time, the shape of the shadow which be measured by distant observers is also gotten from the photons motion around the parametrized black hole. We have calculated the frequency and LE of the unstable circular orbits of photon in order to get the information on QNMs. We have also obtained the frequency and decay rate of ringdown which can be used to construct a waveform model for the KRZ black hole merger.
In the special case when the photon orbits in the equatorial plane, we have found that only two primary parameters δ 1 and δ 2 of the KRZ metric affect on photon trajectory. It has been shown that with the increase of the spin parameter a, δ 1 and δ 2 the radius of photon circular orbits decreases. However, the change of the radius of photon circular orbits are more sensitive to the change of spin parameter a, especially when the spin parameter exceeds the value 0.4. It has been found out that effects of spin parameter a and δ 2 on the frequency of photon orbits is the same order while the effect of δ 2 is weaker. With these fitted formula, one can in principle get the frequency of ringdown waveforms from a perturbed black hole described by the KRZ metric.
The decay rate of QNM amplitude can be obtained from the LE characterizing the rate of divergence of nearby null geodesics. In this work we hve shown that when δ 1 increases from 0 to 1, LE will have two maximum values and there is no obvious downward trend in the overall trend of the dependence. When δ 2 increases from 0 to 1 LE will have one maximum value and the overall trend of the graph is down. This is related to the initial conditions of the equation of motion. Though the angular momentum l increases with the decrease of δ 1 and δ 2 . We have shown that the decay of ringdown is very sensitive to the KRZ parameters.
We have also studied the shadow of the black hole described by the parametrized KRZ spacetime. It has been shown that the parameters δ 1 and δ 4 can change the size of the shadow but with opposite effect. The parameter δ 2 makes the shadow deviate from the center, the parameter δ 5 also can change the size of the shadow but effect is relatively weak. The parameters δ 3 and δ 6 can change the contour shape. The shape of the shadow has different sensitivity to different parameters: the shadow is more sensitive to parameters δ 1 , δ 2 , δ 4 and less sensitive to parameters δ 3 , δ 5 , δ 6 . The effects of the parameters δ 7 and δ 8 to the shape of shadow is very weak and almost negligible (compare with [31] where authors have neglected the effects of δ 7 and δ 8 on the iron line in the X-ray spectrum of black holes).
One of the main results of this paper is the analysis of the dependence of the average radius of the shadow and asymmetry (distortion) parameter from the spin parameter and KRZ parameters. Among the effects of other parameters the effect of parameter δ 1 is dominant in changing the average radius of the shadow. On the other hand the main contribution to the change of the asymmetry parameters comes due to the presence of the parameter δ 2 . In principle, one may see that the two observable quantities (radius of shadow and asymmetry parameter) could provide the rough estimation of the parameters δ 1 and δ 2 . Further analysis of the KRZ prameters and comparison with particular black hole solutions may provide a useful tool to probe the gravity models.
The ringdown and shadow reflect the dynamical process and geometric properties in the strong field of the black hole, respectively. The former can be observed by the ground and space-borne GW detectors and the later can be observed by the EHT. From the GWs and image of EHT, one can reveal the nature of the black holes and test if they are described by the Kerr spacetime which is an exact solution of Einstein field equation in general relativity and assumed to describe the astrophysical black holes. Fortunately, both of these two phenomenon are related to the photon orbits around the black hole at the light-ring. In the present work, by calculating the photon's motion at the light-ring, we qualify the QNMs and shadows of generally axisymmetric black holes under general parametrized metrics. Perturbing the supermassive black hole and radiating the ringdown signals can be expected in our Galaxy [113], and the imaging of this nearest supermassive black hole is a target of the EHT project. Our results may play a role to construct a joint constraint of dynamical and stationary spacetime of black hole with both LISA and EHT observations.