Hawking radiation by spherically-symmetric static black holes for all spins: I -- Teukolsky equations and potentials

In the context of the dynamics and stability of black holes in modified theories of gravity, we derive the Teukolsky equations for massless fields of all spins in general spherically-symmetric and static metrics. We then compute the short-ranged potentials associated with the radial dynamics of spin 1 and spin 1/2 fields, thereby completing the existing literature on spin 0 and 2. These potentials are crucial for the computation of Hawking radiation and quasi-normal modes emitted by black holes. In addition to the Schwarzschild metric, we apply these results and give the explicit formulas for the radial potentials in the case of charged (Reissner--Nordstr\"om) black holes, higher-dimensional black holes, and polymerized black holes arising from loop quantum gravity. These results are in particular relevant and applicable to a large class of regular black hole metrics. The phenomenological applications of these formulas will be the subject of a companion paper.


INTRODUCTION
Black holes (BHs) are fascinating astrophysical objects. As the ultimate stage of the gravitational collapse of stars, they probe the limits of general relativity and our understanding of high energy and high density physics. With the recent rise of experimental gravitational wave detection, they have become the natural arena to seek and test modified theories of gravity. For this reason, it is essential to analyze all facets of their phenomenology. At the theoretical level, the study of physical properties of black holes sets them at the interface between general relativity, thermodynamics and quantum theory.
Since Hawking discovered that black holes emit a quasi-thermal radiation [1], and therefore slowly evaporate away, a vast literature has studied the characteristics of this Hawking radiation. Following Hawking's seminal work, Teukolsky, Press, Page, Chandrasekhar, and Detweiler have worked out the equations governing the perturbations of rotating and charged Kerr-Newman BHs for perturbations with spins 0, 1, 2 and 1/2 in general relativity. From this, they have then deduced the resulting rates of emission of Hawking radiation [2][3][4][5][6][7][8][9][10][11] (see [12] for a complete mathematical review). Since then, it has been understood that general relativity is extremely likely to acquire corrections in both the infrared and ultraviolet regimes. These corrections naturally affect black hole physics. On the one hand, in the context of cosmology, general relativity has been challenged by the discovery of dark matter and dark energy and, for instance, the presence of a cosmological constant in the Einstein equations leads to anti-de Sitter types of BH metrics with modified Hawking radiation [13]. On the other hand, the attempts to reconcile general relativity with quantum theory have led to theories of quantum gravity extending general relativity into the deep quantum regime, such as string theory and loop quantum gravity. Although the purpose of such theories is to propose an ultraviolet completion of general relativity, there is a sense in which they naturally lead to observable effects at large scales. For instance string theory leads to extra spatial dimensions [14,15], and loop quantum gravity leads to effective modifications in turn implying the avoidance of cosmological [16] and black hole singularities [17][18][19][20][21]. From this perspective, black holes act as probes and testbeds, translating the deep quantum corrections to the high curvature regime of general relativity within the horizon into semi-classical corrections to black hole properties as seen from outside the horizon.
Following this logic, various black hole solutions to the corrected Einstein equations in these many modified gravity frameworks have been proposed over the last decades, accompanied by the computation of the corresponding Hawking radiation. Recent work includes e.g. massive gravity [22], cubic gravity [23], hairy BHs [24], Einstein-Gauss-Bonnet BHs [25,26], higher-dimensional BHs [14,15], Kerr-Newman massive scalar emission [27], Hayward BHs [28,29], general "regular" BHs [30] and Kerr-Newman-de Sitter BHs [31]. Beside the Hawking radiation, another important near-horizon property of black hole, which is very sensitive to modifications of general relativity and can be measured from the outside, is the detail of quasi-normal modes. They constitute the ringdown signal of a black hole relaxing towards its equilibrium state. This has become especially relevant in view of the recent gravitational wave detections from black hole mergers by LIGO/VIRGO (see [32] and references therein). Indeed, the increasing sensitivity of the gravitational wave detectors promises an access to the fine structure of the quasi-normal modes resulting from black hole mergers. Through this, we aim to push general relativity to its limits of validity. Indeed, there is (justified) hope that the measure of those quasi-normal mode gravitational waves will give access to the precise characteristics of black hole horizons and thus to their correct metric description. Recent work has focused for example on charged Bardeen BHs [33], Gauss-Bonnet BHs [34,35], Palatini gravity [36,37], f (R) gravity [38], Kerr-de Sitter BHs [39], conformal gravity [40,41], higher derivative gravity [42], and so-called polymerized BHs within loop quantum gravity [43][44][45].
The computation of both Hawking radiation and quasi-normal modes is related to the response of black holes to perturbations. Thus, understanding the physically-measurable consequences of modified gravity on the Hawking radiation and quasi-normal modes requires to work out the equations of motion of the various spin perturbations to black hole metrics. This means generalizing the work of [2][3][4][5][6][7][8][9][10][11][12] to all black hole metrics predicted by the various modified gravity theories. In the present paper, we focus on spherically-symmetric static metrics of the form (1.1), and show how the equations of motion can be written in a form similar to the Regge-Wheeler equation for Schwarzschild BHs, i.e. as a one-dimensional Schrödinger-like radial wave equations with a short-ranged potential. This potential depends on the spin of the perturbation field and we give its explicit expression for each spin 0, 1, 2 and 1/2. This derivation already exists in the literature for fields of spins 2 and 0 (see e.g. respectively [44] and [46]), but here we extend it to spins 1 and 1/2. The potential for spin 1 already appears in [30], but however without the intermediate Teukolsky equation leading to the result (instead, the authors use an argument related to conformal invariance). Hawking radiation rates for spin 1/2 fields in the case of polymerized BH metrics have been studied in [47]. As we will show in the forthcoming companion paper devoted to the detailed study of the Hawking spectra, our results differ from that of [47]. We believe that this is due to the fact that this reference does not use the analytic form of the spin 1/2 short-ranged potential. This is the reason for which in the present paper we derive once and for all the analytic form of the Teukolsky equations and short-ranged potentials in the case of the general metrics (1.1). We also give a general derivation of the intermediate Teukolsky equation for spins 0, 1, 2, 1/2 and 3/2 for these generalized metrics. This is given in formula (2.3). This result will be particularly useful in future work when studying the Hawking emission spectra for metrics of the general type (1.1). In particular, the formulas for the Teukolsky equations and the short ranged potentials can be applied to the case of metrics with independent time and radial components (F = G in the metric ansatz (1.1)), which is particularly interesting since these metrics arise as regular black holes in several effective models of quantum gravity [30,44,47].
We note that Kodama and Ishibashi (see e.g. [48]) develop a formalism that goes straight from the metric to the short range potentials, by means of the stress-energy tensor and without computing the intermediate Teukolsky equations. If it were to be extended to spins 1/2 and 3/2 as well as metrics with different time and radial components, this formalism would provide a complementary way of deriving the potentials. We have checked that our results are similar to those of [49], to which [48] then compares (in the case of Schwarzschild-AdS BHs).
The paper is organized as follows. In section 1, we present the equations of motion using either a direct metric development or the Newman-Penrose formalism. Section 2 shows how to separate these equations to extract the onedimensional radial Teukolsky equation for all spins. Section 3 presents the computation of the short-ranged potentials for all spins. In particular, the calculation for spins 1 and 1/2 requires the use of a Chandrasekhar transform. Finally, section 4 is devoted to the study of some examples of potentials for various black hole metrics, and their comparison with the Schwarzschild case. A mode detailed application of the formalism to various metrics will appear in the companion paper [50].

METRIC AND NEWMAN-PENROSE EQUATIONS
We consider spherically-symmetric static metrics, which constitute a subset of Petrov type D metrics. In fourdimensional Boyer-Lindquist coordinates, the general form of such metrics is where dΩ 2 = dθ 2 + sin θ dϕ 2 is the solid angle in spherical coordinates. Within this family of metrics, we further focus on solutions to the Einstein equations which are asymptotically flat. This means that at spatial infinity the functions F , G, and H must satisfy the asymptotic conditions Many usual metrics fall into this category, such as charged BHs, higher-dimensional BHs or effective BH metrics inspired by (loop) quantum gravity. One particular case that will be especially relevant is to which we refer as tr-symmetric (for time-radius symmetric). For instance, charged and higher-dimensional BHs are tr-symmetric.
We now have to describe the dynamics of matter fields in these types of spacetimes. This can be done either by studying the equations of motion written in terms of the metric, or by using the Newman-Penrose formalism. In the following, we will use the most direct method to obtain the results. Starting with the spin 0 case, we consider a massive scalar field φ. In this case, it is easier to write the Proca equation in curved spacetime where m φ is the mass of the field. For the other types of matter fields, the multiplicity of the vector, spinor or tensor components makes it difficult to obtain a single equation of motion when working directly with the metric. A simple and efficient way to bypass this difficulty is to exploit the Newman-Penrose formalism [12,51], which relies on a reformulation of the equations of motion using a null tetrad field. A choice of null tetrad such that g ab = −l a n b − n a l b + m amb +m a m b is given by where m andm are complex conjugate. This tetrad satisfies l · n = −1 and m ·m = 1, while all other scalar products vanish. Introducing e a i = (e a 1 , e a 2 , e a 3 , e a 4 ) = (l a , n a , m a ,m a ), we define the λ-coefficients as These coefficients enter the definition of the so-called Ricci spin (or rotation) coefficients and some specific linear combinations of these Ricci coefficients are then denoted by For the family of metrics (1.1), the only non-vanishing components are real and given by where X ≡ ∂ r X denotes the derivative in the radial direction. In the tr-symmetric case, these spin coefficients are the same as in [14]. We define the covariant derivatives along the four directions of the tetrad (1.5) as (1.10) These derivatives satisfy the general commutation relation where p and q are arbitrary constants. This identity, which is valid for type D metrics (see equation (2.11) of [2]),is pivotal in what follows. In particular, for the family of spherically-symmetric static metrics (1.1) that we focus on, it reduces to We are now equipped with the necessary material to write down the Newman-Penrose equations of motion for fields of various spins. Massless spin 1. For a massless gauge boson, satisfying the Einstein-Maxwell field equations dF = 0 and d * F = 0, the general form of the Newman-Penrose equations is [2,12,14] where the three Maxwell scalars are (1.14) The cancellation of many of the Ricci coefficients for the family of metrics (1.1) allows to write the first and third equations as a coupled system involving φ 0 and φ 1 only, i.e.
These coupled first order equations can then be turn into a pair of decoupled second order differential equations. One applies δ to the first equation, and applies D − 3ρ to the second one. Adding the two resulting equations, and using the identity (1.12) with p = 0 and q = −2, gives a differential equation involving φ 0 only: This is the equation of motion for a massless spin 1 field.
Massless spin 2. For purely gravitational perturbations, which are equivalent to a massless spin 2 graviton field, the general form of the Newman-Penrose equations is [2,12] where the ψ i are the perturbed components of the Weyl tensor (e.g. ψ 0 ≡ −C abcd l a m b l c m d ), ψ • 2 is the only nonvanishing background component, and the tilde on a spin coefficient indicates a perturbed quantity. We now specialize to the family of metrics (1.1). If we remove the vanishing unperturbed spin coefficients, apply the operator δ − 2β to the first equation, and the operator D − 5ρ to the second one, add the two and make use of identity (1.12) with p = 2 and q = −4, we obtain an equation involving solely ψ 0 , with theσψ • 2 andκψ • 2 contributions replaced by ψ 0 ψ • 2 thanks to the third equation. The resulting equation reads where the background ψ • 2 is given by the Ricci identity as [12] ψ where χ i are the two components of the spinor. We now specialize to the metrics (1.1). We remove the vanishing spin coefficients, apply the operator δ − α to the first equation, apply the operator D − 2ρ to the second one, subtract the two and make use of identity (1.12) with p = −1 and q = −1. This produces a decoupled differential equation for χ 0 only: This is the equation of motion for a massless spin 1/2 field. Massless spin 3/2. Finally, the general form of the Newman-Penrose equations for a Rarita-Schwinger massless spin 3/2 field is [52] where H 000 = (δ − 2β −ᾱ +π)ψ 000 − (D − 2 +¯ −ρ)ψ 001 is a combination of the spinor components, and ψ • 2 is the same background component as in (1.19). Specializing to the metric ansatz (1.1), we remove the vanishing spin coefficients, apply the operator δ + α to the first equation, apply the operator D − 4ρ to the second one, subtract the two and use identity (1.12) with p = 1 and q = −3. This leads to an equation on H 0 ≡ H 000 only, which reads as an ansatz for the wavefunctions. Here S s ,m are the spin-s weighted spherical harmonics for angular modes , m, satisfying the equation where the separation constant is λ s ≡ ( + 1) − s(s + 1). In the spin 0 case, S 0 ,m = Y ,m are just the spherical harmonics. As we are here considering metrics with spherical and not axial symmetry, the dependency on the angular momentum projection m factorizes as S s ,m (θ, ϕ) = S s (θ)e imϕ . Expanding with (2.1) the equations of motion obtained above for all spins will now allow us to decouple the angular and radial equations, just like in the Schwarzschild and Kerr cases [2,4]. Furthermore, the time symmetry replaces time derivatives by the energy ω of the field.   [14], which is where in [14] the notation is ∆(r) ≡ r 2 h(r).

SHORT-RANGED POTENTIALS
The next step towards an applicable formulation of the equations of motion, for the computation of both quasi normal modes and Hawking radiation, is to write the Teukolsky equations (2.3) in the form of a Schrödinger wave equation with short-ranged potentials. Even if the equations can in principle be solved in the form (2.3), precise and stable numerical computations require to work with potentials which fall off at least as 1/r 2 at infinity. Furthermore, working with real-valued potentials also constitutes an appreciable bonus. We therefore need to get rid of the first order radial derivatives and of the complex isω terms in equation (2.3), which have a 1/r behaviour at infinity.
For all of this section, it will be convenient to define a generalized tortoise coordinate r * as [44,46] dr In what follows we will give the expressions of the potentials with both the r * and r coordinates, because the first one is more concise and the second one is better suited for numerical calculations. Furthermore, we will also consider the general redefinition of the wave function as where all quantities are functions of r and we keep track of the spin s. Finally, for each spin our goal will be to find a wave function Z s satisfying the general Schrödinger-like equation with spin-dependent potentials V s , and where ∂ * denotes the derivative with respect to the tortoise coordinate r * . Spins 0 and 2 are already treated in the literature, while spins 1/2 and 1 require more work, and in particular the use of the Chandrasekhar transformation. We now study in detail these aspects.

Spins 0 and 2
For the massive spin 0 field, there is no complex term in (A.2), and all the terms are already decreasing faster than 1/r 2 at infinity because of the fall-offs (1.2). Applying the transformations (3.1) and (3.2), we obtain simply a Schrödinger wave equation for Z 0 ≡ Ψ 0 , with a potential given by [46] V 0 r(r * ) = −Gm 2 φ + This is the short-ranged potential for the massive spin 0 field in the metric (1.1). For the massless spin 2 field, reference [44] follows [12]. They consider clever combinations of the metric components and the vanishing of the Ricci tensor components at first order in the perturbation to obtain directly a decoupled radial equation of the Schrödinger-like form, with the potential This is the short-ranged potential for the massless spin 2 field in the metric (1.1).

Spins 1 and 1/2
We now complete the above results, which are already present in the literature, by deriving the short-ranged potentials for spins 1 and 1/2. These represent the main results of this article. In order to do so, we follow the method which has been used by Chandrasekhar and Detweiler to find the short-ranged potentials for the Kerr metric [8][9][10][11], and perform a Chandrasekhar transformation of the radial Teukolsky equations (A.4) and (A.8).
Let us briefly look at the massless spin 1 and spin 1/2 fields separately before going back to general expressions for spin s. For the massless spin 1 field, applying the transformations (3.1) and (3.2) to (A.4) gives For the massless spin 1/2 field, (A.8) becomes The general form of these equations is The only way to suppress the complex term without re-introducing first order derivatives is to change the unknown function Ψ s by a linear combination of itself and its first order derivative. In this context, this is called the Chandrasekhar transformation. In order to achieve this, we first define the intermediate function Y s by This function is such that equation (3.8) can be written in the form with two functions P s and Q s , and the operators with σ ≡ −ω. When written using (3.9), equation (3.8) becomes Comparing this result with (3.10) then reveals that the two new functions are defined by the requirements One can then show that this gives where the expressions for Q 1 and Q 1/2 can explicitly be checked using (3.6) and (3.7). Note that Q s takes a remarkably simple form, as displayed here, in the case of spin 1/2 and 1. Unfortunately, this is not true for spin 2 and 3/2, in which case the explicit expression is actually much more complicated. The solution for α s , however, is valid for all spins.
In order to continue with a lighter notation, from now on we remove the explicit spin label s from all the various functions involved. We simply need to keep in mind that all the functions encountered below depend on the spin s. Now, let us further decompose Y as a linear combination of the function Z satisfying the Schrödinger wave equation (3.3), by writing where on the right-hand side we have two unknown functions f and W . The Schrödinger equation (3.3) takes the form Λ 2 Z = V Z, where V is the short-ranged potential that we are trying to determine for spin 1 and 1/2. Acting on (3.15) with Λ − and using Λ + = Λ − + 2iσ then leads to where on the right-hand side we have introduced two unknown functions β and R. Acting once again with Λ − on both sides gives Next, we can use Λ + = Λ − + 2iσ once again to rewrite equation (3.10) in the form where P is given in equation (3.13). Matching the Z and Λ + Z terms of these two different expansions for Λ − Λ − Y now tells us that we must have 19) in addition to which we should remember that, because of (3.16), we also have the definitions Now, one can check by a direct substitution that the four previous equations lead to the conservation equation which is a generalization of Chandrasekhar's result [8][9][10][11]. We call this constant K, and we will see later on that it simplifies the calculations neatly. We also define the quantity T ≡ W + 2iσ. Using the identity (3.21) to remove an unwanted derivative of the potential V (which would have caused further difficulties), we finally obtain that (3.19) and (3.20) reduce to the following system of four equations: where (3.22c) has been obtained by combining (3.20) and (3.21). This is the system that we have to solve in order to prove that a solution Z satisfying the Schrödinger wave equation in the potential V does indeed exist. This system follows from the form of the Chandrasekhar transformation, and is valid for all spins. 1 Chandrasekhar and Detweiler have solved it for the Kerr metric and for spins 0, 1, 2 and 1/2. We will now solve it in the general case of the metric (1.1) for spin 1 following [8], and for spin 1/2 following [11]. Spin 1 In the case of spin 1 we look for a simple solution, i.e. we suppose that the unknown quantities are linear in σ and of the form A = A 1 + 2iσA 2 , and that the desired potential V is of course independent of σ (together with Q which is the initial potential without the iω part). Looking at the system (3.22) tells us that the only σ 2 term will come from R 2 , meaning that we need to actually choose R 2 = 0. Then, if we do not wish to carry out the integrations, we further assume that ∂ * T 2 = 0. The only remaining term in iσ then come from f 2 and ∂ * β 2 , so we take f 2 = 0 and β 2 to be constant. Indeed as in [8], both R and f are also independent of σ with these hypotheses. We therefore only need to decompose With all these assumptions, the system (3.22) simply becomes Identifying the no-σ and σ terms in (3.24b) then gives us the two equations while doing the same in (3.24c) leads to We can see from the last equation that the numerical value of the constant T 2 can be absorbed in the other unknown quantities, so we set T 2 = 1 and define κ ≡ K 1 + 4σ 2 β 2 . In order to rewrite the system in an elegant way, we now define the function F ≡ α 2 Q = ( + 1) . (3.27) With this, the second equation in (3.25) gives which can then be injected in the first equation of (3.25) to find We can now use (3.26) to eliminate β 1 from all other equations, and (3.28) to eliminate R. We can then write the previous equation as and (3.24c) can be rewritten in the form Substituting the expression (3.30) for T 1 , we finally obtain an identity on F which reads [8] F The goal is now to find a set of constants β 2 , κ, and K 2 compatible with this identity. Since F is a constant given simply by F = ( + 1), this is actually straightforward. We deduce that β 2 = ± ( + 1) and K 2 = 0, while κ is unconstrained. We then choose K 1 = 0 and obtain κ = 4σ 2 β 2 . Finally, since (3.30) does not constrain T 1 , we choose T 1 = 0, which implies β 1 = 0 thanks to (3.26). We have therefore found a consistent set of constants satisfying the assumptions, and all the remaining functions can be analytically computed. At the end, equations (3.24a) and (3.24d) lead to the very simple result This is the short-ranged potential for a massless spin 1 field in the metric (1.1), which is moreover coherent with the potential obtained in [30].

Spin 1/2
In order to study the case of spin 1/2, we first note that the definition of F gives the simple result In spite of this simple form, using the same hypothesis as in the previous subsection, leading to (3.32), we find that the latter has no solution. We therefore need to make fewer assumptions than above. We will in fact follow [11], and go back to the system of equations (3.22). In this system, integrations can be avoided by assuming which in turn implies thatR ≡ α 2 R is a constant. Thus we have the system where, in order to obtain (3.36a), we have used (3.36d). We see that (3.36a) already gives us the potential as a function of β andR. The goal is therefore to determine these functions. For this, we set T = 2iσ by analogy with the final result of [11] and the result of the spin 1 calculation. Equation (3.36b) then becomes and (3.36c) givesR 2 + 4σ 2 (λ 1/2 + 1) = α 2 (4σ 2R + K) . (3.38) In this equation, everything is a constant except α, we can therefore identify separately We have therefore found a set of constants and a function β satisfying the system (3.36). We can finally use (3.36a) to write the potential as This is the short-ranged potential for the massless spin 1/2 field in the metric (1.1).

Summary
In this section we have obtained the short-ranged massless potentials for all spins in elegant forms using the tortoise coordinate r * . This is summarized as where we have defined for conciseness ν 0 ≡ ( + 1) ≡ ν 1 , ν 2 ≡ ( + 1) − 2 and ν 1/2 ≡ ( + 1) + 1/4. These results are surprisingly compact, and extend the existing literature to the case of spin 1 and 1/2 massless fields in the metric (1.1). We can now focus on specific examples of metrics, and see what the resulting potentials look like when compared to the Schwarzschild ones.

SOME EXAMPLES
There are numerous physically-motivated spherically-symmetric and static metrics of the form (1.1). These examples come from both classical general relativity and modified theories of gravity with e.g. quantum gravity corrections. In this section we will study three examples of potentials, for charged BHs, higher-dimensional BHs, and finally so-called polymerized BHs with corrections from loop quantum gravity. We will have a more extensive discussion of the applications of the formalism presented here in the companion paper [50].
There is already important physical information which can be extracted directly from the form of the potentials given below. For example, BHs which have a higher potential barrier than in the Schwarzschild case will have a lower Hawking emission rate. While this cannot directly be seen from the plots below because they are rescaled, it can easily be seen by looking at the analytic form of the potentials.

tr-symmetric case
Charged and higher-dimensional BHs fall within the family of tr-symmetric metrics (1.3). We therefore first give general results about this case. Using a tr-symmetric ansatz, which depends only on a single function h(r), the massless potentials (3.41) become (4.1d) The first three potentials, which are bosonic, can be written as a single master potential We note that the last term in this master potential is absent from equation (6) of [28], but is coherent with our results and with that of [44]. 2 For the Schwarzschild metric, we recall that where r S = 2M ≡ r H is the Schwarzschild radius of the horizon.

Charged black holes
After the Schwarzschild solution (4.3), the simplest tr-symmetric physically-relevant BH which is solution of classical general relativity equations is the charged BH with where r 2 Q ≡ Q 2 and Q < M is the charge of the BH (since we are working with natural units 4πε 0 = 1, the fine structure constant is α em = 1). The exterior horizon is given by For neutral particles (i.e. with no additional coupling between the charge of the BH and that of the particle), the potentials take the form In figure 1 we show these potentials compared to the Schwarzschild ones for the minimum possible angular momenta = s and for r Q = r S /3 (that is to say Q = 2M/3).

Higher-dimensional black holes
Another simple case of tr-symmetric metrics describes (4 + n)-dimensional BHs [14,15]. In this case the geometry is specified by where the horizon radius is 3

CONCLUSION
In this paper we have studied the dynamics of massless fields of all spins in the general spherically-symmetric and static black hole metrics (1.1), deriving a generic one-dimensional radial Teukolsky equation. For the spin 1 and spin 1/2 cases, we have computed the short-ranged potentials, following the transformation developed by Chandrasekhar, thus completing the existing literature on spins 0 and 2. We have applied our general formalism to three examples, namely charged black holes, higher-dimensional black holes and polymerized black holes (coming from effective models of loop quantum gravity), and compared them to the case of Schwarzschild black holes. We have seen that the resulting potentials can largely deviate from the Schwarzschild case, in particular for spin 2 fields. The two main applications of these potentials is the computation of quasi-normal modes and Hawking radiation. The latter will be the subject of the companion paper [50], which will tackle the in-depth computation of Hawking radiation for modified gravity black holes in order to identify critical differences with the emission by a standard Schwarzschild black hole.
Using the ansatz (2.1), the radial part decouples and becomes This is the radial Teukolsky equation for a massless spin 1 field. Massless spin 2. For the massless spin 2 field, equation (1.18) takes the explicit form Using the ansatz (2.1), the radial part decouples and becomes This is the radial Teukolsky equation for a massless spin 2 field. Massless spin 1/2. For the massless spin 1/2 field, equation (1.21) takes the explicit form Using the ansatz (2.1), the radial part decouples and becomes where we have used the ansatz (2.1) to obtain the final expression. One can observe that this equation is also valid for s = 0, as shows a comparison with equation (A.2).