Shadow and stability of quantum-corrected black holes

Recently the quantum Oppenheimer–Snyder gravitational collapse model has been proposed in loop quantum gravity, providing quantum-corrected Schwarzschild spacetimes as the exterior of the collapsing dust ball. In this paper, the quantum gravity effects on the black hole shadows in this model are studied, and the stability of the quantum-corrected black holes is also analyzed by calculating the quasinormal modes. It turns out that the quantum correction always shrinks the radius of shadows, and the quantum-corrected black holes are stable against the scalar and vector perturbations.


I. INTRODUCTION
The singularity theorems in general relativity (GR) predicts the formation of black holes (BHs) and spacetime singularities [1,2].The appearance of singularities implies the breakdown of GR when spacetime curvature increases unboundedly.Therefore, it is nature to expect that some quantum theory of gravity becomes domain in these regions and can cure the singularities.Up to now, some competitive quantum theories of gravity have been proposed.Among them, loop quantum gravity (LQG), a background-independent and nonperturbative approach to quantum gravity, has received considerable attention and has been widely investigated (see, e.g., [3,4] for books, and [5][6][7][8][9][10] for articles).After thirty years of theoretical research, the canonical (Hamiltonian) and the covariant (Lagrangian) formulations of LQG have achieved individual and remarkable successes.Besides the predictions of discretized geometries and the interpretation of BH entropy in the canonical formulation [11][12][13][14][15][16][17][18], well-defined dynamical descriptions in both formulations were proposed, and the consistency between them was checked in certain senses [19][20][21][22][23][24][25][26][27][28].Gravity coupled to matters were studied and revisited in LQG, leading to the resolution of some long-standing problems [20,29,30].Furthermore, loop quantization program has been successfully extended to high-dimensional theories of gravity as well as alternative gravitational theories [31][32][33][34][35][36][37][38][39][40][41][42][43].Aiming to test the ideas and techniques of full LQG and to quantize the symmetry-reduced phase space of the theory, loop quantum cosmology (LQC) and loop quantum black hole (LQBH) have been studied, and substantial progresses have been made in solving the cosmological big-bang singularity and the BH singularity [44][45][46][47][48][49][50][51][52].Although a systematical derivation of these symmetric-reduced models from LQG is still absent up to now, some progress has been made on the relation between LQG and these models by calculating the expectation value of the Hamiltonian in LQG under certain coherent state peaked at some point in the classical phase space [53][54][55][56][57].
The spherically symmetric gravitational collapse plays an important role in understanding the formation of BH and its singularity in GR.Classically, the gravitational collapse model was constructed by Oppenheimer and Snyder [58] (OS model).In this model, the interior sourced matter is assumed to be a spherically symmetric and homogeneous pressureless dust, and thus it can be modeled by a FRW model.Due to the simplicity, this model is exactly solvable, providing us with a new window to understand more complicated and realistic dynamical processes of gravitational collapse.There is no doubt that to have a complete description of the gravitational collapse, one needs to incorporate the quantum gravity effects.In this direction, the quantum effects of LQG on the gravitational collapse models are being studied [59][60][61][62][63][64][65][66][67][68].In particular, an effective Hamiltonian and an effective metric for the vaccuum exterior solution for the OS collapse model were derived in the effective context of spherically symmetric spacetime [62].An interior to exterior approach for quantum OS model has been proposed recently in [69], where the effective interior spacetime with LQC corrections is carried out to its exterior one by certain matching condition on its boundary surface.It turns out that in preliminary stage the effective exterior spacetime is a quantum-corrected Schwarzschild spacetime.As collapse goes on, a quantum-corrected BH forms with the occurrence of horizon.In the late stage, the collapsing phase will be bounced to an expanding phase, resulting in a transformation from a BH to a white hole.Moreover the global causal structure of the maximal extension of the quantum-corrected Schwarzschild spacetime has been studied in [69].
The successful detection of gravitational-wave from a pair of BHs and the observation of the images of supermassive BH are promoting strongly the observational investigation of BHs [70][71][72].On one hand, the systematical analysis of BH shadow and photon ring has been proposed [73,74], providing a possibility to distinguish various BHs in different gravitational theories.On the other hand, it is widely believed that the study of quasinormal modes (QNMs) of BH plays important roles not only in analyzing the stability of BHs, but also in understanding the gravitational wave signals [75][76][77].These works attract much more attentions in the fields of both arXiv:2211.04263v3[gr-qc] 19 Jul 2023 astrophysics and theoretical physics.In this paper, we will study how the quantum correction to a BH spacetime affects the shadow and QNMs of the BH, by comparing them to those of the classical Schwarzschild BH.We are going to reveal the effects from the effective exterior spacetime, though there are also other quantum gravity effects in BH image [78].
The rest of this paper is organized as follows.In Sec.II, we recall an interior to exterior approach for the OS model to obtain the effective exterior spacetime.In Sec.III, we study the quantum gravitational effects on shadows and QNMs of BHs.Our results are summarized in Sec.IV.The detailed derivation of the junction conditions in the main text will be presented in Appendix A.

II. THE STRATEGY TO GENERATE EXTERIOR SOLUTIONS
In this section, we recall the OS model.In this model, the spacetime manifold M is divided into two regions, the interior region M − and the exterior region M + , by its timelike boundary 3-surface Σ = M − ∩ M + .The former is assumed to be a FRW spacetime region sourced with a spherical dust ball (cloud) with homogenous density with the following line element where τ denotes the proper time of the comoving observer, a(τ) is the scale factor, dΩ 2 ≡ dθ 2 + sin 2 θ dφ 2 , and k = −1, 0, +1 for the open, flat and closed FRW models, respectively.The boundary surface Σ of M − is located at a constant r = r0 , in the comoving coordinates, and can be parametrized by (τ, θ, φ).While the latter is assumed to be a static and spherically symmetric spacetime region with the line element, and can be expressed in Schwarzschild coordinates (t, r, θ, φ) as To match these two regions at Σ, suitable boundary conditions need to be imposed.In the classical theory, the conditions have been studied by Darmois and Israel [79,80], and hence are called as the Darmois-Israel junction conditions.By matching the exterior to the interior of dust collapsing star via the Darmois-Israel junction conditions, the exterior solutions could be generated dynamically by the interior solutions.It turns out that, by matching M + to M − along Σ generated by geodesics, the junction conditions lead to the following relations between the two regions [69,81] (see Appendix A for details) where dτ is the Hubble parameter.Now we impose the dynamical equation for the interior region M − .In the classical theory, the dynamics of the interior spacetime satisfies the Friedmann equation where G is the Newtonian gravitational constant, ρ(τ) = M 4π 3 r(τ) 3   is the energy density in M − .Inserting Eq. (2.5) into Eq.(2.4) and using Eq.(2.3), we have where R s ≡ 2GM denotes Schwarzschild radius.Hence the exterior metric is nothing but Schwarzschild metric as expected, i.e., We assume that the Darmois-Israel junction conditions be still valid for the effective theory, which is widely adopted in the literature [59,65,82].In addition, to achieve the aim that the quantum effects in the interior can be carried into the exterior, one needs also to specify the interior effective dynamical equation.The effective Friedmann equation reads [46,83,84] where ρ c = 3 8πGγ 2 ∆ is the critical energy density.The parameter ∆ = 4 √ 3πγ 2 p is the area gap, where 2 p = G and γ denotes the Immirzi parameter whose value has been fixed as 0.2375 by calculating BH entropy [85,86].Note that the resulting effective Friedmann equation (2.8) is suitable for arbitrary matter sources though it is derived with a massless scalar field.The reason behind this is that the quantum corrections to Friedmann equation arise completely from the quantum modification in the gravitational part.We now consider the dynamics of the interior spacetime determined by the effective equation (2.8).Therefore, repeating the matching procedure, the quantum-corrected exterior metric can be obtained as [69] ds 2 where It is easy to see that the effective metric (2.9) tends to the classical metric (2.7) as ∆ → 0. This behavior also holds for large scale with r >> √ ∆.Hence the effective exterior metric (2.9) goes to the classical Schwarzschild metric (2.7).It is remarkable that the effective exterior metric with k = 0 in Eq. (2.9) is just the one obtained directly by loop quantizing the spherically symmetric exterior spacetime in [87].In contrast to the classical case where the exterior metric always takes the Schwarzschild one for different k, the quantum-corrected exterior metrics for different k are actually different from each other in the effective theory.
The energy density ρ + in M + measured by the static observer can be obtained as where G tt and g tt denote the tt-components corresponding to the Einstein tensor G ab and the metric tensor g ab , respectively.
Here the Hubble parameter H should be understood as a function of ρ = M 4π 3 r 3 , and thus as a function of r, by using the Friedmann equation.In the classical case, the vanishing energy density determined by Eq. (2.11) indicates the the vaccum exterior.A straightforward calculation shows that the effective energy density reads (2.12) As collapse goes on, a quantum-corrected BH will form with the occurrence of horizon for a large M. Note that the exterior metric (2.9) is static.Hence the vector field ξ a = ∂ ∂t a is a Killing vector field.By definition, the Killing horizon K is a null hypersurface determined by ξ a ξ a | K = 0. Therefore, the Killing horizon corresponds to the surface with radius r satisfying f (r) = 0.The functions f (r) of the quantum-corrected BHs are compared to that of Schwarzschild BH in Fig. 1.Interestingly, it is easy to see from Fig. 1 that each of the quantum-corrected BHs has two Killing horizons at f (r) = 0, while there is only one horizon in the BH of classical theory.It turns out that the three horizons, the event, Killing and apparent horizons, are the same for the spacetime with metric (2.9).
In the late stage of the evolution in the effective LQC model, the collapse of the dust cloud will stop whence v(r) ≡ ṙ = 0 is satisfied, which leads to f (r) = 1 − kr 2 0 .At this point, r takes the minimum value where Thus the exterior effective metric is well defined for r r m .In the limit r → r m , the interior (equivalent) energy density ρ − 3 8πG k a 2 of the dust cloud reaches its maximum value ρ c .The functions v(r) are plotted for the Schwarzschild and quantum-corrected cases in Fig. 2. It is shown that, during the collapse with decreasing r, the collapsing velocity of the boundary surface increases to its maximum and then decreases to 0 in the effective theory, while it increases to infinity in the classical case.After that, the collapsing phase will be bounced to an expanding phase in the effective theory, resulting in a transformation from a BH to something similar to a white hole.Therefore, the quantum gravitational effects resolve the classical BH singularity and replace it by a quantum bounce at small scale, and the effective metric agrees quite well with the Schwarzschild metric at large scale.The global causal structure of the maximally extended BH spacetime with k = 0 in the effective theory has been studied in [69].

III. THE QUANTUM GRAVITATIONAL EFFECTS
In this section, we will study the quantum gravity effects on observables such as the BH shadow, the QNMs and the ringdown waveform, and compare them to those of Schwarzschild case.

A. Shadows, rings and lensing rings
Considering the spherically symmetric BHs, the trajectory of a light outside a BH always lies on a plane.Taking into ac-count the conserved energy E and angular momentum J, the orbital equation of light ray approaching to the BHs described by Eq. (2.9) can be obtained as an ordinary differential equation of the radius r in terms of the azimuthal angle φ on the orbital plane as where b = J/E is the impact parameter associated to the light ray.Hence the trajectory of light ray is completely determined by its impact parameter b.The radius r ph of the photon sphere formed by a bounded orbit of light is determined by , the light ray will be binded by the BH and then escapes to spatial infinity.For the Schwarzschild BH, the orbital radius of photon sphere and the corresponding impact parameter read respectively as For the quantum-corrected BHs, these physical quantities depend on ∆ as well as γ.It should be noted that, in Fig. 1 and the following numerical calculations, the parameter R s is fixed as R s = 2, and hence the impact parameter b is scaled by GM.
The parameters r ph and b c for different values of parameters γ and ∆ are plotted in Fig. 3 for the quantum-corrected BHs and compared to those of the Schwarzschild BH.In the left and middle panels of Fig. 3, the Immirzi parameter is chosen as γ = 0.2375 and γ = 1, respectively.It is clear that as ∆ increases, both r ph and b c monotonically decrease for each of the quantum-corrected BHs.The right panel of Fig. 3 shows that both r ph and b c have a similar behavior with fixed ∆ = 0.25 and varying γ.It indicates that for small ∆, both r ph and b c for each of the quantum-corrected BHs are always smaller than those of Schwarzschild BH with the same mass.It should be noted that by choosing ∆ = 0.25, its effect on the shadows of BHs are obviously different for different theories, although ∆ should be a very small number.Also, we will fix γ = 1 for simplicity in the following.
To study the light bending near a BH, it is convenient to introduce a variable u = 1/r.Then the orbital equation (3.1) can be expressed in terms of u as Hence the total change in azimuthal angle outside the horizon of the trajectory for the situation (ii) with b < b c is given by [74,88] where u h := 1/r h with r h being the radius of the (outermost) horizon.In the situation (iii) with b > b c , the total change in angle reads where u max := 1/r min with r min being the light ray's radial minimal distance from its trajectory to the BH.Let n = φ/(2π) be the total number of orbits, which is a function of b, satisfying [74,88] For each given m, there will be two solutions for Eq.(3.10) [74,88], denoted by b ± m with b − m and b + m being the minimum and the maximum solutions, respectively.Then the rays can be classified as follows: (1) Direct: . Note that the orbit equation of the time-like geodesic reads, The radius r isco of the innermost stable circular orbit (isco) is determined by (3.12) In Table I, the involved physical quantities are shown for the classical and quantum-corrected BHs.It indicates that for the fixed parameters R s = 2, γ = 1, ∆ = 0.25, and r0 = 0.5, the quantum effect always shrinks the corresponding quantities.It is worth noting that the universal conjecture [   is satisfied for Schwarzschild BHs as well as quantumcorrected BHs.The behaviors of photons in the effective spacetime of the quantum-corrected BH with k = 0 are plotted in Fig. 4. The left panel of Fig. 4 depicts the total number of orbits, and the right panel shows the trajectories of light rays surrounding the quantum-corrected BH.Now we study the shadows of the two kinds of BHs surrounded by an optically and geometrically thin accretion disk on the equatorial plane of BHs, with an observer located at the north pole.Let us consider the simple case, where the emission originates from the accretion disk near BHs, and the emission intensity I em ν depends only on the radial coordinate r.Here ν denotes the emission frequency in a static frame.It turns out that the observed intensity is related to the emission intensity by [74] where I em (r) := I em ν (r)dν is the integrated intensity, and r m (b) (m = 1, 2, 3, • • • ) is the so-called transfer function describing the radial position of the mth intersection of the light ray and the accretion disk outside the horizon at φ = 2m−1 2 π.Here the absorption and reflection of light by the accretion disk are neglected for simplify.For the case of m > 3, it turns out that the contributions from the corresponding photon rings to the total luminosity can be ignored.The first three transfer functions r m (b) (m = 1, 2, 3) can be expressed as ) where u(φ, b) denotes the solution to the orbit equation (3.7).The first three transfer functions for Schwarzschild spacetime and the effective spacetime are plotted in Fig. 5.According to Eq. (3.14), the local brightness can potentially reach extremely high values as a result of the superposition of contributions from all the intersections.However, as emphasized in Ref. [74], the detector should measure the average brightness, which is proportional to the flux detected by the detector.As shown in Fig. 5, the second and third transfer functions for the quantum corrected BH with k = 0 only have support within the "photon ring" regions b ∈ (4.85475, 6.12553) and b ∈ (5.07672, 5.13117) respectively.Thus, the support of the higher-order transfer function is dramatically narrower than that of the lower-order one, since the light rays would converge due to the effect of the strong gravitational field, and hence those interacting with the accretion disk more times would converge towards a smaller region.Consequently, the images of the disk resulting from the contributions of m > 1 are highly demagnified, with demagnification factors determined by the slope dr m /db.For the case of m = 3, the numerical result yields dr 3 /db ≈ 581.908, indicating that the tertiary image contributes around 1/580 of total flux.Since the images for m > 3 are even greater demagnified compared to m = 3, their contribution to the total flux is negligible.Moreover, as shown in Ref. [74], the negligible contribution of the photon ring to the total flux can be analytically understood by examining the total azimuth angle φ given in Eqs.(3.8) and (3.9), with the corresponding numerical results displayed in Fig. 4. By applying the method of matched asymptotic expansion (see Ref.This logarithmic divergence is consistent with the divergent behavior in the Schwarzschild spacetime.Therefore, following the analysis presented in Ref. [74], it can be concluded that the contributions of photon rings with m > 3 to the total flux decay exponentially.
To study the observational appearance of emission, one needs to specify the intensity of emission I em .Now let us consider the following three specific intensities of emission, where the emission intensities are peaked at r isco , r ph and r h , respectively, and decay sharply for the first two cases and decay slowly for the last case.Here I 0 denote the maximum value of the emitted intensities.It is easy to see from Table I that the difference of the three values of r isco , r ph and r h between the Schwarzschild and effective spacetimes are very tiny, resulting in almost the same observational appearance of emission originating near BHs.The observational appearances of the thin disk near the BHs with these three different profiles are shown in Fig. 6.For comparison, the left and middle panels in each row of Fig. 6 depict the emission intensities and observational intensities for the Schwarzschild BH and its quantum corrections.We observe that for the fixed parameters R s = 2, γ = 1, ∆ = 0.25, and r0 = 0.5, the quantum correction always shrinks the radius of shadows.
Considering the spherical symmetry of the spacetime, the perturbation field Ψ can be expressed as where Y(θ, φ) denotes the spherical harmonics.It turns out that the equation of motion for the perturbation field Ψ can be uniformly written in the following Schrödinger-like wave equation [92]: where r * is the tortoise coordinate as the solution of which maps the region (r h , ∞) into the region (−∞, +∞), and V(r * ) denotes the effective potential with the form here l is the multipole quantum number, and s denotes the spin of the perturbation field with values 0, 1 and 2 for the scalar perturbation, vector perturbation, and axial perturbation, respectively.
Assuming that ψ(t, r * ) can be split into Eq. (3.22) reduces to the time-independent wave equation The complex frequency ω can be obtained by solving Eq. (3.26) under appropriate boundary conditions for the wave at event horizon (r * = −∞), and infinity (r * = +∞).
We now calculate the QNM frequencies ω under the scalar and vector perturbations in the case of l = 2, starting from the time-independent wave equation (3.26) and the time evolution wave equation (3.22), respectively.Notice that both Eq.(3.26) and Eq.(3.22) are completely determined by the effective potential V(r * ).To obtain the precise values of V at r * , we integrate numerically Eq. (3.23) under suitable initial data, e.g., r(r * = 0) = 10R s , and inserting the result into the right-hand side of Eq. (3.24).The effective potentials of the scalar and vector perturbations outside the classical and quantum-corrected BHs are plotted in Fig. 7.As it can be seen, for both perturbations, the effective potential for the effective theory always takes the maximum at each point r * , and the effective potentials in the effective theory have the form of a barrier and take constant values at the event horizon and spatial infinity.Hence, it is convenient to impose the condition such that the wave is purely incoming at the horizon, while it is purely outgoing at spatial infinity, i.e., [94] Moreover, negative imaginary part of ω indicates that ψ is damped and thus is stable, while positive imaginary part means an instability.
To solve ω from Eq. (3.26), we adopt the WKB method, by treating it as the problem of scattering near the peak of the barrier potential in quantum mechanics.This method is initially utilized in [94], and is improved in [117,118].The resulting values for the fundamental QNM frequencies up to the 6thorder approximation are shown in Table II as well as in Fig. 8. Also, one can employ the time domain integration method to calculate the evolution of ψ(t, r * ) at a fixed point r * and obtain the time-domain profile, from which the frequency ω can be extracted by using the Prony method (See, e.g., [119,120] for reference).To do this, we implement the time domain integration by employing the finite difference method.Introducing the light-like coordinates then Eq. (3.22) can be expressed in terms of u and w as To study numerically the differential equation (3.29), one can discretize it on the u-w null grid as [77] ψ i, j =ψ i, j−1 + O(δ 4 ), (3.30)where δ denotes the overall grid scale factor, ψ i, j := ψ(iδ, jδ) and V i, j := V ( j−i)δ 2 .To implement the discretized evolution (3.30), one needs to specify certain initial data.It is shown that the QNMs depend neither on the initial data nor on the small δ.Following Ref. [77], on the null boundary u = 0, ψ is specified as a Gaussian wavepacket where w 0 and σ are the median and width of the wavepacket, and on the null boundary w = 0, ψ is specified as a constant determined by ψ 0,0 .Then one can calculate the values of ψ i, j on the whole grid from Eq. (3.30).In this way, one can generate the ringdown waveform, namely, the time domain profile of perturbations by extracting the values of ψ at constant r * .
To determine the quasinormal frequency corresponding to the profile, we use the Prony method to fit the profile data, and set w 0 = 0, σ = 1, δ = 0.5.The ringdown waveforms |ψ| of the classical and quantum-corrected BHs under the scalar and vector perturbations with l = 2 are plotted in Fig. 9.It shows that the quantum-corrected BHs and the Schwarzschild  BH have almost identical waveforms.All the waveforms have the same power-law tail due to the fact that the effective potentials of these BHs have the same asymptotic behavior.We observe that the oscillation frequency for each of the quantumcorrected BHs is higher than that of the Schwarzschild one in each of the perturbations.Moreover, the damping rate for the Schwarzschild BH is the higher in each of the perturbations.These results are in agreement with those provided in Table II.To further check the consistency between the generated ringdown waveforms in Fig. 9 and the QNMs obtained by the WKB method and presented in Table II, the frequencies ω are also extracted by the Prony method from the ringdown waveforms in Fig. 9 and shown in Table III.By comparison, we find that the data presented in Table II agree quite well with those shown in Table III.The fact that the imaginary frequencies always take negative values indicates that these oscillations will die off with time evolution, and thus the classical and quantum-corrected BHs are all stable against the scalar and vector perturbations.To understand how quantum corrections affect the observables of BHs, the shadows and QNMs of quantum-corrected BHs have also been studied in the previous sections.On one Re(ω),Im(ω)  III.Moreover, we found that the quantum correction increases the real part of the QNMs and decreases the absolute value of the imaginary part.In particularly, the imaginary frequencies of QNMs always take negative values, and thus the quantumcorrected BHs are stable against the scalar and vector perturbations.
) d V(r) dr r=r ph = 0. (3.3) Due to the spherical symmetry, Eq. (3.3) can be reduced to d dr f (r) r 2 r=r ph = 0. (3.4) Inserting the solution r = r ph of Eq. (3.4) into Eq.(3.2), one obtains the critical compact parameter b c corresponding to the photon sphere as b c = r ph f (r ph ) .(3.5)According to the value of b, the light trajectory near a BH can be classified into the following three situations: (i) b = b c , the light ray will surround BHs in the circular orbit; (ii) b < b c , the light ray will approach and fall into the BH; (iii) b > b c

FIG. 3 .
FIG. 3. Behavior of r ph and b c for Schwarzschild BH (black) and quantum-corrected BHs corresponding to k = −1 (red), k = 0 (blue) and k = +1 (green): In the left and middle panels, we fix γ = 0.2375 and γ = 1, respectively, by varying ∆.In the right panel, we fix ∆ = 0.25 and vary γ.The parameters R s and r0 are taken to R s = 2 and r0 = 0.5.

FIG. 4 .
FIG.4.Behavior of photons in the effective spacetime with k = 0 as a function of b: On the left panel, the total number of orbits, n = φ/(2π), is shown.The black, gold and red lines correspond to the direct, lensed and photon rings, respectively.On the right panel, a selection of associated trajectories in the Euclidean polar coordinates (r, φ) is plotted.The impact parameter spacing is 0.1, 0.01 and 0.001 in the direct (black), lensed (gold) and photon rings (red), respectively.The BH is represented by a black disk, while the circular orbit of light is shown as a dashed black circular.The parameters are R s = 2, γ = 1, and ∆ = 0.25.
[78] and references therein), it can be analytically calculated that φ(b) diverges logarithmically as b approaches b c , i.e., φ ∼ −C log |b − b c |, where C is a positive constant depending on the sign of b − b c .

FIG. 5 .
FIG. 5.The first three transfer functions for a face-on thin disk in Schwarzschild spacetime (left) and the effective spacetime with k = 0 (right), representing the radial coordinate of the first (black), second (gold), and third (red) intersections with a face-on thin disk outside BH: The parameters are R s = 2, γ = 1, and ∆ = 0.25.

55 FIG. 6 .
FIG. 6.The observational appearances of the thin disk near the BHs with the three different profiles: In each row, the first two panels show the emission intensity I em /I 0 and observational intensity I obs /I 0 , normalized to the maximum value I 0 , of a thin disk near the quantum-corrected BHs, corresponding to k = −1 (red), k = 0 (blue) and k = +1 (green), compared to those of the Schwarzschild BH (black), and the third panel depicts the density plot of I obs /I 0 of a thin disk near the quantum-corrected BH with k = 0.The parameters are R s = 2, γ = 1, ∆ = 0.25, and r0 = 0.5.

FIG. 8 .FIG. 9 .
FIG.8.The real (top) and imaginary (bottom) parts of the frequencies ω as a function of the order of the WKB method up to 6th order for the scalar (left panel) and vector (right panel) perturbations of l = 2, corresponding to the Schwarzschild BH (black) and the quantum-corrected BHs corresponding to k = −1 (red), k = 0 (blue) and k = +1 (green).The parameters are R s = 2, γ = 1, ∆ = 0.25, and r0 = 0.5.

TABLE III .
The QNMs ω of the scalar and vector perturbations of