## Shadow of a Renormalization Group Improved Rotating Black Hole

Luis A. S´ anchez 1, ∗

1 Departamento de F´ ısica, Universidad Nacional de Colombia, A.A. 3840, Medell´ ın, Colombia

## Abstract

We present a study on quantum gravity effects on the shadow of a rotating black hole (BH) obtained in the setting of the asymptotically safe gravity. The rotating metric, which results from a static regular one recently presented in the literature, is generated by using the generalized Newman-Janis algorithm. The novelty of the static regular metric lies in the fact that it is the outcome of an effective Lagrangian which describes dust whose spherically symmetric collapse is non-singular as a consequence of the antiscreening character of gravity at small distances. The effective Lagrangian includes a multiplicative coupling, denoted as χ , with the Lagrangian of the collapsing fluid. The resulting exterior metric for large radii depends on a free parameter ξ which captures the quantum gravity effects. The form of the coupling χ and its connection with the quantum parameter ξ are determined by the running of the Newton coupling G k ( ) along a renormalization group trajectory that stops at the ultraviolet non-gaussian fixed point of the asymptotic safety theory for quantum gravity. Varying both the spin parameter a /star and the quantum parameter ξ , we explore the quantum gravity effects on several astronomical observables used to describe the morphology of the shadow cast by rotating BHs. In order to obtain constraints on the parameter ξ , we confront our results with the recent Event Horizon Telescope (EHT) observations of the shadows of the supermassive BHs M87 ∗ and Sgr A . ∗ We find that the ranges of variation of all the studied shadow observables fall entirely within the ranges determined by the EHT collaboration. We then conclude that the current astronomical data do not rule out the renormalization group improved rotating BH.

keywords : Quantum aspects of black holes, asymptotic safety, black hole shadow.

∗ lasanche@unal.edu.co

## I. INTRODUCTION

It is well known that, in the classical theory of GR, singularity theorems establish that, as long as certain energy conditions are satisfied, the gravitational collapse of a massive body inevitably gives rise to a curvature singularity that renders geodesic incomplete the space-time geometry [1-3]. Moreover, the cosmic censorship hypothesis affirms, in its weak form, that the singularity arising from the gravitational collapse to a BH is concealed by an event horizon which hidden it from detection by any external observer [4]. The strong version positions this conjecture at the level of a universal principle as asserts that every singularity arising from gravitational collapse is always hidden behind an event horizon. So, the strong version excludes the physical existence of naked singularities. Since gravitational singularities (or singular space-times) imply the loss of predictability of the classical GR, they are considered as unphysical pathologies in the ultraviolet (UV) regime where the theory breaks down. Thus, there is a general consensus that the solution to the so-called 'singularity problem' requires a quantum theory of gravity that, at trans-Planckian energies, should supersede classical GR and that should provide us with non-singular BH solutions (see Ref. [5] for a recent review on regular BH).

Among the several proposals for a quantum theory of gravity, asymptotic safety (AS) is a promising contender that, using the methods of the standard quantum field theory and the techniques of the functional renormalization group (FRG), gives rise to a non-perturbative renormalizable quantum theory of the gravitational interaction [6]. The basic ingredient of the theory is the conjecture of the existence of an UV non-gaussian fixed point which attracts the flow of the scale dependent quantum effective action for gravity. The fixed point determines a finite transition scale, identified as the Planck scale, beyond which quantum scale invariance is realized such that a finite number of relevant dimensionless couplings, as the Newton coupling G k k ( ) 2 , become constant without dependence of the energy scale. In contrast, the corresponding dimensionfull Newton coupling G k ( ) falls down running as the inverse squared of the energy scale and, consequently, quantum fluctuations weaken the gravitational interaction in the UV [7-9]. It is just this antiscreening effect of the quantum fluctuations at high energies that opens the possibility of solving the singularity trouble in the AS scenario. The running Newton coupling leads to the construction of the so-called renormalization group improved (RGI) metrics which allow including the leading quantum gravity effects on gravitational phenomena in the strong field regime [10-14]. In this regard, the problem of the singularity-free gravitational collapse in AS was studied for the first time in [15] where it was concluded that stronger deviations from GR than the provided by AS, are required to prevent the final singularity. In [16], by analyzing the modifications induced by AS gravity on the gravitational collapse of generic Vaidya spacetimes, Misner-Sharp mass functions were obtained that lead to infer a softening of the singularity as compared with GR. The problem was revisited in [17] where, by assuming a RGI regular exterior metric, the effective interior RGI metric is deduced from compatibility with the exterior one, resulting in an regular interior dust collapse. Indications that the weakening of gravity in the UV in the context of the AS theory could solve the singularity problem, have also been discussed in [10, 18-20].

Recently, a more satisfactory top-down approach to the singularity resolution has been adopted by constructing an effective Lagrangian for the gravitational collapse of dust that includes a multiplicative coupling χ with the matter Lagrangian [21]. The form of this coupling is obtained from the running of the Newton coupling G k ( ) along a renormalization group trajectory that ends at the Reuter ultraviolet (UV) fixed point of the AS theory. In the resulting exterior metric, the coupling G k ( ) gives rise to a Misner-Sharp mass function M ( r, ξ ) depending on a free parameter ξ linked to χ by suitable junction conditions. The exterior metric remains everywhere regular

without the existence of singularities and, for large radii, the associated RGI static metric is a Schwarzschild-like one from which the classical spherically symmetric line element is recovered when ξ = 0. Since there is growing observational support that astrophysical BH candidates rotate, it is pertinent to construct the RGI rotating counterpart of the static metric and verify if the RGI rotating version remains regular. For this purpose we generate the rotating metric through the Newman-Janis (NJ) algorithm [22] and, in the light of a theorem that establish the necessary and sufficient conditions for the absence of scalar curvature singularities [23], we discuss the issue of the regularity of the RGI rotating BH.

On the other hand, the recent observation of the shadow images of the supermassive BHs M87 ∗ and Sgr A ∗ by the Event Horizon Telescope (EHT) collaboration [24, 25], has opened new possibilities to investigate gravity theories in the strong field regime. The important fact about these images is that they are determined only by the space-time background metric making visible horizon scale features of the BHs. Immediately after the publication of the EHT results, a considerable number of interesting papers, focused on testing GR, alternative theories of gravity as well as quantum gravity scenarios, have appeared in the literature (see e.g. [26] and references therein for a large and, already, incomplete list of works). In the context of the AS theory, investigations have been done to identify quantum gravity modifications to the shadows of both non-rotating and rotating regular RGI BHs, and to constraint the corresponding quantum parameter [27-30].

In this paper we aim to study quantum gravity effects on the shadow cast by a BH whose specetime is described by the RGI rotating metric that, after using the NJ algorithm, results from the static exterior RGI metric obtained in [21]. In this context, we analyze the quantum gravity effects on several selected astronomical observables used to describe the shadow shape. By confronting our results with the recent EHT observations of the shadows of the supermassive BHs M87 ∗ and Sgr A , ∗ we also search for constraints on the quantum parameter ξ appearing in the mass function M ( r, ξ ).

This paper is organized as follows. In Sect. 2 we present a sketch of the procedure in [21] to obtain the exterior metric for a static RGI regular BH from the collapse of dust in the AS scenario. In Sect. 3 we use the NJ algorithm to construct the metric of the RGI rotating BH and we discuss whether this metric is regular or not. In Sect. 4 we study the photon geodesics, the general formulae describing the shadow contour of a rotating BH whose metric involves the mass function M ( r, ξ ), and analyze the shadow of the RGI rotating BH. In Sect.5 we describe the shadow observables that we will use, confront our results with the EHT data, and search for constraints on the parameter ξ . In the last section we present our conclusions.

## II. THE STATIC REGULAR BLACK HOLE IN ASYMPTOTIC SAFETY

In this section we briefly outline the procedure to obtain the static regular BH in AS from the collapse of dust (for details see Ref. [21]). The starting point is the introduction of an action for an effective lagrangian which includes a gravity-matter coupling χ ( E ) with the lagrangian L m = -E of a non-dissipative fluid with rest-mass density ρ and four-velocity u µ satisfying δρ/ρ = δ E / p ( ( E )+ E ) and u u µ µ = -1. This action reads

$$S = \frac { 1 } { 1 6 \pi G _ { 0 } } \int d ^ { 4 } x \sqrt { - g } [ R + 2 \chi ( \mathcal { E } ) \mathcal { L } _ { m } ],$$

where G 0 is the Newton coupling and R is the Ricci scalar. The field equations that follow from the total variation of S are

$$R _ { \mu \nu } - \frac { 1 } { 2 } g _ { \mu \nu } R = 8 \pi G ( \mathcal { E } ) T _ { \mu \nu } - \Lambda ( \mathcal { E } ) g _ { \mu \nu } \equiv T _ { \mu \nu } ^ { \text{eff} ( - ) },$$

where 8 πG ( E ) = ∂ ( E χ /∂ ) E , Λ( E ) = -E 2 ∂χ/∂ E , and T µν = ( ( p E ) + E ) u u µ ν + pg µν , with G ( E ) and Λ( E ) being the effective Newton coupling and the effective cosmological constant, respectively. Now, for a spherically collapsing fluid, the metric in co-moving coordinates is

$$d s _ { ( - ) } ^ { 2 } = g _ { 0 0 } d t ^ { 2 } + g _ { 1 1 } d r ^ { 2 } + g _ { 2 2 } d \Omega ^ { 2 } = - e ^ { - 2 \nu ( r, t ) } d t ^ { 2 } + e ^ { - 2 \psi ( r, t ) } d r ^ { 2 } + C ^ { 2 } ( r, t ) d \Omega ^ { 2 }.$$

where d Ω 2 is the metric on the unit 2-sphere. Then, the field equations acquire the form

$$\frac { 1 } { C ^ { 2 } } \frac { d \mathcal { M } _ { \text{eff} ( - ) } } { d r } \left ( \frac { d C } { d r } \right ) ^ { - 1 } & = \ 8 \pi G ( \mathcal { E } ) \mathcal { E } + \Lambda ( \mathcal { E } ) = \mathcal { E } _ { \mathcal { X } } ( \mathcal { E } ) \equiv \mathcal { E } _ { \text{eff} ( - ) }, \\ 1 & \quad d \mathcal { M } _ { \text{eff} ( - ) } \quad \, \not \, d C \, \cap ^ { - 1 } &$$

$$\frac { \stackrel { \cdot } { d } } { d r } \left ( \frac { \stackrel { \cdot } { d C } } { d t } \right ) \, = \, \frac { d C } { d t } \frac { d \nu } { d r } + \frac { d C } { d r } \frac { d \psi } { d t }, \\. \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \,.$$

$$\frac { \mu } { C ^ { 2 } } \frac { \stackrel { \cdots } { \mathcal { M } } _ { \text{eff} ( - ) } } { d t } & \left ( \frac { \stackrel { \cdots } { d C } } { d t } \right ) ^ { - 1 } = \, - 8 \pi G ( \mathcal { E } ) \mathcal { E } + \Lambda ( \mathcal { E } ) \equiv - p _ { \text{eff} ( - ) }, \\ & d \, \not \, \not \, d C \, \not \, \quad d C \, d \nu \quad d C \, d \psi \quad \stackrel { \cdots } { \cdots }.$$

where M eff( -) is the interior Misner-Sharp mass function given by

$$\mathcal { M } _ { \text{eff} ( - ) } = C \left ( 1 - \left ( \frac { d C } { d r } \right ) ^ { 2 } e ^ { - 2 \psi } + \left ( \frac { d C } { d t } \right ) ^ { 2 } e ^ { - 2 \nu } \right ),$$

After using the Bianchi identity to integrate Eq. (6), the interior metric (3) can be written as

$$d s _ { ( - ) } ^ { 2 } & = - d t ^ { 2 } + \frac { 1 } { 1 - K r ^ { 2 } } \left ( \frac { d C } { d r } \right ) ^ { 2 } d r ^ { 2 } + C ^ { 2 } d \Omega ^ { 2 }, \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{$$

where the integration constant K is the curvature of the 3-space. Using the dimensionless scale factor a t ( ) to rescaling the area-radius C as C = a t r ( ) and rescaling the Misner-Sharp function as M eff( -) = r 3 m eff( -) , Eqs. (4) and (7) produce

$$m _ { \text{eff} ( - ) } ( a ) \, & = \, \frac { a ^ { 3 } } { 3 } \mathcal { E } _ { \text{eff} ( - ) } = a \left ( \left ( \frac { d a } { d t } \right ) ^ { 2 } + K \right ) \\ & = \, - a V ( a ) = - a \left ( - \frac { a ^ { 2 } } { 3 } \int _ { 0 } ^ { \mathcal { E } ( a ) } G ( \mathcal { E } ^ { \prime } ) d \mathcal { E } ^ { \prime } \right ). \\ \text{int the next two ston to connect with the A S gravity is to assume that the evolution of}$$

At this point, the next key step to connect with the AS gravity is to assume that the evolution of G k ( ) with the energy scale k follows a renormalization group trajectory that ends at the UV fixed point g /star and, accordingly, is given by [31]

$$G ( k ) = \frac { G _ { 0 } } { 1 + G _ { 0 } k ^ { 2 } / g _ { * } },$$

where g /star = 570 π/ 833. Here it must be remarked that it is also assumed that the presence of the collapsing fluid does not perturb the flow enough as to alter the fixed point. That this assumption is robust, it has been shown in [32].

If the collapsing fluid is dust, then p = 0, E ∝ a -3 and ρ = E . Additionally, the link between E and k is achieved by using the identification between the proper distance d and the energy scale k in AS gravity, i.e., d ∼ 1 /k which produces d ∼ E -1 2 / , and the running Newton coupling becomes

$$G ( \mathcal { E } ) = \frac { G _ { 0 } } { 1 + \xi \mathcal { E } }, \quad \quad \quad$$

where the parameter ξ , with units of inverse mass to the fourth power, includes the constants G 0 and g /star but it is a free parameter of the AS theory that must be constrained from the observational data. The link implies

$$\chi ( \mathcal { E } ) = \frac { \L o g ( 1 + \xi \mathcal { E } ) } { \xi \mathcal { E }$$

When the present formalism is brought to the Oppenheimer-Snyder approach to the gravitational collapse of dust [33], it turns out that m eff → m 0 , T µν = E u u µ ν and E = 3 m /a 0 3 , and the scale factor obeys the equation

$$\frac { d a } { d t } = - \sqrt { \frac { \text{Log} (1 + 3 m _ { 0 } \xi / a ^ { 3$$

In the limit t →∞ , a t ( ) behaves as

$$a & \sim e ^ { - t ^ { 2 } / 4 \xi },$$

such that, provided ξ &gt; 0, the value a = 0 is never reached. This means that, due to the antiscreeneing character of gravity at small distances in AS, a collapsing boundary C t, r ( b ) = a t r ( ) b at an interior co-moving boundary r = r b , never collapses to a singularity. The space-time is then geodesically complete.

Following [34-36], the matching with the exterior geometry proceeds as follows. First, a generic exterior spherically symmetric line element is assumed

$$d s _ { ( + ) } ^ { 2 } = - f ( R ) d T ^ { 2 } + f ( R ) ^ { - 1 } d R ^ { 2 }$$

with f ( R ) = 1 -2 M ( R /R ) and where now { T, R, θ, φ } are exterior coordinates. Now, for the collapsing boundary C t, r ( b ) = a t r ( ) b , the induced metric on the matching surface Σ is written, in co-moving coordinates, as

$$d s _ { \Sigma } ^ { 2 } = - d t ^ { 2 } + a ^ { 2 } r _ { b } ^ { 2 } d \Omega ^ { 2 }$$

Parameterizing the collapsing boundary as R = R T b ( ), Eq. (16) takes the form

$$d s _ { \Sigma } ^ { 2 } = - \left ( f ( R _ { b } ) - \frac { 1 } { f ( R _ { b } ) } \left ( \frac { d R _ { b } } { d T } \right ) ^ { 2 } \right$$

For the interior metric, the second fundamental form (which reflects the extrinsic geometry of the interior surface), in co-moving coordinates, is

$$K _ { t t } ^ { ( - ) } = 0, \quad K _ { \theta \theta } ^ { ( - ) } = a r _ { b } \sqrt { 1 - K r _ { b } ^ { 2 } }, \\ \quad \quad \quad$$

and the exterior extrinsic curvature is given by

$$K _ { t t } ^ { ( + ) } = - \frac { 1 } { 2 \Delta ( R _ { b } ) } \left ( 2 \frac { d ^ { 2 } R _ { b } } { d t ^ { 2 } } + \frac { d f ( R _ {$$

where

$$\Delta ( R _ { b } ) & = \sqrt { 1 - 2 \frac { \mathcal { M } ( R _ { b } ) } { R _ { b } } + \left ( \frac { d R _ { b } } { d t } \right ) ^ { 2 }$$

Finally, imposing the junction conditions R b ( T t ( )) = a t r ( ) b and

$$K _ { t t } ^ { ( + ) } = K _ { t t } ^ { ( - ) }, \quad K _ { \theta \theta } ^ { ( + ) } = K _ { \theta \theta } ^ { ( - ) },$$

the functional form of the exterior Misner-Sharp mass M ( R,ξ ) is obtained as

$$\mathcal { M } ( R, \xi ) = \frac { R ^ { 3 } } { 6 \xi } \text{Log} \left ( 1 + \frac { 6 \xi M _ { 0 } } { R ^ { 3 } } \right ) \\ \text{--} \quad \text{--}$$

with M 0 = m r / 0 3 b 2. Then, for R &lt;&lt; 1

$$\mathcal { M } ( R, \xi ) = \frac { R ^ { 3 } } { 6 \xi } \text{Log} \left ( \frac { 6 M _ { 0 } \xi } { R ^ { 3 } } \right ) + O ( R ^ { 6 } ), \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

and given that R ≥ R b = a t r ( ) b and a t ( ) &gt; 0, the solution (22) is regular. Note that again ξ &gt; 0 is required in order for M ( R ) to be well defined in this limit. On the other hand, in the IR limit R &gt;&gt; 1 ( ξ &lt;&lt; 1)

$$\mathcal { M } ( R, \xi ) = M _ { 0 } - \frac { 3 M _ { 0 } ^ { 2 } \xi } { R ^ { 3 } } + O ( \xi ^ { 2 } ),$$

and we see that for ξ = 0 the classical constant Schwarzschild BH mass M 0 is recovered. It is then clear that the free parameter ξ captures the quantum gravity effects of the AS scenario.

## III. THE METRIC OF THE RGI ROTATING BLACK HOLE

In the previous section we saw how the matching conditions imposed between the interior and the exterior geometry, make the running Newton coupling G k ( ) of the AS theory turns into the Misner-Sharp mass function M ( R,ξ ) appearing in the Schwarzschild-like form of the exterior metric. Renaming the exterior coordinates { T, R, θ, φ } to the most commonly used { t, r, θ, φ } , the exterior ds 2 (+) ≡ ds 2 static regular metric is given by

$$d s ^ { 2 } & = - \left ( 1 - \frac { 2 \mathcal { M } ( r, \xi ) } { r } \right ) d t ^ { 2 } + \left ( 1 - \frac { 2 \mathcal { M } ( r, \xi ) } { r } \right ) ^ { - 1 } d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 }. \quad \ \ ( 2 5 ) \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

The associated rotating BH is obtained from the generalized NJ algorithm [22], which is a procedure that includes the following steps (for details see [37-40]). In the first step, a general static spherically symmetric metric is considered

$$d s ^ { 2 } = - f ( r ) d t ^ { 2 } + g ( r$$

and the Boyer-Lindquist coordinates ( t, r, θ, φ ) are transformed into advanced null coordinates ( u, r, θ, φ ). Then, the contravariant form of the line element is written in terms of a null tetrad

Z µ a and, to include rotation, the advanced null coordinates x µ are transformed into a new set of complex coordinates ˜ x µ = x µ + ia δ ( µ r -δ µ u ) cos θ , such that the metric coefficients take the new form: f ( r ) → F ( r, a, θ ), g ( r ) → G r,a,θ ( ), and h r ( ) → Σ = r 2 + a 2 cos 2 θ , and such that the null tetrad vectors Z µ a undergo the transformation Z µ a → ˜ Z µ a (˜ x ν ) with the requirement that this transformation recovers the old tetrad and the old line element when ˜ x µ = x µ . In the last step, the new metric is reversed back to Boyer-Lindquist coordinates. In this way, the rotating counterpart of the metric in Eq. (26) takes the form

$$\begin{array} {$$

Comparing now Eqs. (25) and (26) and replacing into Eq. (27), the RGI rotating metric acquires the form

$$d s ^ { 2 } \ & = \ - \Big ( 1 - \frac { 2 \math$$

where a is interpreted as the BH spin and

$$\Delta = r ^ { 2 } - 2 \mathcal { M } ( r, \xi )$$

$$\Sigma = r ^ { 2 } + a ^ { 2 } \cos \theta.$$

Note that in the limit ξ &lt;&lt; 1, M ( r, ξ ) is given by Eq. (24) from which we retrieve the classical expression for ∆ when ξ = 0.

Let us now discuss the regularity of the RGI rotating metric in Eq. (28). This metric is Petrov type D and Segrev type [(1 , 1)(11)] with an algebraically complete set of four second order invariants [37]. For this type of metrics, it has been proved in [23] that if M ( r, ξ ) is a C 3 function, all the second order scalar curvature invariants of the metric (28) are finite at ( r = 0, θ = π/ 2) if, and only if , M (0 , ξ ) = M ′ (0 , ξ ) = M ′′ (0 , ξ ) = 0 (here, the prime denotes derivative with respect to r ). For the mass function M ( r, ξ ) in Eq. (22) is direct to show that: lim r → 0 M ( r, ξ ) = 0, lim r → 0 M ′ ( r, ξ ) = 0, but lim r → 0 M ′′ ( r, ξ ) →∞ . So, in the light of this theorem, the rotating metric does not satisfy the necessary and sufficient regularity conditions. Notice however that, as mentioned in the previous section, the collapsing dust ball never reaches the singularity and therefore, the theorem, which involves an evaluation at r = 0, does not apply. In other words, as it happens for the static exterior metric, given that R ≥ R b = a t r ( ) b and a t ( ) &gt; 0, the rotating metric in Eq. (28) is regular, that is, this metric inherits the absence of singularities of the collapsing dust ball from which it comes.

The BH horizons are the solutions to the transcendental equation ∆ = 0. In the left panel of Fig. 1 we plot ∆( x ) versus x where x = r/M 0 . The plot shows that for fixed values of the spin parameter a /star = a/M 0 and as long as the dimensionless parameter ˜ = ξ ξM 4 0 is lees than a critical value ˜ , ξ c there are two horizons: one inner Cauchy horizon r -and one outer event horizon r + . When ξ = ξ c the two horizons merge, and for ˜ ξ &gt; ξ ˜ c a scalar remnant appears. The right panel of Fig. 1 shows the parameter space ( a/M ,ξ 0 ˜ ). The green curve corresponds to the critical values ˜ ξ c

✁

Figure 1. Left panel: Plot of ∆( x ) versus x , with x = r/M 0 , for a /star = 0 5 and for . ˜ = ˜ ξ ξ (0 5) . c = 0 275119633455 . (black dashed), ˜ ξ &gt; ξ ˜ (0 5) . c (red), ˜ = 0 15 ξ . &lt; ξ ˜ (0 5) . c (blue), and ˜ = 0 05 ξ . &lt; ξ ˜ (0 5) . c (green). Right panel: the parameter space ( a/M ,ξ 0 ˜ ). The green curve corresponds to the critical values ˜ ξ c for each value of a /star .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-f5d29969-09c3-454d-b2d6-bab32e9f9893.jpg)

for each value of a/M 0 in such a way that for ˜ ξ &lt; ξ ˜ c (shaded region) there are RGI rotating BHs, whilst for ˜ ξ &gt; ξ ˜ c (unshaded region) one has a scalar remnant. Consistent with our discussion in the previous section, we will consider for our study values of ˜ satisfying 0 ξ &lt; ξ ˜ ≤ ˜ . ξ c

At this point, it should be highlighted that the horizon becomes more compact as ˜ grows (see ξ left panel of Fig. 1). Moreover, the line element in Eq. (28) has the form of a Kerr-like metric with the BH mass replaced by the 'effective' mass M ( r, ξ ) which, consequently, contains all the information on the quantum gravity effects on the space-time geometry. Fig. 2 are plots of M ( r, ξ ) for growing values of ˜ : ξ ˜ = ˜ ξ ξ (0 99) . c = 0 003425142584 which is the critical value for . a/M 0 = 0 99 . (blue curve), ˜ = 0 0025 (red curve), and ξ . ˜ = 0 0015 (green curve). ξ . Since all of these values of ˜ are allowed for ξ a/M 0 = 0 99, the blue, red, and green vertical dashed lines locate the position . of the outer event horizon r + for this value of a/M 0 and for the referred values of ˜ , ξ respectively. If we fix our attention on the horizon r + for ˜ = ξ ˜ ξ (0 99) . c , we see that even tough the difference in the 'effective' mass hidden behind the horizon is negligible for different values of ˜ , ξ the quantum gravity effect on the location of the horizon is considerable. For example, when passing from r + for ˜ = 0 0025 to ξ . r + for ˜ = ˜ ξ ξ (0 99) . c , the size of the horizon decreases by ∼ 30%. The shrinking of the horizon is then a quantum and strong-field effect and should gives rise to strong distortions of

Figure 2. Plots of the mass function M ( r, ξ ) for a /star = 0 99 and for the values of . ˜ as shown. ξ The vertical dashed lines locate the position of the outer event horizon r + for the same value of a /star and for each one of the values of ˜ . ξ

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-4cd851cc-32a0-4d10-a0ff-e11d81454305.jpg)

the shadow of the RGI rotating BH. As we will see in the next section, this is just what happens specially for high values of a/M 0 and ˜ ξ

## IV. PHOTON GEODESICS AND SHADOW

The space-time with the metric in Eq. (28) is stationary and axisymmetric. There are, then, two Killing vector fields and, consequently, two conserved quantities. They are the energy E and the axial component of the angular momentum L z . Associated to the Killing-Yano tensor field, there is a third conserved quantity, the Carter constant Q . These three quantities are related to each other through the quadratic integral of motion K in the form [41]

$$\mathcal { Q } = \mathcal { K } - ( L _ { z } - a E ) ^ { 2 }.$$

The null-geodesic equations are

$$\Sigma \frac { \dot { r } } { E } \, & = \, \pm \sqrt { \mathcal { R } ( r ) }, \\ \Sigma \frac { \dot { \theta } } { = } \, & = \, \pm \sqrt { \Theta ( \theta ) }.$$

$$\Sigma \frac { \stackrel { \sim } { \dot { \phi } } } { E } \ = \ \frac { \eta ^ { 2 } } { \sin ^ { 2 } \theta } + \frac { a } { \Delta } \left [ 2 r \mathcal { M } ( r, \xi ) - a \eta \right ],$$

$$\Sigma \frac { \stackrel { R } { \theta } } { E } \, & = \, \pm \sqrt { \Theta ( \theta ) }, \\ \Sigma \frac { \stackrel { \stackrel { \beta } { \phi } } { } } { \, } \, & = \, \frac { \eta ^ { 2 } } { \lambda } + \frac { a } { \pi } \left [ 2 r \mathcal { M } ( r. \xi ) - a \eta \right ].$$

$$\Sigma _ { \frac { i } { E } } ^ { \iota } \, = \, - a ^ { 2 } \sin ^ { 2 } \theta + \frac { 1 } { \Delta } [ \left ( r ^ { 2 } + a ^ { 2 } \right ) ^ { 2 } - 2 a r \mathcal { M } ( r, \xi ) \eta ] \,, \\ \text{dots denote derivative respect to an affine parameter along the geodesic and}$$

where the dots denote derivative respect to an affine parameter along the geodesic and

r

$$\mathcal { R } ( r ) \ & = \ r ^ { 4 } + \left ( a ^ { 2 } - \eta ^ { 2 } - \chi \right ) r ^ { 2 } + 2 r \left [ \chi + ( \eta - a ) ^ { 2 } \right ] \mathcal { M } ( r, \xi ) - a ^ { 2 } \chi, \\ \Theta ( \theta ) \ & = \ \chi + a ^ { 2 } \cos ^ { 2 } \theta - \eta ^ { 2 } \cot ^ { 2 } \theta, \\. \. \. \. \. \. \. \. \. \.$$

χ

+(

η

a

)

2

(

r, ξ

)

a χ,

2

(36)

and where the two conserved parameters χ and η are defined as

$$\chi = \frac { \mathcal { Q } } { E ^ { 2 } }, \quad \eta = \frac { L _ { z } } { E }.$$

As stated above, the shadow of a BH is a pure strong-field effect that probes gravity theories on horizon scales and, therefore, the background spacetime geometry. The boundary of the BH shadow corresponds to the innermost unstable spherical photon orbit with radius r p determined by the conditions R ( r ) = 0 and ∂ R ( r /∂r ) = 0, which provide the following expressions for the conserved impact parameters [42]

$$\chi \, & = \, \frac { \cdot \, \cdot \, \cdot } { a ^ { 2 } ( r - \mathcal { M } - r \mathcal { M } ^ { \prime } ) ^ { 2 } } \left [ 4 a ^ { 2 } \mathcal { M } - r ( r - 3 \mathcal { M } ) ^ { 2 } - 2 r \mathcal { M } ^ { \prime } ( r ^ { 2 } - 3 r \mathcal { M } + 2 a ^ { 2 } ) - r ^ { 3 } ( \mathcal { M } ^ { \prime } ) ^ { 2 } \right ], \, \quad & ( 3 9 ) \\ \eta \, & = \, \frac { ( r ^ { 2 } - a ^ { 2 } ) \mathcal { M } - r \Delta - r ( r ^ { 2 } + a ^ { 2 } ) \mathcal { M } ^ { \prime } } { a ( r - \mathcal { M } - r \mathcal { M } ^ { \prime } ) }, \\ \text{where the prime denotes derivative with respect to $r$}.$$

where the prime denotes derivative with respect to r .

Let us now consider photon sources located at infinity and distributed uniformly in all directions, and an observer also at infinity. The shadow of the BH is obtained by solving the scattering problem of photons falling onto the BH from infinity as seen by the observer. The celestial coordinates of

Figure 3. Shadows cast for the RGI rotating BH for the observer inclination angle ι = π/ 2. Each shadow silhouette is constructed for the values ( a , ξ /star ˜ ( a /star ) c ) as shown. The left cusp effect fades as a /star decreases.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-6592203b-ec16-4237-a0a0-452c478a2e79.jpg)

the observer ( x, y ) are the impact parameters of photons falling just into the innermost unstable photon orbit and correspond to the apparent angular distances of the image on the celestial sphere measured from the direction of the observer's line of sight. In this sense, x and y are measured in the directions perpendicular and parallel to the projected rotation axis onto the celestial sphere, respectively. Let ι be the inclination angle of the observer defined as the angle between its line of sight and the rotation axis of the BH. Then, with ( p ( ) t , p ( r ) , p ( θ ) , p ( φ ) ) being the tetrad components of the photon momentum respect to a locally non-rotating reference frame, the coordinates ( x, y ) can be written in terms of the parameters ( χ, η ) and the angle ι as

$$x \, = \, \lim _ { r \to \infty } \frac { - r p ^ { ( \phi ) } } { p ^ { ( t ) } } = - \eta \csc \iota,$$

$$y \, = \, \lim _ { r \to \infty } \frac { r p ^ { ( \theta ) } } { p ^ { ( t ) } } = \pm \sqrt { \chi + a ^ { 2 } \cos ^ { 2 } \iota - \eta ^ { 2 } \cot ^ { 2 } \iota }, \\ \text{is sign is for prograde orbits and the minus sign is for retrograde orbits. The appear}$$

where the plus sign is for prograde orbits and the minus sign is for retrograde orbits. The apparent shape of the rotating BH (its shadow) is obtained by doing the parametric plot of y r ( ) versus x r ( ) for fixed values of M 0 , a , ι and ξ (in the following we will write M instead of M 0 for simplicity).

## A. The shadow of the RGI rotating black hole

From Eqs. (41) and (42) it follows that the effect of the spin on the shadow contour is more pronounced for an observer inclination angle ι = π/ 2. For this value of ι , Fig. 3 shows the shadow contour for selected values of a /star when ˜ = ˜ ξ ξ c for each one of these values. We see that for ι = π/ 2 and for high values of the spin ( a/M /greaterorsimilar 0 8), . just at the left end of the contour, given by one of the solutions of y r ( ) = 0, appears a kind of cusp-like effect or 'point pressure deformation' from the right that is a smooth dent for spins a /star /lessorsimilar 0 8 but increases as the spin goes to greater . values. The left end of the shadow corresponds to photons in prograde orbits while the right end

Figure 4. The same as in Fig. 3 but including values of ˜ less than ξ ˜ ξ c and for a /star = 0 99 (left panel), . a /star = 0 7 . (central panel), and a /star = 0 3 (right panel). . For each value of a /star , the comparison with the shadow of a Kerr BH ( ˜ = 0, black dashed contour) is done. ξ

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-8739b4b2-b3d8-47f3-8cd3-a61a8603f65a.jpg)

corresponds to photons in retrograde orbits. Consequently, the frame-dragging effect is stronger for prograde photons which leads to a greater distortion of the shadow of the RGI rotating BH than for the shadow of a classical Kerr BH. This is in line with the discussion in the previous section where we show that high values of the BH spin and high values of the parameter ˜ lead to a more ξ compact horizon and, therefore, to prograde photons traveling through a region of higher curvature.

Fig. 4 shows the same as Fig. 3 but in three separate panels in order to make evident two additional shadow features: (i) that the left cusp effect appears only in the immediate vicinity of ˜ = ˜ ξ ξ c (see the left panel), and (ii) that the area enclosed by the shadow contour of the RGI rotating BH is always smaller that the one enclosed by the shadow of a Kerr BH ( ˜ = 0) which is shown as the ξ black-dashed contour in each one of the three panels. Fig. 5 reveals that the cusp effect is not detectable for observer inclination angles ι different from π/ 2, and that the increase in the spin value is accompanied by a displacement of the shadow to the right along the x -axis due also to

Figure 5. Dependence of the shadow of the RGI rotating BH on the observer inclination angle. Left panel: ι = 17 ◦ and different values of ( a , ξ /star ˜ ( a /star ) c ). Right panel: The same as in the left panel but for ι = 30 . ◦

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-5e441e8e-8e51-4810-a1bd-17f835642c7d.jpg)

frame-dragging.

Interestingly, the same cusp effect that we notice is reported in Refs. [27, 29]. There the authors succeed in obtaining a rotating regular metric in AS by using invariant scalars (as the Kretschmann scalar K ) as the identification criterion with the energy scale k in the form k 2 = αK 1 2 / (where α is a dimensionless constant of order one), instead of the energy scale identification with the proper distance k ∼ 1 /d discussed in Sec. 2. This seems to mean that for any RGI regular rotating BH, not only the presence of a smooth dent and a smaller shadow area are robust, but so is the presence of the cusp effect for ι = π/ 2, for high values of a /star and for ˜ ξ ≈ ˜ . ξ c Obviously, as pointed out in [27], if the quantum parameter ξ is strictly tied to the Planck scale ξ ∼ m -4 p , as is done in AS, this effect is too tiny to be detectable with the current observational facilities (sew the following section).

## V. SHADOW OBSERVABLES AND COMPARISON WITH THE EHT DATA

In this section we aim to constraint the quantum parameter ˜ using the recently EHT observaξ tions of the shadow images of the supermassive black holes M87 ∗ and Sgr A ∗ [24, 25].

To characterize the shape of the shadow, several authors have proposed different sets of distortion parameters that can be measured from the observed EHT shadow images [43-50]. Here we use the following three parameters:

- (i) The areal radius r A , which quantifies the size of the shadow and is defined by

$$r _ { A } = \sqrt { \frac { \mathcal { A } } { \pi } },$$

where A is the area enclosed by the contour of the shadow

$$\mathcal { A } = 2 \int y ( r ) d x ( r ) = 2 \int _ { r _ { p } ^ { - } } ^ { r _ { p } ^ { + } } y ( r ) \frac { d x ( r ) } { d r } d r,$$

and where r -p and r + p are the radii of photons, at unstable circular orbits, that co-rotate and counter-rotate, respectively, with the BH. These are the real positive solutions of y 2 ( r ) = 0.

- (ii) The oblateness D , which quantifies the deformation of the shadow, defined as

$$\mathcal { D } = \frac { \Delta y } { \Delta x } = \frac { y _ { t } - y _ { b } } { x _ { R } - x _ { L } } = \frac { 2 y _ { t } } { x _ { R } - x _ { L } }, \\ \cdot \, - \, - \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdots$$

where the subscripts L R , stand for the left and right ends of the shadow silhouette, and t , b stand for the top and bottom ends. The factor 2 in Eq. (44) and the equality y b = -y t in Eq. (45) are due to the reflection symmetry of the shadow respect to the x -axis as shown in Fig. 6.

- (iii) The fractional deviation parameter δ sh of the shadow diameter from that of a Schwarzschild BH, given by

$$\delta _ { s h } = \frac { R _ { s h } } { 3 \sqrt { 3 } } - 1,$$

Figure 6. The ( x, y ) coordinates of the apparent shadow shape (solid black) of a rotating BH used to define the distortion parameters (for details see the text). The red-dashed curve is the so-called fitting circle which is a circle passing through the three points ( x , y y t ), ( x R , 0) and ( x , y b b ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-51b64989-582f-49df-8485-a5c5e8ae0dd9.jpg)

where R sh is the average radius of the shadow

$$R _ { s h } = \sqrt { \frac { 1 } { 2 \pi } \int _ { 0 } ^ { 2 \pi } \left [ ( x - x _ { c } ) ^ { 2 } + ( y - y _ { c } ) ^ { 2 } \right ] d \phi },$$

with

$$\phi & = \tan ^ { - 1 } \left ( \frac { y } { x - x _ { c } } \right ), \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \end{array}$$

and ( x , y c c = 0) is the geometric center of the shadow with x c calculated from x c = ( x R -x L ) / 2. The left x L and right x R points are the celestial coordinate x r ( ) evaluated at the real positive solutions of y 2 = 0, while the upper point y t of the shadow is obtained from the condition dy/dx = 0, that is

$$( 1 + \mathcal { M } ^ { \prime } ) r ^ { 3 } - 3 \mathcal { M } r ^ { 2 } + a ^ { 2 } ( 1 - s ^ { 2 } ) ( 1 + \mathcal { M } ^ { \prime } ) r + a ^ { 2 } ( 1 - s ^ { 2 } ) \mathcal { M } = 0.$$

where s ≡ sin ι .

## A. Constraints from the M87 ∗ data

The 2017 measurements of the shadow images of M87 , obtained by EHT collaboration, deter-∗ mined the areal radius and the oblateness of this BH, respectively, in the ranges [24]

$$4. 3 1 M \leq r _ { A } \leq 6. 0 8,$$

and

$$1 \leq \mathcal { D } \leq 4 / 3,$$

For the value of the observer inclination angle we take ι = 17 ◦ [52].

Figure 7. Left panel: the areal radius per unit mass r A /M vs log( ˜ ) of the RGI rotating BH for ξ ι = 17 ◦ and for values of a /star in the range 0 1 . &lt; a /star &lt; . 0 99. Right panel: the oblateness D vs log( ˜ ) for the same value of ξ ι and for the same range of values of a /star .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-e6dfd479-f64f-494e-bb6c-1b1b2adc675b.jpg)

For the RGI rotating BH, the left panel of Fig. 7 shows plots of r A /M vs. log ˜ for different values ξ of a /star , while the right panel shows plots of D vs. log ˜ for the same values of the spin parameter. ξ The right end of each curve corresponds to the value of ˜ ξ c for each value of a/M and the parameter ˜ runs in the rage 10 ξ -4 ≤ ˜ ξ ≤ ˜ . ξ c From the left panel we see that the quantum gravity effects render the size of the shadow smaller as a /star increases and that, for fixed value of a /star , the size of the shadow decreases a ˜ grows. The change in size is more noticeable as ξ ˜ ξ → ˜ . ξ c The right panel shows that, in correspondence with the behavior of r A , the shadow silhouette becomes more distorted as a /star becomes larger and, again, that the distortion is greater as ˜ ξ → ˜ . ξ c However, from these plots we also observe that both r A and D cover ranges entirely included within those established by the EHT collaboration, indeed: 4 840390 . ≤ r A /M ≤ 5 191641 and 1 000047 . . ≤ D ≤ 1 014533 for . 0 1 . ≤ a/M ≤ 0 99. . This implies that, even tough the EHT measurements of r A and D are both satisfied, we can not determine constraints on the quantum parameter ˜ of the AS theory. ξ In this sense, the present EHT measurements of M87 ∗ do not rule out the RGI regular rotating BH we are considering.

## B. Constraints from the Sgr A ∗ data

To contrast our results with the EHT measurements for Sgr A ∗ we use the data obtained by the Very Large Telescope Interferometer (VLTI) and Keck observatories for the areal radius and for the fractional deviation parameter δ sh , respectively, as follows [25]

$$4. 5 M & \leq r _ { A } \leq 5. 5, \quad ( \text{VLI} ) \\ 4. 3 M & \leq r _ { A } \leq 5. 3, \quad ( \text{Keck} )$$

$$( 5 2 )$$

$$- 0. 1 7 & \leq \delta _ { s h } \leq 0. 0 1, \quad ( \text{VLTI} ) \\ - 0. 1 4 & \leq \delta _ { s h } \leq 0. 0 5, \quad ( \text{Keck} )$$

The EHT collaboration has determined the observer inclination angle of Sgr A ∗ to be ι &lt; 50 ◦ for two class of simulation models, one with a /star = 0 5 and the other with . a /star = 0 94, . that fully pass a set of EHT constraints for all the general relativistic magnetohydrodynamic (GRMHD) simulations, including the variability constraint. However, in order to fix a value of ι for our calculations, we take ι = 30 ◦ which is the value for the same two class of models passing all the

/a115

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00226v4-figure-cefa6ecc-7b4f-43e7-af3d-a340898eb44f.jpg)

✠

Figure 8. Left panel: the areal radius per unit mass r A /M vs log( ˜ ) of the RGI rotating BH for ξ ι = 30 ◦ and for values of a /star in the range 0 1 . &lt; a /star &lt; 0 99. . Right panel: The fractional deviation parameter δ sh for the same value of ι and for the same range of values of the spin parametr a /star .

EHT constraints but the variability constraints [25]. Our results for r A /M and δ sh are shown in the left and right panels of Fig. 8, respectively, when 10 -4 ≤ ˜ ξ ≤ ˜ ξ c and for different values of a/M . First we see that the set of curves for δ sh are almost a rigid displacement along the y -axis of the set of curves for r A . This is because the values of r A and R sh calculated form the Eqs. (43) and (47), respectively, are such that r A /approxorequal R sh . The behavior of the size of the shadow ( r A ) for different values of a /star and increasing values of ˜ is ξ the same as the already discussed for the M87 ∗ BH. The absolute value of fractional deviation parameter δ sh increases with increasing values of a /star , with the increasing being more conspicuous when ˜ ξ → ˜ . ξ c The values of the areal radius run in the range 4 852527 . ≤ r A /M ≤ 5 193447 . and the values of δ sh cover the range -0 0661306 . ≤ δ sh ≤ -0 0000217. . Consequently, also in this case these ranges are subsets of the ones determined by the EHT collaboration, and we arrive again to the conclusion that the current EHT data for Sgr A ∗ do not allow us to impose constraints on the parameter ˜ and that Sgr A ξ ∗ could be a RGI regular rotating BH.

On the other hand, a recent measurement of the spin of the M87 ∗ BH from its observed twisted light, establishes that this BH rotates clockwise with a spin parameter a /star = 0 9 . ± 0 05 [51]. . Also, in a recent analysis of the spin of Sgr A ∗ using the outflow method, results obtained with the preferred data sets were combined and indicate a practically identical spin parameter a /star = 0 9 . ± 0 06 [52]. . Then if, without pretending to be rigorous, we assume the central value a /star = 0 9 as the spin of both . M87 ∗ and Sgr A , the quantum parameter ∗ ˜ is constrained to be in the range 0 ξ ≤ ˜ ξ ≤ ˜ ξ (0 9) . c , that is 0 ≤ ˜ ξ ≤ 0 040967851. . This a very weak constraint if the value of the real physical parameter ξ is conditioned by the Planck scale since, using for M the mass of the BH Sgr A , we get ∗ ξ = ˜ ξ (0 9) . M /m 4 4 P /lessorsimilar 1 3 . × 10 177 , where m P is the Planck mass. Clearly, the constraint is even weaker if we use the mass of M87 . ∗

## VI. SUMMARY AND CONCLUSIONS

The horizon-scale images of the supermassive BHs M87 ∗ and Sgr A ∗ by the EHT collaboration, have opened an unprecedented window to test, in the strong-field regime, general relativity, alternative theories of gravity, as well as theories of quantum gravity. In this work we have studied quantum gravity effects on the shadow shape of a rotating BH in the framework of the AS theory for quantum gravity. We start from the RGI static regular exterior metric obtained in [21], which

results from the spherically symmetric collapse of dust. The collapse is regular due the to the antiscreening character of gravity at small distances in AS that leads to a ball of dust that never ends in a singularity at any finite time. Next, we have used the generalized NJ algorithm to construct the RGI rotating metric. This metric differs from the Kerr solution by the replacement of the constant BH mass by a Misner-Sharp mass function M ( r, ξ ˜ ) depending only on the radial coordinate and on a dimensionless free parameter ˜ which captures the quantum gravity effects. ξ

The study of the causal structure of the RGI rotating BH space-time has revealed that, associated to each value of the spin parameter a /star (0 &lt; a /star ≤ 1), there is a critical value ˜ ξ ( a /star ) c such that for 0 &lt; ξ &lt; ξ ˜ ˜ c there are two horizons, one outer horizon and one inner Cauchy horizon; for ˜ = ˜ ξ ξ c the two horizons coincide, and for ˜ ξ &gt; ξ ˜ c a horizonless scalar remnant arises. It is important to note that for growing values of the BH spin and growing values of ˜ ( ˜ ξ ξ ≤ ˜ ) the horizon shrinks, that ξ c is, becomes more compact. We have also argued that, being the outcome of a ball of dust that never reaches the singularity, the metric describing the RGI rotating BH is also free of singularities.

When analyzing the shadow cast by the RGI regular rotating BH, we have found the following main features: (i) the area enclosed by its shadow contour is always smaller that the area enclosed by the shadow of a Kerr BH ( ˜ = 0), (ii) ξ only for an observer inclination angle ι = π/ 2 and for values of the spin parameter a /star /greaterorsimilar 0 8 and . ˜ ξ ≈ ˜ , a kind of cusp-like effect appears which increases as ξ c a /star increases (see Fig. 3). Since the left end of the shadow corresponds to photons in prograde orbits, the frame-dragging effect is stronger for these photons than in the case the clasical Kerr BH. This is in accordance with the fact that high values of the BH spin and high values of the parameter ˜ lead ξ to a more compact horizon. The cusp effect bears great similarity to the found in [27, 29] and lead us to argue that the two referred features are robust for shadows cast by RGI regular rotating BHs.

To confront our findings with the EHT images of M87 ∗ and Sgr A , we have used the following ∗ shadow observables: for comparison with the M87 ∗ data we use the areal radius r A , that inform us on the size of the shadow, and the oblateness D , which quantifies the deformation in the shadow, and for comparison with the Sgr A ∗ data we use the areal radius r A and the fractional deviation parameter δ sh of the shadow diameter with respect to the shadow of a Schwarzschild BH ( δ sh,S = 0). Our results for the RGI rotating BH are shown in Figs. 7 and 8. We have found that the range of values of all these observables are proper subsets (in the mathematical sense) of the ranges established by the EHT collaboration for them. Our numerical results for r A are essentially the same as the reported in [53]. We then are forced to conclude that the current EHT data for M87 ∗ and Sgr A ∗ do not allow us to impose constraints on the quantum parameter ˜ and, ξ for the same reason, that the RGI regular rotating BH is not ruled out by these data. We note, however that, as discussed above, if we take for a /star the central value of the recently established spin parameters: a /star = 0 9 . ± 0 05 for M87 . ∗ [51], and a /star = 0 9 . ± 0 06 for Sgr A . ∗ [52], and we use the mass of the BH Sgr A , we can impose on the physical parameter ∗ ξ the very weak constraint ξ = ˜ ξ (0 9) . M /m 4 4 P /lessorsimilar 1 3 . × 10 177 .

It is to be expected that future improvement of the image fidelity of the EHT array will allow to better constrain the image morphology of the BH shadows and, hopefully, to identify quantum gravity effects from them.

## ACKNOWLEDGEMENTS

The author is thankful for useful and insightful discussions with Oleksandr Stashko. The author acknowledge financial support from Universidad Nacional de Colombia.

- [1] R. Penrose, Phys. Rev. Lett. 14 , 57 (1965)
- [2] S.W. Hawking, Proc. R. Soc. Lond. Ser. A, Math. Phys. Sci. 300 (1461), 187 (1967)
- [3] S.W. Hawking, R. Penrose, Proc. R. Soc. Lond. Ser. A 314 , 529 (1970)
- [4] R. Penrose, Nuovo Cimento Numero Speziale 1 , 252 (1969)
- [5] C. Bambi (ed.), Regular Black Holes: Towards a New Paradigm of Gravitational Collapse , Springer, Singapore (2023)
- [6] M. Reuter, Phys. Rev. D 57 , 971, (1998)
- [7] M. Reuter and F. Saueressig, New J. Phys. 14 , 055022 (2012)
- [8] D.F. Litim and F. Sannino, JHEP 1412 , 178 (2014)
- [9] K. Falls, D.F. Litim, K. Nikolakopoulos and C. Rahmede, Phys. Rev. D 93 , 104022 (2016)
- [10] A. Bonanno and M. Reuter, Phys. Rev. D 62 , 043008 (2000)
- [11] M. Reuter and E. Tuiran, Phys. Rev. D 83 , 044041 (2011)
- [12] B. Koch and F. Saueressig, Class. Quan. Grav. 31 , No. 1, 015006 (2014)
- [13] J.M. Pawlowski and D. Stock, Phys. Rev. D 98 , 106008 (2018)
- [14] C. Gonz´ alez and B. Koch, Int. J. Mod. Phys. A 31 , No. 26, 1650141 (2016)
- [15] R. Casadio, S.D.H. Hsu and B. Mirza, Phys. Lett. B 695 , 317 (2011)
- [16] A. Bonanno, B. Koch and A. Platania, Class. Quan. Grav. 34 (9), 095012 (2017)
- [17] R. Torres, Phys. Lett. B 733 , 169 (2014)
- [18] A. Ad´ e` ıf´oba, A. Eichhorn and A. Platania, Class. Quan. Grav., e 35 (22), 225007 (2018)
- [19] L. Bosma, B. Knorr and F. Saueressig, Phys. Rev. Lett. 123 , 101301 (2019)
- [20] R. Torres, Gen. Rel. Grav. 49 (6), 74 (2017)
- [21] A. Bonanno, D. Malafarina, and A. Panassiti, Phys. Rev. Lett. 132 , 031401 (2024)
- [22] E.T. Newman and A.I. Janis, J. Math. Phys. 6 , 915 (1965)
- [23] R. Torres and F. Fayos, Gen. Rel. and Grav. 49 , 2 (2017)
- [24] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 875 , L1 (2019)
- [25] K. Akiyama et al. (Event Horizon Telescope Collaboration), Astrophys. J. Lett. 930 , L12 (2022)
- [26] S. Vagnozzi et al., Class. Quant. Grav. 40 , 165007 (2023)
- [27] A. Held, R. Gold and A. Eichhorn, JCAP 06 , 029 (2019)
- [28] R. Kumar, B.P. Singh and S.G. Ghosh, Annals Phys. 420 , 168252 (2020)
- [29] A. Eichhorn and A. Held, JCAP 05 , 073 (2021)
- [30] A. Eichhorn and A. Held, JCAP 01 , 032 (2023)
- [31] A. Bonanno, T. Denz, J.M. Pawlowski and M. Reichert, SciPost Phys. 12 , 001 (2022)
- [32] A. Baldazzi, K. Falls , Y. Kluth and B. Knorr, arXiv: 2312.03831 [hep-th]
- [33] J.R. Oppenheimer and H. Snyder, Physical Review 56 , 455 (1939)
- [34] W. Israel, Nuovo Cimento B Serie 44 , 1 (1966)
- [35] F. Fayos, X. Ja´n, E. Llanta, and J.M.M. Senovilla, Phys. Rev. D e 45 , 2732 (1992)
- [36] F. Fayos, J.M.M. Senovilla, and R. Torres, Phys. Rev. D 54 , 4862 (1996)
- [37] S.P. Drake and P. Szekeres, Gen. Rel. and Grav. 32 , 445 (2000)
- [38] C. Bambi and L. Modesto, Phys. Lett. B 721 , 329 (2013)
- [39] M. Azreg-A¨ ınou, Phys. Rev. D 90 064041, (2014)
- [40] E. Zakhary and C.B.G. McIntosh, Gen. Rel. and Grav. 29 , 539 (1997)
- [41] S. Chandrasekhar, The Mathematical Theory of Black Holes , Oxford Univ. Press, New York (1992)
- [42] I. Dymnikova and K. Kraav, Universe 5 , 163 (2019)
- [43] K. Hioki and K. Maeda, Phis. Rev. D, 80 , 024042 (2009)
- [44] T. Johannsen and D. Psaltis, Astrophys. J. 718 , 446 (2010)
- [45] N. Tsukamoto, Z. Li and C. Bambi, JCAP 06 , 043 (2014)
- [46] A. A. Abdujabbarov, L. Rezzolla and B.J. Ahmedov, MNRAS 454 , 2423 (2015)

- [47] O.Y. Tsupko, Phys. Rev. D, 95 , 104058 (2017)
- [48] R. Kumar and S.G. Ghosh, Astrophys. J., 892 , 78 (2020)
- [49] M. Afrin, R. Kumar and S.G. Ghosh, MNRAS 504 , 5927 (2021)
- [50] D. Ayzenberg and N. Yunes, Class. Quan. Grav. 35 , 235002 (2018)
- [51] F. Tamburini, B. Thid´ e and M. Della Valle, MNRAS Lett. 492 , L22 (2020)
- [52] R.A. Daly, M. Donahue, C.P. O'Dea, B. Sebastian, D. Haggard and A. Lu, MNRAS 527 , 428 (2024)
- [53] O. Stashko, arXiv: 2407.07892 [gr-qc]
Luis A. Sánchez
Shadow of a Renormalization Group Improved Rotating Black Hole
We present a study on quantum gravity effects on the shadow of a rotating
black hole (BH) obtained in the setting of the asymptotically safe gravity. The
rotating metric, which results from a static regular one recently presented in
the literature, is generated by using the generalized Newman-Janis algorithm.
The novelty of the static regular metric lies in the fact that it is the
outcome of an effective Lagrangian which describes dust whose spherically
symmetric collapse is non-singular as a consequence of the antiscreening
character of gravity at small distances. The effective Lagrangian includes a
multiplicative coupling, denoted as chi, with the Lagrangian of the collapsing
fluid. The resulting exterior metric for large radii depends on a free
parameter xi which captures the quantum gravity effects. The form of the
coupling chi and its connection with the quantum parameter xi are determined by
the running of the Newton coupling G(k) along a renormalization group
trajectory that stops at the ultraviolet non-gaussian fixed point of the
asymptotic safety theory for quantum gravity. Varying both the spin parameter
a/M and the quantum parameter xi, we explore the quantum gravity effects on
several astronomical observables used to describe the morphology of the shadow
cast by rotating BHs. In order to obtain constraints on the parameter xi, we
confront our results with the recent Event Horizon Telescope (EHT) observations
of the shadows of the supermassive BHs M87 and Sgr A. We find that the ranges
of variation of all the studied shadow observables fall entirely within the
ranges determined by the EHT collaboration. We then conclude that the current
astronomical data do not rule out the renormalization group improved rotating
BH.