Parametrized black hole quasinormal ringdown. I. Decoupled equations for nonrotating black holes

Black hole solutions in general relativity are simple. The frequency spectrum of linear perturbations around these solutions (i.e., the quasinormal modes) is also simple, and therefore it is a prime target for fundamental tests of black hole spacetimes and of the underlying theory of gravity. The following technical calculations must be performed to understand the imprints of any modified gravity theory on the spectrum: 1. Identify a healthy theory; 2. Find black hole solutions within the theory; 3. Compute the equations governing linearized perturbations around the black hole spacetime; 4. Solve these equations to compute the characteristic quasinormal modes. In this work (the first of a series) we assume that the background spacetime has spherical symmetry, that the relevant physics is always close to general relativity, and that there is no coupling between the perturbation equations. Under these assumptions, we provide the general numerical solution to step 4. We provide publicly available data files such that the quasinormal modes of {\em any} spherically symmetric spacetime can be computed (in principle) to arbitrary precision once the linearized perturbation equations are known. We show that the isospectrality between the even- and odd-parity quasinormal mode spectra is fragile, and we identify the necessary conditions to preserve it. Finally, we point out that new modes can appear in the spectrum even in setups that are perturbatively close to general relativity.


I. INTRODUCTION
The coalescence of two black holes (BHs) is a fascinating process from a theoretical perspective: two regions devoid of matter interact violently to produce enormous amounts of gravitational radiation, leaving behind a stationary, vacuum solution of the field equations.In general relativity (GR), this process is particularly interesting since both the initial BH spacetimes and the final state are characterized by only two parameters: their mass and angular momentum [1,2].
The transition from the initial to the final state is described by the so-called "ringdown" phase, during which the distorted remnant sheds away the nontrivial multipolar structure through gravitational waves, relaxing to the final stationary solution.This exponential approach to equilibrium is a property of many physical systems.It can be mathematically described as a small deviation away from the final state, or as the time window in which the spacetime response to generic perturbations can be decomposed into a set of quasinormal modes (QNMs) characterized by complex frequencies [3,4].The ringdown is of special relevance in the context of BH physics: if GR is correct, then QNM frequencies must depend only on the BH mass and spin, and therefore ringdown observations can provide simple and clean tests of GR [5][6][7][8][9][10].
It is generally accepted that GR does not offer the ultimate description of gravity.Rather, discrepancies in its strong-field or large-scale descriptions will eventually become apparent in a manner consistent with some of its known conceptual problems [8,11].Thus, one expects that BH solutions will differ, if only slightly, from those of GR, and that the dynamics of BHs will also be described by slightly different equations.The correct description of gravity at all scales and energies is presently unknown.Traditionally, the strategy has been to identify modified theories of gravity which contain new features in their dynamics that could be smoking guns of new physics.Once the theory is identified, one computes its BH solutions, studies linearized perturbations around these BH solutions, and finally solves the relevant equations to compute QNM frequencies.
Recent efforts have focused instead on understanding the possible functional form of such corrections, using resummation methods.For example, effective field theory (EFT) arguments can produce the general form of the corrections in the absence of extra fundamental fields in the problem [12].Such corrections can then be used to study BH solutions and their dynamics, most notably their QNM spectra [13].Once new fundamental fieldsfor example scalars -are allowed, the dynamics become intrinsically more complicated and richer phenomenology arises from, for instance, the coupling between the scalar and the gravitational degrees of freedom [14][15][16].
The (technical) burden of computing the QNMs cannot be avoided.In the majority of the literature, arXiv:1901.01265v2[gr-qc] 8 Feb 2024 tools yielding approximate answers, such as WKB techniques [17,18], have been used.The errors involved in the WKB approximation can be rather large even in GR (see e.g.Fig. 3 of Ref. [4]), and they are poorly known in modified theories of gravity [19].Clearly, when testing GR precision physics will be involved.In addition, it quickly becomes apparent that one is repeatedly solving the same equation (or systems of equations) with the same boundary conditions and with results (the QNM frequencies) which are but small modifications to the GR predictions.This work is intended as a first step in the general program to addressing the repetitive task of finding the eigenvalues (i.e., the QNM frequencies) of these ordinary differential equations.Once the linearized equations of motion around a BH spacetime are known, our results can be used to find the corresponding QNM frequencies.In this first work, for simplicity, we assume that there are no couplings between the perturbation equations.These complications will be addressed in a forthcoming publication [20].
We also make the rather restrictive working hypothesis that the underlying equations are separable, and that the background GR solution is nonspinning.Many known modified theories of gravity indeed yield separable equations for nonspinning background spacetimes, but it is unclear whether this feature is generic.Certainly BHs in our universe are spinning, generically breaking this hypothesis.Nevertheless, our results are interesting on their own and provide a better understanding of GR itself.Our formalism is theory-agnostic and backgroundagnostic.While in principle it only applies to nonspinning backgrounds, in Sec.V C we will show that, in fact, it can be used to accommodate (in a perturbative sense) the QNM frequencies of slowly spinning Kerr BHs.Therefore our results seem to apply more broadly than these restrictions might imply.
The plan of the paper is as follows.In Sec.II we present the general perturbative framework underlying our analysis.In Sec.III we present our main results and their implications for isospectrality.In Sec.IV we discuss analytically and numerically the convergence of the perturbative expansion, and we present a simple example of a perturbation to the known GR potential that may lead to divergences.In Sec.V we verify our formalism in some interesting cases where QNM solutions are known: an EFT calculation of deviations from GR, the Reissner-Nordström (RN) solution, and scalar perturbations are slowly rotating BHs.In Sec.VI we present a preliminary estimate of the conditions under which perturbative deviations from the GR frequencies could be detectable.We conclude in Sec.VII with a discussion of the limitations of this analysis and directions for future work.In Appendix B we prove that our assumed form for the perturbative potentials is general enough for our present purposes, and in Appendix C we study the large-ℓ (eikonal) limit of our approximation.

II. FRAMEWORK
In GR, gravitational fluctuations have two degrees of freedom, which can be described by a linear combination Φ of the metric perturbations.In a nonspinning background, such modes are described by an "axial" (or odd, or Regge-Wheeler) and a "polar" (or even, or Zerilli) equation of the form [21,22], where f = 1 − r H /r, and r H is the horizon radius in standard Schwarzschild coordinates.The effective potential V ± is for odd perturbations, and for even perturbations, where λ = ℓ 2 + ℓ − 2, and ℓ is an angular harmonic index labeling the tensorial spherical harmonics used to separate the angular coordinates from the problem (the azimuthal index can be set to m = 0 for spherically symmetric backgrounds).The dynamics of test scalar or vector fields in GR is also described by similar master equations, with potentials [4] V s = ℓ(ℓ + 1) where s = 0, 1 for scalar and vector perturbations, respectively.
Given the strong observational constraints on deviations from GR [23], it might be expected that a more fundamental description of gravity will add small corrections to the GR description of astrophysical BHs.Such corrections will slightly disturb both the equilibrium configurations (the Schwarzschild geometry) and the equations governing dynamical fluctuations [24].Our guiding light is asymptotically flat spacetime and, therefore, we parametrize the modified master equations by adding power-law corrections to the effective potential: for gravitational fluctuations and for scalars and vectors.Here, α ± j , β s j are constant coefficients.Since the maximum value of the j-th term f (r)α ± j (r H /r) j is α ± j (1 + 1/j) −j /(j + 1), the smallness of the corrections translates to the criterion: Note that α ± 0 and β s 0 in Eqs. ( 6) and ( 8) are effectively mass-square terms, hence they should be positive.Similar constraints may apply to the j = 1 terms, which dominate the potential at large distances and which are therefore constrained by weak-field physics.For example, a 1/r term appears when dealing with charged scalars in a charged background (see, e.g., Eq. (4) in Ref. [25] or Eq.(4.54) in Ref. [26]).As we show in Sec.V, the leading-order coefficients (α 0 , β 0 and β 1 , in the notation of this section) are nonzero for charged (RN) or slowly rotating Kerr BHs.Although these terms are dominant in the power-law expansion -falling off even slower than the GR terms -they should still be considered small when compared to the GR potential V ±,s , due to the smallness of the parameters α ± j , β s j .In general, perturbation equations can be reduced to the form (1), but with f ̸ = 1 − r H /r. In Appendix A we prove this for scalar and vector perturbations of static, spherically symmetric spacetimes, which indeed obey a generalized wave equation [Eq.(A4)].However in Appendix B we show that, at linear order in perturbations, Eq. ( A4) is in fact equivalent to Eq. ( 1).Gravitational perturbations depend on the chosen modified gravity theory, so it is nontrivial to provide a general proof that they always reduce to Eq. (A4), but in Appendix A we provide arguments supporting this conclusion.Note that, as we show below, even perturbations of slowly rotating BHs can be accommodated within our formalism, so its regime of validity may be larger than anticipated.
The parametrization used here embodies both the change in the background spacetime and the change in the dynamical equations.It assumes no couplings to other fundamental degrees of freedom, and it also assumes no coupling between the odd and even modes; these issues will be addressed in a forthcoming study [20].
perturbations, must be computed numerically.Below, when the context is clear, we omit the superscripts from the coefficients (e j , d j ) for simplicity.
We should note that, generically, there are frequencydependent terms in the corrections to the potential.In our approach, any ω-dependent term in the perturbing potential should be set to ω 0 , thereby corresponding to a redefinition of α ± j .We have indeed verified numerically that the induced corrections are linear in the perturbations to the potential, and that ω-dependent terms can be handled by setting them equal to their GR value.
Corrections to the QNM frequencies scale linearly with corrections to the potential appearing in the master equation.In a forthcoming paper we will show that, when couplings between fluctuations are allowed, corrections in the frequency appear at quadratic order in the couplings.In some modified gravity theories both corrections are of the same order of magnitude, and therefore couplings cannot be neglected.These issues will be addressed in a separate study [20].

III. RESULTS
The most important results of our analysis are the numerical values of the complex factors e ± j and d s j , which were obtained via high-precision direct integration of the relevant equations.Figure 1 and Table I show representative results for the coefficients e ± j and d s j in ( 10) and ( 11) for odd gravitational and scalar modes, respectively.Results for all fields and the lowest multipole numbers ℓ = |s|, |s + 1|, ..., 10 are available online [30].Our calculated e ± j (and d s j ) have a 5-digit accuracy for j ≲ 30.

A. Isospectrality
A remarkable result in GR concerns the isospectrality of potentials (2) and (3).These potentials are sometimes referred to as "superpartners," since they can both be on the completeness of this basis.Re [e obtained from a "superpotential" W 0 [4,31]: When the potentials are perturbed from their GR values as in Eq. ( 6), isospectrality will in general be broken, unless the coefficients α + j are related in some way to α − j .The linearized superpartner relations now yield with δW the induced change in the superpotential W with respect to its GR value W 0 [Eq.( 13)].
Let us admit that the theory predicts a set of α − j .The relations above give us a unique possibility for α + j , namely: These relations show that isospectrality is very fragile, and that in general it will be broken.

IV. ASYMPTOTIC BEHAVIOR AND CONVERGENCE OF SERIES
Having defined the "basis vectors" e j , it is crucial to investigate their behavior, and the overall convergence of the QNM frequency expansion (10).For large j, the contribution to the potential from the j-th term can be described by a Gaussian profile depending on the tortoise coordinate2 r ⋆ , i.e.
showing that the location of the peak is proportional to − ln j (in terms of the tortoise coordinate r * ).Moreover, since the amplitude of f δV j is proportional to 1/j for large j, we can expect that e j → 0 for j → ∞.Following [32], a perturbative expansion of the QNM frequency yields the following correction: where Φ 0 is the unperturbed mode itself.
Let us assume that f δV j has support only in −r H ln j − r H σ ≤ r * ≤ −r H ln j + r H σ, with a constant σ = O(1) (the particular choice of σ does not affect our results.)Then the integral can be estimated as follows: In conclusion we find Re [e where ω 0 = ω R + iω I .Eq. (21) shows that the basis frequency e j decays as 1/j 1+2r H ω I , and oscillates as sin(2r H ω R ln j).For tensor perturbations3 with ℓ = 2, we have e j ≃ e 1.5i ln j j 0.64 .This result is also confirmed by a numerical study of e j for large j.We fitted our large-j numerical results to the following functional form: with (κ, β, γ, δ) numerical coefficients to be determined.We find β ≃ 0.66 and γ ≃ 1.5, in very good agreement with the analytical estimate of Eq. ( 21). Figure 2 compares the fit with the actual data for e j .The agreement is extremely good, and it is strong evidence in support of the asymptotic behavior (21).
Having assessed the asymptotic properties of the basis e j , we can now study in detail the convergence of the frequencies ω QNM .We assume that the effective potential of the specific theory under consideration is C ∞ , such that we can expand it for large distances as δV = ∞ j=0 aj r j , where the a j 's are constant coefficients.Moreover, following our formalism developed in Sec.II, we expand the corrections with respect to the Schwarzschild term according to Eq. ( 6), such that a mapping exists between the α ± j and the coefficients of the theory a j , i.e. α ± j = r 2 H aj r j H .To study the convergence of the QNM expansion (10), we can compute the ratio The series converges when Υ < 1. Replacing the explicit form of the expansion coefficients α ± j , and using the fact that the frequency "basis vectors" e j behave as in Eq. ( 22) -or equivalently Eq. ( 21) -we find The previous equation implies that, in our approach, the convergence of the QNM frequencies depends on the behavior of the coefficients of the (expanded) potential for the specific theory of gravity.As a nontrivial example to test the formalism, consider the case where ρ 0 is a constant, and α ± j = (−1) j+1 10 −2 (ρ 0 /r H ) kj .For this potential the convergence criterion reduces to and therefore the perturbative formalism is valid if ρ 0 < 1.

V. EXAMPLES
In this section we will provide some specific examples of theories of gravity whose gravitational perturbations can be described in terms of a master equation with the same functional form of Eq. ( 1).In order to determine the domain of validity and the accuracy of our semianalytical approach, we will compare the results obtained by solving the exact equations of the theory on a case-by-case basis with the frequencies computed through Eq. (10).

A. Effective Field Theory
An interesting scenario to test against our formalism is given by the BH solutions derived within an EFT approach in Refs.[12,13].Deviations from GR in this theory are characterized by three coupling parameters, which define the scale of the strong curvature modifications.For the sake of simplicity, we focus on a single coupling parameter, namely Λ (ϵ 2 in the notation of Ref. [13], which we adopt).In this setup, nonspinning BHs are described by the Schwarzschild geometry, and its axial perturbations lead to the following wave equation: where V − is given by Eq. ( 2), while the extra term in the effective potential reads, Note that in this case f = 1 − r H /r and dr/dr ⋆ = f .According the formalism developed in Sec.II, the only coefficient of the expansion ( 10) is then From Table I we read off r H ω = r H [ω 0 + (0.0663354 + 0.117439i)ϵ 2 (ℓ − 1)(ℓ + 1)(ℓ + 2)].For ℓ = 2 this leads to a fractional deviation from GR of ϵ 2 (2.13, −15.84) in the real and imaginary components, respectively, in excellent agreement with published results: cf.Eqs.(20) of Ref. [13].We find a similar level of agreement for higher-ℓ modes.

B. Reissner-Nordström black holes
The odd-mode gravitational perturbations of RN BHs are driven by the master equation with We will consider the charge Q -or, equivalently, the Cauchy horizon radius r − ∼ Q 2 /(2M ) -to be a small perturbation.Note, however, that r H is not to be assumed to be equal to 2M , as it would introduce an error of the same order of the leading order correction.According to the procedure in Appendix B, introducing ϕ = √ ZΦ = (1 − r − /r) 1/2 Φ, we can write the master equation in the form where f = 1 − r H /r and we ignored O(r 2 − ) terms.We can now apply our formalism to RN perturbations.The coefficients of the effective potential expansion (10) read: Thus, the QNMs for ℓ = 2 can be written as where QNM frequency, so Ω 0 = 0.3736716844180418... − i0.0889623156889357... for ℓ = 2. 4 Using the notation of Ref. [28], these results imply , in excellent agreement with reported values in the literature [28,33,34].
We can quantify the accuracy of our approach and estimate the uncertainty on the non-GR modifications of the QNM frequencies by defining where ω full is computed by numerically solving the exact master equation without any approximation, and ω QNM is given by Eq. (34).For Q/M = 0.2, for example, the QNMs frequencies derived using Eq. ( 34) agree to within 0.05% with those in Ref. [33].The errors obtained for different values of Q/M are listed in Table II.
The relative uncertainties ∆ R,I grow with Q/M up to ∼ 10% (4%) for the real (imaginary) component of the mode's frequency when Q/M = 0.5.This behavior is consistent with the assumptions made to obtain Eq. (34), in which we neglected terms that are order O((α ± j ) 2 ), the latter corresponding to ignoring O(Q 4 ) corrections to the RN perturbations.Therefore, for large values of Q, a better accuracy would require to compute ω QNM by including second-order terms in both α + j and α − j .

C. Scalars around a slowly-spinning black hole
The massless Klein-Gordon equation □Ψ = 0 around a slowly rotating Kerr BH can be written as at O(a), where f = 1 − 2M/r, a/M is the BH angular momentum, and we assumed Ψ = e −iωt Y ℓm (θ, ϕ)Φ(r)/r. 4Because ω QNM in our formula is written as a function of r H and α j , our formalism requires that one first writes down ω QNM as a function of r H and α j .We note that we should formally use ω 0 in Eq. ( 10) as 2Ω 0 /r H .
a/M ∆R ∆I 0 0% 0% 10 −4 0.0050% 0.83% 10 −3 0.049% 5.1% 10 −2 0.49% 34% We can rewrite this equation in the following form: where we used r H = 2M + O(a 2 ).Thus and we find: Comparing these results with numerical data for ℓ = m = 2 [4,30], we find the relative percentage errors (∆ R , ∆ I ) listed in Table III.Our formula is a good approximation for very small a/M ∼ 10 −4 .For a/M ∼ 10 −2 the agreement gets worse, especially in the imaginary component.We find that a/M ≪ 10 −2 is a necessary condition for our approximation to be valid.

VI. STATISTICAL ERRORS
In this section we analyze the detectability of the modifications of the QNM spectrum by space and terrestrial GW interferometers.We follow the approach described in Ref. [5], and we refer the reader to this paper and references therein for further details.Since we are interested only in an order-of-magnitude estimate of the associated observational errors, we assume that the mass of the BH is known (and therefore that the fundamental GR frequencies are known).For a more sophisticated analysis, we refer the reader to Refs.[7,35].
The gravitational waveform measured by the interferometers is a linear superposition of two polarization states of the form h = h + F + + h × F × , where F + and F × denotes the standard pattern functions (which depend on the source orientation with respect to the detector and on a "polarization angle").In the frequency domain, the two GW components are simply given by where the amplitude coefficients A +,× ℓmn and the phase coefficients ϕ +,× ℓmn are real, S ℓmn represent the (complex) spin-weighted spheroidal harmonics of spin weight 2, which depend on the polar and azimuthal angles, and b ± (f ) are the Breit-Wigner functions: For simplicity, we consider A × ℓmn = N × A + ℓmn and ϕ × ℓmn = ϕ + ℓmn + ϕ 0 , where N × and ϕ 0 are a scale factor and a phase shift, which are assumed to be known.
Moreover, consistently with the approach introduced in Sec.II, we express the modes' frequency and damping time as (small) deviations of the corresponding Schwarzschild values: ω ℓmn = ω 0 ℓmn + δω, where δω = δω 1 + iδω 2 .The GW signal is then completely specified by 4 parameters, (A + ℓmn , ϕ + ℓmn , δω 1 , δω 2 ).In order to estimate the errors on the coefficients δω 1,2 , we consider a Fisher matrix approach, which is a reliable approximation of a full Bayesian analysis for signals with high signal-to-noise ratios [36].For this reason we focus on possible detections by the space interferometer LISA [37], an ideal tool to exploit the full potential of QNM spectroscopy [38][39][40].Finally, we set N × = 1, and ϕ + ℓmn = ϕ 0 = 0. We compute the errors on the fundamental ℓ = 2 mode of nonspinning BHs, averaging over sky orientation through the following identities: Note that the amplitude of the signal depends in general on the mass, the source distance, and the energy released during the ringdown [5].For simplicity, here we fix A + ℓmn by requiring that for δω 1 = δω 2 = 0 the signalto-noise ratio is equal to 10 4 (a typical value for LISA sources [40]).
The contour plot of Fig. 3 shows the relative errors σ δω1 /δω 1 as a function of various combinations of (δω 1 , δω 2 ).Dashed white curves identify all the possible GR modifications that lead to uncertainties of 20% and 1%.We immediately note from the plot that the values of σ δω1 /δω 1 are almost completely degenerate with respect to δω 2 .More interestingly, all the parameter space sampled seems to indicate that LISA would be able to provide high-precision measurements of the frequency shifts, with a relative accuracy better than 1%.This is extremely promising for some of the GR modifications which are currently under investigation.For the EFT described in Sec.V A, and identified by a red star in Fig. 3, the frequency change δω 1 could be constrained with an error of ∼ 20%.
We obtain similar results for the second parameter δω 2 , with relative uncertainties being always below ∼ 50% for all the configurations analyzed.Note finally that a different choice of the waveform parameters, like the quality factor Q ℓmn = Re[ω ℓmn ]/Im[ω ℓmn ]/2, may provide better results.We defer more detailed investigations of the statistical errors, as well a multiband analysis including third-generation detectors [6,41], to future work.

VII. DISCUSSION
Our results yield (modulo an algebraic multiplication and sum) the QNM frequencies of any BH spacetime, provided that it describes a small deviation away from the Schwarzschild geometry.Thus, step 4 in the abstract can be avoided altogether by using our results, which are available online.
There are some important open questions, in particular the extension to coupled master equations [20].The expansion (6) describes many known theories, and is "reasonable."However, our formalism requires Υ in eq. ( 23) to be smaller than unity.This restricts the possible space of theories, as we discussed in the example of Eq. ( 24).In fact, we suspect that potentials localized far away from the horizon may be hard to accommodate: perturbation theory projects the perturbing potential onto the undisturbed QNMs, but these diverge exponentially at large distances.Thus, there are known, exponentially sensitive corrections when the potentials have structure at large distances [28,32,42].
Alternative expansions are possible, and have been used.Since the GR QNMs are strongly connected to the light ring [4,43,44], a tempting alternative is to expand the perturbing potential around the light ring [16,19].Such an approach is adequate for approximate, WKBlike expansions, but it may not be suitable for accurate numerical results.
It is intrinsic to our approach (and to many others) that the calculations are performed in a very special coordinate system, and the gauge invariance of the final results is not obvious.On the other hand, background solutions in generic theories are also usually written in very special coordinate systems.This weakness comes hand in hand with the power of the present approach: it handles any potential, but it is blind to the underlying theory or background.This blindness prevents us from understanding the connection between the expansion parameters in the potentials and the metric multipoles, or to explore the possible connection between geodesic motion and QNMs [44].Finally, we should highlight the fact that our results apply only to those modes which differ slightly from those of GR.The addition of a small perturbation to a known potential causes the appearance of both a perturbatively small correction to the modes, and new nonperturbative solutions.Take for instance the term β 0 for scalar perturbations, which is equivalent to a mass parameter µ 2 = β 0 /r 2 H .It is well known that for small µ, the GR QNMs are slightly displaced, and a new quasibound state with frequency ω 2 ∼ µ 2 − r 2 H µ 4 /(4(ℓ + 1) 2 ) branches off from a zero-mode in the zero-mass limit [26,45,46].We find, numerically and analytically, that a similar mode exists when α 0 is turned on.Such phenomena occur in electromagnetic waveguides as well [47,48].Our analysis is unable to cover systematically the emergence of these new families of modes.
thankfully acknowledge the computer resources, technical expertise and assistance provided by CENTRA/IST.Computations were performed at the cluster "Baltasar-Sete-Sóis" and supported by the H2020 ERC Consolidator Grant "Matter and strong field gravity: New frontiers in Einstein's theory" grant agreement no.MaGRaTh-646597.
our formalism to this case.From our QNM formula, the relation A numerical calculation shows that |1 − 2e − 2 /(3e − 3 )| is 0.016 for ℓ = 10, and 0.0082 for ℓ = 20.This is consistent with the fact that the error is of order O(ℓ −2 ).We should note that since our QNM formula contains an error of order O((α 2 ) 2 ) = O(ℓ −4 ), we cannot obtain relations for e j with j ≥ 4.

FIG. 2 .
FIG.2.Values of the frequency components e − j (red dots) for the odd tensor modes with ℓ = 2, compared against the numerical fit of Eq. (22) (black dashed curve).

FIG. 3 .
FIG. 3. Contour plot for the relative errors on the parameters δω1 which modifies the frequency of the ℓ = 2 fundamental BH mode.Vertical dashed lines represent configurations that lead to 20% and 1% fixed uncertainty.The red star identifies the values on (δω1, δω2) for the EFT model described in Sec.V A.

TABLE II .
Relative percentage errors on the real and imaginary parts of the QNMs for RN BHs, as a function of the charge-to-mass ratio Q/M .
and in the last line we have expanded the equation up to order O(Q 2 ).Here, Ω 0 is the dimensionless Schwarzschild

TABLE III .
Relative percentage errors in the real and imaginary parts of the QNM frequencies for scalar perturbations around a slowly spinning black hole, as a function of the BH angular momentum a/M .