Quasinormal modes of a regular black hole with sub-Planckian curvature

This paper explores the properties of the quasinormal modes (QNMs) of a regular black hole(BH) characterized by a Minkowski core and sub-Planckian curvature. When focusing on a special case, this regular BH exhibits identical large-scale behavior with the Hayward BH and some loop quantum gravity corrected (LQG-corrected) BH. A notable characteristic of the QNMs in this regular BH is the pronounced outburst of overtones when compared to the Schwarzschild BH (SS-BH). This outburst can be attributed to the deviation from the SS-BH in the near-horizon geometry region due to the quantum gravity effect. Furthermore, we compare the QNM properties of the regular BH with those of the Hayward BH and the LQG-corrected BH. A similar phenomenon of overtone outburst is observed in the modes of the overtone. As a conclusion, the QNMs may be a powerful tool for detecting the quantum gravity effect and distinguishing different BH models.


I. INTRODUCTION
General relativity (GR) is plagued by the well-known singularity issue, wherein the scalar curvature becomes infinitely large at the core of the black hole (BH) [1][2][3][4][5].Singularities signify GR's failure.Therefore, it makes sense to anticipate that a quantum theory of gravity would eventually take hold in these regions and solve the singularity issue.Prior to the establishment of a full theory of quantum gravity, physicists intensively investigated various non-singular BHs at the phenomenological level, commonly known as regular BHs.
These BHs can be constructed by incorporating quantum gravity effects or exotic matter that typically violates the standard energy conditions.Typically, these regular BHs usually do not fulfill the vacuum Einstein equations.For comprehensive reviews on regular BHs, we can refer to [6][7][8].
In terms of the asymptotic behavior near the center of the regular BHs, we can classify them into two types: one featuring a de-Sitter (dS) core, and the other featuring a Minkowski core.The regular BH with dS core includes the well-known Bardeen BH [9], Hayward BH [10] 1 and Frolov BH [12].The key characteristic of these regular BHs is that their Kretschmann scalar curvature remains finite everywhere.
The BH with Minkowski core was initially postulated in [13], featuring a Newton potential that is exponentially suppressed.Subsequently, this solution has been extended to a class of regular BHs with various forms of exponential potentials in [14][15][16][17][18][19][20][21].Additionally, as the Hawking temperature vanishes during the final phase of evaporation for this type of regular BH, it provides a more realistic depiction of the final phase of BH evaporation and a more reasonable behavior of its remnant.However, it is important to note that these regular BHs featuring a Minkowski core are meaningful only for BHs with relatively small mass at the Planck scale or during the last phase of BH evaporation.Furthermore, based on the generalized uncertainty principle over curved spacetime, a novel regular BH with Minkowski core is proposed in [22].This novel regular BH overcomes the drawback that the Kretschmann scalar curvature is unbounded by the Planck mass density.Consequently, the metric is applicable to BHs with arbitrarily large masses, making it suitable throughout the process of BH evaporation.Especially, these BHs can reproduce the geometry of Bardeen BH, Hayward BH, or Frolov BH at large scales, depending on the chosen potential form.Furthermore, the authors in [23][24][25] explore the features of the photon sphere and the extremal stable circular orbit (ESCO) for massive particles, the accretion disk, and the shadow over this novel regular BH.These analyses contribute to the differentiation of this novel BH from others through astronomical observation.In this paper, we aim to study the properties of the quasinormal modes (QNMs) over this novel regular BH.
QNMs are the oscillatory patterns that emerge when a BH is perturbed, offering insights into the stability and spacetime properties of the BH [26].The frequencies and damping times associated with these oscillations serve as key characteristics of QNMs, providing essential information about the internal structure and properties of the BH.Departures from GR are anticipated to leave a discernible mark on the QNM spectra, serving as a distinct means of investigating quantum gravity effects [27][28][29][30][31][32][33][34][35][36].
Recently, the analysis of gravitational wave (GW) data from LIGO/Virgo has provided evidence supporting the existence of overtone patterns [37,38].Theoretically, overtones offer more precise information about the mass and spin of the remnant BH compared to results derived solely from the fundamental mode [39].Furthermore, a number of researches have also highlighted the significance of overtones in effectively modeling the ringdown phase [39][40][41][42].These findings challenge the dominant notion that the fundamental mode predominantly dictates the signal and suggest that the quasinormal ringing may start earlier than previously expected.
Moreover, studies have discovered that the characteristics of high overtones are linked to the geometric properties around the event horizon of a BH [43].It is widely recognized that modifications to Einstein's theory will definitely result in significant deformations of the BH's geometry near the event horizon.While the fundamental mode remains relatively unaffected by alterations in the event horizon's geometry, even minor changes in this region can exert a substantial influence on the first few overtones, as observed in [43].This highlights the potential to explore the geometry of the BH's event horizon through the investigation of overtones.Such findings motivate further exploration of overtones in various BH scenarios [44][45][46][47].
In this paper, our objective is to explore the properties of the QNMs, including both the fundamental mode and the high overtones, for regular BHs with a Minkowski core and sub-Planck curvature.Then, we will compare the QNM properties with those of regular BHs featuring a dS core, as well as loop quantum gravity corrected (LQG-corrected) BH.The paper is structured as follows.Section II provides a concise overview of the regular BH with a Minkowski core and sub-Planck curvature, along with the perturbations of the scalar field over this regular BH.In particular, we delve into the properties of the effective potential, which is intricately connected to the characteristics of the QNMs.In Section III, we study the properties of quasinormal frequencies (QNFs), including both the fundamental mode and the overtones.Additionally, a comparison is made among the QNFs of the regular BH we studied, Hayward BH, and LQG-corrected BH in Section IV.We summarize our findings in Section V.The numerical approaches and error analysis of QNMs are briefly elaborated in Appendix A.

II. SCALAR FIELD OVER THE REGULAR BLACK HOLE
The metric of the regular BH with Minkowski core and sub-Planckian Kretschmann scalar curvature has the following form [22] Here, the function f (r) is given by: where P (r) represents the modified Newton potential defined as: In this equation, M represents the mass of the BH, while α 0 , x, and c are dimensionless parameters.Specifically, α 0 represents the degree of deviation from the Newton potential.
We would like to emphasize that the aforementioned modified Newton potential exhibits an exponentially suppressing form, playing a crucial role in the construction of a regular black hole with an asymptotically Minkowski core.When α 0 = 0, the modified Newton potential given by Eq.( 3) becomes equivalent to the standard Newton potential.As a consequence, the regular BH is transformed into the Schwarzschild BH (SS-BH).While the parameters x and c characterize the unique features of the regular BH with sub-Planckian curvature.To guarantee that the existence of a BH horizon and the scalar curvature to be sub-Planckian, both the parameters x and c must satisfy the conditions c ≥ x ≥ c/3 and c ≥ 2. Without loss of generality, we set x = 1 and c = 3 throughout the paper.
Subsequently, we can compute the Hawking temperature of the regular BH, which is expressed as follows: where r h is the event horizon of the BH.In the extremal case, i.e., when T = 0, the deviation parameter α 0 should satisfy the following condition: Note that the e in the above equation is just the natural constant, also known as the Euler's number.To illustrate the above relationship, we present the Hawking temperature as a function of the deviation parameter α 0 in Fig. 1 for visualization.We proceed to study how this regular BH responds to the perturbations of a probing massless scalar field Φ.We can describe the dynamics of the scalar field using the Klein-Gordon (KG) equation: Due to the spherical symmetry of spacetime we are investigating, we may utilize the technique of separating variables using spherical harmonics: Here, Y l,m (θ, ϕ) is the spherical harmonics, with l and m representing the angular and azimuthal quantum numbers, respectively.This transformation reformulates the Klein-Gordon equation into a unified Schrödinger-like form: where r * is the tortoise coordinate related to r as dr * /dr = 1/f (r).For a given l and m, we have simplified the notation by removing the subscript from Ψ l,m (t, r) and denoted it as Ψ.
In addition, the effective potential is expressed as follows in the equation above: The angular quantum number l takes the values of 0, 1, . ... for the scalar field, showcasing various deviation parameters α 0 .We can find the effective potential is always positive, indicating that the system is stable under the scalar field perturbation.In the case of l = 0, the first term disappears, leaving only the second term.In this scenario, an increase in the deviation parameter α 0 leads to a decrease in the peak of the potential, as evident in the inset shown in Fig. 2. When l > 0, the first term of the effective potential becomes non-zero.In this scenario, it is observed that as the deviation parameter α 0 increases, the peak of the potential rises (the right panel in Fig. 2), which demonstrates a contrasting trend compared to the l = 0 case.

III. QUASINORMAL MODES
In this section, we will thoroughly investigate the properties of the QNMs associated with this regular BH.To determine the QNM spectra, it is imperative to apply the following boundary conditions:  Here, the minus sign and the plus sign correspond to purely ingoing waves near the event horizon and purely outgoing waves at the boundary, respectively.Consequently, there are no waves emanating from either of the boundaries.This implies that a response in perturbations is only detected when the source of perturbations ceases its action.In other words, proper oscillation frequencies are observed once the perturbation source is no longer active.
Several methods have been developed for determining QNMs, including the Wentzel-Kramers-Brillouin(WKB) method [48][49][50][51][52][53], the Horowitz-Hubeny method [54], the continued fraction method [55], the asymptotic iteration method [56][57][58], the pseudo-spectral (PS) method [59,60], and others.Among these, the WKB method stands out as an extremely efficient approach.The WKB method has undergone significant development from 1st to 13th order to determine the QNFs [49][50][51][52][53]61].It is important to highlight that, when employing WKB methods, a higher order does not automatically ensure higher accuracy in the results.Through meticulous error analysis, we have observed that the results obtained using the ninth-order WKB approximation in this study exhibit a minimum error (refer to the Appendix A).Consequently, for our investigation, we will utilize the ninth-order WKB approximation to determine the QNFs.
However, the effectiveness of WKB method loses when seeking higher overtones, particularly in the case of n > l.Table I displays the QNFs for the fundamental modes with various deviation parameter α 0 at l = 0 and l = 1, using WKB and PS methods.It is noteworthy that the difference between the WKB and PS methods is consistently less than a thousandth.This trend holds true even for l > n, as demonstrated in Tables II, III and IV.
But with an increase in the value of n, this difference becomes more evident.This difference can be evaluated by the relative error, defined as: Here, we assume that the results produced from the PS approach are precise.Clearly, in the case of l = 0, the relative error has reached two percent even for the fundamental mode (refer to Table II).However, for l = 1, the relative error is one percent for n = 3.Meanwhile, for l = 2, we can calculate the QNF with n = 5 with a relative error of one percent.
Therefore, in this paper, we primarily utilize the WKB method as a cross-check, specifically for determining fundamental modes.For higher overtones, we commonly resort to the PS method-a fully numerical approach renowned for its power in determining QNFs [29,30,60,[62][63][64][65][66][67][68][69].Nevertheless, it is essential to note that, like many numerical methods, the PS method may be time-consuming compared to the semi-analytic WKB method.
Therefore, when employing the PS method, a thoughtful balance between computational time and desired accuracy is crucial.For a concise introduction to the WKB method and the PS method, please refer to the Appendix A.    red dots represent various deviation parameters α 0 .As the deviation parameter α 0 grows, we see that the real part ω R shows a distinct non-monotonic pattern (left plot in Fig. 3).
More precisely, when the parameter α 0 grows, indicating a larger departure from the SS-BH, ω R first rises, indicating a more pronounced oscillation of the system, and subsequently declines, suggesting a reduction in the amplitude of the oscillation.The turning point of this non-monotonic behavior occurs at approximately α 0 ≈ 13/20.It is important to mention that when the deviation parameter α 0 is increased to a value near to the extremal case, the magnitude of ω R becomes less than that of the SS-BH.Simultaneously, when α 0 increases, the magnitude of the imaginary part of the QNFs, |ω I | − ω R , consistently drops, suggesting a reduced damping rate.This finding suggests that the presence of quantum gravity effects leads to decelerated decay modes.
However, it is clearly observed in the right plot in Fig. 3 that the non-monotonic behavior disappears for the case of l = 1.More precisely, when the parameter α 0 grows, both ω R and ω I increase in a monotonic manner.This indicates that the system exhibits stronger oscillations and slower decay.In addition, we also analyze the fundamental mode for values of l > 1 and reach a comparable conclusion to that of the situation where l = 1.A similar observation has been noted in the fundamental modes associated with the LQG-corrected BH [31].This suggests that for the fundamental modes, the influence of the angular number surpasses that of quantum gravity corrections.In the future, it deserves further investigation into the university of the non-monotonic behavior for l = 0 and exploration of the underlying reasons behind this phenomenon.As revealed above, the fundamental modes exhibit interesting properties in the context of this regular BH pattern.However, it is essential to note that upon closer analysis, we observe that the absolute values of the fundamental modes change only by a few percent compared to those of the SS-BH as the deviation parameter α 0 increases (refer to Table I).In Fig. 4, the blue curves corresponding to n = 0 clearly highlight this discrepancy.However, while examining the overtone modes, a notable discrepancy arises in QNFs between the regular BH and the SS-BH.
The most crucial feature of the overtone modes is the phenomenon known as the "outburst of the overtone", indicating a relatively significant change in the spectra of the overtone modes compared to that of the SS-BH [30, 31, 43-45, 47, 70, 71].In our current model, such an outburst can also be distinctly observed for l = 1 beginning from the third overtone (refer to Fig. 4).Furthermore, from Fig. 4, we also note that as the deviation parameter α 0 increases, this regular BH approaches the extremal case and an oscillatory behavior can be clearly observed.This oscillation behavior becomes more pronounced for the higher overtones.Especially, we would like to point out that such oscillation behavior has also been observed in LQG-corrected BH [30,31].This oscillatory pattern may be associated with the extremal effect.But the reason behind this phenomenon deserves further exploring in the future.
Recall that the non-monotonic behavior of QNFs appears in the case of l = 0 for the fundamental modes, but this non-monotonic pattern disappears for l ≥ 1. Surprisingly, for the overtones, we observe the reemergence of this non-monotonic behavior (refer to Fig. 5).
Specially, for l = 1 and n = 2, as α 0 increases, the magnitude of ω R becomes smaller than that of the SS-BH, similar to the behavior observed in the fundamental modes with l = 0.
However, we observe that the critical point of α 0 at which the magnitude of ω R becomes smaller than that of the SS-BH is smaller than that of the fundamental modes with l = 0.  Upon further examination of the phase diagram |ω I | − ω R for l = 1 and higher overtone modes, a spiral-like structure is identified in Fig. 6.As the overtone number increases, this spiral-like structure becomes more pronounced.Notably, for l = 0, a clearly visible spirallike structure is observed for the first overtone mode (see Fig. 7).This spiral-like pattern is, in fact, a manifestation of the oscillatory behaviors observed in Fig. 4.Such a spiral-like structure has been observed in the Hayward BH [44] and the LQG-corrected BHs [30,31,72].
Previous studies have suggested that the outburst in overtones may be associated with subtle differences between the SS-BH and the modified effective SS-BH [30, 31, 43-45, 47, 70, 71].This assertion is further supported in our current model.Fig. 8 illustrates a comparison of the effective potential between this regular BH and the SS-BH for different values of α 0 , considering l = 0 and l = 1.This comparison is measured by the quantity ∆V , which quantifies the difference in the effective potentials (see the inset in Fig. 8).Evidently, there is a minor change in the region surrounding the event horizon.Increasing the deviation parameter α 0 causes a more significant change in the effective potential.This change in the effective potential may act as a trigger for the outburst of overtones.Therefore, the deviation parameter α 0 significantly influences the outburst of overtones.

IV. A COMPARISON AMONG THE QNFS OF THE REGULAR BH, HAYWARD BH, AND LQG-CORRECTED BH
Throughout this paper, we have employed the values x = 1 and c = 3.In this setting, the redshift factor (1) of the regular BH can be expanded at large scales, i.e., as r → ∞, If we also expand the redshift function of the Hayward BH at large scales, we observe: Remarkably, we find that the large-scale behavior of the Hayward BH is identical to that of the regular BH with x = 1 and c = 3.This similarity has been pointed out in [22].Interestingly, a LQG-corrected BH within the entire spacetime has the following form [73]: which is just the large-scale behavior of the Hayward BH and the regular BH with x = 1 and c = 3.
Next, we will compare the properties of the QNFs of the regular BH studied in this paper with those of a Hayward BH and a LQG-corrected BH.This comparative analysis seeks to evaluate the viability of QNFs as a reliable probe for discerning among these distinct models.
To begin with, we illustrate the difference in the effective potentials between the Hayward BH and the regular BH in Fig. 9, and likewise, the difference between the effective potentials of the LQG-corrected BH and the regular BH in Fig. 10.For a small deviation parameter α 0 , it appears that all three models have identical effective potentials, as seen in the left plots of Figs.9 and 10.However, when doing meticulous investigation, we observe nuanced distinctions among them.Particularly, the difference is obvious in the vicinity of the event horizon (refer to the insets in the left plots in Figs.9 and 10).The disparity becomes more evident with a larger deviation parameter α 0 , as see the right plots in Figs.9 and 10 for It is expected that subtle differences in the vicinity of the event horizon can be distinguished by the QNFs, as discussed earlier.To investigate this, we calculate the QNFs for these three black hole models with the same parameter α 0 .Furthermore, we quantify the difference between the regular BH and either the Hayward BH or the LQG-corrected BH using the following definition: where ω Regular represents the QNF of the regular BH, and ω # denotes the QNF of either the  From Table V, we observe that for α 0 = 10/52, the discrepancy in the fundamental modes among the three BH models is merely 0.01 percent.However, as the overtone number increases, the disparity grows, reaching approximately 0.15 percent for n = 5.Increasing the deviation parameter to α 0 = 30/52, we note that the difference among these models reaches approximately 0.15 percent even for the fundamental modes, particularly escalating to around 4.3 percent for some higher overtones.These trends can be characterized as the outburst of overtones in the regular BH when compared to the Hayward BH and the LQGcorrected BH.In conclusion, the QNFs undeniably serve as a reliable tool for distinguishing these distinct models.

V. CONCLUSION AND DISCUSSION
In this paper, we investigate the properties of the QNMs of a regular BH characterized by a Minkowski core and sub-Planckian curvature.Our specific focus is on the instance of this regular BH with fixed parameters x = 1 and c = 3, where the only free parameter is the deviation parameter α 0 .This regular BH shares identical large-scale behavior with the Hayward BH, and notably, the large-scale behavior of both these BHs corresponds to the geometry of the LQG-corrected BH in [73].
Owing to the influence of quantum gravity, the fundamental modes display a nonmonotonic trend concerning the deviation parameter α 0 when l = 0.However, this nonmonotonic pattern diminishes for l ≥ 1, indicating that the impact of the angular quantum number supersedes that of the quantum gravity effect.This non-monotonic behavior reemerges in the modes of higher overtones.
Furthermore, we investigate the properties of the overtones.The most crucial feature of the overtone modes is the outburst of the overtones compared to that of SS-BH.results.In this subsection, we present a succinct overview of the WKB method.We will also perform an error analysis for various orders of the WKB method when determining the QNFs for our current model, aiming to identify the most suitable order for handling our specific model.
Typically, the higher-order WKB formula can be formulated as follows: In this context, V 0 denotes the maximum of the effective potential, V 2 stands for the secondorder derivative of the effective potential with respect to r at its maximum position, and K assumes a value that is a half-integer.A k (K 2 ) is the correction term of order k, comprising polynomials in K 2 with coefficients as rational numbers.This correction term depends on the higher derivative of the effective potential V eff at its maximum position.For more details, please refer to [61,74,75].
To assess the most suitable order of the WKB method for our model, we adopt the error estimation function ∆ k in the work [61]: The left plot in Fig. 11 illustrates the error estimation function ∆ k as a function of the WKB order for the fundamental mode, considering different quantum-corrected parameters α 0 .Our analysis shows that the ninth-order WKB approximation produces a comparatively tiny ∆ k in the model.In addition, we illustrate the relationship between ∆ k and α 0 for the ninth-order WKB approximation, clearly showing that it is consistently less than 0.105%.
These findings lead us to the conclusion that the ninth-order WKB approximation provides a superior degree of accuracy when compared to WKB approximations of other orders in our model.Therefore, we will follow the ninth-order WKB technique to calculate the QNFs in this case.

Pseudo-spectral method
The PS method, which is a fully numerical technique, is highly effective in determining QNFs.It has been widely utilized in various studies [29,30,60,[62][63][64][65][66][67][68][69], particularly for identifying high overtone modes [30].The fundamental principle of the PS approach illustrating its dependence on the WKB order for various deviation parameters α 0 , is presented in the left plot.Simultaneously, the right plot illustrates its dependence on α 0 for the ninth-order WKB approximation.
is to discretize the differential equations and subsequently solve the resultant generalized eigenvalue equations.Specifically, we substitute the continuous variables with a discrete collection of collocation points known as the grid points, and we expand the functions using specific basis functions referred to as cardinal functions.Typically, the Chebyshev grids and Lagrange cardinal functions are employed: x − x i x j − x i , i = 0 , ... , N .

(A3)
Another crucial aspect of the PS approach is to operate in the Eddington-Finkelstein coordinates.In this coordinate system, Eq.( 8) linearizes with respect to the frequency ω.
In particular, working in the Eddington-Finkelstein coordinates offers enhanced simplicity in enforcing the boundary conditions.In order to accomplish this objective, we implement the below transformations: r → r h u and Ψ = e iωr * (u) ψ .(A4) Then, the wave equation ( 8) transforms into the subsequent expression: It is evident that Eq.(A5) satisfies the ingoing boundary condition at the horizon.Moreover, we implement the subsequent transformation to ensure that the Eq.(A5) adheres to the

FIG. 1 :
FIG. 1: The Hawking temperature as a function of the deviation parameter α 0 .Here, we have set r h = 1 and M = 1.

Fig. 2
Fig.2 depicts the behavior of the effective potential V eff as a function of the coordinate r

FIG. 2 :
FIG.2:The effective potential V eff (r) of the scalar field with different deviation parameters α 0 .The left panel is for l = 0, while the right panel is for l = 1.The black dotted line represents the maximum of the effective potential V eff (r).

FIG. 3 :
FIG.3:The fundamental modes of the scalar field perturbation over the regular BH are depicted in the left plot for l = 0 and the right plot for l = 1.The black dot corresponds to the SS-BH (α 0 = 0), while the red dots correspond to various deviation parameters α 0 .

FIG. 4 :
FIG.4:The real part ω R is depicted on the left plot, while the right plot illustrates the imaginary part ω I .Both plots represent the fundamental mode and the first five overtones from top to bottom, with α 0 as the varying parameter.

FIG. 8 :
FIG.8:The difference in the effective potentials between the regular BH and the SS-BH for various deviation parameters α 0 and angular quantum numbers l.

FIG. 9 :FIG. 10 :
FIG.9:The difference in the effective potentials between the regular BH and the Hayward BH for various deviation parameters α 0 .Here we have fixed the angular quantum number l = 1.
This effect has been noticed in several geometric situations other than the SS-BH, such as the Reissner-Nordström BH, Bardeen BH, higher-derivative gravity model, and certain quantum-corrected BHs.The outburst of the overtones is expected to detect the effects deviating from the SS-BH, particularly the impact of quantum gravity.In addition, we also analyze and contrast the characteristics of the QNFs of the regular BH with those of the Hayward BH and the LQG-corrected BH.We observe that the disparity in the QNFs between the regular BH and the Hayward BH or the LQG-corrected BH is roughly 4.3 percent for certain high overtones.Such a relatively large change of the overtone modes in the regular BH, as compared to the Hayward BH and the LQG-corrected BH, may also be referred to as the outburst of the overtones.Hence, the QNFs unquestionably function as a dependable instrument for discerning various BH types.The discovery of GWs has provided a novel means of detecting the effects of quantum gravity and distinguishing between different theories of BHs.The next crucial step is the investigation of gravitational perturbations and the analysis of their QNMs.By employing QNMs, encompassing both the fundamental mode and high overtone modes, one can achieve an accurate alignment with gravitational waves emitted during the ringdown phase.As a result, by analyzing QNMs and GW signals during the ringdown phase, we may evaluate quantum gravity effects and differentiate between various BH models.The investigation of these theories and experimental findings has the potential to enhance our comprehension and study of quantum gravity.It advances our understanding of BH physics and GW astronomy while also uncovering fundamental natural laws.

FIG. 11 :
FIG.11:The error evaluation for the fundamental modes.The error estimation function ∆ k ,

TABLE I :
The fundamental modes of the regular BH are presented for various α 0 with l = 0 and l = 1.Left column is for the WKB method, whereas the right column is for the PS method.

TABLE III
: A comparison of the QNFs between the results obtained through the WKB method and the PS method.Here we have fixed α 0 = 1 2 for l = 1 with various n.

TABLE IV :
A comparison of the QNFs between the results obtained through the WKB method and the PS method.Here we have fixed α 0 = 1 2 for l = 2 with various n.
Hayward BH or the LQG-corrected BH.The results are presented in Table V for α 0 = 10/52 and Table VI for α 0 = 30/52.

TABLE VI :
A comparison of the QNFs between the regular BH and either the Hayward BH or the LQG-corrected BH.Here we have fixed α 0 =3052 for l = 1 with various n.