Quasi-normal modes of static spherically symmetric black holes in f(R) theory

We study the quasi-normal modes (QNMs) of static, spherically symmetric black holes in f(R) theories. We show how these modes in theories with non-trivial f(R) are fundamentally different from those in general relativity. In the special case of f(R)=αR2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(R) = \alpha R^2$$\end{document} theories, it has been recently argued that iso-spectrality between scalar and vector modes breaks down. Here, we show that such a break down is quite general across all f(R) theories, as long as they satisfy f′′(0)/(1+f′′(0))≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f''(0)/(1+f''(0)) \ne 0$$\end{document}, where a prime denotes derivative of the function with respect to its argument. We specifically discuss the origin of the breaking of isospectrality. We also show that along with this breaking the QNMs receive a correction that arises when f′′(0)/(1+f′(0))≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f''(0)/(1+f'(0)) \ne 0$$\end{document} owing to the inhomogeneous term that it introduces in the mode equation. We discuss how these differences affect the “ringdown” phase of binary black hole mergers and the possibility of constraining f(R) models with gravitational-wave observations. We also find that even though the iso-spectrality is broken in f(R) theories, in general, nevertheless in the corresponding scalar-tensor theories in the Einstein frame it is unbroken.


Introduction
Recent observations of gravitational wave (GW) signals in LIGO and Virgo from compact object mergers have opened up a new chapter in the history of physics [1][2][3][4]. So far these detectors have already observed several binary black hole mergers and a binary neutron star merger [5]. The latter observation has put strong constraints on the equation of state of the matter inside neutron stars [6][7][8][9][10]. It has also provided a strong bound on the graviton mass [11] and the deviation of the velocity of GWs from that of the light [5,12]. As a result it has already helped us understand not only the merga e-mail: skdatta@iucaa.in b e-mail: sukanta@iucaa.in ers themselves but also fundamental aspects of the nature of gravity [12][13][14][15][16][17]. This scenario is expected to get even more interesting when LISA is launched in several years from now since the space detector will give us the opportunity to test gravity in a different frequency band [18][19][20][21][22].
To date general relativity (GR) has been very successful in observational tests of its predictions (see, e.g., Refs. [23,24]). However, certain aspects of GR, e.g., black hole and cosmological spacetime singularities, have invited proposals for higher-curvature correction terms to the Einstein-Hilbert action [25][26][27][28][29][30][31][32][33] in attempts toward resolving these pathologies. These corrections come in various forms, such as functions of Ricci scalars and different combinations of the contractions of the Riemann tensor. These suggestions have gained importance due to other reasons as well. It has been shown that several actions with correction terms can drive inflation and others have been successful in explaining the late time acceleration of the universe [34][35][36][37]. Such actions have also arisen from the low-energy limit of quantum corrections or string theory [38][39][40] and in Loop Quantum Gravity [43,44]. These results have initiated extensive research in such alternative theories of gravity in the last few decades. Several studies have already investigated modifications in QNMs for other types of corrections (quadratic in curvature) [41,42]. Now the availability of GW observations helps us examine if these theories can be subjected to useful tests.
As has been shown in numerical relativity, when the two black holes in a binary system merge, a distorted black hole is created. This remnant radiates GWs as it settles down into a Kerr black hole [45][46][47][48][49]. Asymptotically, these GWs can be expressed as superpositions of damped sinusoidal modes, termed as quasi normal modes (QNMs). These QNMs in GR depend only on a couple of parameters characterizing the black hole, namely, its mass and spin, for astrophysical black holes. This is a manifestation of the No-Hair theorem [50][51][52]. Therefore, observation of these QNMs in the future is an anticipated testing ground for the No-Hair theorem and alternative theories of gravity.
It is well known that in GR there are two modes that contribute to GW observables. These are the odd (vector) and even (scalar) modes, and they share the net emitted gravitational energy equally. It has been argued in Refs. [53,54] that in the case of the (R + α R 2 ) action, where R is the Ricci scalar, QNMs emitted from black holes will not have the same spectrum as in GR. In this work, we show that such a break down in iso-spectrality is quite general across all f (R) theories, as long as they satisfy f (0)/(1+ f (0)) = 0, where a prime denotes derivative of the function with respect to its argument. We specifically discuss the origin of the breaking of isospectrality. We also show that along with this breaking the QNMs receive a correction that arises when f (0)/(1 + f (0) = 0 owing to the inhomogeneous term that it introduces in the mode equation. We also discuss the resulting structural modification of the QNMs.
In an important piece of work, Tattersall et al. [55,56] demonstrated that in Horndeski theory of gravity the QNM equation acquires a source term, which results in the breaking of iso-spectrality. Going by these results alone, one might expect that when an f (R) theory is related to an STT through a conformal transformation, and both theories have the Schwarzschild solution, then iso-spectrality should break in that f (R) theory as well for QNMS in that solution. Interestingly, and perhaps somewhat counterintuitively, we show here that when an f (R) theory and its conformally related STT both have the Schwarzschild metric as solution, and the value of the conformal factor relating these respective solutions is unity, then isopectrality is not broken for the Schwarzschild QNMs in STT even though it is broken for the Schwarzschild QNMs in f (R).
In Sect. 2 general results of an arbitrary f (R) theory and the modified Einstein equations in such theories are discussed. In Sect. 3 we present the perturbation equations in the context of GR around the Schwarzschild metric. In Sect. 4 we discuss the perturbation of the modified Einstein equation in general f (R) theories. In Sect. 4.1 we deduce the different modes and the final perturbation equations in f (R) theories. Then in Sect. 5 we deduce the changes in the QNMs in comparison to GR. Finally, in Sect. 7 we summarize our results. Throughout the work we have adopted the {−, +, +, +} signature and set G = c = 1.

General equations for f (R) theories
As mentioned above, here we focus our attention on f (R) theories of gravity. Our primary motivation for doing so is that they offer resolutions for certain limitations in the standard cosmological model and the fact that f (R) terms arise as correction terms from loop quantum gravity. Another rea-son is that these theories are not affected by the Ostrogradsky instability [57]. Among all the f (R) theories we only focus on those that allow R = 0 solutions, for the obvious reason that such solutions exist in GR and are observationally relevant. We discuss this aspect in more detail later.
Let us begin by considering the gravity model that has a Lagrangian of the following form: where κ 2 = 8π . Varying this action with respect to the metric g μν yields the equation of motion [58], It is known that Eq.
(2) has a constant curvature solution, R =R. In that case, the above equation takes the form Taking its trace gives Since our objective here is to study perturbations of the Schwarzschild black hole, we focus only on such models of f (R) that haveR = 0 as a root of Eq. (4). It is straightforward to verify that GR is but one example of such theories. Since our main goal is to look for deviations from GR, we will separately track the GR part of the f (R) action, namely, the R term. Therefore, we will express the Lagrangian as where ϕ(R) denotes terms in the gravity action beyond GR. Henceforth, we rename ϕ(R) as f (R). Thus, our working Lagrangian becomes R + f (R).

QNMS in general relativity
For the purpose of our calculations we adopt the gauge invariant formalism discussed in Refs. [59][60][61][62]. Our focus remains the Schwarzschild spacetime. In the case of a spherically symmetric spacetime, it is possible to split the manifold into an orbit space and a unit sphere. We will be focusing on the four-dimensional non-spinning black hole background.
In that case the manifold splits into a two-dimensional orbit space and a two-sphere. We will use the co-ordinates x a for the orbit space and z A for the two-sphere. The covariant derivatives on the orbit space and the two-sphere are represented by D a and D A respectively. The d'Alembertian operators on the orbit space and the two-sphere are represented by˜ andˆ , respectively.
With proper choice of co-ordinates it is possible to cast the metric in the following form, where the overbar represents quantities of the unperturbed background, g(r ) = 1 − 2M r , and d 2 represents the metric on the unit two-sphere. The metric perturbation on the background can be expressed as follows [59][60][61][62]: where ds 2 p represents the perturbed part of the line element. To simplify the problem, the metric functions are usually expanded in spherical harmonic basis. After linearizing the Einstein equations, they can be cast in terms of the coordinates of the orbit space alone. This problem has been well studied. We follow the formalism defined by Kodama, Ishibashi and Seto (KIS). With proper identification of the scalar and vector mode it is possible to show that for each multipole l ≥ 2 they satisfy the following equations [59][60][61]: where λ = (l−1)(l+2) 2 and the time dependence has been taken to be e −iωt . The subscripts S and V represent the scalar mode (Zerilli mode) and the vector mode (Regge-Wheeler mode), respectively. Moreover, "GR" in the superscript distinguishes these modes from those in f (R).
These QNMs have been studied extensively in the literature. It has been shown that the transmission and the reflection coefficients of the V S and V V are equal. They make equal contributions to GWs asymptotically. It has also been demonstrated that they share the same frequency spectrum [49]. We will show below that this does not hold, in general, in f (R) theories.

Perturbation in general f (R)
Our main goal is to investigate the QNM structure of the Schwarzschild black hole. Therefore, we take the background metric to be that given in Eq. (6). Now the perturbed metric can be written as whereḡ μν and h μν represent the background metric and its perturbation, respectively. Using the perturbed metric from Eq. (10) in the equation of motion Eq. (2), we find the equation for the perturbation as follows, where δ R μν and δ R are the perturbed Ricci tensor and Ricci scalar, respectively. Equation (11) can be re-expressed as where , Taking the trace of Eq. (11) it can be shown that where h is the trace of the perturbation. As we show more explicitly below this equation contains a massive scalar degree of freedom that is generically present in f (R) theories, unlike general relativity. Further details are available in Refs. [63][64][65]. This equation is equivalent to the equation of a massive scalar field with a source term that depends on the perturbation.

QNMs in f (R)
In the present work our main goal is to find the perturbed equation of Schwarzschild like solutions in f (R). For that reason we can fixR = 0 along side Eq. (4). This translates to f (0) = 0. Under these conditions Eq. (14) becomes This equation is a source-free massive scalar equation, with δ R identified as the scalar field and the mass-squared being m 2 ≡ 1+ f 3 f . Henceforth, whenever there is no scope for ambiguity, we will omit the 0 from the argument value of f and all its derivatives. This massive scalar longitudinal mode [66,67] is absent in GR. Now the perturbed Einstein equation takes the following form: accompanied by the simplified where β = f 1+ f . Heretofore, we will use both expressions interchangeably. The important thing to note is that there are possible f (R) theories, such as α R 3 , for which f = 0 (at R = 0); in that case β becomes zero. Another key point is that the information of a particular f (R) theory enters only through the β. So, apart from β, dynamically speaking, all f (R) theories are the same. Now we use the separation of variables to isolate the angular dependence. For this reason we write the perturbed Ricci scalar as are the scalar spherical harmonics. For radial coordinate we sometimes utilize the tortoise coordinate r * , defined The time dependence is taken to be e −iωt . After simplification, Eq. (14) takes the following form: wherẽ This result for f (R) = α R 2 theory was discussed in Ref. [53]. It is well known that each component of a tensor behaves differently under rotation group. As a result, the effective energy momentum tensor can be separated into scalar, vector and the tensor modes [59][60][61], as follows: where S, V A and T AB are the pure scalar, vector and tensor spherical harmonics, respectively. Rest of the tensors are defined from the pure spherical harmonic tensors. Further details, including the expression for T can be found in Refs. [59][60][61] and in Appendix A.
We find the various components of the energy momentum tensor from Eq. (20) to be The scalar and vector master equations get modified as follows: where and It should be noted that as β → 0 these equations, appropriately, reduce to their GR counterparts even if the theory is We also find that the tensor mode gives rise to the following equation [61], where˜ represents the (r * − t) part of the d'Alembertian operator and Above, λ L is the eigenvalue of the Lichnerowicz operator [61]. Therefore we see that the tensor and the vector modes do not undergo any modification. This is understandable because the effective energy-momentum tensor arises from the different combinations of the derivatives of the Ricci scalar. It has been shown in Ref. [68] that the tensor harmonic functions are identically zero in transverse-traceless (TT) gauge. For this reason we will not pursue the study of the tensor mode any further. Since the background in our studies is a R = 0 solution, its metric does not have any extra hairs. However, we notice from Eq. (22) that even though the background black hole metric has no additional hair and is similar to the black holes in GR, at the level of perturbations the Zerilli mode gets modified from GR due to an extra source term. This source term originates because at the level of perturbations the extra massive mode that was absent in the background now gets excited. The excited massive mode "generates" hair via β. This tells us that the perturbations will therefore get modified from GR. As a result, we can expect that the detection of QNMs through gravitational wave observations will, in principle, provide us with the opportunity to distinguish between GR and f (R) theories of gravity.

Solution in radiation zone
As we can see from the structure of the equations, the scalar equation gets modified by a source term. Therefore, it is understandable that the solution of that equation will have two parts: (a) One of these will be for the homogeneous part of the equation; this is identical to the GR solution; (b) The second will be the solution of the inhomogeneous part. Solution of the vector equation will be identical to the GR solution. Therefore, we can write This implies that from the observational perspective the nonzero inhomogeneous part will be the signature of any deviation from GR. If we take the r → ∞ limit, then we find that where H = k 2 − 2. However, due to Eq. (18), The source term becomes Equation (30) has solution of the form e iχr * , with From Eq. (22) we find, This gives the form of the inhomogeneous term arising due to the f (R) theory: Therefore, one can see that there is a deviation from GR in this f (R) theory, howsoever tiny its effect might be on observables.
In Eq. (33) we derived the asymptotic deviation of the f (R) QNMs from GR. This provides the opportunity to constrain f (R) theories by using it. In principle, there exists the possibility of bounding β by observationally tracking how much energy is lost by QNMs and if that deviates from the prediction of GR. Solar system tests have already constrained f (R) theories (details can be found in Appendix E). In comparison, the results found in the current work lend themselves to possible tests using GW observations. The inhomogeneous term arises due to the existence of the massive polarization mode, as seen from Eqs. (15) and (17). This occurs because the scalar mode and the massive mode are coupled, as was shown in Eq. (22). Through this inhomogeneous term, energy transfer can occur between the massive mode and the scalar mode. Therefore, tests for non-GR polarizations and energy content in each polarization may, in principle, constrain these f (R) theories.
Another observational aspect of these findings is the breaking of iso-spectrality. It has been shown in Ref. [49] that the scalar and vector modes have identical spectra in GR. But now due to the source term it becomes evident that in general f (R) theories this iso-spectrality will not be satisfied. The homogeneous equation will have an iso-spectral solution to the vector mode but the inhomogeneous part will not be iso-spectral. This point has been discussed in Sect. 5.2. Reference [53] claimed that the breaking of iso-spectrality will be there for all f (R) theories. We prove above that in general the iso-spectrality will be violated if the theory has nonzero β. The details of the breakdown of the iso-spectrality has been discussed in 5.2.

Broken iso-spectrality
Breaking of iso-spectrality between the scalar and vector modes in the context of α R 2 theory has already been discussed in Ref. [53]. Even though they claim it to be a universal phenomenon for f (R) theories, we showed above that this is not true unless β = 0. For this reason we will investigate for the first time in detail exactly how iso-spectrality breaks down when β = 0.
Iso-spectrality in GR was first discovered by Chandrasekhar [49]. We will follow his methods for this investigation. Say, Z 1,2 satisfy the following equations w.r.t. r * : where V 1 and V 2 are the corresponding potentials. Using the ansatz that Z 1 = pZ 2 + q Z 2 and taking its derivative twice it is possible to show where the prime denotes derivative w.r.t. r * . Comparing it with Eq. (35) Chandrasekhar showed that Now, finding a solution for p and q establishes a relation between Z 1 and Z 2 . For the QNMs of the Schwarzschild solution in GR it is possible to identify, It is also possible to find p and q. As a result, one can show that where λ = (l−1)(l+2) 2 . Due to this relationship it becomes evident that if S depends on the radial coordinate as e iωr * then so does V . Therefore, they share the same spectrum [69].
By inspecting Eq. (35) it is clear that the reason for such a simplification is its homogeneous nature. If we modify Eq. (35) with a source term, as in then an ansatz of the previous form, Z 1 = pZ 2 + q Z 2 , again leads to But comparison with Eq. (39) does not lead to equations of p, q that are independent of Z 1 . So, in the sourced case it is not possible to have Z 1 = pZ 2 + q Z 2 as a solution. This eventually results in the violation of iso-spectrality.
One key point is worth mentioning. It is always possible to separate Z 1 as and Therefore, it is possible to have Z G R 1 = pZ 2 + q Z 2 and, consequently, there will exist an iso-spectral solution. Hence, S in f (R) theories will have a part that is iso-spectral to V and another part that arises solely due to the source term that does not obey iso-spectrality.

Massive mode
As already shown above in Eq. (14) there exists a massive mode in f (R) theories. The equation for this scalar mode is exactly equivalent to the equation of a scalar field around a Schwarzschild black hole in GR. The QNM of a massive scalar field around a Schwarzschild black hole in GR has been studied extensively [70,71]. Because of the similarity between the two problems, the results for the QNM frequency of the massive mode relevant for our current work will be exactly the same as the ones found in those works. Naturally, the QNM frequency values found in those works depend on the mass of the scalar field. In our current work that mass depends on the specific f (R) theory chosen. Even though there are stringent restrictions on the value of this mass from solar system observations (see Appendix E), it is possible that independent restriction may be found from GW observations. Current constraints on the mass imply that if an f (R) theory is the correct theory of gravity then, m ≥ 1.25 × 10 4 m −1 .
One key point is worth mentioning here. Usually ω is a complex number, which implies that χ should be such a number as well. Consequently, the QNMs are damped. On top of that, the large value of the mass of the massive mode will give rise to a strong damping pattern. As a result, the contribution of the inhomogeneous part in Eq. (33) will be very less. But massive scalar fields in Schwarzschild backgrounds can have purely real ω [72]. This originates mainly from the sub-dominant asymptotic contribution arising due to the irregular singularity at infinity [72]. These modes are called quasi-resonance frequency modes [71]. The requirement for the existence of such modes is ω Q RM < m [72], where m is the mass of the massive mode. Further numerical study revealed that mass has a crucial role in the damping of that mode. They showed that the greater the mass the lesser is the damping rate. Therefore, purely real modes originate corresponding to the non-damping oscillation. They also showed that for a given mass of the field a certain number of lower overtones disappear. However, this disappearance happens only for the lower overtones, while the remaining overtones are still damped. Hence, owing to the non-zero mass of the massive mode the inhomogeneous term will contain some quasi-resonance contribution that will not decay as fast as the other damped overtones.

Connecting with scalar-tensor theories
It has been shown in several works that there is an equivalence between f (R) theories and Brans-Dicke scalar-tensor theories (STT) [30,[73][74][75][76][77]. In this section we investigate what implication can be found from such equivalence for a perturbed Schwarzschild black hole. For this reason first we will discuss the perturbation equations in STT around a Schwarzschild black hole in Einstein frame. Then we discuss the connection between the QNMs in f (R) theory and STT.
The Einstein equations in STT can be written as, where, and the scalar field φ satisfies, In this work we focus on the perturbation of a hairless Schwarzschild BH. Therefore, the LHS of Eq. (43) is zero for the background. This implies that for the background one has φ = φ C , where φ C is some constant value of the scalar field φ that satisfies V (φ C ) = 0. The perturbed part of Eq. (43) takes the following form, where δφ is the perturbation in scalar field φ.
For simplicity we will write V (φ)| φ=φ C as V (φ C ). To isolate the dynamic part of the perturbation we separate δφ in orbit space and angular space, as has been described earlier. For that purpose we write δφ = ST T (y a ) r S(z A ) ≡ ST T (y a )S(z A ), where S(z A ) are the scalar spherical harmonics. For radial coordinate we use the tortoise coordinate (r * ). The time dependence is taken to be e −iωt . After simpli-fication, Eq. (45) takes the following form: wherẽ The scalar and vector master equations get modified as follows: where Details of the calculation can be found in the Appendix C.

Radiation zone solution and broken iso-spectrality
As we can see from the structure of the equations, the scalar equation gets modified by a source term. Therefore, it is understandable that the solution of that equation will have two parts: (a) One of these will be for the homogeneous part of the equation, and will be identical to the GR solution; (b) The second part will be the solution of the inhomogeneous part. The solution of the vector equation will be the same as the corresponding case in GR. Therefore, we can write From the observational perspective this implies that the nonzero inhomogeneous part is the signature of deviation from GR. If we take the r → ∞ limit then we find that the source term becomes where H = k 2 − 2. Asymptotically Eq. (47) has solution of the form e iχr * with Hence from Eq. (22) we find, Therefore, the inhomogeneous term arising in the STT: Therefore, we find that there is a deviation from GR in STT, as is manifest when comparing with Eq. (28), howsoever tiny the effect of this term might be on observables. Therefore, in STT also the iso-spectrality is broken as long as V (φ C ) = 0.

Equivalence with f (R)
In this section we will discuss the correspondence between the QNMs of a Schwarzschild BH in STT and those of such a BH in f (R) theory. Let us begin by considering the f (R) gravity model (known as Jordan frame) that has a Lagrangian of the following form: where f (R) is the same as that used in Eq. (1); we will retain this definition until Eq. (61) below. Conformal transformation of the Jordan frame metric yields the new metric g μν = 2 g μν , where is the conformal factor. Under the identification the action reduces to which is known as the Einstein-frame action; here, the potential V (φ) can be written as In our convention, all the metric dependent quantities that originate from a conformal transformation of the Jordan frame are boldfaced. The field equation for g μν can be found by varying the Einstein frame action, and is given by [73,85] In the current work our objective is to compare the QNMs in these two theories. The background metrics are taken to be identical. In general, however, it is not possible to have identical background in both the theories that are related to each other through Eq. (58). We now return to the usage f (R) → R + f (R), as has been described in Sect.

For
Schwarzschild solution in f (R) we find 2 = 1+ f (0) from Eq. (58). For both of the theories to have the same metric as conformally related solutions we need 2 = 1, in case of Swarzschild solution it is achievable only if f (0) = 0. This is not true in general. Hence, in general, a Schwarzschild solution in f (R) theory does not conformally transform to a Schwarzschild solution in STT, even though it is equivalent to some solution in STT. Therefore, QNMs in f (R) theory do not conformally transform to QNMs of the same background in STT in general. This property is not limited to R = 0 solutions alone. For a general background, one has 2 = 1 + f (R), and the two conformally related theories can have identical metric solutions if and only if f (R) = 0. If this condition is not satisfied then even though the corresponding theories are conformally equivalent to each other, they can not have an identical metric solution that are conformally related.
As we are discussing QNMs of Schwarzschild BHs, we will focus on the case 2 = 1 + f (0). When f (0) = 0 a Schwarzschild solution in f (R) theory corresponds to a Schwarzschild solution in the corresponding STT. For this reason we will focus on this particular category. In this case 2 = 1 for the background. Therefore, for perturbation we can write where δφ and δ R are perturbed scalar field in STT and perturbed Ricci scalar in f (R) theory, respectively. If we separate out the angular dependence as δφ = ST T (y a )S(Z A )/r , then we can identify From Eq. (62), we see that in STT the background scalar field,φ = 0. Similarly in f (R) theory extra scalar mode is not present. Therefore, the extra scalar mode gets excited in the presence of a perturbation; otherwise it remains at zero. Vector QNMs are identical to GR in both the theories. Therefore, the vector mode in f (R) theory is equivalent to the vector mode in STT; in fact, they are identical. As has already been discussed, the scalar QNM in both the theories have a homogeneous part that is iso-spectral to the vector QNM and an inhomogeneous part that is not isospectral to the vector QNM. The homogeneous part in both the theories is identically equal to the GR QNM. Hence, the homogeneous part of the scalar mode in f (R) theory is equivalent to the homogeneous part of the scalar mode in STT: they are identical.
The difference arises in the inhomogeneous part. This difference can be understood by computing δV (φ), where δ R can be represented in terms of δφ from Eq. (62). ForR = 0 this gives dV (φ) dφ = 0. Hence, the source term in the Eq. (49) vanishes for STT, leading to a vanishing inhomogeneous solution. The Regge-Wheeler and Zerilli mode in these STTs are identical to GR.
This can be understood further with an example, namely that of the f (R) = α R 2 theory. The conformal transformation can be expressed as follows: The potential of the corresponding STT can be written as On the other hand, if we take f (R) = α R 2 in Eq. (22), we find that β = 2α, which leads to a non-vanishing source term. Therefore, the QNMs in that theory are different from those in GR, and isospectrality is violated in this example of f (R). But in the corresponding STT the Zerilli QNM is exactly like GR and the iso-spectrality stays unbroken. This result indicates some interesting features of the correspondence between f (R) theory and STT. We have shown explicitly that the iso-spectrality breaking is not present in the QNMs after one transforms to the Einstein frame using a conformal transformation. A particular metric solution of an f (R) theory gets mapped to some metric solution of the conformally related scalar-tensor theory. This by itself does not, however, imply that the two spacetime solutions are identical; therefore, it is not necessary for them to share all of the same characteristics. Indeed, in general, a particular solution in f (R) will not get mapped to exactly the same solution in the conformally related scalar-tensor theory. This is often the case when the conformal factor is not unity. In that case there is no obvious reason for their QNMs to share similar properties. What we have shown here is that when an f (R) there and an STT both have the Schwarzschild metric as solution and conformal factor relating these respective solutions is unity, then isopectrality is not broken for the Schwarzschild QNMs in STT even though it is broken for the Schwarzschild QNMs in f (R).

Discussion
In this work we have focused on the QNMs of a Schwarzschild black hole in a general f (R) theory that allows the R = 0 solution. We came to the conclusion that in general f (R) theories the iso-spectrality between the scalar and the vector mode will not be valid. Breaking of iso-spectrality was discussed in the work [53] for the specific case of α R 2 theory. They claimed this result will be true for all other theories also. We found that this claim is not entirely true. In general, there will be breaking of isospectrality, but it is possible to have theories where it is not true, such as when f (R) = α R 3 . Indeed, it is possible to conclude that the isospectrality will be violated if f /(1 + f ) = 0.
Secondly, we discussed that the structure of QNMs gets modified at the asymptote. This modification arises as an inhomogeneous term alongside the GR contribution. Therefore, not only the iso-spectrality will be violated but also the QNM structure will get modified. But for f /(1 + f ) = 0, this contribution along with the breaking of the isospectrality will be absent, resulting in zero deviation from GR in first order perturbation. This does not imply that the ringdown spectrum will be exactly like GR. If we take second order perturbation then we find, ForR = 0 solution S reduces to, As it depends on f , the ringdown mode will contain information of the theory but the first order QNM will be exactly like GR. Therefore a detailed numerical relativity evolution will capture these deviations. As a result one can in principle constrain β and d n f (R) d R n R =0 solely from GW observations.
One key-point is worth mentioning here. In our current work we have focused only on the perturbation around a vac-uumR = 0 solution. We have used these similar assumptions while interpreting the solar system test. Owing to this choice we have not been able to pursue the study of Chameleon screening in f (R) theories [82]. This screening behaviour is the origin of the different behaviour of the theory at different scales. Due to the screening effect the mass of the massive scalar mode depends on the mass density of energy momentum tensor (T μν ) in the space-time. Effective mass of the massive scalar mode becomes heavier in the region with high mass density and lighter in the low density region. Therefore, in different environments the mass of the massive scalar mode will have different values depending on the density profile. This gives rise to the screening mechanism. Therefore, interpreting and using the solar system constraints should be done carefully. But unfortunately due to vacuum assumption (T μν = 0), we have not been able to explore this sector in our current work. Hence, a rigorous analysis with non-vacuum condition is an important endeavour that may be taken up in the future.
QNMs in scalar-tensor theories have recently been calculated in Refs. [55,56,[79][80][81]. It has been shown before that a class of scalar-tensor theories are equivalent to f (R) theories. Tattersall et al. have also found results similar the ones noted here in the context of scalar-tensor theories. Due to the equivalence of scalar-tensor theories and f (R) theories this is, of course, expected. The equivalence between Brans-Dicke with ω 0 = 0 and f (R) theory is valid under the condition f (R) = 0 (related to the sufficient condition for invertibility between the two theories). When f is not defined, or vanishes, the equality φ = f (R) and the equivalence between the two theories cannot be guaranteed, although this is not a priori excluded by f = 0 [73]. Therefore, there are families of f (R) that in principle are not equivalent to scalar-tensor theories. Interestingly, in our work we see that this is precisely the sector where f (R) does not bring any modification w.r.t. GR, mainly due to the vanishing coupling between the GR mode and the extra massive scalar mode.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permit-ted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Spherical harmonics
Any field on a spherically symmetric background can be classified as a scalar, vector and tensor depending on their transformation under the rotation group. For this reason three kinds of spherical harmonics are used to separate out the angular dependence.
One subset of these functions are the scalar spherical harmonics. Any scalar function on a spherically symmetric spacetime is expressible as a sum of the scalar spherical harmonics, Y lm (θ, φ). We call them S(Z A ).
A vector field U A can be expressed as, U A = V A + D A S, where S is the scalar spherical harmonic defined above. V A are called the vector spherical harmonics.
A tensor field X AB on spherically symmetric spacetime can be expressed as, where γ AB is the metric on the two sphere.
where k 2 = l(l + 1), and l is a non-negative integer. From the pure spherical harmonics several useful vector and tensor fields can be defined. We only mention a few among them, namely, (A3)

Appendix B: Separation of effective energy momentum tensor
Let us assume that the metric of a manifold is as follows, ds 2 = g ab (x)dx a dx b + r 2 (x)d 2 , In this case the Christoffel symbols take particular forms as discussed below. A hat represents connection on two sphere and 2 a bc represents the connection in the orbit space. The nonzero components of the connection are as follows, (B2) Using these results for connection we can calculate various covariant derivatives. Various components of double covariant derivatives of a scalar function F(x a , z A ) can be calculated to be as follows, It has been shown that the Einstein equation gets modified in our case. As result we get an effective energy momentum tensor of the following form, f (R) theory of gravity. To study the propagation of GW