Greybody factor and quasinormal modes of Regular Black Holes

In this work, we investigate the quasinormal frequencies of a class of regular black hole solutions which generalize Bardeen and Hayward spacetimes. In particular, we analyze scalar, vector and gravitational perturbations of the black hole both with the semianalytic WKB method and their time-domain profiles. We analyze in detail the behavior of the spectrum depending on the parameter $p/q$ of the black hole, the quantum number of angular momentum and the $s$ number. In addition, we compare our results with the classical solution valid for $p = q = 1$.


Introduction
One of the most striking predictions of General Relativity (GR) [1] is the existence of Black Holes (BHs), objects which produce a region where not even light can escape. They can be formed in the extreme final stages of the gravitational collapse of stars. Such astrophysical objects are remarkably simple, being characterized by three parameters: mass, charge and angular momentum, by virtue of the so-called "no-hair" theorems [2][3][4]. Since gravitational radiation emitted by an oscillating black hole can carry information about its inner properties like mass and charge [5,6], this enables to use BHs as a laboratory for studying gravity in strong regimes, where quantum phenomena might take uttermost importance. Also, after Hawking's seminal papers where the radiation from the black hole horizon was explained [5,6], black boles become an excellent laboratory to study and also enhance our comprehension about quantum gravity. Naively, the Hawking radiation is taken as black body radiation a e-mail: victor_santos@fisica.ufc.br b e-mail: angel.rincon@pucv.cl arXiv:2009.04386v1 [gr-qc] 9 Sep 2020 parametrized by the hawking temperature T H . However, the latter is just an approximated picture because emitted particles feel an effective potential barrier in the exterior region. Such barrier backscatters a percentage of the outgoing radiation back into the black hole [7]. Thus, the spectrum of the Hawking radiation as seen by an asymptotic observer has not a complete blackbody distribution: it is better described by a greybody distribution. The greybody factor can obtained from the transmission amplitude as the field modes pass from near horizon region to an asymptotic observer through the effective potential induced by the spacetime geometry. Estimation of this greybody factor is usually a difficult task and often one has to resort to approximations, usually in low/high frequency limits. There are monodromy methods [8,9] and computations in a variety of scenarios [10][11][12][13], where one can also employ numerical approaches to estimate them [14]. In particular, it should be noticed that a considerable part of the literature is dedicated to the case of black holes with singularities. In order to fix that problem, it is expected that quantum ingredients play a dominant role. Thus, although a complete theory of quantum gravity is still under construction, regular black hole solutions can be obtained utilizing additional matter, for instance, taking advantage of non-linear electrodynamics, as was previously pointed out by Ayon-Beato and Garcia in [15]. Additionally, an interesting feature of the gravitational collapse is that during the formation phase there is the emission of gravitational radiation, strong enough to be detected by the gravitational wave detectors [16]. This connection was strongly provided by the historical direct detection of GWs by LIGO, that two binary merger objects can coalesce into a super-massive black hole [17]. Moreover, recently was reported the first image of the super-massive black hole at the center of the giant elliptical galaxy Messier 87 (M87) by the Event Horizon Telescope [18,19] which increased the interest on the physics behind this class of intriguing objects.
In the framework of GR, gravitational radiation arises as perturbations of spacetime itself. Due to the nonlinear nature of Einsteins equations, it is often very hard to find closed solutions, and usually one has to resort to perturbation theory. such perturbations give rise to a set of damped vibrations called quasinormal modes [20], complex numbers whose real part represents the actual frequency of the oscillation and the imaginary part represents the damping. A simple example of such oscillations are the oscillations of stars which are damped by internal friction [21]. It is known that in general relativity damping might occur even in frictionless scenarios. This effect arises because energy may be radiated away towards infinity by gravitational waves. As it was previously pointed in Ref. [22], even linearized perturbations of black holes exhibit quasinormal modes. QNMs are also important in the investigation of the post-merger remnant of a binary black hole (BBH), as it coalescence settles to a Kerr black hole at a sufficiently late time after the merging phase. The stationary state is reached when the perturbed BH remnant emits gravitational waves (GWs) during a process known as the ringdown (RD), and The GW corresponding to this phase is described by the linear perturbation theory of a Kerr BH [23,24].
As was previously said, the study of the quasinormal frequencies is quite relevant, because they encode information on how a black hole relaxes after it has been perturbed. Even more, such frequencies depend on i) the geometry and ii) the type of perturbations [25]. Along the years, the study of QNM becomes more essential than ever, and certain seminal works have been performed up to now, for instance see [26-28] and more recent works [29][30][31][32][33][34][35][36][37][38][39][40][41][42][43].
One issue of the classical description of gravity is the existence of spacetime singularities [44]. Although the known solutions of Schwarzschild, Kerr and Reissner-Nordstrom have singularities protected by an event horizon [45], the prediction that BHs emit radiation cause them to shrink until the singularity if reached. This is still a conundrum, as singular-ities are generally regarded as an indicative of the breakdown in the theory, requiring new Physics for a proper description. It is a common belief that only a consistent quantum theory of gravity could solve it properly [46]. Since there are no consistent models yet, phenomenological models have been proposed, and most of them are based on the avoidance of the central singularity. Such non-singular solutions are called regular black holes [25, 47,48].
Given the increasing interest in gravitational wave astronomy and on QNMs of black holes, it would be interesting to investigate the QNM spectra from regular BHs. In previous works quasinormal modes of regular black holes were computed by several authors, see e.g. [49][50][51][52]. In this paper we propose to investigate the QNMs of a class of regular black holes which generalize Bardeen and Hayward spacetimes. In particular, we analyze scalar, vector and gravitational perturbations of the black hole both with the semianalytic WKB method and compute the greybody factor. Furthermore, it is essential to note that there is a vast collection of works where QNMs are calculated using the WKB approximation (see [31,35,[53][54][55][56] and references therein).
The work is structured as follows. In section 2 we introduce the spacetime background of regular black holes and how to describe the field perturbations. In section 3 we compute the quasinormal frequencies emplying the 6th order WKB approximation, and in section 4 we compute the greybody factor. Section 5 is devoted to the conclusions.

Regular Black Holes
The class of regular BH solutions of interest in this work have the line element of a static, spherically symmetric black hole where dΩ 2 is the metric of the 2-sphere and the lapse function f (r) = 1 −2 m(r)/r explicitly depends on the matter distribution. The particular mass function [57] m(r) = M 0 ensures an asymptotically flat spacetime for positive p and q. M 0 and r 0 can be interpreted as mass and length parameters. The regular BH solutions proposed by Bardeen [58] and Hayward [47] correspond to the choices (p, q) = (3, 2) and (p, q) = (3, 3) in equation (2). The limit of large r is the only restriction on (3) to obtain an asymptotically Schwarzschild solution is to have q > 0. To compare our results with results found for Bardeen and Hayward solutions, we will restrict ourselves to the case p = 3 and q > 0. The behaviour of some typical cases is shown in Fig. (1). Thus, should be noticed that when the parameter q increases the mass function m(r) tends to its constant value M 0 .  (2) with p = 3 and q > 0. The parameter r 0 is typically assumed to be microscopic (r 0 M 0 ) and hence the exterior region of the BH can be very close to the Schwarzschild space-time. The q → ∞ case, as evidenced by the case q = 90, corresponds to the usual matching between de Sitter and Schwarzchild solutions in the interior region of the BH. The grey shaded region indicates the interior region r < 2M 0 . In this graphic we have used r 0 /M 0 = 10 −1 .

Field perturbations
This subsection is devoted to introduce the basic formalism as well as the theoretical grounds regarding perturbations in black holes in four dimensional spacetime. Thus, we will focus on scalar, electromagnetic and gravitational field perturbations in the fixed background given by Eq. (1) and mass function (2). Due to spherical symmetry, the fields can be decomposed into spherical harmonics, and the equations of motion can be reduced to the form where the tortoise coordinate r * is defined according to the differential equation and the effective potential is where denotes the multipole number of the spherical harmonics decomposition and s = 0, 1, 2 is the spin of the perturbation. Its behavior is shown in Fig. 2. From the behaviour showed in 2 we can see that in the exterior region r > 2M 0 the effective potential does not varies sensibly for q > 1. Therefore, one should not expect large differences in the perturbation oscillation frequencies in the limit of large q. This can be further inferred by observing the behavior of the effective potential with respect to the tortoise coordinate (5), as shown in Fig. 3. Employing the stationary ansatz Ψ (t, r * ) ∼ e iωt ψ(r * ), equation (6) becomes which precisely takes the form of a Schrödinger-like equation. The frequencies of the temporal decomposition take the form ω = ω R + iω I , where ω R is the oscillation frequency of the QNM and ω I is the damping time. Therefore, any mode with ω I < 0 is unstable.

QNMs of regular BHs in the WKB approximation
The QNMs can be computed by imposing proper boundary conditions on (7), where the fields are purely ingoing at the BH horizon and purely outgoing at the spatial infinity. With such boundary conditions, the resulting frequencies are complex. If the potential (6) is peaked and falls to a constant in the asymptotic region, one can compute the QNM frequencies from (7) employing the WKB approximation described by Schutz, Iyer and Will [59][60][61] and posteriorly improved by Konoplya [62]. With this last improvement the QNMs can be computed by where the correction term Λ i can be obtained for different orders of approximation. n is the overtone number and Q (i) 0 is the i-th derivative of Q = ω 2 − V computed at the maximum of the potential. It is worth to mention that the accuracy of the WKB method is dependent of the multipole number and overtone: in general the approximation is appropriate for > n and it is not applicable for < n.
The results for the quasinormal modes are presented in tables 1, 2 and 3. We can notice that all QNMs possess a negative imaginary part, conferring therefore stability for the black hole. Also, both real and imaginary parts of the QNMs do not seem to be wildly sensible to the parameter q in the limit q → ∞, indicating that the damping time is, with good approximation, independent of it. However, we can observe a change the case q ≤ 1, when the BH is smaller damping time than the classical case.

Absorption cross section
Another interesting aspect for investigate the field perturbations around a black hole spacetime is the amount of plunging field which is absorbed by the black hole, the absorption cross section. It embodies the likelihood of a particle to be scattered/deflected by the black hole. Based on the quantum mechanics analogy, the cross-section can be written as where we defined the partial absorption cross section σ ( ) where T ,s (ω) is the transmission coefficient. To compute the transmission coefficient we employ the transfer matrix method, as described in [63]. The partial absorption cross sections versus M 0 ω for = 0 to = 3 are plotted in Fig. 7. One can notice that the = 0 contribution is the largest, responsible for the non-vanishing of the cross section in the small energy limit. We observe how σ /M 2 0 decreases when the value of increases, due to the fact that the effective potential peak increases for larger values of .

Conclusion
In the present paper we have investigated the quasinormal frequencies and the greybody factors of gravitational perturbations around a class of regular black hole solutions. To obtain the corresponding quasinormal frequencies, we perturbed the background and taking advantage of the symmetry we are able to decompose the original function into spherical harmonics, writing down the differential equation of the problem involved. Then we consider the sixth order WKB approximation to obtain the QN frequencies. Further, as can be observed in tables, our solutions reveal that the family of black holes here analysed recover the classical behaviour for q → ∞, and possess a stronger absorption section in the microscopic regime. That means that BH regularity amends the dissipative effect of the black hole on its neighborhood. Finally, as can be observed in tables, for q > 1, the imaginary part of the quasinormal frequencies are also negative. Thus, our results indicate that all modes are found to be unstable.

Acknowledgments
The author V. S. would like to thank the Fundação Cearense de apoio ao Desenvolvimento Científico e Tecnológico (FUNCAP) for financial support. The author A. R. acknowledges DI-VRIEA for financial support through Proyecto Postdoctorado 2019 VRIEA-PUCV.