Observational redshift from general spherically symmetric black holes

,


I. INTRODUCTION
Black holes are one of the most interesting compact objects obtained from Einstein's general relativity theory and they are simple objects, completely described by only a handful of parameters, namely the black hole mass M and angular momentum J.The existence of black holes has been identified to date through the pieces of evidence provided by observations of the orbital motion of stars at the center of our galaxy [1][2][3][4], captured shadow images of supermassive black holes hosted at the core of M87 and Milky Way galaxies [5,6], and gravitational wave detections by LIGO-Virgo collaborations [7,8].
Recent discoveries in identifying black holes in nature have motivated astrophysicists to invent and develop new methods in order to obtain black hole parameters from different types of observations and gain a better understanding of the underlying physics.In this direction, a robust general relativistic method has been provided in [9][10][11] which consists of obtaining the black hole parameters, such as mass and angular momentum, in terms of the observational redshift/blueshift of photons emitted by geodesic particles circularly orbiting in the vicinity of the Kerr black hole.Such a procedure allows one to obtain analytical expressions for the black hole's mass and spin in terms of directly measurable quantities [10].
The primary work [9] has been introduced based on the kinematic frequency shift which is not a directly oba Electronic address: diego.martinezval@alumno.buap.mxb Electronic address: mmomennia@ifuap.buap.mx,momennia1988@gmail.comc Electronic address: aherrera@ifuap.buap.mxservable quantity and has been applied to a wide variety of black hole spacetimes, such as regular black holes in a spherically symmetric background [12], Kerr-Newmann black holes in de Sitter geometry [13], and higher dimensional black holes [14].Then, this general relativistic method has been extended to incorporate observational total frequency shift in the mathematical modeling of black hole rotation curves [10].The new extended method has been employed to express the parameters of static polymerized black holes in terms of the total redshift [15].More recently, it was shown that by employing a similar procedure, it is possible to extract the Hubble law from the Kerr black hole in asymptotically de Sitter spacetime and add a new independent general relativistic approach to measure the late-time Hubble constant [11].
In addition, one can apply this general relativistic formalism to real astrophysical systems with available data on the frequency shift and positions of orbiting particles.This method has been used so far to estimate the mass-to-distance ratio of 17 supermassive black holes in the center of active galactic nuclei (AGNs) which possess an accretion disk of water vapor clouds [16][17][18][19].It was demonstrated that this relativistic approach allows us to quantify the gravitational redshift, hence identifying the relativistic effects in such astrophysical systems.Recently, by introducing the "redshift rapidity" which is an observable element, we disentangled M and D in the Schwarzschild black hole spacetime and expressed mass and distance to the black hole just in terms of observational frequency shifts [20].The redshift rapidity has also been employed to express the mass, charge, and distance to the Reissner-Nordström black hole in terms of a few directly observable quantities, such as the total frequency shift and aperture angle of the telescope [21].
On the other hand, by considering recent electromagnetic and gravitational wave detections from black holes, we find that there are still large uncertainties in estimating the background test fields and the black hole parameters, such as electric charge and angular momentum [1][2][3][4][5][6][7][8].This uncertainty in obtaining parameters leads to the opportunity to explore their strong field regime and allows for the existence of black hole solutions beyond Einstein theory of relativity.The new black hole solutions of gravitational models beyond Einstein theory of relativity, like modified theories of gravitation and loop quantum gravity, as well as Einstein gravity coupled to various matter fields will contain new free parameters (characterizing the gravitational model) in addition to the black hole mass and angular momentum.The signature of these new free parameters describing the underlying modified gravity is encoded in the frequency shift of photons that we measure here on the Earth.Therefore, extending the general relativistic formalism [9][10][11] to general spherically symmetric black hole spacetimes allows us to decode the information of the free parameters of new black hole solutions in modified theories of gravitation and probe new physics beyond Einstein gravity with the help of observational redshift.
In the present study, we aim to extend the general relativistic formalism [9][10][11] to general spherically symmetric black hole spacetimes with the flat asymptote that would help to estimate or constrain the free parameters of the extended theories of gravitation.Then, we explore regular black holes in conformal gravity as a special example to demonstrate the application of this extension.This paper is outlined as follows.In the next section, we briefly introduce the general spherically symmetric and static black hole background in asymptotically flat spacetime.In section III, we consider the geodesic motion of the massive and massless particles in this background and express the nonvanishing components of the 4-velocity of the geodesic massive particles in terms of the metric functions that characterize the black hole parameters.Then, we continue to find the light bending parameter in terms of the metric functions with the help of the geodesic motion of emitted photons.In section IV, we obtain a general expression for the total observational frequency shift of photons emitted by test particles circularly orbiting black holes in general spherically symmetric backgrounds in terms of the metric components by using the results obtained in Sec.III.In Sec.V, we employ the relations obtained in Secs.III-IV in order to analyze the effects of free parameters of nonsingular black holes in conformal gravity on the observational frequency shift, hence showing an application of our general results.Finally, some concluding remarks are presented.

II. BLACK HOLE SPACETIME CONFIGURATION
Here, we take into account a general static and spherically symmetric black hole background with a flat asymptote.The line element takes the following form in the Schwarzschild coordinates (t, r, θ, ϕ) [we use c = 1 = G units] ds 2 = g µν dx µ dx ν = g tt dt 2 + g rr dr 2 + g θθ dθ 2 + g ϕϕ dϕ 2 , (1) where the metric components g tt (r), g rr (r), and g θθ (r) are arbitrary functions of the radial coordinate r as well as g ϕϕ (r, θ) = g θθ (r) sin 2 θ.Hereafter, we remove the arguments r and θ from the metric functions for notational simplicity [22].
Since we are interested in asymptotically flat spacetimes, the metric components of the line element (1) are required to behave as lim r→∞ at spatial infinity.In this paper, we are interested in investigating the effects of spacetime curvature, produced by a general static and spherically symmetric black hole described by the metric (1), on the frequency shift of photons emitted by geodesic massive particles orbiting the black hole.For instance, the photon sources could be stars, gas or water masers in the case of AGNs.

III. GEODESICS OF TIMELIKE AND NULL PARTICLES
The trajectory of particles revolving around the spherically symmetric background (1) is governed by the following equation of motion for both timelike (κ = −1) and null (κ = 0) particles.
Besides, ẋµ = dx µ /dλ denotes the derivative of the particle's position where λ is the affine parameter for photons and it represents the proper time λ = τ in the case of massive particles.Moreover, the background spacetime (1) possesses two Killing vector fields as follows ψ µ = (0, 0, 0, 1) rotational Killing vector field, (5) which characterize the symmetries of the spacetime and can be employed to obtain the conserved quantities of the geodesic motion of the test particles orbiting the black hole.
Due to the spherical symmetry of the line element (1), we restrict the particle's motion to the equatorial plane θ = π/2 without loss of generality.For inclined orbits of geodesic massive particles as seen from the Earth, one can apply the rotation to the frame by Euler angles.However, in the case of the accretion disks which can be only observed from an edge-on view from the Earth [23][24][25][26], this equatorial plane has a great interest.

A. Geodesics of massive particles
The geodesic equation (3) for massive particles reads where U µ is the 4-velocity of the particle.On the other hand, the conserved quantities associated with the Killing vector fields ( 4) and ( 5) are given by where E and L are, respectively, the total energy and angular momentum of the test particle per unit mass.Now, by considering the particle motion confined in the equatorial plane where U θ = 0, and substituting these relations into Eq.( 6), we have which has the energy conservation law structure such that the first term is the kinetic energy of the particle moving in an effective potential with the following form From Eq. ( 9), the radial velocity of the massive particle reads that is a function of the radial coordinate r only.
Taking into account the real astrophysical systems and spherical symmetry of the spacetime background, investigating circular equatorial motion is relevant.Besides, taking θ = π/2 allows us to express some relations in a simpler way without loss of generality.The following conditions govern the massive particles in circular orbits where the prime symbol denotes the derivative with respect to the radial coordinate r.From these conditions, one can find E and L in terms of the metric functions g tt and g ϕϕ and their derivatives as below Now, with the aid of Eqs. ( 7) and ( 8), t− and ϕ− components of 4-velocity can be expressed as follows where subscript e refers to the emitter radius r e .For stability of the circular orbits, it is required for the second derivative of V r (r) to satisfy the following condition where the equality indicates the innermost stable circular orbit (ISCO) and characterizes the inner edge of the accretion disk.It is possible to find an explicit expression for the second derivative of V r (r) in terms of the metric functions, with the aid of Eqs. ( 14) and (15), which reads

B. Geodesics of null particles
In the case of massless particles, which follow null geodesics, the equation of motion (3) reads where k µ = (k t , k r , k θ , k ϕ ) is the 4-wave vector of the photons.These photons move outside the event horizon described by the metric (1), which depends on the explicit form of the functions g tt , g rr , and g ϕϕ .By considering the motion of photons on the equatorial plane, the component k θ of the 4-wave vector vanishes, hence Eq. ( 20) reduces to Due to the symmetries of the spacetime described by the Killing vector fields (4)-( 5), the following conserved quantities also hold for the photon motion where E γ and L γ are the energy and angular momentum of the photons, respectively.Now, by substituting these relations into the equation of motion of photons (21), we get Note that for the points where the radial component of the 4-wave vector k r vanishes (which are the diametrically opposite points whose joint line is perpendicular to the line of sight), Eq. ( 24) simplifies in the following way By introducing the following definition of the light bending parameter in the latter equation we can find this parameter in terms of the properties of the spacetime as follows We recall that this relation shows the deflection of light sources when they are located on either side of the black hole at the midline characterized by ± sign.

IV. FREQUENCY SHIFT
The frequency of photons at some position x µ p = x t , x r , x θ , x ϕ | p reads where the index p refers to either the point of emission x µ e or detection x µ d of the photon.
The most general expression for shifts in the frequency ω p in static and spherically symmetric backgrounds of the form (1) can be written as [9][10][11] In addition, since the motion of the massive test particles is restricted to the equatorial plane, implying that U θ = 0, and from the circular orbit condition U r = 0, the frequency shift formula (29) reduces to By considering the observer at a very large distance from the black hole, its 4-velocity is given by U µ d = (1, 0, 0, 0), and therefore the frequency shift simplifies to where we used the light deflection parameter defined in Eq. ( 26).
One may note that U ϕ e presented in Eq. ( 30) also posses ± sings according to Eqs. ( 15) and ( 17) which amounts to the rotation direction of massive particles (+/− for counterclockwise/clockwise).However, due to the spherical symmetry of the background, it does not matter whether the particle is co-rotating or counterrotating, hence we choose the plus sign without loss of generality.
Finally, by substituting Eqs. ( 16), ( 17) and ( 27) into Eq.( 31), one obtains the following general expression for the shift in the frequency of photon sources circularly orbiting a general spherically symmetric black hole spacetime at the radius r = r e .In this relation, z BH 1 /z BH 2 corresponds to the +/− sign that represents the redshift/blueshift at the midline and on either side of the line of sight.It is important to remark here that the aforementioned expression is the total frequency shift where the first term corresponds to the gravitational frequency shift and the second term is the kinematic frequency shift such that

V. BLACK HOLES IN CONFORMAL GRAVITY AS A SPECIAL EXAMPLE
Here, in order to show the application of the results obtained in the previous sections regarding the frequency shift of test particles revolving spherically symmetric spacetimes, we take into account nonsingular black holes in conformal gravity as a special example.
The curvature singularity at the center of black holes is a crucial and outstanding problem and we need a quantum version of gravity to describe the origin.Since there is no self-consistent and complete quantum gravity theory so far, one can extend the gravitational theory beyond general relativity in order to cure this essential singularity.
In this context, conformal gravity is an extended theory of gravitation that provides a way to get rid of singularities by means of a conformal transformation applied to the Schwarzschild metric g µν , giving rise to a new metric ĝµν which is conformally related to g µν by the conformal factor Ω(x) as follows ĝµν (x) = Ω (x) g µν (x) , (34) where Ω (x) is a function of the spacetime coordinates that satisfies In this regard, one assigns the conformal symmetry to the Schwarzschild background, hence the spacetime singularity can be removed after a suitable conformal transformation.Nonsingular black holes conformally related to the Schwarzschild solutions and the Kerr solutions have been proposed in [27,28] in which their spacetime is geodesically complete.In this paper, we employ the following conformal factor [28] where l is a length factor and N is a positive integer.Therefore, the new line element dŝ 2 conformally related to the Schwarzschild black holes ds 2 Schw = g µν dx µ dx ν with the conformal factor Ω(r) is given by dŝ 2 = Ω(r)ds 2 Schw = Ω(r)g µν dx µ dx ν = ĝµν dx µ dx ν , (37) where the Schwarzschild metric components are as follows in which M is the total mass of the Schwarzschild black hole.Now, by taking into account Eqs. ( 36)-( 38), the explicit form of the metric components of the conformal nonsingular black holes reads ĝtt As the next step, we can consider massive geodesic test particles circularly orbiting the conformal black hole described by (37) and apply all the formalism developed previously.Besides, by substituting the expressions for the metric components (39)-(42) into Eq.( 19), r ISCO can be obtained by solving and we have stable orbits for the radii r ≥ r ISCO .Note that this equation is a 5th order polynomial which cannot be solved analytically and we shall use its numerical solutions in our future analysis.A numerical investigation of the solutions to Eq. ( 43) in black hole mass unit for N = 1, 2 and 0 ≤ l ≤ 1 indicates that there is only one positive definite root larger than the photon sphere radius r ph and the rest of the roots are either imaginary/complex-valued or smaller than r ph [r ph can be obtained through the condition (g ′ tt g ϕϕ ) r=r ph = g tt g ′ ϕϕ r=r ph and this constraint is deduced from Eq. ( 24) by considering the conditions k r = 0 and (k r ) ′ = 0].Hence, the positive definite root of Eq. ( 43) satisfying r ISCO > r ph indicates the ISCO radius.
However, the black curves in Figs.1-2 show the dependency of r ISCO on the free parameters N and l of the theory.From these figures, we see that for nonvanishing values of the free parameters, r ISCO in conformal gravity is always less than r (Schw) ISCO = 6 for the Schwarzschild black hole case, r ISCO .It is also important to remark that although it is possible to have circular orbits for r e < r ISCO , these orbits are not stable.Hence, the region of interest corresponds to values of r greater than r ISCO for which bounded orbits are stable.
On the other hand, from Eqs. ( 32) and ( 39)-(42), the total frequency shift in the conformal black hole background reads where one can identify and as, respectively, the gravitational frequency shift and kinematic frequency shift satisfying z 1,2 = z g + z kin± .Note that these formulas reduce to the Schwarzschild black hole case [16] for either N = 0 or l = 0, as it should be.In addition, since z 1,2 are observable quantities, the aforementioned relations can be employed alongside the real astrophysical data in order to estimate (or constraint) the free parameters N and l appearing in this model of conformal gravity.
Note that the length scale parameter l has been constrained with the help of x-ray observational data of the supermassive black hole in 1H0707-495 to be l ≤ 0.6 for N = 2 [29] and l ≤ 0.225 for N = 1 [30].Therefore, we concentrate our attention on the restricted values N = 1, 2 and 0 ≤ l ≤ 1 for the free parameters of the conformal gravity in our analysis.Figs.1-2 illustrate the density plots of the redshift and blueshift in the r e − l parameter space for N = 1, 2. As one can see from these figures, r e and l have opposite effects on the redshift/blueshift and as r e (l) increases, the shift in frequency decreases (increases).Therefore, it is expected that these parameters balance the frequency shift at some points, denoted by the continuous gray curves in Figs. 1  and 2, which represent constant level sets for z 1,2 .
In addition, from Fig. 3, one can see that the redshift, blueshift, and gravitational redshift increase when the free parameter N increases.Moreover, this figure shows that as the parameter N takes greater values, the frequency shift curves increase faster as a function of l since N appears as an exponent in the conformal factor Ω. As the free parameters take higher values, z 1,2 significantly deviates from the Schwarzschild redshift/blueshift depicted in Fig. 3 as N = 0.It is worth mentioning that although z g is positive and increases as the free parameters increase (see the right panel of Fig. 3), z 2 is still negative and acquires lower values for nonvanishing N and l since z kin− -term in Eq. ( 44) is dominant (see the middle panel of Fig. 3).Furthermore, by defining R = 1 + z 1 and B = 1 + z 2 , one can readily verify which leads to an expression for the mass of the nonsingular black hole in terms of the observable frequency shifts z 1 and z 2 as well as the free parameters l and N of the conformal gravity theory in the following way where reduces to the Schwarzschild black hole mass formula for either vanishing N or l, as it should be.It is worth mentioning that the maximum values of z 1 and z 2 for revolving particles in stable circular orbits occur at the radius r ISCO .In this regard, the frequency shift values on the black curves in Figs.1-2 indicate the possible maximum values of z 1 and z 2 for each black hole configuration.In addition, figure 3 is evaluated for the emitter orbiting circularly the black hole at r e = r ISCO , hence the curves indicate the maximum values of the frequency shift.In this picture, the frequency shift vanishes as the emitter radius tends to infinity.
From an observational point of view, we can refer to the measured values of z 1 and z 2 from H 2 O megamaser systems on the accretion disks of supermassive black holes hosted at the core of AGNs.For instance, observations from the megamaser system at the center of NGC 4258 galaxy show that the redshift lies in the range 4 × 10 −3 < z 1 < 6 × 10 −3 , the blueshift is within −1.7 × 10 −3 < z 2 < −1.0 × 10 −3 , and the emitter radius r e extends in the sub-parsec region from 0.04pc to 0.5pc from the center of the disk [31].Similarly, from the data set of the megamaser systems NGC 1194, NGC 2273, NGC 2960, NGC 6264, NGC 6323, and UGC 3789, we find the ranges 1.1 × 10 −3 < z 1 < 3.3 × 10 −3 for the redshift, −2.6 × 10 −3 < z 2 < −0.9 × 10 −3 for the blueshift, and 0.028 pc < r e < 1.33 pc for the emitter radius [18,32,33].One may note that applying our general relativistic formalism to such astrophysical systems, including the ones mentioned above, allows us to quantify the gravitational redshift z g (45), hence identifying a general relativistic effect in these systems [16][17][18][19].In the case of megamaser systems, it was shown that the gravitational redshift ranges 1.6 × 10 −7 < z g < 3.1 × 10 −5 , depending on the mass of the supermassive black hole and its distance to the emitter.
In addition, it is worth noting that by taking into account the frequency shift formulas (44) and the conformal factor (36) for the special case N = 1, one can show that Now, by employing the trigonometric solution of cubic equations for the irreducible case [34], we solve these relations for M and l to find the nonsingular black hole mass and length scale only in terms of the observational redshift/blueshift and orbital parameter of the emitter as follows where the parameter σ is given by and the solutions M and l + for p = π are valid for M re < 1 6 , whereas M and l − for vanishing p = 0 are valid for M re ≥ 1 6 .

VI. DISCUSSION AND FINAL REMARKS
In this work, we have extended a general relativistic method [9][10][11] for measuring the black hole parameters to general spherically symmetric spacetimes by considering massive geodesic particles circularly orbiting a static black hole.Then, we have analytically obtained a general formula for the observational frequency shift of photons emitted by the test particles in terms of the metric functions and their derivatives that characterize the black hole parameters.
In addition, we have studied a special case by applying the former results to a nonsingular black hole conformally related to the Schwarzchild solutions as a concrete example of this general relativistic approach.The nonsingular black hole mass was expressed in terms of the observational redshift/blueshift.Furthermore, we have investigated the effects of the free parameters of the conformal gravity theory on the observational redshift and compared results with those of the standard Schwarzschild black hole.Specifically, we have seen that for nonvanishing values of the free parameters N and l of the conformal gravity theory, the total observational redshift/blueshift increases.
Furthermore, it would be interesting to estimate (or constraint) the free parameters N and l with the help of real astrophysical data of supermassive black holes hosted at the core of AGNs by making use of Bayesian fitting methods and comparing results to the previous estimations based on x-ray data [29,30].We leave this study for future work.
Within this general relativistic formalism, the detector receives the information of spacetime background encoded in the frequency shift of photons emitted by the photon source.If the spacetime background is constructed based on gravitational theories beyond Einstein gravity, it would characterized by some new free parameters of the theory, hence the information of these free parameters encoded in the observational redshift.In this regard, if we develop the mathematical modeling of the general relativistic formalism [9][10][11] to general spherically symmetric spacetimes, this would help us to extract the information of the free parameters of the underlying theory from observations.Therefore, the generalization of this formalism to general spherically symmetric spacetimes, that contain an arbitrary number of the free parameters and the metric functions depend on the modified gravity, is a crucial extension of this approach and could be useful to extract information on black hole parameters and test extended theories of gravity.

4 FIG. 1 : 26 FIG. 2 :
FIG. 1: The redshift in the re − l plane for M = 1.The gray curves stand for constant values of z1 and the black curve denotes re = rISCO.

FIG. 3 :
FIG.3:The redshift z1 (left), blueshift z2 (middle) and gravitational redshift zg (right) versus the length scale l for various values of N and M = 1.The emitter is orbiting circularly at re = rISCO and the continuous green curves refer to the frequency shift in the Schwarzschild spacetime where N = 0.