Quasinormal modes of a massive scalar field nonminimally coupled to gravity in the spacetime of Self-Dual Black Hole

In this work, we investigate the quasinormal modes for a massive scalar field with a nonminimal coupling with gravity in the spacetime of a loop quantum black hole, known as the Self-Dual Black Hole. In this way, we have calculated the characteristic frequencies using the 3rd order WKB approach, where we can verify a strong dependence with the mass of scalar field, the parameter of nonminimal coupling with gravity, and parameters of the Loop Quantum Gravity. From our results, we can check that the Self-Dual Black Hole is stable under the scalar perturbations when assuming small values for the parameters. Also, such results tell us that the quasinormal modes assume different values for the cases where the mass of field is null and the nonminimal coupling assumes $\xi=0$ and $\xi=1/6$, i.e., a possible breaking of the conformal invariance can be seen in the context of loop quantum black holes.


I. INTRODUCTION
General Relativity (GR) is accepted as the best description of gravitation physics. One of the most striking predictions of Einstein's theory is the prediction of the black holes (BHs), which are objects from which nothing (even light signals) can escape after crossing the event horizon.
The interest in BHs goes beyond astrophysics because they have been pointed out as possible objects that can help us to understand one of the most intriguing problems in theoretical physics nowadays: the conciliation of quantum physics and gravitation through a quantum gravity theory (QG). Because it is expected that in the presence of a very strong gravitational field, the quantum nature of spacetime becomes revealed.
One of the main candidates for a theory of quantum gravity is the Loop Quantum Gravity (LQG) [2], which competes with string theory [3]. In the context of LQG, it is possible to get interesting theoretical models that provide an insight into the quantum characteristics of spacetimes revealed by BHs. An important scenario corresponds to the quantum version of the Schwarzschild black hole (SchBH), which is called self-dual black hole (SDBH) [4]. The SDBH solution has a very interesting property, self-duality. Wherefrom such property, the physical singularity of BHs can be replaced by an asymptotically flat region, which is an expected effect in a regime of QG.
In the last years, it has increased the interest in BHs physics because of gravitational waves (GWs) observations, originated from a binary BHs merging and neutron stars. This class of events has been observed with great precision through the LIGO and Virgo collaborations [5].
These GWs are due to the perturbations of BHs, once we cannot find BHs completely isolated in the Universe. Thus, the BHs are always in perturbed states due to the interaction with other compact objects in their neighborhood, for instance, other BHs or neutron stars. Also, the perturbations are characterized by a set of complex eigenvalues (frequencies) of the wave equations, called quasinormal modes (QNMs), and, therefore, they can be observed through experiments with gravitational interferometers.
As aforementioned the QNMs are complex values, where the real parts give the oscillation frequencies, while the imaginary parts determine the damping rates. The QNMs of BHs depend only on BHs parameters and not on how they were perturbed. Thus, the QNMs are known, as the "fingerprint" of a BH. The studies of the QNMs are of great interest and importance in different contexts [6,43].
In recent years, it has been suggested that the QNMs can play a very important role in understanding quantum aspects of gravity theories. Especially, in the LQG context, it has been also suggested that the QNMs can be used to fix the value of the Immirzi parameter, a parameter measuring the quantum of the spacetime [8], a fundamental issue that remains open in this field.
In this context, the QNMs have been analyzed for the SDBH solution in [9], where the authors have considered the perturbations of the massless scalar field. Also, considering gravitational perturbations the analyses have been made in [10,11], where the results showed for the first time in the context of BHs in LQG, a breaking of the isospectrality. Still, in the context of LQG, it is made an analysis of the QNMs for a rotating black hole solution in [12].
The present work has the main goal to perform an analysis of the QNMs spectrum for the metric of SDBH. Thus, we will consider the perturbations through a massive scalar field coupled nonminimally with gravity [13,14]. Some efforts in that sense were considered [15][16][17][18][19].
For instance, in [17] it was analyzed the long-lived quasinormal modes and the instability in a Reissner-Nordström spacetime. The motivations for considering the nonminimal coupling with gravity can be found in different contexts: cosmology [20][21][22][23][24], general relativity [25][26][27], superstring theory [28], induced gravity [29] and loop quantum gravity [30]. To compute the QNMs we will use the WKB (Wentzel -Kramers -Brillouin) approach, which has been introduced by Schutz and Will [31], and after improved by Iyer and Will [32] and more recently by Konoplya [33,36], Matyjasek and Opala [34] and Hatsuda [35]. In this context, we will study the influence of the parameters of LQG, the mass of the scalar field, and the nonminimal coupling with gravity in the QNMs, and consequently in the stability of the SDBH solution.
The paper is organized as follows. In Section II, we briefly review the SDBH solution and discuss their self-duality property. In Section III, we derive a Schrödinger-like equation, where is considered a massive scalar field nonminimally coupled with gravity. In Section IV, we calculate the QNMs through the WKB method, and finally, we summarize our results and draw concluding remarks in Section V. Throughout this work, we use natural units = c = G = 1 and metric signature (−, +, +, +).

II. SELF-DUAL BLACK HOLE
In this section, we will briefly introduce the SDBH solution that arises from a simplified model of LQG. This solution is consisting of an asymmetry-reduced model corresponding to homogeneous spacetimes, and for a more detailed study of the model, see Ref. [4].
The structure of SDBH corresponds to a quantum version of the SchBH and is described by the metric In the Eq. (II.1), we have the presence of an external horizon localized in r + = 2m, an intermediate in r * = √ r + r − and a Cauchy horizon localized at r − = 2mP 2 . Here, the polymeric function P is given by where = γδ b , γ is the Barbero-Immirzi parameter, and δ b is polymeric parameter used for the quantization in LQG. Also, in the Eq. (II.1) appears the parameter a 0 defined by with A min being the minimal area in the context of LQG.
It is important to notice, that the Eq. (II.1) is written in terms of the SDBH mass m that is associated with the ADM mass as follows Also, using the multiplicative factor of the angular part of Eq. (II.1), we can define a new radial Here, R measures the circumference distance and is equal to the r coordinate only in the asymptotic limit. Furthermore, from the Eq. (II.5), we can see an important characteristic of the internal structure of the SDBH. When r decreases from infinity to zero, the R coordinate decreases from infinity to R = √ 2a 0 in r = √ a 0 , and then increases again to infinity. Considering the Eq. (II.5) in the external event horizon, i.e., in r = r + , we get A very interesting characteristic of this scenario is the self-duality of the metric in Eq. (II.1).
The self-duality means that, if we introduce new coordinates,r = a 0 /r andt = tr 2 * /a 0 , the form of metric is preserved. The dual radial coordinate is given byr = √ a 0 and corresponds to a minimal element of surface. Furthermore, the Eq. (II.5) can be written in the form R = √ r 2 +r 2 that clearly shows an asymptotically flat space, that is, a Schwarzschild region in the place of singularity in the limit as r tends to zero. This region corresponds to a wormhole with the size of the order of the Planck length. The Carter-Penrose diagram for the SDBH is shown in Fig. 1.

III. SCALAR PERTURBATIONS AND EFFECTIVE POTENTIAL
In this section, we will study a simple way to introduce the black holes perturbations, called scalar perturbations. Here, we shall consider a massive scalar field coupled to the Ricci scalar that is associated with background geometry.
The scalar field dynamics is described by the modified Klein-Gordon equation, given by Here, µ is the mass of the scalar field Φ and R is the Ricci scalar. The nonminimal coupling, ξ, have a similar form to the models used in cosmological context [37][38][39][40]. A conformal or Weyl transformation is the rescaling of the metric and the scalar field given bỹ where Ω(x) is a real, continuous, finite and non-vanishing function. In special, two values of ξ are interesting: the ξ = 0 called the minimally coupled case, and the conformally coupled case, ξ = 1/6. It is possible to show that, if µ = 0 and ξ = 1/6 the field equation given by Eq. (III.1) is invariant under conformal transformations (conformal symmetry) of Eq. (III.2) [41]. Also, in the Schwarzschild limit (P = 0 and a 0 = 0) the Ricci scalar vanishes, and, consequently, we have the minimally coupled scalar field propagating in the Schwarzschild background.
So, we will start by developing the Eq. (III.1) considering the background metric given by where Ricci scalar, is given by Now, we should introduce the scalar field through standard Ansatz in the following form: Here, the Y m l (θ, φ) denotes the spherical harmonics functions [42] that satisfy the relation: In this point, we can redefine the function ψ as: where the ω parameter denotes the quasinormal modes. Also, we should introduce the tortoise coordinate, x, through the relation: that performing the integration of Eq. (III.9), we get x =r − a 2 0 rr − r + + a 2 0 (r − + r + ) log(r) .
(III.10) like wave equation for the scalar perturbation of SDBH spacetime in the following form: where the effective potential is given by: a 2 0 + r 4 4 (r + r * ) 2 l(l + 1)r 2 a 2 0 + r 4 2 (r + r * ) 4 + a 2 0 r 6 4ξ −r 2 * 3r 2 + 3r − 2ξr 10 r * r 2 (3(r − + r + ) + 4r * ) + 2r(r − (r * − 2r + ) + r * (r + + 2r * )) − 3r − r + r * + r 3 * + µ 2 r 12 (r + r * ) 4 . i.e. r −→ ∞, the effective potential exhibits the following behavior V eff (r) ∼ µ 2 . Also, it is observed that for small masses of the scalar field, the effective potential still has the form of barrier potential. However, by increasing values of µ, the form of the effective potential is turned a potential barrier into a step potential. In Fig. 3, we can observe that for large couplings values (orange curve) we have an instability phase, where the effective potential is negative. Also, the limit case when ξ −→ 0 denoted by the black curve, is nearly coincident with the conformal gravity value (ξ = 1/6). Due to these facts, we should consider parameters with small values so we can apply the WKB approximation correctly for calculation of QNMs.

IV. QUASINORMAL MODES
In this section, we will focus on the computation of the QNMs for the quantum black hole described through metric of the Eq. (II.1). As we saw in the last section, after considering the scalar perturbations to SDBH spacetime we can find a Schrödinger-like equation given by Eq. As earlier mentioned, the QNMs are complex frequencies and can be expressed in the following form: where the real part (ω R ) determines the normal frequency of the oscillations, while the imaginary part (ω I ) represents the damping time of vibration modes. Also, we can get information about Due to behavior of V eff (x) showed in Fig. 4, we can make a direct analogy with the problem of scattering near the peak of the barrier potential of quantum mechanics, where ω 2 in Eq. (III.11) plays the role of the energy. Several methods to compute the QNMs have been developed [43][44][45][46].
However, we chose to apply an approximated method, which is the well-known WKB approach introduced by Schutz and Will [31]. This treatment was later improved to the 3rd order by Iyer and Will [32], and is applied in our calculus. Thus, the QNMs (ω = ω n ) that appear in Eq.
(III.11) are determined by the following equation: Here, we have α = n + 1/2 and V (n) 0 denotes the n-order derivative of the effective potential on the maximum point x 0 . Other developments recent have extended the WKB method of 3th to 6th order [33], well as to higher-orders [34,36]. However, these approaches of higher orders should reproduce results very nearly of results reproduced in the third order, in cases that n l.
So, as we have interesting in small overtones number, n, the 3rd order WKB method has good accuracy to QNMs.
Thus, using the Eq. (III.12) and Eq. (IV.2), we can calculate the QNMs for the self-dual black hole considering different values of the parameters. In Tab. I, we show the fundamental QNMs ( i.e., n = 0) for the case l = 0 and considering different values of the other parameters.
We can note that the QNMs values considering µ = 0, for the case ξ = 0 are different from the conformal coupling case (ξ = 1/6). Furthermore, for the case when polymeric P and minimal area a 0 parameters tend to zero, the results converge to SchBH [47,48]. Also, the QNMs for the SDBH considering the first overtone numbers (n) are shown in Tab. II for the case with l = 1.
We have adopted some typical values to the other parameters, for instance, P = 0.1, 0.3 and 0.5.
However, for the case ξ = 0 and µ = 0 the results agree with the values obtained in [9]. Finally, for better visualization of the effects due to the new parameters in the QNMs spectrum, we have shown the behavior by graphs in Figs. 6, 7 and 8, where has been plotted the real and imaginary parts of ω for the case l = 2.  In plot (a) is shown the real part, while the imaginary part is shown in (b).   So, analyzing our results, we can verify that the scalar quasinormal modes depend strongly on the field mass, the parameter associated with the nonminimal coupling with gravity, and the loop quantum gravity parameters. Our results show that as the parameter P grows, the real part of the QNMs suffers an initial increase and then starts to decrease, while the magnitude of the imaginary part decreases, considering the fixed-parameter a 0 . This behavior is also verified in [9][10][11]. This characteristic reveals that the damping of massive scalar perturbations coupling nonminimally with gravity for the self-dual black hole is slower and the oscillations are faster or slower depending on the value of P .
Another new phenomenon found here is that for fixed l, n, and ξ, the real and imaginary parts of the quasinormal frequencies increase and decrease, respectively, as the mass of the scalar field, µ, increases. As we know, the µ has a maximum value, so the QNMs also have a limited value, and due to this fact can be easier to detect by experiments. Also, the introduction of the scalar field mass can help us to understand different aspects of the self-dual black hole, for instance, the aspects of thermodynamics [30].
We have verified, that by increasing the value of the nonminimal coupling with gravity, ξ, the values of QNMs are modified. Also, for large values of ξ, the number of stable modes will significantly decrease. An interesting fact is that for the case when we assume µ = 0 and ξ = 1/6, the obtained results differ from the case ξ = 0. Therefore, the conformal symmetry can be broken in the presence of LQG corrections. Based on these results, we can conclude that the self-dual black hole has a stable behavior under perturbations provided by a massive scalar field with a nonminimal coupling with gravity.
So, the present study together with [9][10][11] can help us to understand the stability of the selfdual black hole, and also opens a discussion about the interesting issue of conformal symmetry violation in the context of loop quantum black holes. Further analysis considering charged and rotating extensions of SDBH, as well as, other black hole solutions in LQG [49][50][51][52][53][54][55][56], must improve our understanding of this issue. Also, is very important to consider future analysis applying the extended WKB method to higher orders to the self-dual black hole, well as future investigations must address the issue of quasinormal modes considering the nonminimal coupling with gravity in the context of more recent BHs in LQG.