Quasinormal modes and shadow of noncommutative black hole.

In this paper we investigate quasinormal modes (QNM) for a scalar field around a noncommutative Schwarzschild black hole. We verify the effect of noncommutativity on quasinormal frequencies by applying two procedures widely used in the literature. The first is the Wentzel-Kramers-Brillouin (WKB) approximation up to sixth order. In the second case we use the continuous fraction method developed by Leaver. Besides, we also show that due to noncommutativity, the shadow radius is reduced when we increase the noncommutative parameter. In addition, we find that the shadow radius is nonzero even at the zero mass limit for finite noncommutative parameter.

Introduction

Initial studies for black hole perturbations were done by Regge and Wheeler[1] and Zerilli[2] for Schwarzschild geometry, as well as for the Kerr black hole[3]. In 1970 Vishveshwara[4] identified a type of disturbance subject to special conditions such as outgoing waves in the spatial infinity and ingoing waves in the vicinity of the event horizon. These disturbances were called quasinormal modes valid only for a group of complex frequencies[5]. These quasinormal frequencies present a real part that provides the oscillation frequency while the imaginary part determines the damping rate of the modes. The dominant quasinormal modes can be seen in gravitational wave signals and in this case, the emitted waves are related to many physical processes such as astrophysical phenomena involving the evolution of binary systems and stellar oscillations or other highly dense objects in the early universe. In this way these quasinormal modes have been observed experimentally by LIGO/VIRGO[6,7].

The analysis of quasinormal modes has been widely explored in the literature[8-18] by using different mechanisms, such as the WKB approximation and numerical methods. The first works using the WKB approximation to find quasinormal modes were done by Schutz and Will[19]. Improvements in the method were made by Iyer and Will by adding corrections up to third order[20], and so the results for Schwarzschild black hole are very close to those obtained by the numerical method of Leaver[21] in regime. Aiming at a new improvement in the WKB approach, studies made by Konoplya[22] extended the method up to sixth order leading to more accurate results. Currently, we can find extensions of the WKB approximation up to thirteenth order[23]. As aforementioned, the numerical method is another way to obtain quasinormal frequency modes. Therefore, the first numerical approach to calculate quasinormal frequencies was described by Leaver, and the applied mathematical procedure is called a continuous fraction[24]. In[21], Leaver has obtained quasinormal modes for Schwarzschild and Kerr black holes and also for Reissner–Nordström black hole in[25]. Several works[26-29] have applied this numerical method that has presented a good precision.

In addition to quasinormal modes, in recent years several authors have devoted themselves to study the shadow of the black hole[30-34]. This shadow requires information about the geometry around the black hole, which makes its study a very important way to understand the properties near the event horizon. Moreover with the advancement and improvements of the experimental techniques allowed us the first image of a supermassive black hole in the center of the M87 galaxy by the Event Horizon Telescope[35,36] by using the properties of the shadow. These experimental results have been studied by several authors[37-39] stimulated by the possibility of understanding phenomena in the regime close to the event horizon.

In this work we will use the shadow ray to better understand the proximity of the noncommutative Schwarzschild black hole horizon. The noncommutative gravity has been extensively investigated, particularly in black hole physics, mainly due to the possibility of better understanding the final stage of the black hole—see[40-43] for further details. We know that there is a relationship between quasinormal modes and of the black hole shadow and that several investigations contributed to this understanding. One of the first studies that certainly served as a basis for structuring this relationship was made by Mashhoon[44], which describes an alternative method to calculate quasinormal modes at the eikonal limit. The most detailed geodesics study is shown in Cardoso et al.[45], which shows that the real part of the quasinormal modes is related to the angular velocity of the null circular orbit and the imaginary part is associated with the Lyapunov exponent. Stefanov et al.[46] in the eikonal regime established a connection between black hole quasinormal modes and lensing in the strong deflection limit. Currently, important results have been obtained at the eikonal limit, such as the relation between the real part of quasinormal frequencies and the black hole shadow radius[47-49].

Studies related to quasinormal modes of noncommutative black holes have been extensively carried out by several authors[50-55]. In this paper, we aim to determine the quasinormal modes of noncommutative Schwarzschild black hole via Lorentzian mass distribution in order to verify the changes caused by the noncommutative parameter. Moreover, we show that contrary to the case of the Schwarzschild black hole, the shadow radius presents a non-zero result at the zero mass limit. Therefore, at this limit, the shadow radius is proportional to a minimum mass. This result has not been obtained analytically by using the Gaussian distribution. Thus, considering the Lorentzian distribution, some results in an analytical way are more easily explored than in the Gaussian case where this is done numerically. Furthermore, in[56,57], we have also found a similar result when investigating the zero mass limit in the black hole absorption process. In that case, we have obtained a non-zero absorption at the zero mass limit. By considering a Lorentzian mass distribution to introduce the noncommutativity, in[58], we have explored the process of scattering and absorption of scalar waves through a noncommutative Schwarzschild black hole. Moreover, in[59-74], the thermodynamics of the BTZ and Schwarzschild black holes in the noncommutative background has been investigated by using the WKB approach in tunneling formalism[75-77]. An advantage of using the Lorentzian distribution in analytical calculus has been investigated in[59] and[60] where logarithmic corrections for entropy and the condition for black hole remnant formation were obtained.

We organize the paper as follows: in section “Noncommutative black hole with Lorentzian smeared mass distribution” we implemented the effect of noncommutativity in the Schwarzschild black hole metric by a Lorentzian smeared mass distribution, and we analyze the results for quasinormal frequencies. In section “Null geodesic and Shadow of a noncommutative black hole” we apply the null geodetic method to determine the shadow of the noncommutative black hole. In section “Conclusions” we make our final considerations.

Noncommutative black hole with Lorentzian smeared mass distribution

On this section we begin by considering a Lorentzian distribution[43,78] given bywhere is the noncommutative parameter of dimension and M is the total mass diffused throughout the region of linear size . Thus, the smeared mass distribution function becomes[58]

Hence, the line element of the Schwarzschild black hole in the noncommutative background is now given bywithwherewhich represent the radius of the event horizon and the Cauchy horizon, respectively.

The next step we consider the case of the massless scalar field described by the Klein–Gordon equation in the background (7)

Thus, we apply the separation of variables method in the above equation by using the following Ansatzwhere is the frequency and are the spherical harmonics.

Now, we can obtain a radial equation for :

We can reduce the radial equation (12) into a Schrödinger like equation by introducing a new coordinate (called tortoise coordinate) given by andso that the radial equation becomeswhere

In the following sections, we will obtain the quasinormal frequencies by two methods that are widely used in the literature. The first uses a sixth order WKB approximation, and the second method introduced by Leaver and improved by Nollert[82], consists of using the continuous fraction method to find numerically quasinormal modes.

WKB approximation

Quasinormal modes correspond to solutions of the wave equation (14) that satisfy the conditions of the purely outgoing waves at infinity and purely incoming waves at the event horizon, i.e.,

In this section we will use the WKB approximation to find the quasinormal modes. The first works using the WKB approximation to evaluate the quasinormal modes were done by Schutz and Will[19]. Improvements in the method were made using corrections up to third order[20,79] and up to sixth order by Konoplya[22]. The quasinormal modes are obtained by using the sixth order corrections for the WKB approximation as followswhere are the correction terms of the model. We have that is the maximum effective potential at point and refers to the second derivative with respect to the tortoise coordinate. We can obtain the values of by making . In Fig. 1 we show the curves of the effective potential for and , where we define for . In the Tables 1, 2 and 3 we present the tabulated quasinormal frequencies using sixth order WKB method.

Figure 1

The effective potential as function of the tortoise coordinate (a) and (b) .

Table 1

QN frequencies for .

Θ	ω_0		ω_1		ω_2
	6th order WKB	Numerical	6th order WKB	Numerical	6th order WKB	Numerical
0.00	0.292910–0.097762i	0.292936–0.0976600i	0.264471–0.306518i	0.264449–0.306257i	0.231014–0.542166i	0.229539–0.540133i
0.05	0.316239–0.099311i	0.316243–0.0992441i	0.290671–0.309586i	0.290591–0.309431i	0.260297–0.543374i	0.258427–0.542162i
0.10	0.351481–0.097397i	0.351435–0.0973776i	0.330103–0.300074i	0.329892–0.300209i	0.300773–0.519039i	0.299190–0.519345i
0.12	0.371967–0.091934i	0.371932–0.0919646i	0.346286–0.282086i	0.346358–0.282153i	0.300083–0.491675i	0.301082–0.491385i

Table 2

QN frequencies for .

Θ	ω_0		ω_1		ω_2
	6th order WKB	Numerical	6th order WKB	Numerical	6th order WKB	Numerical
0.00	0.483642–0.096766i	0.483644–0.0967588i	0.463847–0.295627i	0.463851–0.295604i	0.430386–0.508700i	0.430544–0.508558i
0.05	0.521842–0.098436i	0.521844–0.0984288i	0.504088–0.300022i	0.504087–0.299997i	0.474038–0.514099i	0.474067–0.513965i
0.10	0.580028–0.096808i	0.580028–0.0968027i	0.565520–0.293676i	0.565511–0.293660i	0.540100–0.499152i	0.540039–0.499083i
0.12	0.615658–0.091534i	0.615659–0.0915307i	0.599809–0.276911i	0.599817–0.276896i	0.569259–0.469194i	0.569307–0.469133i

Table 3

QN frequencies for .

Θ	ω_0		ω_1		ω_2
	6th order WKB	Numerical	6th order WKB	Numerical	6th order WKB	Numerical
0.00	0.675366–0.096501i	0.675366–0.0964996i	0.660671–0.292288i	0.660671–0.292285i	0.633591–0.496011i	0.633626–0.496008i
0.05	0.728594–0.098197i	0.728594–0.0981957i	0.715414–0.297060i	0.715414–0.297056i	0.691124–0.502924i	0.691137–0.502911i
0.10	0.809869–0.096642i	0.809869–0.0966413i	0.799161–0.291636i	0.799159–0.291633i	0.779124–0.491455i	0.779117–0.491442i
0.12	0.860294–0.091410i	0.860294–0.0914089i	0.848871–0.275415i	0.848871–0.275412i	0.826463–0.462957i	0.826463–0.462950i

Another possibility is to study scattering by using the WKB method done in[20]. In order to develop this investigation, we use the boundary conditions for Eq. (14) in the form:

Notice that as we want to obtain the reflection and transmission coefficients, as done for tunneling in quantum mechanics, we need the condition . The quantity in (14) is assumed to be purely real, and with these imposed conditions we can findwhere is purely real and are coefficients that depend on the effective potential and is a purely imaginary quantity. This way of studying scattering by using the WKB approximation can be also found in[23,80]. Thus, using the relationship between and the reflection and transmission coefficients obtained in[20] we get the following:

Now to find the coefficients, we just calculate the value that can be obtained by solving the Eq. (19). This method has a good approximation for as we can see in Fig. 2, where we have a comparison between the numerical results and the WKB approximation for the transmission coefficient for and . Notice that the curves obtained by the WKB approximation are very close to those obtained by the numerical method used in the paper[58], showing that the method presents excellent results. The sixth-order WKB approximation does not show good results for , improving only when we take large . This problem is also mentioned in[23].

Figure 2

Transmission coefficients for three multipole (from left to right) and . The approximation between the two methods is very good.

Leaver’s continued fraction

The numerical method is another way to obtain the quasinormal frequency, and this procedure has been described by Leaver[21,25], and which is also found in other works[26,27] showing that it is a method with good precision.

We start analyzing the radial equation (12) which is subject to boundary conditions at infinity and near the event horizon , such that one obtains the asymptotic solutionsWe can obtain a solution that has the desired behavior on the horizon (), and that can be written in the form

By replacing the solution (23) in the Eq. (12), we obtain the recurrence relationand

The recurrence relation coefficients , and are simple functions of k and the parameters , l and the radius and :

See that the solution (23) applied directly to the radial equation returns a three-term recurrence relation which makes it easier to use continued fractions. This is because our metric for a non-commutative black hole results in a generalized spheroidal wave equation whose solutions are connected by three-term recurrence relations—see “Appendix” for further details. However, the Leaver’s method is not necessarily limited to recurrence relations of this type as shown in[25], where it was considered an equation that describes odd-parity perturbations of a charged black hole, which equation has series solutions whose coefficients are connected by four-term recurrence relations.

The boundary condition at infinity will be satisfied for quasinormal frequency values , so that the series in (23) is absolutely convergent. Hence, we have a recurrence relation with three terms to determine the coefficient , and we can write in terms of a continuous fraction[81]which can also be found as follows

We can obtain the characteristic equation for quasinormal frequencies by assigning in (30) and comparing with reason obtained from (24),

With the equation above, we can obtain the quasinormal frequencies , by just calculating its roots numerically. However, the Eq. (31) is more used to find the fundamental frequency in the case of the more stable root, and another way of finding the modes is to invert this equation to a large number of k as follows

To complement the analysis, we will check the behavior for very large k, as done in[82]. Now, we have to reorganize the Eq. (25) dividing , and by to obtain

We can see that, , and we can get a more complete expression by expanding in power series in terms of ,

By considering the series up to the third term and admitting we have

Now, to obtain the values of and , we replace these two expressions above into Eq. (33) in the very large k regime by making the multiplications and restricting up to terms of the order to getsuch that we find the following:

Now we have,

Here, we see that at the limit and , we get the result for the Schwarzschild case initially found by Leaver.

Results

In the Tables 1, 2 and 3, we present some results for the quasinormal modes computed using the sixth-order WKB approximation, and the Leaver’s continues fraction method described in the previous section by admitting for various values of , l and n. We see that the results between the methods approach when , this is, due to the instability of the WKB method for small multipole numbers mainly close to zero this instability can also be seen in the graphs of Fig. 3. An important detail is in the sign of the imaginary part that is always negative when the frequency is associated with the scalar field. A justification for this is due to the exponential drop of the quasinormal modes over time by losing energy in the form of scalar waves.

Figure 3

Real (top) and imaginary (bottom) parts of the quasinormal frequencies as function of the n. We see that for the frequency curves are very dispersed with the increase of n, while for the curves of the real part of the frequency are more constant.

We can see the influence in the quasinormal modes for the noncommutative case by admitting values for where returns to the Schwarzschild case. With the increase of the non-commutative parameter , we have an increase in the real part of the quasinormal frequency, while the imaginary part begins to grow and then decreases. Another way of visualizing the effects of the noncommutative parameter is through the graphs of Fig. 3, these plots were obtained using the WKB method. We have quasinormal modes where in the plot we depicted the real part (top) and the imaginary part (bottom). The modes are based on n for the following multipoles numbers and 4, for which we see that the results become more linear with the results varying and and 4. Thus, in Fig. 3 we observe that, by varying l and the parameter , the imaginary part of the frequency does not cross the horizontal axis or change sign, thus indicating that the black hole remains stable due to scalar perturbation. In addition, it is interesting to make a plot for the complex plane as in Fig. 4 where we consider three families of multipoles and varying as follows 0 (black), 0.05 (red), 0.10 (blue), 0.12 (green). The left panel was obtained by the WKB approximation and the right panel by the continuous fraction method (numerical). We can see that the frequency curves incline more closely when we use the WKB approximation. We also see that for the extreme case the curve tilts more to the left in both methods.

Figure 4

Complex plane of the QNMs. In plot (a) we have the results obtained by the WKB approximation, while in (b) we use continuous fraction. The markers denote the multipole number as: (circle), (square) and (diamond), while the colors denote the value of the noncommutative parameter (black), (red), (blue) and (green).

Null geodesic and Shadow of a noncommutative black hole

We know that in the vicinity of a black hole all the photons are absorbed so that a distant observer looking at the black hole, in absence of any other source, will see a spot created by this absorption, and this spot is usually called the shadow of the black hole. These shadows have been studied long ago by Synge[83] and Luminet[84] who started studies for Schwarzschild black hole while Kerr black holes were studied by Bardeen[85].

As these shadows correspond to the apparent shape of the photon capture orbits, the space-time metric itself is enough to determine them and thereby better understand the geometry of the near horizon. One method of determining the apparent shape of the black hole is through the shadow boundary that can be studied by the equations of null geodesics. Similar studies have been done in different contexts[87-92].

Null geodesic

We can find the geodesics from Eq. (7) by taking a Lagrangian in the form

Thus, we havewhere the “” is the derivative with respect to an affine parameter.

We are interested in the path of a ray of light in the described metric, which is spherically symmetrical, so if we analyze in a plane, any ray of light that begins with a certain angle must remain with the same angle. We will then consider an equatorial plane by setting the angle to .

Thus, two equations are enough to describe the movement of a beam of light. We can put together a system with these equations that give rise to two geodesic motion constants E and L, which correspond to energy and angular momentum respectively:

Now, as our goal is to study the null geodesics we have to , and using the Eq. (43) we can write

Introducing a new variable we can write the orbit equation as followswhere is the impact parameter defined as the perpendicular distance (measured at infinity) between the geodesic and a parallel line that passes through the origin. So differentiating (45) we have,

By solving the Eqs. (45) and (46) numerically, we can obtain the behavior of the geodesic lines for different values of the impact parameter b. In the Fig. 5, we verify the change of the geodesic lines for different impact parameters b and also by varying the values of the noncommutative parameter. In the figures, we have a black disk that represents the limit of the event horizon, the internal dotted circle is the radius for the photon sphere (critical radius), and the external dashed circle is the critical impact parameter (shadow). Hence, we see that the noncommutative parameter decreases the effect of the black hole on the light beams. For a similar effect see also[93].

Figure 5

Geodesics surrounding a noncommutative black hole. The impact parameters defined as and 6.5 are the same for all graphs, assuming . We can clearly see the influence of the noncommutative parameter on the geodesic curves from (a) (Schwarzschild case) to (d) .

Critical orbit and shadows

It is known that the shadow of the black hole is directly related to the impact parameter for the photon orbit, as we will see below. So to determine the shadow limit, we will start by studying the effective potential that satisfies the equation of the null geodesic as followswhere using (44), we have

In this case, we can obtain a critical radius or critical circular orbit for a photon and critical impact parameter , by using the following conditions: and . So we find

Let us now compute the size of the black hole shadow that can be expressed via celestial coordinates as follows[85]where is the observer position at infinity.

As our study is restricted to the equatorial plane, the radius that delimits the size of the shadow is equivalent to the critical impact parameter, and so we havewhere

In a semiclassical description of the scattering[86], the impact parameter is associated with each partial wave in the large l regime. As shown in[45], the real part of the quasinormal frequencies at the eikonal limit corresponds to the angular velocity for the last null circular orbit and the imaginary part is associated with the Lyapunov exponent which determines the unstable timescale of the orbitbeing the angular velocity given by

This suggests a relationship between quasinormal frequencies and shadow radius as done by Jusufi[47] at the eikonal limit, showing that the real part of quasinormal modes is inversely proportional to the radius as follows

Note that this expression is valid only for large values of l in most of the cases although fails for Einstein–Lovelock theory as shown by Konoplya and Stuchlik[94]. We can see in Fig. 3 that the higher the value of l the smaller the contribution of n. Thus by using the results obtained by the WKB approximation we compare with the shadow radius .

In the Table 4 we see that by increasing the value of l, the results between the real part of the quasinormal frequencies and the black hole shadow radius approach each other. Now, we can express by considering small, and so we get the following approximate expressionNote that for , we have the shadow radius for the Schwarzschild black hole case. Therefore, we notice that the shadow radius is reduced when we change the parameter . In Fig. 6, we see the circles that represent the shadow boundaries of the noncommutative black hole for different values of . Note that we have a reduction in the circles when we vary . Furthermore, taking in (55), we obtain a non-zero result for the shadow radius, that iswhere is the minimal mass[60]. Therefore, at the limit of the shadow radius is proportional to the minimum mass and the black hole becomes a black hole remnant. We have shown this behavior in Fig. 7. In Fig. 7, we show the behavior of the shadow radius by keeping fixed and assuming small values of M. In Fig. 8, We show the behavior of the shadow radius by keeping M fixed and assuming small values of .

Table 4

The real part of quasinormal frequencies for large l compared with shadow radius .

M=1	Re(ω)/l+1/2		R_s^-1
Θ	l=100	l=1500	–
0.00	0.192450707	0.192450092	0.192450090
0.05	0.207583654	0.207582949	0.207582946
0.10	0.230750607	0.230749836	0.230749832
0.12	0.245314933	0.245314347	0.245314347

Figure 6

We see the influence of non-commutativity in the shadow admitting and and 0.12.

Figure 7

We see the influence of non-commutativity in the shadow admitting (a) and . (b) and .

Figure 8

We see the influence of non-commutativity in the shadow admitting (a) and . (b) and .

Conclusions

In summary, in this work, we investigate the quasinormal frequencies for a noncommutative Schwarzschild black hole by two different methods in order to investigate and compare the results. Using the sixth-order WKB approximation and Leaver’s continuous fraction, we found that there is a small difference between the quasinormal frequencies obtained by each method mainly for small multipoles. The effects of the noncommutative parameter cause an increase in the real part of the quasinormal frequencies, while the magnitude of the imaginary part begins to grow and then decreases. For the black hole shadow we use the results obtained by the WKB method to verify that in large l regimes the real part of quasinormal modes is inversely proportional to the shadow radius. In addition, we have shown that the shadow radius is non-zero at the zero mass limit. Therefore being proportional to a minimum mass. However for , we recover the shadow radius for the Schwarzschild black hole case. Finally, we also notice that the shadow radius is reduced when we increase the noncommutative parameter.