An overview of quasinormal modes in modified and extended gravity

TL;DR

This paper reviews how quasinormal modes in various gravity models are affected by modifications, providing qualitative insights and deriving the Regge-Wheeler equation for static, spherically symmetric spacetimes.

Contribution

It offers a qualitative analysis of quasinormal mode frequencies in extended gravity theories and derives the Regge-Wheeler equation for general static, spherically symmetric metrics.

Findings

Qualitative tendencies of quasinormal mode frequencies in modified gravity.

Derivation of the Regge-Wheeler equation for general static, spherically symmetric metrics.

Insights into how gravitational wave observations can probe new physics.

Abstract

As gravitational waves are now being nearly routinely measured with interferometers, the question of using them to probe new physics becomes increasingly legitimate. In this article, we rely on a well established framework to investigate how the complex frequencies of quasinormal modes are affected by different models. The tendencies are explicitly shown, for both the pulsation and the damping rate. The goal is, at this stage, purely qualitative. This opportunity is also taken to derive the Regge-Wheeler equation for general static and spherically symmetric metrics.

Figures

Equations

$$e^{-i\omega t}=e^{-i({\omega }_R+i{\omega }_I)t},$$

$$V_{\ell }^{\text{RG}}(r)=\left(1-\frac{2M}{r}\right)\left[\frac{\ell (\ell +1)}{r^2}-\frac{6M}{r^3}\right],$$

$$V_{\ell }^{\text{Z}}(r)=\frac{2}{r^3}\left(1-\frac{2M}{r}\right)\times$$

$$\frac{d^{2}Z}{d{r}^{\ast 2}}+V(r)Z=0,$$

$$d{s}^2=f(r)d{t}^2-f(r{)}^{-1}d{r}^2-r^{2}d{\theta }^2-r^2{\sin }^2\theta d{\varphi }^2,$$

$$d{r}^{\ast }=\frac{1}{f(r)}dr.$$

$$V(r)=f(r)\left(\frac{\lambda +2(f(r)-1)}{r^2}-\frac{f^{\prime }(r)}{r}\right).$$

$$Z\sim e^{-i\omega r^{\ast }}{r}^{\ast }\to -\infty ,$$

$$Z\sim e^{i\omega r^{\ast }}{r}^{\ast }\to +\infty .$$

$$d{s}^2=A(r)d{t}^2-B(r{)}^{-1}d{r}^2-H(r)d{\theta }^2-H(r){\sin }^2\theta d{\varphi }^2.$$

$$\frac{d}{d{r}^{\ast }}=\sqrt{AB}\frac{d}{dr}.$$

$$d{s}^2=$$

$$e^{2\nu }=A(r),e^{-2{\mu }_2}=B(r),$$

$$e^{2{\mu }_3}=H(r),e^{2\psi }=H(r){\sin }^2\theta ,$$

$$\sigma =q_2=q_3=0.$$

$$R_{12}=R_{13}=0.$$

$$R_{12}=\frac{1}{2}{e}^{-2\psi -\nu -{\mu }_3}\times$$

$$[(e^{3\psi -\nu -{\mu }_2+{\mu }_3}{Q}_{02}{)}_{,0}-(e^{3\psi +\nu -{\mu }_2-{\mu }_3}{Q}_{32}{)}_{,3}],$$

$$Q_{ab}=q_{a,b}-q_{b,a}\mathrm{and}{Q}_{a0}=q_{a,0}-{\sigma }_{,a}\mathrm{for}a,b=2,3.$$

$$(H{\sin }^3\theta \sqrt{AB}{Q}_{23}){,}_3=-H^2{\sin }^3\theta \sqrt{\frac{B}{A}}{Q}_{02,0}.$$

$$Q=\sqrt{AB}H{Q}_{23}{\sin }^3\theta ,$$

$$\sqrt{\frac{A}{B}}\frac{1}{H^2{\sin }^3\theta }\frac{\partial Q}{\partial \theta }=Q_{20,0}.$$

$$\frac{\sqrt{AB}}{H{\sin }^3\theta }\frac{\partial Q}{\partial r}=-Q_{30,0}.$$

$$\sqrt{\frac{A}{B}}\frac{1}{H^2{\sin }^3\theta }\frac{\partial Q}{\partial \theta }=-{\omega }^2{q}_2-i\omega {\sigma }_{,2},$$

$$\frac{\sqrt{AB}}{H{\sin }^3\theta }\frac{\partial Q}{\partial r}={\omega }^2{q}_3+i\omega {\sigma }_{,3}.$$

$$\sin^{3}\theta\frac{\partial}{\partial\theta}\bigg{(}\frac{1}{\sin^{3}\theta}\frac{\partial Q}{\partial\theta}\bigg{)}+H^{2}\sqrt{\frac{B}{A}}\frac{\partial}{\partial r}\bigg{(}\frac{\sqrt{AB}}{H}\frac{\partial Q}{\partial r}\bigg{)}+\omega^{2}\frac{QH}{A}=0.$$

$$Q(r,\theta )=R(r)C_{l+2}^{-3/2}(\theta )$$

$$\bigg{(}\frac{d}{d\theta}\sin^{2m}\theta\frac{d}{d\theta}+n(n+2m)\sin^{2m}\theta\bigg{)}C^{m}_{n}(\theta)=0.$$

$$H^{2}\sqrt{\frac{B}{A}}\frac{\partial}{\partial r}\bigg{(}\frac{\sqrt{AB}}{H}\frac{\partial R(r)}{\partial r}\bigg{)}+\bigg{(}\frac{H}{A}\omega^{2}-\mu^{2}\bigg{)}R(r)=0,$$

$$\frac{d^{2}Z}{d{r}^{\ast 2}}+({\omega }^2-V(r))Z=0,$$

$$V(r)=\frac{1}{2H^{2}}\bigg{(}\frac{dH}{dr^{*}}\bigg{)}^{2}+\frac{\mu^{2}A}{H}-\frac{1}{\sqrt{H}}\frac{d^{2}}{dr^{*2}}\bigg{(}\sqrt{H}\bigg{)}.$$

$$Z\sim exp\bigg{[}\frac{1}{\epsilon}\sum_{n=0}^{\infty}\epsilon^{n}S_{n}(x)\bigg{]},\quad\quad\epsilon\rightarrow 0.$$

Full text

I Introduction

General relativity (GR) is our best theory of spacetime. Although the Lovelock theorem Lovelock (1971) ensures that it cannot be easily modified there are quite a lot of attempts to relax some hypotheses and build a deeper model to describe the gravitational field. From effective quantum gravity to improved infrared properties, the motivations to go beyond GR are countless. So are the situations, both in astrophysics and cosmology, where extended gravity theories can, in principle be tested. In practice, reaching the level of accuracy useful to probe the relevant range of parameters is obviously far from trivial. In this article we focus on a specific aspect of gravitational waves that would be emitted during the relaxation phase of a deformed black hole (BH).

We will consider quasinormal modes associated with the ringdown phase of a BH merger. The modes are not strictly normal due to energy losses of the system through gravitational waves. The boundary conditions for the equation of motion are unusual as the wave has to be purely outgoing at infinity and purely ingoing at the event horizon. The time component of the radial part reads (an introductory review can be found in Chirenti (2018))

[TABLE]

the complex pulsation $\omega$ being split in a real part $\omega_{R}$, which corresponds to the frequency, and an imaginary one $\omega_{I}$, which is the inverse timescale of the damping. Stability requires $\omega_{I}<0$. Although real-life BHs are spinning, we focus on Schwarzschild solutions in this article. The details of these predictions can not be used to directly compare with observations. We, however, expect the general tendencies and orders of magnitudes to remain correct, as it can be checked for the general relativistic case in Krivan et al. (1997).

The linearized Einstein equations lead to wave equations with different potentials whether one considers “axial” or “polar” perturbations. In GR, the (so called Regge-Wheeler) potential for axial perturbations is

[TABLE]

while the (so-called Zerilli) one for polar perturbations is

[TABLE]

where $a=\ell(\ell+1)/2-1$. Throughout all the paper we use Planck units. In the purely gravitational sector, one needs $\ell\geq 2$. Interestingly, both those equations have the very same spectrum of quasinormal modes (QNMs). This property, called isospectrality Chandrasekhar (1985) is not always true in modified gravity (see Moulin and Barrau (2019) for an extension and a discussion of the original proof). Basically, quasinomal modes are described by their multipole number $\ell$ and their overtone number $n$. The fundamental quadrupolar mode ($n=0$ and $\ell=2$) for a Schwarzschild BH in GR is given by $M\omega\approx 0.374-0.0890i$.

There are many different ways to calculate the QNMs: continued fractions, Frobenius series, Mashhoon’s method, confluent Heun’s equation, characteristic integration, shooting, WKB approximations, etc. In this article we focus on the last approach. For most models considered here, the QNMs have already been calculated in previous studies. However, this has most of the time been done for $s=0$ or $s=1$, not for $s=2$ as we have done it here. More importantly, it is in addition very useful to rely on the very same method to investigate all models so that the differences underlined are actually due to physical effects and not to numerical issues. Even when the same approach is considered, the way it is implemented is often different enough, between articles, so that it is hard to directly compare the results. This is why we have here tried to consider methodically several modified gravity models with a well controlled WKB approximation scheme used in the same way in all cases so as to compare the tendencies between modified gravity proposals. This is not mandatory for this qualitative step but this will become useful in future quantitative studies.

The determination of the complex frequencies of QNMs is difficult (see Kokkotas and Schmidt (1999); Nollert (1999) for historical reviews and Berti et al. (2004); Dorband et al. (2006) for results based on numerical approaches). This work is based on the WKB approach described in Konoplya (2003). Following the pioneering work Mashhoon (1982), the WKB method for QNMs was developed in Schutz and Will (1985); Iyer and Will (1987); Iyer (1987); Kokkotas and Schutz (1988). This formalism leads to fairly good approximations, especially for high multipole and low overtone numbers. In the following, we restrict ourselves to $n<l$ and use the 6th order WKB method developed by Konoplya Konoplya (2003) (see also Konoplya (2009); Konoplya and Zhidenko (2011); Konoplya et al. (2019)). This allows one to recast the potential appearing in the effective Schrödinger equation felt by gravitational perturbations in a complexe but tractable form.

The aim of this introductory paper is to investigate how several modified gravity theories impact the QNMs at the qualitative level. There are several ways to go beyond GR: extra dimensions, weak equivalence principe violations, extra fields, diffeomorphism-invariance violations, etc. Beyond those technicalities, there are strong conceptual motivations to consider extended gravity approaches, from the building of an effective quantum gravity theory to the improvements of the renormalisation properties, through the implementation of a dynamical cosmological constant. Among many others, examples of recent relevant works on QNMs can be found in Blázquez-Salcedo et al. (2016, 2017); Chen et al. (2019); Chen and Chen (2019).

II Perturbation dynamics

The QNMs are solutions of a perturbation equation with the specific boundary conditions given in the previous section. The radial and angular parts can be separated. The radial part is governed by a Schrödinger-like equation:

[TABLE]

where Z is the radial part of the “perturbation” variable, assumed to have a time-dependance $e^{i\omega t}$, and $r^{*}$ is the tortoise coordinate. For a metric such that

[TABLE]

the tortoise coordinate is defined by

[TABLE]

It tends to $-\infty$ at the event horizon and to $+\infty$ at spatial infinity.

As explained previously, BH gravitational perturbations can be of two different types distinguished by their behavior under a parity transformation. For an angular momentum $l$, axial perturbations transform as $(-1)^{l}$ under parity, while polar perturbations transform as $(-1)^{l+1}$. This leads to the two different potentials in Eq.(4). The potentiel for the gravitational axial perturbations reads in full generality (see Chirenti (2018) and references therein) for the metric given by Eq. (5):

[TABLE]

In this work we will not consider the isospectrality-violation issues and we will focus only on such perturbations. It should anyway be kept in mind that, in principle, isospectrality might not hold.

The boundary conditions can be expressed as

[TABLE]

We shall now derive the Regge-Wheeler equation for the more general (spherical and static) metric:

[TABLE]

For this metric, the tortoise coordinate is defined by

[TABLE]

The general form of an axisymmetric metric can be written as Chandrasekhar (1985):

[TABLE]

where $t=x^{0}$, $\phi=x^{1}$, $r=x^{2}$ and $\theta=x^{3}$. For the metric given by Eq. (10), the correspondance is:

[TABLE]

A perturbation of this kind of spacetime is described by $\sigma$, $q_{2}$ and $q_{3}$, assumed to be first order quantities, and by infinitesimal increments, $\delta\nu$, $\delta\mu_{2}$, $\delta\mu_{3}$, of the other quantities. We focus here on axial perturbations. The point is to linearize the field equations about the solution given by Eq. (10), considering components where $\sigma$, $q_{2}$ and $q_{3}$ are only function of $t$, $x^{2}$ and $x^{3}$. The equations governing $\sigma$, $q_{2}$ and $q_{3}$ are described by the vanishing of the Ricci tensor components:

[TABLE]

For Eq. (12), one has Chandrasekhar (1985)

[TABLE]

with

[TABLE]

The comma indicates the derivative. The notation $Q_{0a}$ is used to mean $-Q_{a0}$. The component $R_{13}$ is also given by Eq. (15) by switching indices $2$ and $3$.

The perturbed field equation are obtain by $\delta R_{\alpha\beta}=0$. After replacing $\nu$, $\mu_{2}$, $\mu_{3}$ and $\psi$ by their expressions, $\delta R_{12}=0$ leads to

[TABLE]

By defining

[TABLE]

one obtains

[TABLE]

For $\delta R_{13}=0$, one is led to

[TABLE]

We assume that perturbation have a time dependance given by $e^{i\omega t}.$ This implies that Eqs. (19) and (20) read

[TABLE]

[TABLE]

Taking the derivative of Eq.(21) with respect to $\theta$, the derivative of Eq. (22) with respect to $r$, and combining the results leads to:

[TABLE]

As suggested in Chandrasekhar (1985), one can then separate the variables $r$ and $\theta$ using

[TABLE]

with $C^{m}_{n}$ the Gegenbauer function satisfying

[TABLE]

Inserting Eq. (24) into Eq. (23), one is led to following radial equation:

[TABLE]

where $\mu^{2}=(l-1)(l+2)$. Defining $Z$ so that $R=\sqrt{H}Z$ and using the tortoise coordinate, we are led to a Schrödinger-like equation:

[TABLE]

where the potential is

[TABLE]

The potential reduces to Eq. (7) for $A(r)=B(r)$ and $H(r)=r^{2}$. This derivation is useful to calculate QNMs for general static and spherically symmetric metrics.

III The WKB approximation

The WKB approximation Schutz and Will (1985); Iyer and Will (1987); Iyer (1987) is known for leading to good approximations (compared to numerical results) for the QNMs. The potential is written using the tortoise coordinate so as to be constant at $r^{*}\rightarrow 0$ (which represent the horizon of the BH) and at $r^{*}\rightarrow+\infty$ (which represents spatial infinity). The maximum of the potential is reached at $r_{0}^{*}$. Three regions can be identified: region $I$ from $-\infty$ to $r_{1}$, the first turning point (where the potentiel vanishes), region $II$ from $r_{1}$ to $r_{2}$, the second turning point, and region $III$ from $r_{2}$ to $+\infty$. In region $II$, a Taylor expansions is performed around $r_{0}^{*}$. In regions $I$ and $III$, the solution is approximated by an exponential function:

[TABLE]

This expression can be inserted into Eq. (4) so as to obtain $S_{j}$ as a function of the potential and its derivative. We then impose the boundary conditions given by Eq. (9) and match the solutions of regions $I$ and $III$ with the solution for region $II$ at the turning points $r_{1}$ and $r_{2}$ (respectively). The WKB approximation has been usefully extended from the third to sixth order in (Konoplya, 2003).

This allows one to derive the complex frequencies as a function of the potential and its derivatives evaluated at the maximum. For the sixth order treatment, one is led to:

[TABLE]

where the expressions of the $\Lambda_{j}$s can be found in (Konoplya, 2003). In the following, we use this scheme to compare different modified gravity models and we present results only in the range of validity of the WKB approximations.

Interesting recent considerations on the convergence on the WKB series are given in Hatsuda (2019). Details on the expansion parameter used in this work can be found in Konoplya and Zhidenko (2011). The consistency of the WKB approximation has been checked for the presented results.

IV Modified gravity models and results

Throughout all this section we investigate some properties of the QNMs for several extended gravity approches. We pretend, in no way, to do justice to the subtleties of those models and, when necessary, we explicitly choose a specific of simplified setting to make the calculations easily tractable.

As we focus on phenomenological aspects, the more interesting mode is the fundamental one: $n=0$ and $l=2$. We therefore focus on a few points around this one (keeping in mind that the accuracy is better for higher values of $l$). In all the figures, the lower overtone $n$ is the one with the smallest imaginary part.

We first consider models with a metric of the form:

[TABLE]

and, then, investigate a model with two different metric functions, using the result obtained in Eq. (28).

IV.1 1. Massive gravity

In GR, the graviton is a massless spin-2 particle. One of the first motivations for modern massive gravity – which can be seen as a generalization of GR – was the hope to account for the accelerated expansion of the Universe by generating a kind of Yukawa-like potential for gravitation D’Amico et al. (2011). The initial linear approach to massive gravity was containing a Boulware-Deser ghost, which was cured in the dRGT version de Rham and Gabadadze (2010); de Rham et al. (2011); Hassan and Rosen (2011); de Rham (2014). Massive gravity also features interesting propertied for holography (see, e.g. Camara dS et al. (2013)).

Starting from the action

[TABLE]

where $R$ is the Ricci scalar and $\mathcal{U}$ is the potential for the graviton, the following black hole solution can be derived Ghosh et al. (2016); Eslam Panah et al. (2018):

[TABLE]

where $\Lambda$, $\gamma$ and $\epsilon$ are, respectively

[TABLE]

$a$ and $b$ being two dimensionless constants and $c$ is positive. It should also be pointed out that a positive value of $\gamma$ might raise consistency issues Eslam Panah et al. (2018).

The results are presented in Fig. 1. The values chosen for the constants do of course change the amplitude of the displacement of the QNMs. The global trend, which is the point of this study, however remains the same. Increasing of the graviton mass $m$ tends to increase the real part of QNMs, that is the frequency of the oscillations. The difference in frequency between the fundamental and the first overtone also increases with $m$. The effect on the imaginary part is hardly noticeable on the plot even though a slight increase should be noticed, which is actually 50% less important, in relative variation, than the shift in frequency. The values considered here for the mass are, of course, way out of the known bounds but this is clearly not the point. As a specific feature, one can notice that the frequency shift due to massive corrections decreases for higher overtones. The shift patterns are mostly the same whatever the multipole number considered.

IV.2 2. Modified STV gravity

The Scalar-Tensor-Vector modified gravitational theory (MOG) allows the gravitational constant, a vector field coupling, and the vector field mass to vary with space and time Moffat (2006). The equations of motion lead to an effective modified acceleration law that can account for galaxy rotation curves and cluster observation without dark matter. Although it has recently been much debated and put under pressure, the theory is still worth being considered seriously. We consider the field equation for the metric tensor Moffat (2015) :

[TABLE]

where the gravitational coupling is $G=G_{N}(1+\alpha)$, with $G_{N}$ the Newton’s constant. The gravitational strength of the vector field $\phi_{\mu}$ (spin 1 graviton) is $Q_{g}=\sqrt{\alpha G_{N}}M$. With $B_{\mu\nu}=\partial_{\mu}\phi_{\nu}-\partial_{\nu}\phi_{\mu}$, the energy-momentum tensor for the vector field is :

[TABLE]

the constant $\omega$ of Moffat (2006) being set to one. Solving the vacuum field equations

[TABLE]

and

[TABLE]

with the appropriate symmetry leads to the metric

[TABLE]

We focus on the case where the field equations for $B_{\mu\nu}$ are non-linear, as the phenomenology is then richer, and we consider the relevant choice $\alpha<\alpha_{c}=0.67$ where there are two horizons and an appropriate potential behavior for the WKB approximation to hold. An up-to-date investigation of QNMs in MOG can be found in Manfredi et al. (2018).

The results are given in Fig. 2. The imaginary part of the QNMs is nearly the same whatever the value of $\alpha$: the modified metric has no effect on the damping rate. However, increasing $\alpha$ does increase of the real part, that is the frequency. The effect is important for values near the critical value $\alpha_{c}$. The slope of the Imaginary part versus the real one, at a given $l$ for different values of $n$, is nearly independent of $\alpha$. This slope is not directly observable but it shows how the structure of the QNMs changes with the overtone number. The curves remain here parallel one to the other: this means that increasing the deformation parameter does not change the frequency shift between overtones.

IV.3 3. Hořava-Lifshitz gravity

Hořava-Lifshitz gravity bets on the fundamental nature of the quantum theory instead of relying on GR principles. It is a renormalizable UV-complete gravitational theory which is not Lorentz invariant in 3 + 1 dimensions Horava (2009). The relativistic time with its Lorentz invariance emerges only at large distances. Black hole solutions have been found O Colgain and Yavartanoo (2009); Kehagias and Sfetsos (2009); Lu et al. (2009) and QNMs were studied Chen and Jing (2010).

Using the ansatz

[TABLE]

in the action, one is led to the Lagrangian

[TABLE]

wherei $w=8\mu^{2}(3\lambda-1)/\kappa^{2}$. For $\lambda=1$, the solution is

[TABLE]

with $\beta=1/(2w)$, $w$ being the deformation parameter enterring the action given in Kehagias and Sfetsos (2009). There are two horizons for $M^{2}>\beta$.

The results are given in Fig. 3. The frequency increases with an increase of $\beta$. Interestingly, the imaginary part of the overtones is highly sensitive to $\beta$. This remains true for higher multipoles. The relative variation of the imaginary part is nearly the same whatever the overtone number. It therefore becomes large in absolute value for high $n$ values.

IV.4 4. $\hbar$ correction

It has been known for a long time that quantum corrections to the Newtonian gravitational potential can be rigorously derived without having a full quantum theory of gravity at disposal (see, e.g., Donoghue (1994a, b); Bjerrum-Bohr et al. (2003a, b); Bjerrum-Bohr et al. (2015) to cite only a few works from a very long list). Recently, a quite similar approach was developed Bargueño et al. (2017) requiring that the quantum mechanically-corrected metric reproduces the corrected Newtonian limit, reproduces the standard result for the entropy of black holes including the known corrections, and fulfills some consistency conditions regarding the geodesic motion.

The resulting metric is

[TABLE]

We use, as previously, natural units and the coefficients of the last term, $\gamma$, is proportional to $\hbar$ in these models. It is worth noticing that there has been a long controversy about the value and the sign of the $\gamma$ factor. From the phenomenological perspective, we do not fix it to a particular value but we keep it negative, in agreement with the latest expectations.

The results are given in Fig. 4. For large values of $\gamma$, the effects are noticeable on the frequency. It is remarkable that, from our analysis, the real part of the complex frequency is only decreased, which is not the case for the other models that have been considered in this study. The higher the absolute value of $\gamma$, the larger the difference of frequency between the fundamental and the overtones. This effect however remains quite subtle.

IV.5 5. LQG polymeric BH

Loop Quantum Gravity (LQG) is a non-perturbative and background-independant quantum theory of gravity Ashtekar (2013). In the covariant formulation, space is described by a spin network Rovelli (2011). Each edge carries a “quantum of area”, labelled by an half integer $j$, associated with an irreducible representations of $SU(2)$. Each node carries a “quantum of space” associated with an intertwiner. A key result is that area is quantized according to

[TABLE]

with $\gamma_{BI}$ the Barbero-Immirzi parameter. Black holes are usually described in LQG through an isolated horizon puncturing a spin network Perez (2017) and the phenomenology is very rich, depending on the precise setting chosen Barrau et al. (2018). We focus here on the model developed in (Alesci and Modesto, 2012), as this is the one leading to metric modifications outside the horizon, where a regular lattice with edges of lengths $\delta_{b}$ and $\delta_{c}$ is considered. Requiring the minimal area to be one derived in LQG, one is left with only one free parameter $\delta$. From this minisuperspace approximation, a static spherical solution can be derived and is given by

[TABLE]

where $d\Omega^{2}=d\theta^{2}+\sin^{2}\theta d\phi^{2}$, $r_{+}=2m$ and $r_{-}=2mP^{2}$ are the two horizons, and $r_{*}=\sqrt{r_{+}r_{-}}=2mP$, $P$ being the polymeric function defined by $P=(\sqrt{1+\epsilon^{2}}-1)/(\sqrt{1+\epsilon^{2}}+1)$, with $\epsilon=\gamma_{BI}\delta$, and the area parameter $a_{0}$ is given by $a_{0}=A_{min}/8\pi$, $A_{min}$ being the minimum area appearing in LQG. The parameter $m$ in the solution is related to the ADM mass $M$ by $M=m(1+P)^{2}$.

The results are given in Fig. 5. The damping rate does not depend at all on the polymerization parameter. The real part of the complex frequency does, however, first decrease with $\delta$. Noticeably, the slop is unchanged and varying the deformation parameter just lead to a horizontal translation of the QNM frequency in the complex plane. This means that the frequency shift between the fundamental and the overtones does not depend on the amplitude of the quantum gravity corrections, as in modified gravity. Interestingly, for higher values of $\delta$, the frequency begins to increase. This is the only model considered in this study with a non-monotonic behavior. For $\delta\approx 10^{-0.7}$ the “polymerization” effect nearly exactly compensates the “area discretization” effect and one recovers the GR frequencies (and damping rates).

V Conclusion

This study shows the evolution of the complex frequency of quasinormal modes of a Schwarzschild black hole for the fundamental and the first overtones for a few multipole numbers. We have considered massive gravity, STV gravity, Hořava-Lifshitz gravity, quantum corrected gravity, and loop quantum gravity. All the results were derived using the very same WKB approximation scheme which makes a meaningful comparison possible. It will be especially useful for future quantitative studies.

Obviously, distinguishing between those models with observations is more than challenging. First, because there exist degeneracies, for given overtone and multipole numbers, between the models – when taking into account that the values of the parameters controlling the deformation are unknown. Second, because the intrinsic characteristics of the observed black holes are also unknown, which induces other degeneracies. In addition, this study should be extended to Kerr black hole, which also adds some degeneracies in addition to the complexity.

Some interesting trends can however be underlined. For all models, the effect of modifying the gravitational theory are more important for the real part than for the imaginary part of the complex frequency of the QNMs. Otherwise stated, the frequency shift is more important than the change in the damping rate. Obviously, it does not make sense to quantitatively compare the results from various models as the deformation parameters are different. However, the “trends” are clearly specific to each studied theory and there is no need to define comparable “steps” in the deformation parameters (which do not have the same units anyway) to draw significant conclusions about the directions in which the different models considered deviate from GR. In addition, the sign of the frequency shift, and its dependance upon the overtone and multipole numbers is characteristic of a given extension of GR. The accurate patterns are never the same, which is an excellent point for phenomenology. It can basically be concluded that a meaningful use of QNMs to investigate efficiently modified gravity requires the measurement of several relaxation modes. This is in principle possible Sathyaprakash et al. (2012) but way beyond the sensitivity of current interferometers. If features beyond GR were to be observed, the direction of the frequency shift in the complex plane would already allow to exclude models, as this article shows. The goal of this study was not to perform a detailed analysis of the discrimination capabilities of gravitational wave experiments: it simply aimed at exhibiting the main tendencies for currently considered extended gravity models, as an introduction to this special issue on “probing new physics with black holes”.

VI Acknowledgments

K.M. is supported by a grant from the CFM foundation.