Greybody factor for black holes in dRGT massive gravity

In general relativity, greybody factor is a quantity related to the quantum nature of a black hole. A high value of greybody factor indicates a high probability that Hawking radiation can reach infinity. Although general relativity is correct and has been successful in describing many phenomena, there are some questions that general relativity cannot answer. Therefore, general relativity is often modified to attain answers. One of the modifications is the ‘massive gravity’. The viable model of the massive gravity theory belongs to de Rham, Gabadadze and Tolley (dRGT). In this paper, we calculate the gravitational potential for the de Sitter black hole and for the dRGT black hole. We also derive the rigorous bound on the greybody factor for the de Sitter black hole and the dRGT black hole. It is found that the structure of potentials determines how much the rigorous bound on the greybody factor should be. That is, the higher the potential, the lesser the bound on the greybody factor will be. Moreover, we compare the greybody factor derived from the rigorous bound with the greybody factor derived from the matching technique. The result shows that the rigorous bound is a true lower bound because it is less than the greybody factor obtained from the matching technique.


Introduction
General relativity was formulated in 1915. It offers profound insights into the concept of gravity. General relativity has succeeded in describing many gravitational phenomena such as the gravitational deflection of light, gravitational radiation, the anomalous perihelion of Mercury, and the behavior of black holes.
According to cosmological observations, todays the universe is expanding with acceleration [1,2]. There are a numa e-mail: petarpa.boonserm@gmail.com b e-mail: tritos.ngampitipan@gmail.com c e-mail: pitbaa@gmail.com ber of cosmological models accounting for this current acceleration. Among these models, the simplest one is the Lambda cold dark matter ( CDM) model in which the cosmological constant drives the acceleration of the universe. The model assumes that general relativity is the correct theory of gravity on cosmological scales. Although the CDM model is in excellent agreement with observations, the cosmological constant itself encounters a theoretical problem that the observational value of the cosmological constant is much smaller than the theoretical one. It may be possible that general relativity may not be the most suitable theory to describe the universe on a large scale, at least on a scale larger than the solar system. The dark side of the universe is evidence that some modifications of general relativity at a large scale is required [3].
Such modifications should be correct on the largest scale and then reduce to general relativity on a smaller scale. One of the candidates for such modifications is massive gravity. Massive gravity is a theory of gravity in which a graviton has a mass. The most successful theory of massive gravity is popularly known by the de Rham-Gabadadze-Tolley (dRGT) model [4,5]. Because of the Vainshtein mechanism, adding mass to graviton keeps physics on a small scale equivalent to general relativity, with some small corrections [6]. However, this leads to the modification of gravity on a larger scale.
Regarding the cosmological solutions in the dRGT massive gravity theory, even though all the solutions cannot provide a viable cosmological model, for example, the solutions do not admit flat-FLRW metric [7,8] or the model encounters instabilities [9][10][11], a class of solutions can provide a viable cosmological model [12][13][14]. The solutions for the dRGT massive gravity are not only investigated in a cosmological background, but also in a spherically symmetric background [15]. For spherically symmetric solutions, the black hole solutions have been investigated in both analytical [16][17][18][19][20][21][22][23][24][25] and non-analytical [26,27] forms, depending on the fiducial metric form. However, they still share the same property, which is represented as an asymptotic AdS/dS behavior.
A black hole can emit thermal radiation if the quantum effects are considered. This thermal radiation is known as Hawking radiation [28]. Hawking radiation propagates on spacetime, which is curved by the black hole. The curvature of spacetime acts as a gravitational potential. Therefore, Hawking radiation is scattered from this potential. One part of the Hawking radiation is reflected back into the black hole, while the other part is transmitted to spatial infinity. The transmission probability in this context is also known as the greybody factor.
There are many methods to calculate the greybody factor. For example, one can obtain an approximate greybody factor using the matching technique [29][30][31]. If the gravitational potential is high enough, one can use the WKB approximation to derive the greybody factor [32][33][34]. Other than approximation, the greybody factor can also be obtained using the rigorous bound [35][36][37]. The bound can give a qualitative description of a black hole.
In this work, we investigate the greybody factor using the analytical black hole solution in dRGT massive gravity. In Sect. 2, the structure of the horizons of the solution is analyzed in order to generate a suitable form for the analysis of the properties of the greybody factor. In Sect. 3, the properties of the gravitational potential are investigated for both the de Sitter black hole and dRGT black hole. The height of their potentials are determined by the parameters of the model. In Sect. 4, we derive the rigorous bound on the greybody factor, and the reflection probability for the de Sitter black hole and the dRGT black hole. The value of the rigorous bound on the greybody factor corresponds to the structure of potentials. In addition, the effects of the graviton mass, cosmological constant and angular momentum quantum number on the greybody factors will also be explored. Finally, concluding remarks are provided in Sect. 5.

dRGT black hole background
In this section, we will review the main concept of the dRGT massive gravity theory in the following manner [18]; the analytical solution of the modified Einstein equation due to the graviton mass is first presented, and then the structure of the horizons of the black hole in the dRGT massive gravity theory is investigated. The theory of the dRGT massive gravity is a covariant non-linear theory of massive gravity, which is ghost free in the decoupling limit to all orders. The action of the dRGT massive gravity model in four-dimensional spacetime can be expressed as where R is a Ricci scalar corresponding to a physical metric g μν , m 2 g to the square of the graviton mass, with L i s representing the interactions of the ith order of the massive graviton. In particular, those interactions of the massive graviton are constructed from two kinds of metrics and can be expressed as follows, where the rectangular brackets denote the traces. The tensor K μν is constructed from the physical metric g μν and the fiducial metric f μν as where the square roots of those tensors are defined so that The fiducial metric is chosen as [16,38,39] f μν = diag(0, 0, c 2 , c 2 sin 2 θ), where c is a constant. It can be shown in [38] that the theory of massive gravity with the choice of the singular fiducial metric such as one in equation (6) is absent of Boulware-Deser (BD) ghost. The static and spherically symmetric black hole solution satisfying this theory can be written as [18] where dΩ 2 = dθ 2 + sin 2 θ dφ 2 , and M is an integration constant related to the mass of the black hole. The parameters above can be written in terms of the original parameters as and This solution contains various kinds of black hole solutions found in literature. If m g = 0, the Schwarzschild solution is recovered. In the case of c = 0, the solution reduces to the de Sitter solution for 1+α+β < 0 and reduces to the anti-de Sitter solution for 1+α+β > 0. Moreover, the global monopole solution can be obtained by setting 1 + 2α + 3β = 0. Note that the last term, the constant potential ζ , corresponds to the global monopole term. A global monopole usually comes from a topological defect in high energy physics of the early universe resulting from a gauge-symmetry breaking [40][41][42]. However, in this solution, the global monopole is contributed via the graviton mass. Note that the linear term γ r is a characteristic term of this solution, distinguished from other solutions found in literature. Next, we will consider the structure of the horizons of this solution. Since the solution is an asymptotical AdS/dS solution, we first consider the structure of the AdS/dS solution and then investigate the structure of the horizons of the solution in the dRGT massive gravity theory.

Horizon structure for AdS/dS-like solutions
It is important to note that one can choose c = 0. This corresponds to trivial solutions since the interacting terms (or graviton mass) become constant, which is inferred from K μ ν = δ μ ν . Therefore, the action in Eq. (1) becomes the Einstein-Hilbert action with cosmological constant. In order to investigate the structure of the horizon, let us first consider a simple case where c = 0. As a result, the function in the metric solution becomes From this function, one can see that f (r ) → −∞ where r → 0. In order to have two horizons, f must be increased and then decreased where r is increased. This means that f (r ) → −∞ again when r → ∞. Therefore, in order to obtain two horizons, must be negative. This corresponds to the de Sitter (dS) spacetime, while in the case of anti-de Sitter (AdS) spacetime, > 0, there exists only one horizon. Now let us find the conditions for which there are two horizons for the de Sitter spacetime, where < 0. If two horizons exist, the maximum value of f must be positive. The maximum point of f can be found by solving f = 0. As a result, the maximum point is Substituting this radius into f (r ) in Eq. (13), the maximum value of f can be written as By requiring f (r m ) > 0, the condition for having two horizons can be written as In order to parameterize the solution properly, let us define a dimensionless parameter as where 0 < α 2 m < 1. By using the dimensionless variablē r = r/M, function f can be rewritten as In order to find the horizon, one has to solve the cubic equation; This cubic equation is known as the depressed cubic equation, and the solution can be expressed as where k = 0, 1, 2 for the three distinguished solutions. Since 0 < α m < 1, one can expand the sinusoidal function and then keep only the significant contributions. As a result, for k = 2, r is negative, and for k = 1 and k = 0, the solutions can be respectively approximated as Note that these solutions are well approximated when α m 1. Actually, this approximation can be realized to satisfy the cosmological solution in which the universe expands with acceleration, since the observed value of is very small compared to the black hole mass. The behavior of the horizon with various values of α m is shown in the left panel of Fig.  1.
For the AdS solution, is positive. Therefore, one can find the solution of the horizon by rewriting α m as α m = 9 M 2 . By following the same steps from the de Sitter case, the horizon in the AdS case can be written as As we have discussed above, there exists only one horizon for the AdS solution. This is valid for all values of α m as show in the right panel of Fig. 1.

Horizon structure for the dRGT massive gravity solutions
For the complete massive gravity solution, it is significantly difficult and complicated to find the horizon analytically. One of the conditions for having three horizons is that > 0. Therefore, we can separate our consideration into two classes; the asymptotic AdS solutions for > 0 and the asymptotic de Sitter solution for < 0 [43].
For the general solutions of the dRGT massive gravity, the dimension-length parameter c is not set to be zero. This means that we have to introduce a scale to the theory. It is useful to work out the solution using dimensionless variable, r = r/c, and then find out what scale c would assume. As a result, function f can be written in terms of a dimensionless variable as wherẽ From this equation, it is sufficient to figure out that the scale of c takes place atM ∼ α g . Therefore, one can choose the parameter c as This radius is well known as the Vainshtein radius [44,45]. The theory in which r < r V will approach GR, while the theory in which r > r V , the modification of GR will be active. The horizons can be found by solving for the solution of r through the equation In order to find the conditions for having three horizons for the AdS case and two horizons for the de Sitter case, let us consider the extremum points of the function f using the equation f = 0 or This is the cubic equation. We can solve it by changing the variables to obtain the depressed cubic equation and then analyze the general solution to find the condition for having two real positive roots for the AdS case and one real positive root for the de Sitter case. As a result, one can constrain our consideration to case c 1 = 3(4c 2 2 ) 1/3 . Choosing this condition will guarantee one real positive root for the de Sitter case asr d S = (−2c 2 ) −1/3 and two real positive roots for the AdS case asr Ad S1 = (2c 2 ) −1/3 andr Ad S2 = (1 + √ 3)(2c 2 ) −1/3 . For the asymptotic de Sitter solutions, f at the extremum point can be written as In order to have two horizons, f (r d S ) > 0. We can parameterize the parameter c 0 such that Therefore, the condition for having two horizons in the case of de Sitter-like spacetime is By changing the cubic equation into the depressed cubic equation, one can find the two real positive horizons as the real roots of the depressed cubic equation as follows where Our analysis can be checked using the numerical method as shown in the left panel of Fig. 2. From this figure, one can see that β m can parameterize the existence of two horizons. Now we consider the asymptotic AdS solutions. The f at the extremum points can be written as In order to have three horizons, we must have f (r Ad S1 ) > 0 and f (r Ad S2 ) < 0. By using the parameter of α 0 in Eq. (29), while changing c 2 to −c 2 , the condition for having three horizons can be written as By using the same step as done in the asymptotic de Sitter case, the three real positive horizons for the AdS case can be written as The numerical plot for these horizons is shown in the right panel of Fig. 2. From this figure, one can see that the existence of three horizons satisfy the condition 1 < β m < 2/ √ 3 as we have analyzed. In the next section, we will use the expression for the horizons derived in this section to analyze the properties of the gravitational potential and the greybody factor of the black hole.

Equations of motion of massless scalar field
Classically, nothing can escape a black hole when approaching it. However, when a quantum effect is considered, the black hole can radiate. This radiation is known as Hawking radiation. It is a blackbody spectrum of temperature where r s is the Schwarzschild radius. In this paper, we assume that Hawking radiation is a massless scalar field. The massless scalar field satisfies the Klein-Gordon equation We use the spherical coordinates. The solutions to the wave equation in spherical coordinates are of the form where Y m (Ω) are spherical harmonics. The Klein-Gordon equation becomes The angular part satisfies where is the angular momentum quantum number. Therefore, the Klein-Gordon equation (Eq. 44) is left with the radial part where r * is the tortoise coordinate defined by (47) and V (r ) is the potential given by It can be expressed in terms ofr as Substituting the function f (r ) from Eq. (23), we obtain From the above equation, we can see that the potential is high when the angular momentum quantum number is large. Qualitatively, the leading contribution to the transmission amplitude comes from the mode = 0. Therefore, it is sufficient to qualitatively analyze the potential for the case of = 0. When = 0, the first term in Eq. (49) vanishes and the second term is proportional to f f . Therefore, there are threer -intercepts resulting from f = 0 and f = 0. This behavior significantly differs from the Schwarzschild case, which has only oner -intercept as shown in the right panel of Fig. 3.
For the de Sitter case, the potential can be obtained by setting c 1 = c 0 = 0. As we have analyzed in Sect. 2.2, it is convenient to change parameter as α g c 2 = − α 2 m /(27M 2 ). For this setting, the potential depends only on the parameter α m . In the same strategy as the quantum theory, the shape of the potential controls the transmission amplitude. Therefore, it is worthwhile to consider the maximum value of the de Sitter potential compared to the Schwarzschild potential. As a result, the de Sitter potential for the = 0 case can be written as By solvingr max via V d S = 0 and then substituting the solution back into the above equation, one can find the maximum value of the potential depending on only two parameters,M and α m . The expression is significantly lengthy; we do not present it in the current paper. In order to see the effect of the graviton mass or the cosmological constant, one can fix M and then plot this expression via α m as shown in the left panel of Fig. 3. Note that we also show the result for the = 1 case in this figure. From this figure, one can see that V max contributed from the de Sitter black hole is always less than one from the Schwarzschild black hole. Clearly, the cosmological constant plays a role in reducing the local maximum of the potential. The explicit form of the de Sitter potential is plotted with various values of the cosmological constant as shown in the right panel of Fig. 3. Note that we used the parameter c 2 instead of α 2 m . This is convenient for comparing the results with one from the case of the dRGT massive gravity. Consequently, it might be expected that the transmission amplitude due to the de Sitter black hole should be greater than one in the Schwarzschild black hole. We will clarify this issue explicitly in the next section.
The dRGT potential for = 0 can explicitly be written as By employing the same strategy as used in the de Sitter case, we found that the leading term of V max is proportional to V max ∝ 1/β m . By fixing c 2 = − 1, we have illustrated the explicit behavior of the peak of the potential in the left panel of Fig. 4. We can explicitly see that V max increases as β m decreases. Moreover, the de Sitter potential and the dRGT potential with various values of β m are plotted as shown in the right panel of Fig. 4. It shows that both the de Sitter potential and the dRGT potential increase with radial distance from the black hole. After that, they decrease with radial distance to reach the relative lowest point and then turn to increase again.
In [46], the maximum points of the potentials are not chosen to be equal. We consider this point here. The equality of the maximum points allows us to draw conclusions regarding how high the rigorous bounds on the greybody factors for different types of black holes are. At the highest point, The solution of V (r ) = 0 is not shown here. We find that the equality of the peak of the potentials occur at β m = 0.565375. Moreover, the effects of parameter c 2 are shown in Fig. 5.
To see the effect of the parameter c 2 on the potential, let us fix β m = β mc . The peak of the potential is plotted as shown in the left panel of Fig. 5. The potential is also plotted with various values of c 2 as shown in the right panel of Fig. 5. The parameter c 2 characterizes the strength of the graviton mass. Therefore, the graviton mass will enhance the potential in contrast to the effect of the cosmological constant in the de Sitter black hole.
In this section, we explore the behavior of the gravitational potential for both the de Sitter black hole and the dRGT black hole of a massless scalar field. By making a comparison with the potential in the Schwarzschild black hole, we found that the local maximum of the de Sitter potential is always less than one of the Schwarzschild potential. For the dRGT black hole, the local maximum of the potential depend on the model parameters; β m characterizing the existence of two horizon (0 < β m < 1) and c 2 characterizing the strength of the graviton mass. In contrast to the de Sitter potential, we found that the local maximum of the dRGT potential will be larger than ones for the Schwarzschild and the de Sitter potential by setting parameter β m 1 or c 2 − 1. In the same fashion as quantum theory, the shape of the potential has an effect on the transmission amplitude or the the greybody factor in this context. We will use the information of the potential to analyze the behavior of the greybody factor in the next section.

The rigorous bounds on the greybody factors
From a classical point of view, a black hole is believed to be black because nothing, when having entered the black hole, can escape, not even light. From a quantum point of view, however, a black hole is no longer considered 'black' since it has been proven to emit a type of thermal radiation known as Hawking radiation. At the event horizon of a black hole, Hawking radiation is exactly a blackbody spectrum. While Hawking radiation propagates out from the event horizon, it is, however, modified by the spacetime curvature generated by its black hole source. Thus, an observer at an infinite distance observes the modified form of Hawking radiation, which is different from the original Hawking radiation at the event horizon. This difference can be measured by the so-called greybody factor.
In this section, a greybody factor will be obtained using the rigorous bound [35][36][37]47]. The bound can give a qualitative description of a black hole. It is applied to various types of black holes such as a Schwarzschild black hole [48], a non-rotating black hole [49], a dirty black hole [50], a Kerr-Newman black hole [51], a Myers-Perry black hole [52], and a dRGT black hole [46]. The rigorous bounds on the greybody factors are given by where where h(r * ) is a positive function satisfying h(−∞) = h(∞) = ω. See [35] for more details. We select h = ω. Therefore, The rigorous bounds on the reflection probabilities are given by where, from Eqs. (31) and (32),r H =r d S1 andR H =r d S2 are given bỹ From Eq. (64), we find that the rigorous bound on the greybody factor in massive gravity crucially depends on two parameters, c 1 and c 2 , which determines how the structure of the graviton mass affects the bound. As we have discussed in Sect. 3, the parameter c 1 must be positive in order to have two horizons. Therefore, the last term in Eq. (64) always provides the negative correction to the bound, so that, for the potentials with the same height, the bound from the dRGT black hole is always larger than the bound from the de Sitter black hole. Moreover, this behavior can be qualitatively expressed by analyzing the potential for both cases. From Fig. 4, for the potentials with the same height, the potential from the dRGT case is always thinner than one from the de Sitter case. Therefore, the transmission amplitude for the dRGT case is always greater than one for the de Sitter case as seen in the left panel of Fig. 7. As we have analyzed earlier, the height of the potential can be controlled by two parameters, β m and c 2 . Now let us figure out how the parameters affect the dRGT bound compared to the de Sitter bound. By fixing c 2 , one can see that the bound crucially depends on |R H −r H |, which is proportional to 1/β m . Therefore, one finds that the larger the value of β m , the higher is the value of the bound. This can be seen explicitly by numerically plotting T b with various values of β m as illustrated in the right panel of Fig. 7. Moreover, this behavior can also be seen by analyzing the potential. From Fig. 4, we found that the larger the value of β m , the lower is the peak of the potential. Therefore, one finds that the larger the value of β m , the higher is the value of the bound.
In terms of fixing β m , one can see that the maximum value of the potential decreases when |c 2 | decreases as such that the bound will increase when |c 2 | decreases as shown in the left panel of Fig. 8. Moreover, to compare the bound from the dRGT black hole to one from the de Sitter and the Schwarzschild black hole, we also plot the bound by fixing ω as seen in the right panel of Fig. 8. From this figure, we found that the bound from the dRGT black hole can be larger or smaller than ones from both the de Sitter and the Schwarzschild black holes, depending on c 2 . On the other hand, the bound from the de Sitter black hole is always larger than the bound from the Schwarzschild black hole. Therefore, it is found that there is more room for the dRGT black hole to increase or decrease the greybody factor.

Conclusion
In this paper, we obtain the gravitational potential from Schwarzschild black holes, de Sitter black holes, and dRGT black holes. We also derive the rigorous bound on the grey- Fig. 7 The left panel shows a comparison of the rigorous bound of greybody factor for the dRGT and the de Sitter black holes with = 0,M = α g = 0.1, while the parameter β m is chosen to have the same height of potential. The right panel shows the rigorous bound for c 2 = − 1 body factor for the de Sitter black hole and the dRGT black hole. It is found that the structure of potentials determines how much the rigorous bound on the greybody factor should be. Since V max contributed from a de Sitter black hole is always less than one in a Schwarzschild black hole, the bound for a de Sitter black hole is greater than one for a Schwarzschild black hole. In case of potentials with the same height, the result shows that the bound from a dRGT black hole is always larger than the bound from a de Sitter black hole. Otherwise, the bound from a dRGT black hole can be larger or smaller than ones from both de Sitter and Schwarzschild black holes due to different effects of the parameter c 2 on de Sitter and dRGT spacetimes. Furthermore, we compare the greybody factor derived from the rigorous bound with the greybody factor derived from the matching technique. The results show that the greybody factor obtained from the rigorous bound is less than the one from the matching technique, which means that the rigorous bound is a true lower bound. additional supported by the Naresuan University Research Fund through Grant no. R2559C235 and the ICTP through Grant no. OEA-NET-01 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .