Perturbing microscopic black holes inspired by noncommutativity

We probe into the instabilities of microscopic quantum black holes. For this purpose, we study the quasinormal modes (QNMs) for a massless scalar perturbation of the noncommutative geometry inspired Schwarzschild black hole. By means of a sixth order Wentzel-Kramers-Brillouin (WKB) approximation we show that the widely used WKB method does not converge in the critical cases where instabilities show up at the third order. By employing the inverted potential method, we demonstrate that the instabilities are an artifact of the WKB method. Finally, we discuss the usefulness of the asymptotic iteration method to find the QNMs.


Introduction
The problem of understanding the mechanisms underlying the final evolution of a black hole, which, in turn, is closely related to the problem of the central singularity and the emergence of a possible black hole remnant, triggered several studies in the last two decades [1][2][3][4][5]. Even though we do not have at the moment a final theory of quantum gravity able to offer a definite answer to the aforementioned problem, the existing candidate theories such as string theory, loop quantum gravity, and noncommutative geometry seem to share global features like the noncommutativity of the coordinates at a typical length √ θ less than 10 −18 m [6-9], a new uncertainty principle including gravity effects [10][11][12], the avoidance of physical singularities [13][14][15], and black hole remnants [16,17]. In particular, [18][19][20][21] showed that noncommutativity can provide a cure for the central singularities afflicting the Schwarzschild and Reissner-Nordström metrics. The corresponding noncommutative counterparts of a e-mail: davide.batic@ku.ac.ae b e-mail: nkelkar@uniandes.edu.co c e-mail: mnowakos@uniandes.edu.co d e-mail: kj_redway@hotmail.com these metrics are derived by incorporating noncommutativity only in the matter source while the Einstein action is kept unchanged. As a result of this approach, the central singularity is replaced by a regular region (deSitter core) represented by a self-gravitating, droplet of anisotropic fluid. The price to be paid is that the radial and tangential pressures are always negative, a fact which is difficult to explain. This difficulty was overcome in [22] where by means of a nonlocal equation of state a noncommutative mini black hole model was developed in which the pressure is positive in the interior of the droplet. Furthermore, [23] gave the maximal singularity-free atlas for the noncommutative geometry inspired Schwarzschild metric [18]. This atlas describes an infinite lattice of asymptotically flat universes connected by black hole tunnels. The stability problem of the noncommutative Schwarzschild interior under massless scalar perturbations was attacked by [24,25] leading to conclusions opposite to ours. It is worth mentioning that [26,27] claimed to have derived the noncommutative counterpart of the Kerr metric. Despite the educated guess leading to this result, the authors did not perform any consistency check regarding the question whether or not such a solution satisfies all Einstein field equations coupled to an anisotropic source. In this context, we mention the work of [28] where it was shown that any noncommutative Kerr candidate metric, be it an educated guess or a particular solution, must satisfy a highly complicated system of partial differential equations (PDEs) represented by Eqs. (12)- (19) therein. Moreover, Theorem 1 and 2 in [28] give two independent proofs indicating that the metric found by [26,27] can never satisfy all Einstein field equations and it may lead to contradictions.
It is of general interest to study the stability of microscopic quantum black holes. Moreover, as it will become apparent in the text, there exists indeed a hint of possible instabilities when the third order WKB approximation is used to determine the QNMs. The main question we address here is if these instabilities persist by increasing the order of the WKB or by employing different algorithmic methods.
While there is a vast literature devoted to the study of QNMs for a commutative Schwarzschild spacetime [29][30][31][32][33][34][35][36][37][38], the same cannot be said for the noncommutative Schwarzschild manifold for which only some partial results are available in the literature. For instance, [39] analysed the asymptotic QNMs of a noncommutative geometry inspired Schwarzschild black hole under the assumption that the black hole mass M is much larger than M e , the mass of the extremal noncommutative Schwarzschild solution. At this point a comment is in order. From our Table 2 we can evince that the QNMs of a noncommutative Schwarzschild black hole already coincide with the corresponding ones of its classic counterpart when M is of the same order of magnitude as M e . Therefore, it is not surprising that [39] reached the conclusion that the asymptotic QNMs of a noncommutative Schwarzschild black hole in the regime of large masses remain proportional to ln 3 as is the case for the classic Schwarzschild solution. Recently, [40] studied the QNMs of massless scalar field perturbations in a noncommutativegeometry-inspired Schwarzschild black hole spacetime by means of the third-order WKB approximation. The author considered the cases with = 1, 2, 3 only and computed the QNMs for different values of the noncommutative parameter θ in the interval 0.01 < θ < 0.2758 using the third order WKB approximation. Since the parameter √ θ acts as a quantum of length, it is natural to assume √ θ = P where P = hG/c 3 ≈ 10 −35 m is the Planck length. Then, the interval chosen in [40] over which the noncommutative parameter varies corresponds to a choice of the black hole mass in the interval 3.6M P < M < 19M P with M P the Planck mass. In this particular range studied by the author and within the third order WKB, results for the QNMs were found to be stable. The misleading instabilities occurring in the third order WKB, occur in the range 1.91M P ≤ M ≤ 2.3897M P as can be evinced from our Table 2. Extending the WKB calculations up to the sixth order reveals that the method is not convergent exactly when the presumed instabilities occur. This forces us to turn to another reliable method to pin down the nature of these instabilities. We choose the inverted potential method by which we can assure that the instabilities are due to inefficiency of the WKB in such cases. This case teaches us an important lesson for the determination of QNMs: a single algorithm is not always appropriate to uncover all QNMs. We strengthen this conclusion by investigating the asymptotic iterative method (AIM) and show its advantage over other methods commonly used.
The paper is organised as follows. In Sect. 2, we introduce a suitable rescaling of the noncommutative Schwarzschild manifold and derive the effective potential for a massless scalar field in the aforementioned geometry. Exploiting the short range property of this potential, we obtain a novel inequality linking the event horizon of a noncommutative Schwarzschild black hole with its mass parameter and the angular momentum quantum number of the scalar field. In Sect. 3, we use a third and sixth order WKB approximation to compute the QNMs. Section 4 discusses the inverted potential method followed by Sect. 5 devoted to the asymptotic iteration method.

The noncommutative Schwarzschild black hole
We consider a massless scalar field φ immersed in the noncommutative geometry inspired Schwarzschild background whose line element is [18] with the mass function where M is the total mass in the space-time manifold, θ is a parameter encoding noncommutativity and having the dimension of a length squared, while γ (·, ·) is the incomplete lower gamma function. Taking into account that the incomplete and lower gamma functions are related by the formula [41] the function f (r ) appearing in (1) can be written as the usual Schwarzschild counterpart plus a perturbation due to noncommutativity, i.e.
Observe that the classic Schwarzschild metric can be recovered in the regime r 2 /4θ 1 because in that case (3/2, r 2 /4θ) → 0. Due to the spherical symmetry of the manifold described by (1), the Klein-Gordon equation can be separated into spherical harmonics using the ansatz with = 0, 1, 2, . . .. According to [42], the radial part ψ ω satisfies the following second-order, linear differential equation , where U e f f is the effective potential. If we express the radial variable in units of the black hole mass and the black hole mass in units of the square root of noncommutative parameter, i.e Equation (7) can be cast into the form of an ordinary differential equation (ODE) characterized by two finite singularities at the inner and outer horizon and one singularity at infinity. More precisely, we find , where U e f f is the effective potential, σ = Mω is the dimensionless frequency defined as in [43], and with erf(·) denoting the error function. The above expression for the function f has been obtained by using the following relations in [41], namely (Figs. 1, 2, 3) It can be easily verified that an extreme noncommutative geometry inspired black hole will occur whenever the parameter μ has the critical value μ e = 1.904. In this case, the Cauchy horizon x c , and the event horizon x h coincide with x c = x h = 1.587. For μ < μ e , the line element (1) describes a naked gravitational object regular at x = 0, while there are two distinct horizons if μ > μ e . Furthermore, at short distances, i.e. r √ θ , it turns out that the typical central singularity of a classic Schwarzschild black hole has been replaced by a regular deSitter core [18]. In the present work, we consider both the extreme and non extreme cases. The position of the event horizon x h is obtained by solving the transcendental equation f (x h ) = 0, i.e.
Since the expression of the effective potential has an overall multiplicative term f (x), we expect that U e f f vanishes at the event horizon. Furthermore, in the non extreme case the event horizon is a simple zero of the above equation, while The dashed line patterns represent the asymptotic approximations (19) and (20). Approximation (19)   The plot has been generated by expressing (21) in terms of the variable x with the help of (16). The dotted lines represent the effective potential in the extreme case. Note the change in concavity at the event horizons for the effective potentials in the extreme and nonextreme cases. This behaviour signalizes that the approximation (21) breaks down in the extreme case and it should be replaced by (22) as in Fig. 5. The label on the y-axis refers to the effective potential . The plot has been generated by expressing (22) in terms of the variable x with the help of (16) it becomes a zero of order two in the extreme case [23]. The differential equation (9) can be further simplified by introducing the tortoise coordinate

loses precision as increases
We observe that the function is positive on I . Moreover, with the help of 7.1.23 [41] it is not difficult to verify that asymptotically for x → +∞ the tortoise coordinate exhibits the behaviour Since in the proximity of the event horizon, i.e. for x → x h , the tortoise coordinate will behave there as An inspection of (16) shows that denoting the minimum of f (x) (see dotted line in Fig. 2). Moreover, (16) yields In order to introduce the quasinormal modes boundary conditions, it is convenient to rewrite (9) in terms of the tortoise coordinate. We get As in [42] solutions which are ingoing at the horizon (as Quasi-normal modes are represented by the (complex) frequencies σ n , or equivalently ω n in virtue of the relation σ = Mω, at which A (in) σ = 0. According to [44] the radial functions associated to the quasi-normal modes must diverge exponentially for x * → ±∞ and the wave function φ must exhibit an exponential decrease in the time variable. This requires that Im(σ ) < 0. Note that if we would have taken e iωt in the separation ansatz for the Klein-Gordon equation, then one had to require Im(σ ) > 0 together with the boundary conditions At this point a couple of remarks on the effective potential appearing in (18) are in order. First of all, if we rewrite the lower incomplete Gamma function in terms of the upper incomplete Gamma function with the help of (3), and then, we use (2.24) in [45] which provides an asymptotic expansion for the upper incomplete Gamma function, we discover that the asymptotic behaviour of the effective potential asymptotically away from the black hole is In the case = 0 the corresponding asymptotic expansion of the effective potential is given by In Fig. 2, we compared (19) and (20) with the exact shape of the effective potential. Moreover, when x * → −∞, we find that at the event horizon of a nonextreme Schwarzschild black hole the effective potential can be approximated as follows x h (21) with α given by (16). As discussed above, α is always positive on the interval (x m , ∞), and therefore, β > 0. The comparison of (21) with the exact shape of the effective potential in a neighbourhood of the event horizon can be found in Fig. 4. In the extreme case, the approximation (21) is not sharp enough due to a change in concavity of the effective potential near the event horizon as it can be evinced from Figs. 4 and 5, and therefore, it should be replaced by in order to achieve a better fit with the behaviour of the effective potential in a neighbourhood of the event horizon. The asymptotic expansions (19)- (22) show that U e f f is a shortrange potential. As such, U e f f must satisfy the condition [46] +∞ −∞ Since the above integral can be computed analytically and the effective potential as a function of the tortoise coordinate is everywhere positive on the whole real line, we end up with the following inequality linking the event horizon of a noncommutative geometry inspired Schwarzschild black hole with its mass and the angular momentum quantum number of the scalar field, namely The derivation of (24) can be found in the Appendix A. It is interesting to observe that there is no counterpart of (24) for a classic Schwarzschild black hole.

A WKB approach
We compute the quasinormal mode frequencies of a noncommutative geometry inspired Schwarzschild black hole in the presence of a massless scalar field by means of the WKB (Wentzel-Kramers-Brillouin) approximation carried to sixth order beyond the eikonal approximation [47]. On the way, we discover some shortcomings of the method, especially when applied to the noncommutative Schwarzschild case. In principle, the WKB approximation already at the third order Table 1 Normal modes for scalar perturbations of the noncommutative geometry inspired Schwarzschild metric. The third column represents the numerical values found by [50] in the case of a classic Schwarzschild black hole while the fourth columns reports the corresponding WKB results found by [43]. Here, σ = Mω is the dimensionless frequency emerging after having expressed the radial coordinate in units of the black hole mass M, and σ e denotes the QNM in the case of an extreme non-commutative geometry inspired Schwarzschild black hole  Table 2 Transition from unstable to stable modes in the cases = 0, n = 0, 1, and = 1, n = 2, 3, using the WKB approximation at the third order and make the ansatz with A and φ functions to be determined. Substituting (26) into (25) gives where a prime denotes differentiation with respect to x * . Suppose that A ∼ 0. Recalling that in the case |x * | → ∞ the function Q(σ, x * ) tends to a constant, the solution of (25) is It is clear that this will be the case when Q(σ, x * ) and V e f f (x * ) vary slowly with respect to the wave length λ =  2π/σ . Under this assumption, it is reasonable to expect that ψ σ (x * ) is approximately given by (29) even when A(σ, x * ) changes gradually, i.e. its variation is appreciable only on intervals of total length of the order of some wavelengths. Solving (27) and (28) under the aforementioned assumption gives Observe that the expression above will diverge when Q(σ, x * ) ∼ 0. This aspect turns out to be important in the determination of the QNMs. Since Q(σ, x * ) tends to a constant as |x * | → ∞, the behaviour of the radial wave function will be captured by (29). Furthermore, for a wave coming from space-like infinity and hitting the potential barrier over some interval of finite length only a fraction of the wave will be transmitted towards the event horizon and the transmitted wave will be reduced with respect to the incoming wave by a certain factor. In the case σ V e f f (x * ) we expect instead full transmission. If x * ,m denotes the maximum of −Q(σ, x * ), the first scenario described above will apply for Q(σ, x * ,m ) < 0 while the second one for Q(σ, x * ,m ) ∼ 0 [48]. There are some boundary conditions to be considered. First of all, we need to verify that for x * → −∞ there are only incoming waves towards the event horizon. Furthermore, there must be no waves coming from x * → +∞ For this reason, the interaction of the field with the potential peak must generate a reflected and a transmitted wave, both having the same order of magnitude. In light of the above discussion, we expect that the waves will exhibit characteristic frequen-cies such that Q(σ, x * ,m ) ∼ 0. This implies that (30) is not valid in some neighbourhood (x * ,1 , x * ,2 ) around x * ,m and hence, it should be replaced by At this point we still need to match (31) with the solution on the interval (x * ,1 , x * ,2 ). The problem of matching simultaneously the two exterior WKB solutions across both turning points can be approached as in [43,49]. To this purpose, we consider a Taylor expansion of Q(σ, x * ) at the sixth order in a neighbourhood of the maximum x * ,m in the inner region. This procedure allows to derive an asymptotic approximation to the interior solution which in turn can be used to connect the two WKB solutions in (31). As a result, two connection formulas emerge that relate the amplitude of the incoming and outgoing wave on either side of the potential wall. The boundary condition of a purely outgoing wave leads to the following semi-analytical formula for the QNMs [43,49] with where = d/dx * . Since the formulae for the higher order derivatives of the effective potential can be easily evaluated using Maple, they will not be presented here. Table 1 shows that the extreme noncommutative geometry inspired Schwarzschild black hole is unstable, that is a positive imaginary part of the QNM leads to exponential growth of a scalar perturbation. Using different boundary conditions in the context of black plus mirror, [51] also finds similar instability behaviour for small Kerr-AdS black holes. Furthermore, the non extreme case also exhibits instabilities when the mass parameter μ slightly exceed the corresponding mass parameter of an extreme Schwarzschild black hole. Table 2 also shows that for fixed but increasing n the corresponding QNM of the noncommutative geometry inspired Schwarzschild black hole approaches its classical counterpart already for small values of the mass parameter. These results at third order WKB suggest that the effects due to noncommutative geometry are relevant for microscopic black holes but they can be neglected when dealing with black holes of astrophysical interest.
Extending the above calculations up to the sixth order, as mentioned above, we notice that the instabilities are not a real effect but rather an artifact or limitation of the WKB method. To demonstrate this limitation, in Table 3, we show the QNMs for two values of μ which represent the extreme and non-extreme cases. We observe that the results do not converge when l ≤ n and in fact the convergence gets better when l n. Though the table displays the results up to the sixth order, for the cases where l ≤ n, there is no convergence found even up to the thirteenth order of WKB. To conclude this section, we note that though the WKB approximation has in general been found to lead to precise results when extended up to the 13th order, one should be cautious in the interpretation of new effects such as instabilities while using this method.

The inverted potential method
It has been shown in [52][53][54][55] that the QNMs of a black hole can be related to the bound states of the so-called inverted black hole potentials. Relying on this property, we approximate the effective potential expressed in the tortoise coordinate by means of the Pöschl-Teller potential [56] where U 0 and −2U 0 κ 2 denote the height of the maximum and the curvature of the effective potential at the maximum, more precisely with f (x) given by (10). Taking into account that the bound states of the Pöschl-Teller potential are given by [52][53][54][55] = κ − n + 1 2 the corresponding QNMs associated to the black hole potential barrier are obtained by applying the transformation Table 4 Normal modes for scalar perturbations of the noncommutative geometry inspired Schwarzschild metric. The third and fourth columns represent the numerical values found by [50] and [43], respectively, for a classic Schwarzschild black hole. Here, σ PT is the dimensionless fre-quency computed according to (36), and σ e,PT denotes the QNM in the case of an extreme non-commutative geometry inspired Schwarzschild black hole  (35). In particular, we have where the dependence on enters through the terms U 0 and κ. Fig. 6 shows how this Pöschl-Teller (PT) potential compares with the actual effective potential for different values of l and μ = 1.91. The PT potential reproduces the effective potential quite well in the peak region in all 3 cases.
In Table 4, we present some of the QNMs of a noncommutative geometry inspired Schwarzschild black hole computed according to the method described above. According to the method used, the imaginary part of the QNMs do not exhibit any sign flip. Furthermore, it is gratifying to observe that for large values of the parameter μ, the numerical values of the QNMs in the tenth column of Table 4 coincide with the quasinormal frequencies computed by [52][53][54][55] for the case of the scalar perturbations of the classic Schwarzschild black hole.

QNMs by the asymptotic iteration method (AIM)
From the above exposition it becomes evident that it is important to choose an appropriate method to determine the QNMs.
It is worthwhile to look for other similar cases where conventionally used methods fail.
In this context, it is interesting to observe that the new branch of QNMs for a massless scalar field in the classic Schwarzschild metric derived analytically in [57][58][59] can also be obtained numerically by using the improved AIM [60][61][62]. In order to apply this method, we must consider an ad hoc initial ansatz for a massless scalar field by incorporating its behaviour at the horizon and at infinity as we did in equation (12) in [57]. This gives rise to the homogeneous linear second order differential equation where The implementation of the improved AIM requires first to differentiate (37) n times with respect to the independent variable, which produces the following equation with λ n (x) = λ n−1 (x) + s n−1 (x) + λ 0 (x)λ n−1 (x), (40) s n (x) = s n−1 (x) + s 0 (x)λ n−1 (x).
Then, λ n and s n undergo a Taylor series expansion (42) around some point, x 0 . By replacing the above expansions in (40) and (41) a set of recursion relations can be derived (see details in [60][61][62]) from which one obtain the quantization condition We solved (43) to find the QNMs. In Table 5, we show the fundamental quasi-normal frequencies for a massless scalar field in the classic Schwarzschild metric with = 0, 1, 2. It is worth mentioning that a number of 25 iterations was employed for the improved AIM method. Last but not least, the same method was adopted in [63] to obtain the QNMs computed by [50] with the continued fraction method.

Conclusions
Instabilities of microscopic black holes are an interesting area of research as certain examples in the literature show.
The existence of such instabilities could possibly hint toward the inappropriateness of models of microscopic black holes. We choose here, as an example, the case of a Schwarzschild black hole inspired by noncommutative geometry. Indeed one of our motivation to probe into this case is a hint of instability revealed at the third order of the WKB method. It is then obligatory to pursue further the nature of such a possibility, either by assuring the convergence of the WKB method or, in case the convergence fails, choosing another appropriate algorithm to determine the QNMs. We show by going up to the sixth order of WKB, that the convergence of the WKB method is not guaranteed for the cases under discussion. The inverted potential method is then a more suitable alternative which reveals that these instabilities were only an artifact of the third order WKB. Furthermore, we show how the asymptotic iterative method confirms the existence of a new branch of analytical QNMs recently discovered in [57].