# Gravitational Quasinormal Modes of Regular Phantom Black Hole.

1. IntroductionAs a major topic in cosmology, the accelerated expansion of our universe has caused widespread concern in the scientific community. Since the effect of gravity causes the expansion speed to slow down, the accelerated expansion of the universe implies the existence of an unknown form of energy in the universe. The latter provides a repulsive force to push the expansion of the universe. Such unknown energy is called dark energy (DE). Subsequently, a large number of DE models have been proposed, among which the one with cosmological constant is the most famous. Even though the model of DE with the cosmological constant is reasonable in physical theory and consistent with most observations, two difficulties still remain unsolved, namely, how to derive "vacuum energy" from quantum field theory and why the magnitudes of present DE and dark matter are of the same order.

Many modern astrophysics observations indicated the possibility of pressure to density ratio w < -1. For example, a model-free data analysis from 172 type Ia supernovae (SNIa) resulted in a range of -1.2 < w < -1 for our present epoch [1]. According to the WMAP data during 7 years, w = -[1.10.sup.+0.14.sub.-0.14] (1[sigma]) [2]. By using the data from Chandra telescope, an analysis of the hot gas in 26 X-ray luminous dynamically relaxed galaxy clusters gives w = -[1.20.sup.+028.sub.-0.28] [3]. The data on SNIa from the SNLS3 sample estimates w = -[1.069.sup.+0.091.sub.-0.092] [4]. In fact, several DE models with a supernegative equation of state provide better fits to the above data [5-8]. And all these approaches are in favor of phantom DE scenario [9-13], in which a constant equation of state parameter is used [14, 15]. This implies that the phantom model might be meaningful for in-depth understanding of DE.

In the phantom model, the signature of the metric is +2, and the action of the model reads

S = [integral] x [square root of -g] [d.sup.4] x [R + [epsilon][g.sup.[mu]v] [[partial derivative].sub.[mu]][phi][[partial derivative].sub.v][phi] - 2V ([phi])], (1)

where R is the scalar curvature, V([phi]) is the potential of the scalar field, and [epsilon] = -1 corresponds to a phantom scalar field while [epsilon] = +1 is for a normal canonical scalar field.

Bronnikov and Fabris first investigated the properties of BH with phantom scalar field in vacuum and derived a phantom regular BH solution 10 years ago [16]. Inside the event horizon of such phantom BH there is no singularity similar to the case of regular BHs with nonlinear electrodynamics sources [17]. Outside the event horizon, the properties of a phantom BH are similar to those of a Schwarzschild BH. Due to the absence of the singularity, such phantom regular BH solution has attracted much attention from researchers.

On the other hand, the research of BH perturbation has always been an important issue in BH physics. The first work on QNM in AdS spacetime was about scalar wave in Schwarzschild-AdS spacetime [18], which is then followed by a study on scalar wave in topological AdS spacetime [19]. There are a large number of works on regular BH's QNMs [17, 20-25]. Among various types of perturbation, gravitational perturbation is generally considered to be the most important form due to its practical significance. The intrinsic properties and the stability of a BH can be unfolded through its corresponding gravitational perturbation. In the fifties of last century, Regge and Wheeler began to study the gravitational perturbations of static spherically symmetric BHs. It was pointed out later that [26] the higher dimensional gravitational perturbations can be classified into three types, namely, scalar-gravitational, vector-gravitational, and tensor-gravitational perturbations. The first two types are associated with odd (vector-gravitational) and even (scalar-gravitational) parity in accordance with the spatial inversion symmetry of the perturbations and are of great physical interest [27]. These findings significantly simplify the study of gravitational perturbation of BH. Subsequently, people developed many new methods, and further studies on the gravitational perturbations result in a large number of master equations for various forms of BHs in 4-dimensions [27-30], in higher dimensions [26], and for stationary BHs [31, 32]. In fact, gravitational perturbations of a BH may generate relatively strong gravitational waves (GWs). Recently, the GWs from a binary BH system have been detected by LIGO [33], so BH is proven to be the most probable source of GWs by modern technology. Meanwhile since many alternative theories of gravity can produce the same GW signal within the present accuracy in far field, the reported GW detection still leaves a window for alternative gravity theories [34], which included the theory of phantom BHs. Therefore, the properties of QNMs of gravitational perturbation near the horizon of phantom BH may provide us essential information on the underlying physics of gravity theory. This is the main purpose of the present study.

The thermodynamics of BHs is also an important subject in BH physics. Some works indicated that Hawking radiation can be considered as an effective quantum thermal radiation around the horizon [35, 36], where the corresponding Hawking temperature can be derived from the tunneling rate [35-40]. Furthermore, a natural correspondence between Hawking radiation and QNM has been established recently [35-37,39,41]. Therefore, in this work, we will also investigate the Hawking radiation of regular phantom BH.

The paper is organized as follows. In Section 2, we briefly review the regular phantom BH solutions and discuss their properties in three different spacetimes, namely, asymptotically flat, de Sitter (dS), and anti-de Sitter (AdS). In this work, we focus on the odd parity and even parity gravitational perturbations. As the main component of this paper, Section 3 includes two subsections. In Section 3.1, we derive the master equation for odd parity gravitational perturbation and analyze the corresponding temporal evolution of the perturbed metric; in Section 3.2, corresponding studies are carried out for the even parity gravitational perturbation. In Section 4, we calculate the Hawking radiation of the regular phantom BH. We summarize our results and draw concluding remarks in Section 5.

2. The General Metric for Regular Phantom Black Holes

In this section, we discuss the phantom (e = -1) regular BH solution by considering the following static metric with spherical symmetry:

d[s.sup.2] = -f (r) d[t.sup.2] + d[r.sup.2]/f(r) + p [(r).sup.2] (d[[theta].sup.2] + [sin.sup.2] [theta]d[[phi].sup.2]) . (2)

According to the action, (1), the field equation for a self-gravitating minimally coupled scalar field with an arbitrary potential V([phi]) can be expressed as

[[??].sub.[mu]v] = [R.sub.[mu]v] - [g.sub.[mu]v]/2 (R - [epsilon][[phi].sup.;[alpha]] [[phi].sub.;[alpha]] - 2V ([phi])) - [epsilon][[phi].sub.;[mu]][[phi].sub.;v] = 0. (3)

By combining the scalar field equation,

[epsilon][[phi].sup.;[alpha].sub.;[alpha]] - dV ([phi])/d[phi] = 0, (4)

a regular phantom BH solution can be obtained as

f(r) = {c/[b.sup.2] + 1/[p.sup.2](r) + 3m/[b.sup.2] [br/[p.sup.2](r) + arctan (r/b)]} x [p.sup.2](r), (5)

where

p (r) = [square root of [b.sup.2] + [r.sup.2]], (6)

[mathematical expression not reproducible], (7)

m is the Schwarzschild mass defined in the usual way, and c and b are integration constant and scale parameter, respectively. Then it is necessary to determine the possible kinds of spacetime for such phantom BH, which can be classified as a regular infinity (r [right arrow] [infinity]) to be flat, de Sitter (dS), or anti-de Sitter (AdS). The corresponding parameters c, b, m should be restricted in each spacetime.

For the asymptotically flat spacetime, in accordance with (5), one has c = -3[pi]m/2b and

m = 2[b.sup.3]/3[[pi][b.sup.2] - 2b[r.sub.h] + [pi][r.sup.2.sub.h] - 2[b.sup.2] arctan([r.sub.h]/b) - 2[r.sup.2.sub.h] arctan ([r.sub.h]/b)]. (8)

In this case, the spacetime is asymptotically flat, namely, r [right arrow] [infinity], f(r) [right arrow] 1. And [r.sub.h] is the event horizon of the phantom BH. We note when b [right arrow] 0 that (5) becomes Schwarzschild flat spacetime.

For the de Sitter spacetime, one has

[mathematical expression not reproducible], (9)

where [r.sub.c], [r.sub.h] are the cosmological horizon and event horizon, respectively. We note when b [right arrow] 0 that (5) becomes Schwarzschild dS metric.

For the anti-de Sitter spacetime, we choose f(r) [right arrow] [r.sup.2] with r [right arrow] [infinity] without loss of generality. By expanding (5) around infinity, one finds

[mathematical expression not reproducible], (10)

where [r.sub.h] is the event horizon of AdS spacetime. We note when b [right arrow] 0 that (5) can be returned to Schwarzschild-AdS spacetime. In this context, the parameter b measures the coupling strength between phantom scalar field and the gravity for all three spacetimes.

Since the parameters c, m can be expressed in terms of b, [r.sub.h], and [r.sub.c] (dS), the structures of the regular phantom BH spacetime are completely determined by b, [r.sub.h], and [r.sub.c] (dS) (cf. Figure 1). One can readily verify that all the spacetimes are indeed nonsingular even at r = 0. As for any asymptotically flat spacetime, such flat regular phantom BH has a Schwarzschild-like structure. However, its tendency of approaching flat spacetime at infinity becomes slower with increasing b. In the dS case, the spacetime is bounded by two horizons, that is, [r.sub.h] < r < [r.sub.c]. There is a maximum for f(r), and it decreases with increasing b. For the AdS phantom BH, f(r) [right arrow] [r.sub.2] when r [right arrow] [infinity], and for larger b, the approach to the asymptotic solution becomes slower.

3. Gravitational Quasinormal Frequencies for Regular Phantom Black Holes

As proposed by Regge-Wheeler, two important gravitational perturbations are of odd and even parity. The perturbation gauge for each type has its own definition. In this section, we choose the Regge-Wheeler-Zerilli gauge to discuss the master equation for each perturbation type. Here we take the magnetic quantum M = 0 to make [phi] disappear completely because all values of M lead to the same radial equation [27]. Since the total number of equations in the even parity case is bigger than that in odd parity case, the derivation of the master equation for even parity perturbation is thus more complicated. Once the master equation is derived, the effective potential and the corresponding quasinormal modes can be obtained. We will first discuss the odd parity gravitational quasinormal frequencies in asymptotically flat, dS, and AdS spacetimes and then study the case for even parity. In our work, we consider the metric perturbations not only of the Ricci curvature tensor and scalar curvature (the l.h.s. of the Einstein field equation), but also of the energy momentum tensor (the r.h.s. of the Einstein field equation). On the other hand, we will not consider the perturbations of the phantom scalar field. This is because such perturbations can be canceled out through an appropriate choice of V in the action (see Appendix for details).

In order to discuss gravitational perturbation, one may write down

[g.sub.[mu]v] = [[bar.g].sub.[mu]v] + [h.sub.[mu]v], (11)

where the small perturbation [h.sub.[mu]v] will be divided into odd and even modes in the subsequent sections.

3.1. Master Equation and Quasinormal Modes for Odd Parity Perturbation. The odd parity perturbation [h.sub.[mu]v] has the form as [27]

[mathematical expression not reproducible], (12)

where [Q.sub.p] ([theta]) = sin [theta]d[P.sub.l] (cos [theta])/d[theta] ([P.sub.l] (cos [theta]) is the Legendre function), which satisfies

[Q".sub.p] - cot ([theta]) [Q'.sub.p] ([theta]) = -k[Q.sub.p] ([theta]), (13)

where k = l(l + 1) and l is the angular quantum number.

Then the separation of variables can be carried out by writing [h.sub.0] (r, t) = exp(-i[omega]t) [h.sub.0] (r), [h.sub.1] (r, t) = exp(-i[omega]t) [h.sub.1] (r).

By substituting (5)-(7), (11)-(13) into the field equation (3) and onlykeeping the first order perturbation terms, we obtain the independent perturbation equations as follows:

[mathematical expression not reproducible], (14)

[mathematical expression not reproducible]. (15)

Equation (15) implies [h.sub.0] = -([f.sup.2]/i[omega])([h.sub.1] f'/f + [h'.sub.1]). Substituting [h.sub.0] into (14), we get the master equation for odd parity perturbation:

[mathematical expression not reproducible]. (16)

Finally, we renormalize [h.sub.1] by

[h.sub.1] (r) = B(r) [PHI] (r); B(r) = p(r)/f(r). (17)

By substituting (17) into the master equation, (16), and using a tortoise coordinate [r.sub.*], the Schrodinger-type wave equation for this case can be expressed as

[d.sup.2][PHI]/d[r.sup.2.sub.*] + ([[omega].sup.2] - [V.sub.o] (r)) [PHI] = 0, (18)

where the effective potential for odd parity perturbations [V.sub.o] (r) is

[mathematical expression not reproducible]. (19)

Equation (19) can be used to describe the effective potential of "odd"-type perturbation [V.sub.o] in different spacetimes and be utilized to discuss the relationship between the effective potential and model parameters such as the angular harmonic index I and the parameter b.

(I) Figures 2, 3, and 4 show the potential functions, temporal evolution of the gravitational perturbation, and quasinormal frequency obtained by WKB method in asymptotically flat spacetime.

(i) The form of the effective potential as a function of r for different values of l and b is shown in Figure 2. From the left plot, one sees that as b increases, the shape of the effective potential becomes smoother. The maximum of the effective potential decreases with increasing b and the position of the peak shifts to the right. From the right plot, we see that, with the increase of angular quantum number l, the maximum of the effective potential also increases. V(r) is always found to be positive definite outside the event horizon, which indicates that the corresponding QNMs are likely to be stable.

(ii) We adopt the finite difference method to analyze the stability of such BH. By applying the coordinate transformation (t, r) [right arrow] ([mu], v) with [mu] = t - [r.sub.*], v = t + [r.sub.*] to (16) and integrating numerically using the finite difference method [24, 42-44], we obtain the differential equation for [h.sub.1] ([mu], v). Figure 3 shows the stability of a regular phantom BH in asymptotically flat spacetime with [r.sub.h] = 1. The temporal evolution of each mode in Figure 3 corresponds to a corresponding case in Figure 2. Since V(r) decreases significantly with the growth of b, the decay rate ([absolute value of (Im([omega]))]) becomes smaller and oscillation frequency (i.e., Re([omega])) also drops; as I increases, the values of V(r) are raised correspondingly, so that the oscillation frequency (i.e., Re([omega])) slightly increases together with the decay rate ([absolute value of (Im(w))]).

(iii) We employ the WKB approximation [45, 46] to evaluate the quasinormal frequencies.

The complex frequency [omega] is determined by [47]

[mathematical expression not reproducible], (20)

where [mathematical expression not reproducible] are presented in [47, 48]. By making use of (19), we evaluate the QNM frequencies by employing the 6th order WKB method (see Figure 4). Figure 4 shows that the fundamental quasinormal modes (n = 0) have the smallest imaginary parts, as the modes decay the slowest. As n increases, the imaginary part of the corresponding quasinormal mode becomes bigger for given l and b. For a given principal quantum number n, both the real and the imaginary parts of the frequency decrease with increasing b; on the other hand, the real part of the frequency increases significantly with the angular quantum number b, while the imaginary part also increases slightly with increasing b.

(II) Figures 5, 6, and 7 show the potential function between the event horizon [r.sub.h] and cosmological horizon [r.sub.c], temporal evolution of the gravitational field, and quasinormal frequency obtained by WKB method in dS spacetime.

(i) The form of the effective potential in dS spacetime is similar to that in asymptotically flat spacetime. Figure 5 indicates that, for given l and r, the effective potential decreases with increasing b; and for given b and r, it increases with increasing angular quantum number l. In the range [r.sub.h] < r < [r.sub.c], V(r) is also positive definite, which implies that the corresponding QNM is likely to be stable.

(ii) Figure 6 studies the stability of a regular phantom with [r.sub.h] = 1, [r.sub.c] = 10. Similar to the case in asymptotically flat spacetime, it is found that, with increasing b, the oscillation frequency (i.e., Re([omega])) decreases, while the decay rate ([absolute value of (Im([omega]))]) becomes smaller. Therefore, for smaller value of b, the BH returns to its stable state more quickly as small perturbation dies out faster.

(iii) In accordance with the results of finite difference method, the frequencies presented in Figure 7 by WKB method also show that, for given n and l, the decay rate of perturbation (i.e., [absolute value of (Im([omega]))]) decreases with increasing b. Moreover, it illustrates that, for a given b, the fundamental mode can be found at n = 0, l = 2. Due to the smallness of the real parts of the frequencies, these modes of oscillation persist for longer time.

(III) Figures 8 and 9 show the potential function beyond the event horizon r > [r.sub.h] and the temporal evolution of the gravitational perturbation in AdS spacetime.

(i) The form of the effective potential in AdS spacetime is quite different from that in asymptotically flat and dS spacetime. Figure 8 indicates that the value of V(r) is divergent as r [right arrow] [infinity]. The effects of b and I on the effective potential are similar to previous cases; for a given r, smaller b or larger l leads to bigger V.

(ii) Figure 9 studies the stability of AdS regular phantom BH spacetime with [r.sub.h] = 1. The results are consistent with the above calculated potential function. It is again inferred that the fundamental mode of gravitational perturbation in odd parity occurs for l = 2 and larger b, since such kind of QNM will take a longer time to be stable.

3.2. Master Equation and Quasinormal Modes for Even Parity Perturbation. Another canonical form for the gravitational perturbations is of even parity. After applying the separation of variables, it can be expressed as [27, 28]

[mathematical expression not reproducible], (21)

where [H.sub.0] (r), [H.sub.1] (r), [H.sub.2] (r), K(r) are unknown functions for the even parity perturbation. It is noted that these functions are not independent.

Now we derive the first order perturbation equations by substituting (5)-(7), (11), and (21) into (3) and find the relationships among [H.sub.0] (r), [H.sub.1](r), [H.sub.2](r), K(r):

[mathematical expression not reproducible] (22)

[mathematical expression not reproducible], (23)

[mathematical expression not reproducible], (24)

[mathematical expression not reproducible], (25)

[mathematical expression not reproducible], (26)

where

[mathematical expression not reproducible]. (27)

Solving (23), [H.sub.0] can be expressed as

[mathematical expression not reproducible]. (28)

We define a function [PSI](r) satisfying

[H.sub.1] = -i[omega] r/f ([PSI] + K). (29)

By observing (22), (28), and (29), it turns out that we need to express K by [PSI] in order to express all perturbation functions [H.sub.0] (r), [H.sub.1] (r), [H.sub.2] (r), and K(r) in terms of [PSI]. This can be achieved by substituting (28) and (29) into (25) and (26) and evaluating the subtraction (25)-(26), and one eventually obtains the following expression:

[mathematical expression not reproducible], (30)

where

[mathematical expression not reproducible]. (31)

Substituting (28)-(31) into (24), K can be solved as

K(r) = [A.sub.1] (r) [PSI]' (r) + [A.sub.2] (r) [PSI] (r), (32)

where [A.sub.1] (r), [A.sub.2] (r) are

[mathematical expression not reproducible]. (33)

The resulting master equation can be derived by evaluating

[mathematical expression not reproducible]. (34)

By substituting (28), (29), and (32) into the above equation, the corresponding master equation is given by

[mathematical expression not reproducible], (35)

where [P.sub.v] (r) = -2 + l + [l.sup.2] + 2[p.sup.2] V([phi]) + 3pf' p'. Finally, one defines [PSI](r) = [B.sub.1] (r) [B.sub.2] (r) [PHI] (r) and f(r) = Q(r) F(r) (where Q(r) is the coefficient of [[omega].sup.2]); the master equation can be simplified into the following equation:

[d.sup.2][PHI](r)/d[r.sup.2.sub.*] + {[[omega].sup.2] - [V.sub.e] (r)} [PHI] (r) = 0, (36)

where d[r.sub.*] = dr/F(r), [V.sub.e] = F(r) (2Q(r) (f'(r) Q'(r) + F(r) Q" (r)) - F(r) Q'[(r).sup.2]) + 4[??](r)/4Q[(r).sup.2], (37)

where

[mathematical expression not reproducible]. (38)

Here we study the QNMs for even parity perturbation in asymptotically flat, dS, and AdS spacetimes. Figures 10, 11, and 12 show the effective potential functions and corresponding temporal evolution of even parity perturbation in asymptotically flat, dS, and AdS spacetimes, respectively. We consider the potential in the range r > [r.sub.h] in asymptotically flat and AdS spacetimes and [r.sub.h] < r < [r.sub.c] in dS spacetime. Here, the potential functions are all positive definite so that the corresponding regular phantom BHs are likely to be stable. Furthermore, in asymptotically flat and dS spacetimes, the relationships between Re([omega]), Im([omega]) and b, l are similar to those of odd parity; the differences are in details. However, the potential function in AdS spacetime is mostly a convex function, and therefore it approaches infinity faster than that of odd parity. Nevertheless, the resulting QNMs in AdS spacetime are also found to be stable.

4. Hawking Radiation of the Regular Phantom Black Hole

As mentioned before in Section 2, the interior structures of regular BHs are quite different from those of traditional ones. In order to achieve a better understanding of the properties of such regular phantom BHs, it is meaningful to investigate the corresponding Hawking radiation, which is considered to be a promising method to study the thermodynamics of BHs.

From a quantum mechanics viewpoint, QNM carries information on the quantization of the system around BH's horizon [35-37, 39, 41]. The Hawking radiation spectrum provides an interpretation of such quantization at an effective temperature [25, 35, 36]. In this paper, we adopt a simple and effective method (i.e., Hamilton-Jacobi equation) to analyze the Hawking radiation of the regular phantom BH. Based on the tunneling theory of Hawking radiation [49, 50], this method was first put forward by Srinivasan in 1999 [51] and plays an important role in recent studies [52, 53].

Up to this point, we have been using the natural units in the calculations of QNM. However, since the Planck constant [??] is very important in the study of tunneling radiation, we will explicitly write it down here in the field functions. In other words, [h.sub.0] and hi of (12) are now rewritten as [h.sub.0] (r, t) = exp(-i[omega]t/[??]) [h.sub.0] (r), [h.sub.1] (r, t) = exp(-i[omega]t/[??]) [h.sub.1] (r), and (21) becomes

[mathematical expression not reproducible]. (39)

The corresponding Schrodinger-type equation can be expressed as

[d.sup.2][PHI]/d[r.sup.2.sub.*] + ([[omega].sup.2]/[[??].sup.2] - V(r)) [PHI] = 0. (40)

Since the Hawking radiation reflects BH's radial properties in the vicinity of the horizon, it is worthwhile to explicitly discuss radial field equations. According to (40), near the horizons r [right arrow] [r.sub.h], one has [V.sub.o] [right arrow] 0 and [V.sub.e] [right arrow] 0, so that the field equation can be simplified to

[d.sup.2][PHI](r)/d[r.sup.2.sub.*] + [[omega].sup.2]/[[??].sup.2] [PHI](r) = 0. (41)

By introducing [PHI](r) = and using the semiclassical approximation to neglect the terms of O([??]), as h is small, one obtains

-[f.sup.2] [R'.sup.2] + [[omega].sup.2] = 0. (42)

Therefore, R'(r) = [+ or -] [omega]/f, where [+ or -] represents ingoing or outgoing mode. The above equation is no other than the Hamilton-Jacobi equation at the horizon. We note that the above argument also applies to the case of asymptotically dS spacetime, where the event horizon [r.sub.h] and cosmological horizon [r.sub.c] should be both considered. In fact, according to the Hamilton-Jacobi equation [g.sup.[mu]v] ([partial derivative]S/[partial derivative][x.sup.[mu]])([partial derivative]S/[partial derivative][x.sup.v]) + [m.sup.2] = 0, the radial Hamilton-Jacobi equation can be transformed into (42) at the horizon, so we can use this result to study the Hawking tunneling radiation under semiclassical approximation [54-56].

Owing to the coordinate singularity at the horizon, the integration from inner surface of horizon to outer surface shall be carried out by using Residue Theorem:

R([r.sub.h]) = [+ or -] i[pi][omega]/f'([r.sub.h]). (43)

Then, the tunneling rate of Hawking radiation is [57, 58]

[GAMMA] = exp (-2 Im ([R.sub.+]))/exp(-2 Im ([R.sub.-])) = exp (- 4[pi][omega]/f'([r.sub.h])). (44)

According to the relationship between tunneling rate T and the Hawking temperature, the temperature of the BH in the vicinity of the event horizon [r.sub.h] is

[T.sub.H] = f'([r.sub.h])/4[pi]. (45)

From the inside of the cosmological horizon [r.sub.c] in dS spacetime, we consider the incident wave solution of Hawking radiation. Therefore, there is an extra minus sign in the resulting expression of the Hawking temperature in the neighborhood of the cosmological horizon [59] as follows:

[T.sub.c] = -f'([r.sub.c])/4[pi]. (46)

Figure 13 shows the impact of the parameter b on the thermodynamics of a regular phantom BH for different spacetimes. In asymptotically flat spacetime, the Hawking temperature [T.sub.H] decreases monotonously with increasing b and a given [r.sub.h]; and for a given b, [T.sub.H] becomes smaller with increasing [r.sub.h]. In asymptotically AdS spacetime, however, the Hawking temperature [T.sub.H] increases monotonously with increasing b and a given [r.sub.h], and for a given b, [T.sub.H] becomes larger with increasing [r.sub.h]. In asymptotically dS spacetime, the relationship between [T.sub.H] and the event horizon [r.sub.h] is similar to that in asymptotically flat case, while the Hawking temperature [T.sub.C] in the vicinity of the cosmological horizon [h.sub.c] also decreases monotonously with increasing b, and it decays faster than other cases. The different dependence of Hawking temperature on b probably results from the distinctive spacetime structure of AdS phantom BH from others.

5. Conclusions and Remarks

We studied static spherically symmetric solutions of regular phantom BH in asymptotically flat, dS, and AdS spacetimes and then investigated the QNMs of gravitational perturbations as well as Hawking radiations for these BH spacetimes. In our calculations, besides the metric perturbation [delta][g.sub.[mu]v] of the Ricci curvature tensor and scalar curvature, its effect on the energy momentum tensor of the matter field is also considered. In the derivation of the master equation for even parity, we made use of the method proposed in [26]. However, the BH metric in this paper is a self-gravitating solution of a minimally coupled scalar field with an arbitrary potential rather than the Lovelock equations, so that it does not satisfy the relation of Eq. (5.4) in [26]. As a result, the obtained master equation is more complicated. It is found that the calculated effective potential of AdS spacetime in the limit r [right arrow] [infinity] is very different from those of asymptotically flat and dS spacetimes. For the asymptotical AdS spacetime, the effective potentials diverge at infinity, which implies that the wave function [PHI] should vanish in this limit [60-64]. At the outside of the event horizon, on the other hand, the corresponding wave function must be an incoming wave. The distinct nature of AdS BH spacetime leads to the fact that some traditional numerical methods such as WKB approximation, continuous fraction method cease to be valid. In our calculations, we therefore employed the finite difference method to numerically calculate the temporal evolution of small gravitational perturbations. The relationships between Hawking temperature and parameters such as b, [r.sub.h] in different spacetime are also studied.

Owing to the importance of the parameter b, which carries the physical content of the coupling strength between phantom scalar field and the gravity, we studied the dependence of various physical quantities on this specific parameter, among others. From the above calculated results, we draw the following conclusions. For odd parity, in asymptotic flat and dS spacetime, the perturbation frequency (Re(o>)) becomes smaller and the decay rate ([absolute value of (Im([omega]))]) decreases when b increases. The perturbation frequency (Re([omega])) and the decay rate ([absolute value of (Im([omega]))]) both increase with increasing l. For even parity, the shape of the effective potential in the vicinity of the event horizon is very different from that for odd parity. In particular, [V.sub.e] (r) decreases monotonically in the range r > [r.sub.h], which is caused by the metric transformation F(r) = f(r)/Q(r). However, our numerical calculations show that the resultant gravitational perturbations are stable in all three different spacetime backgrounds. The Hawking temperature [T.sub.H] (and [T.sub.c]) of the phantom BH decreases with increasing b in asymptotically flat and dS spacetimes. On the contrary, in AdS spacetime, the temperature [T.sub.H] increases with increasing b. These results reflect distinct natures of the horizons of different spacetimes.

Recently, another type of regular phantom BH, which is a solution of Einstein-Maxwell equations, was proposed by Lemos and Zanchin [65]. The properties of interior as well as exterior structures were investigated. It is meaningful to investigate the stability of such BH solution and compare it with the phantom BH in the present paper.

Appendix

We note that V([phi](r)) in (7) is just the 0-order solution. However, when we investigate the perturbation in the BH background, the potential could, in principle, be expanded to higher order terms. Since the higher order terms of V([phi]) have not been determined yet, we can freely choose the form of such higher order terms to cancel out the contributions of the perturbation of phantom scalar field. In the Appendix, we carry out explicit calculations to show that this is indeed possible.

We consider the following metric and phantom scalar field:

[g.sub.[mu]v] = [[bar.g].sub.[mu]v] + [h.sub.[mu]v], [phi] = [bar.[phi]] + [[phi].sub.1], (A.1)

where [[bar.g].sub.[mu]v], [bar.[phi]] are background metric and phantom scalar field, respectively, and represents the first order perturbation of phantom scalar field. Then

[g.sup.[mu]v] = [[bar.g].sup.[mu]v] - [h.sup.[mu]v] + [h.sup.[mu].sub.[alpha]] [h.sup.v[alpha]]. (A.2)

The corresponding metric and potential of phantom scalar field can be expressed as

[mathematical expression not reproducible], (A.3)

where the superscript - and subscripts 1, 2 represent the background quantities and the first and the second order perturbations, respectively. Therefore, the action of the gravity with the phantom scalar field, (1), can be expanded as

S = [bar.S] + [S.sub.1] + [S.sub.2] + ..., (A.4)

where

[mathematical expression not reproducible]. (A.5)

In the background equation, there is no effect of [V.sub.1] and [V.sub.2], while in the perturbed equation we can always choose the form of the perturbed potential [V.sub.1], [V.sub.2] as

[mathematical expression not reproducible], (A.6)

to cancel out the contributions from the phantom scalar field perturbation [[phi].sub.1]. Therefore, in the text, we do not explicitly consider the perturbation of the phantom field.

http://dx.doi.org/10.1155/2017/5234214

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by Natural Science Foundation of China Grants no. 11205254, no. 11178018, no. 11573022, no. 11605015, and no. 11375279, the Fundamental Research Funds for the Central Universities 106112016CDJXY300002 and 106112015CDJRC131216, Chongqing Postdoctoral Science Foundation (no. Xm2015027), and the Open Project Program of State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, China (no. Y5KF181CJ1). The authors also acknowledge the financial support by Brazilian foundations Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP), Fundacao de Amparo a Pesquisa do Estado de Minas Gerais (FAPEMIG), Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq), and Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior (CAPES) and the Chinese State Scholarship.

References

[1] U. Alam, V. Sahni, T. D. Saini, and A. A. Starobinsky, "Is there supernova evidence for dark energy metamorphosis?" Monthly NoticesoftheRoyal Astronomical Society, vol. 354, no. 1, pp. 275-291, 2004.

[2] E. Komatsu, "Seven-year wilkinson microwave anisotropy probe (WMAP) observations: cosmological interpretation," The Astrophysical Journal Supplement Series, vol. 192, article 18, 2011.

[3] S. W. Allen, R. W. Schmidt, H. Ebeling, A. C. Fabian, and L. Van Speybroeck, "Constraints on dark energy from Chandra observations of the largest relaxed galaxy clusters," Monthly Notices of the Royal Astronomical Society, vol. 353, no. 2, pp. 457-467, 2004.

[4] M. Sullivan, J. Guy, A. Conleyet al., "SNLS3: constraints on dark energy combining the supernova legacy survey three-year data with other probes," The Astrophysical Journal, vol. 737, no. 2, article 102, 2011.

[5] S. Hannestad and E. Mortsell, "Probing the dark side: constraints on the dark energy equation of state from CMB, large scale structure, and type Ia supernovae," Physical Review D, vol. 66, no. 6, Article ID 063508, 2002.

[6] R. R. Caldwell, "A phantom menace? Cosmological consequences of a dark energy component with super-negative equation of state," Physics Letters B, vol. 545, no. 1-2, pp. 23-29, 2002.

[7] A. Melchiorri, L. Mersini, C. J. Odman, and M. Trodden, "The state of the dark energy equation of state," Physical Review D, vol. 68, no. 4, Article ID 043509, 2003.

[8] J. Weller and A. M. Lewis, "Large-scale cosmic microwave background anisotropies and dark energy," Monthly Notices of the Royal Astronomical Society, vol. 346, no. 3, pp. 987-993, 2003.

[9] R. Gannouji, D. Polarski, A. Ranquet, and A. A. Starobinsky, "Scalar-tensor models of normal and phantom dark energy," Journal of Cosmology and Astroparticle Physics, vol. 2006, no. 9, p. 16, 2006.

[10] A. Vikman, "Can dark energy evolve to the phantom?" Physical Review D, vol. 71, no. 2, Article ID 023515, 14 pages, 2005.

[11] R. R. Caldwell, M. Kamionkowski, and N. N. Weinberg, "Phantom energy: dark energy with w < - 1 causes a cosmic doomsday," Physical Review Letters, vol. 91, no. 7, Article ID 071301, 2003.

[12] Z.-Y. Sun and Y.-G. Shen, "Phantom cosmology with non-minimally coupled real scalar field," General Relativity and Gravitation, vol. 37, no. 1, pp. 243-251, 2005.

[13] L. P. Chimento and R. Lazkoz, "Constructing phantom cosmologies from standard scalar field universes," Physical Review Letters, vol. 91, no. 21, p. 211301, 2003.

[14] I. Maor, R. Brustein, J. McMahon, and P. J. Steinhardt, "Measuring the equation of state of the universe: pitfalls and prospects," Physical Review D, vol. 65, no. 12, Article ID 123003, 2002.

[15] J.-M. Virey, P. Taxil, A. Tilquin, A. Ealet, D. Fouchez, and C. Tao, "Probing dark energy with supernovae: bias from the time evolution of the equation of state," Physical Review D, vol. 70, no. 4, Article ID 043514, 2004.

[16] K. A. Bronnikov and J. C. Fabris, "Regular phantom black holes," Physical Review Letters, vol. 96, no. 25, Article ID 251101, 2006.

[17] A. Flachi and J. P. S. Lemos, "Quasinormal modes of regular black holes," Physical Review D, vol. 87, Article ID 024034, 2013.

[18] J. S. F. Chan and R. B. Mann, "Scalar wave falloff in asymptotically anti-de Sitter backgrounds," Physical Review D, vol. 55, no. 12, pp. 7546-7562, 1997.

[19] J. S. F. Chan and R. B. Mann, "Scalar wave falloff in topological black hole backgrounds," Physical Review D, vol. 59, no. 6, Article ID 064025, 23 pages, 1999.

[20] K. A. Bronnikov, R. A. Konoplya, and A. Zhidenko, "Instabilities of wormholes and regular black holes supported by a phantom scalar field," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 86, no. 2, Article ID 024028, 2012.

[21] K. Lin, J. Li, and S. Z. Yang, "Quasinormal modes of Hayward regular black hole," International Journal of Theoretical Physics, vol. 52, no. 10, pp. 3771-3778, 2013.

[22] S. Fernando and T. Clark, "Blackholes in massive gravity: quasinormal modes of scalar perturbations," General Relativity and Gravitation, vol. 46, no. 12, article, 2014.

[23] C. F. B. Macedo and L. C. B. Crispino, "Absorption of planar massless scalar waves by Bardeen regular black holes," Physical Review D, vol. 90, no. 6, Article ID 064001, 2014.

[24] J. Li, H. Ma, and K. Lin, "Dirac quasinormal modes in spherically symmetric regular black holes," Physical Review D --Particles, Fields, Gravitation and Cosmology, vol. 88, no. 6, Article ID 064001, 2013.

[25] J. Li, K. Lin, and N. Yang, "Nonlinear electromagnetic quasinormal modes and Hawking radiation of a regular black hole with magnetic charge," The European Physical Journal C, vol. 75, no. 3, article 131, 2015.

[26] T. Takahashi and J. Soda, "Master equations for gravitational perturbations of static Lovelock black holes in higher dimensions," Progress of Theoretical Physics, vol. 124, no. 5, pp. 911-924, 2010.

[27] T. Regge and J. A. Wheeler, "Stability of a Schwarzschild singularity," Physical Review, vol. 108, no. 4, pp. 1063-1069, 1957.

[28] F. J. Zerilli, "Effective potential for even-parity Regge-Wheeler gravitational perturbation equations," Physical Review Letters, vol. 24, no. 13, pp. 737-738, 1970.

[29] V. Moncrief, "Gravitational perturbations of spherically symmetric systems. II. Perfect fluid interiors," Annals of Physics, vol. 88, pp. 343-370, 1974.

[30] T. Takahashi and J. Soda, "Hawking radiation from fluctuating black holes," Classical and Quantum Gravity, vol. 27, no. 17, Article ID 175008, 39 pages, 2010.

[31] S. Chandrasekar, The Mathematical Theory of Black Holes, Oxford University Press, Oxford, UK, 1983.

[32] S. A. Teukolsky, "Perturbations of a rotating black hole. I. Fundamental equations for gravitational, electromagnetic, and neutrino-field perturbations," The Astrophysical Journal, vol. 185, pp. 635-647, 1973.

[33] B. P. Abbott, R. Abbott, T. D. Abbott et al., "Observation of gravitational waves from a binary black hole merger," Physical Review Letters, vol. 116, no. 6, Article ID 061102, 16 pages, 2016.

[34] R. Konoplya and A. Zhidenko, "Detection of gravitational waves from black holes: is there a window for alternative theories?" Physics Letters B, vol. 756, pp. 350-353, 2016.

[35] C. Corda, "Effective temperature, Hawking radiation and quasi-normal modes," International Journal of Modern Physics D, vol. 21, no. 11, Article ID 1242023, 11 pages, 2012.

[36] C. Corda, "Effective temperature for black holes," Journal of High Energy Physics, vol. 2011, no. 8, article 101, 2011.

[37] C. Corda, S. H. Hendi, R. Katebi, and N. O. Schmidt, "Effective state, Hawking radiation and quasi-normal modes for Kerr blackholes," Journal of High Energy Physics, vol. 2013, article 8, 2013.

[38] C. Corda, S. H. Hendi, R. Katebi, and N. O. Schmidt, "Hawking radiation-quasi-normal modes correspondence and effective states for nonextremal Reissner-Nordstrom black holes," Advances in High Energy Physics, vol. 2014, Article ID 527874, 9 pages, 2014.

[39] C. Corda, "Black hole quantum spectrum," The European Physical Journal C, vol. 73, article 2665, 2013.

[40] C. Corda, "Quantum transitions of minimum energy for Hawking quanta in highly excited black holes: problems for loop quantum gravity?" Electronic Journal of Theoretical Physics, vol. 11, no. 30, pp. 27-34, 2014.

[41] C. Corda, "Non-strictly black body spectrum from the tunnelling mechanism," Annals of Physics, vol. 337, pp. 49-54, 2013.

[42] C. Gundlach, R. H. Price, and J. Pullin, "Late-time behavior of stellar collapse and explosions. I. Linearized perturbations," Physical Review D, vol. 49, no. 2, pp. 883-889, 1994.

[43] K. Lin, J. Li, and N. Yang, "Dynamical behavior and nonminimal derivative coupling scalar field of Reissner-Nordstroom black hole with a global monopole," General Relativity and Gravitation, vol. 43, no. 6, pp. 1889-1899, 2011.

[44] J. Li and Y. H. Zhong, "Quasinormal modes for electromagnetic field perturbation of the asymptotic safe black hole," International Journal of Theoretical Physics, vol. 52, no. 5, pp. 1583-1587, 2013.

[45] B. F. Schutz and C. M. Will, "Black hole normal modes--a semianalytic approach," The Astrophysical Journal, vol. 291, p. L33, 1985.

[46] S. Iyer and C. M. Will, "Black-hole normal modes: a WKB approach. I. Foundations and application of a higher-order WKB analysis of potential-barrier scattering," Physical Review D, vol. 35, no. 12, pp. 3621-3631, 1987.

[47] Y. Zhang, E.-K. Li, and J.-L. Geng, "Quasinormal modes of massless scalar field perturbation of a Reissner-Nordstroom de Sitter blackhole with a global monopole," General Relativity and Gravitation, vol. 46, no. 5, article 1728, pp. 1-11, 2014.

[48] R. A. Konoplya, "Gravitational quasinormal radiation of higher-dimensional black holes," Physical Review. D. Third Series, vol. 68, no. 12, 124017, 7 pages, 2003.

[49] P. Kraus and F. Wilczek, "Self-interaction correction to black hole radiance," Nuclear Physics B, vol. 433, no. 2, pp. 403-420, 1995.

[50] M. K. Parikh and F. Wilczek, "Hawking radiation as tunneling," Physical Review Letters, vol. 85, no. 24, pp. 5042-5045, 2000.

[51] K. Srinivasan and T. Padmanabhan, "Particle production and complex path analysis," Physical Review D, vol. 60, no. 2, Article ID 024007, 20 pages, 1999.

[52] K. Lin and S.-Z. Yang, "Quantum tunnelling in charged black holes beyond the semi-classical approximation," Europhysics Letters, vol. 86, no. 2, Article ID 20006, 2009.

[53] K. Lin and S. Z. Yang, "Fermion tunnels of cylindrical symmetric black hole and the corrected entropy," International Journal of Theoretical Physics, vol. 48, no. 10, pp. 2920-2927, 2009.

[54] K. Lin and S.-Z. Yang, "Fermion tunneling from higher-dimensional black holes," Physical Review D, vol. 79, no. 6, Article ID 064035, 2009.

[55] K. Lin and S. Z. Yang, "Fermions tunneling of higher-dimensional Kerr-anti-de Sitter black hole with one rotational parameter," Physics Letters B, vol. 674, no. 2, pp. 127-130, 2009.

[56] K. Lin and S. Yang, "Fermion tunnels of higher-dimensional anti-de Sitter SCHwarzschild black hole and its corrected entropy," Physics Letters. B, vol. 680, no. 5, pp. 506-509, 2009.

[57] R. Di Criscienzo, L. Vanzo, and S. Zerbini, "Applications of the tunneling method to particle decay and radiation from naked singularities," Journal of High Energy Physics, vol. 2010, article 92, 2010.

[58] L. Vanzo, G. Acquaviva, and R. Di Criscienzo, "Tunnelling methods and Hawking's radiation: achievements and prospects," Classical and Quantum Gravity, vol. 28, no. 18, Article ID 183001, 2011.

[59] S. Z. Yang and K. Lin, "Tunneling radiation and entropy correction of arbitrary spin particles in Kerr-de Sitter black hole," The Chinese Physical Society, vol. 59, no. 8, 2010.

[60] J. Maldacena, "The large N limit of superconformal field theories and supergravity," Advances in Theoretical and Mathematical Physics, vol. 2, no. 2, pp. 231-252, 1998.

[61] E. Witten, "Anti de Sitter space and holography," Advances in Theoretical and Mathematical Physics, vol. 2, no. 2, pp. 253-291, 1998.

[62] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, "Gauge theory correlators from non-critical string theory," Physics Letters B, vol. 428, no. 1-2, pp. 105-114, 1998.

[63] O. Aharony, S. S. Gubser, J. Maldacena, H. Ooguri, and Y. Oz, "Large N field theories, string theory and gravity," Physics Reports, vol. 323, no. 3-4, pp. 183-386, 2000.

[64] G. T. Horowitz and V. E. Hubeny, "Quasinormal modes of AdS black holes and the approach to thermal equilibrium," Physical Review D, vol. 62, no. 2, Article ID 024027, 11 pages, 2000.

[65] J. P. S. Lemos and V. T. Zanchin, "Regular black holes: Guilfoyle's electrically charged solutions with a perfect fluid phantom core," Physical Review D, vol. 93, no. 12, Article ID 124012, 2016.

Jin Li, (1, 2) Kai Lin, (3, 4) Hao Wen, (1) and Wei-Liang Qian (5, 6)

(1) College of Physics, Chongqing University, Chongqing 401331, China

(2) State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China

(3) Instituto de Fisica e Quimica, Universidade Federal de Itajuba, 37500-903 Itajuba, MG, Brazil

(4) Department of Astronomy China West Normal University, Nanchong, Sichuan 637002, China

(5) Escola de Engenharia de Lorena, Universidade de Sao Paulo, Sao Paulo, SP, Brazil

(6) Faculdade de Engenharia de Guaratingueta, Universidade Estadual Paulista, Guaratingueta, SP, Brazil

Correspondence should be addressed to Jin Li; cqstarv@hotmail.com

Received 24 August 2016; Accepted 9 October 2016; Published 16 February 2017

Academic Editor: Christian Corda

Caption: FIGURE 1: The structures of the regular phantom BH's metric function f(r) for different values of b: (a) asymptotically flat spacetime with [r.sub.h] = 1. (b) de Sitter spacetime with [r.sub.h] = 1, [r.sub.c] = 10; (c) anti-de Sitter spacetime with [r.sub.h] = 1.

Caption: FIGURE 2: The effective potential for odd parity gravitational perturbation in asymptotically flat spacetime, with [r.sub.h] = 1, l = 2 (a) and with [r.sub.h] = 1, b = 5 (b).

Caption: FIGURE 3: The temporal evolution of odd parity gravitational perturbation in asymptotically flat phantom BH, with [r.sub.h] = 1, l = 2 (a) and with [r.sub.h] =1, b = 5 (b).

Caption: FIGURE 4: Calculated quasinormal frequencies of odd parity gravitational perturbation for [r.sub.h] = 1 in asymptotically flat spacetime. Each group of dots, from top-right to bottom-left, corresponds to QNM frequencies obtained by assuming different values of b = 0.3,0.4,0.5,0.7,0.8,0.9,1.0,1.2, respectively.

Caption: FIGURE 5: The effective potential for odd parity gravitational perturbation in dS spacetime with [r.sub.h] = 1, [r.sub.c] = 10, l = 2 (a) and with [r.sub.h] = 1, [r.sub.c] = 10, b = 5 (b).

Caption: FIGURE 6: The temporal evolution of odd parity gravitational perturbation for dS phantom BH with [r.sub.h] = 1, l = 2 (a), and with [r.sub.h] = 1, b = 5 (b).

Caption: FIGURE 7: Calculated quasinormal frequencies of odd parity gravitational perturbation with [r.sub.h] = 1, [r.sub.c] = 10 in dS spacetime. Each group of dots, from top-right to bottom-left, corresponds to QNM frequencies obtained by assuming different values of b = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8, respectively.

Caption: FIGURE 8: The effective potential for odd parity gravitational perturbation in AdS spacetime with [r.sub.h] =1, l = 2 (a) and with [r.sub.h] = 1, b = 1 (b).

Caption: FIGURE 9: The temporal evolution of odd parity gravitational perturbation in AdS phantom BH spacetime with [r.sub.h] = 1, l = 2 (a) and with [r.sub.h] = 1, b = 1 (b).

Caption: FIGURE 10: (a) The effective potential of even parity perturbation in asymptotically flat spacetime. (b) The corresponding temporal evolution of even parity gravitational perturbation field in asymptotically flat spacetime. Calculations are carried out by using [r.sub.h] = 1.

Caption: FIGURE 11: (a) The effective potential of even parity perturbation in dS spacetime. (b) The corresponding temporal evolution of even parity gravitational perturbation. In the calculations, we choose [r.sub.h] = 1, [r.sub.c] = 10.

Caption: FIGURE 12: (a) The effective potential of even parity perturbation in AdS spacetime. (b) The corresponding temporal evolution of even parity gravitational perturbation in AdS spacetime. In the calculations, we choose [r.sub.h] = 1.

Caption: FIGURE 13: The Hawking temperature as a function of b in asymptotically flat and dS spacetimes.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Li, Jin; Lin, Kai; Wen, Hao; Qian, Wei-Liang |

Publication: | Advances in High Energy Physics |

Date: | Jan 1, 2017 |

Words: | 8024 |

Previous Article: | Neutrino Emission from Magnetized Microquasar Jets. |

Next Article: | Top Partner Production at [e.sup.+][e.sup.-] Collider in the Littlest Higgs Model with T-Parity. |

Topics: |