Dirac Quasinormal Modes of New Type Black Holes in New Massive Gravity

We study new type black holes in three-dimensional New Massive Gravity and we calculate analytically the quasinormal modes for fermionic perturbations for some special cases. Then, we show that for these cases the new type black holes are stable under fermionic field perturbations.


I. INTRODUCTION
In recent years, there has been a remarkable interest in the study of three-dimensional models of gravity. Apart of BTZ black hole [1], which is a solution to the Einstein equations with a negative cosmological constant, remarkable attention was addressed to Topologically Massive Gravity (TMG), which is a generalization of three-dimensional GR that amounts to augment the Einstein-Hilbert action adding a Chern-Simons gravitational term, [2]. Here, the propagating degree of freedom is a massive graviton. TMG also admits the BTZ (and other) black holes as exact solutions. The renewed interest on TMG relies on the possibility of constructing a chiral theory of gravity at a special point of the space of parameters, [3].
On the other hand, Bergshoeff, Hohm and Townsend (BHT) introduced another three-dimensional massive gravity theory, which is known as New Massive Gravity (NMG), being the action the standard Einstein-Hilbert term with a specific combination of scalar curvature square term and Ricci tensor square one [4][5][6][7][8][9], and it is equivalent at the linearized level to the (unitary) Fierz-Pauli action for a massive spin-2 field [4]. The model in three dimensions is indeed unitary in the tree-level, but the corresponding model in higher dimensions is not so due to the appearance of non-unitary massless spin-2 modes [10]. Also, NMG admits warped AdS black holes [11], AdS waves [12,13], asymptotically Lifshitz black holes [14], gravitational solitons, kinks and wormholes [15], for further aspects of the BHT theory see [16][17][18][19][20]. Besides, asymptotically AdS and Lifshitz black holes in NMG dressed by a (non)minimally coupled scalar field have been studied recently in [21]. It is worth mention that TMG and NMG share common features, however there are different aspects, one of them is the existence of the new type of black holes, for a specific combination of parameters in the NMG Lagrangian, which was discovered by BHT, and which are known as new type black holes.
The particular motivation of this work is to calculate the quasinormal modes (QNMs) for fermionic field perturbations in the background of the new type black holes in three-dimensional New Massive Gravity and study the stability of these black holes under fermionic perturbations. The QNMs and their quasinormal frequencies (QNFs) are an important property of black holes and have a long history, [22][23][24][25][26][27]. It is known that the presence of event horizons dampens the vibration modes of a matter field that evolves perturbatively in the exterior region. In this way, the system is intrinsically dissipative, i.e., there is no temporary symmetry. In general, the oscillation frequencies are complex, therefore the system is not Hermitian. Nevertheless, the oscillation frequency of these modes is independent of the initial conditions and it only depends on the parameters (mass, charge and angular momentum) and the fundamental constants (Newton constant and cosmological constant) that describe a black hole just like the parameters that define the test field. In three-dimensional spacetime, the QNMs of the BTZ black hole have been studied in [28]- [30], and the QNMs for scalar field perturbations in the background of new type black holes in NMG was studied in [31].
The QNMs give information about the stability of black holes under matter fields that evolves perturbatively in the exterior region of them, without backreacting on the metric. Also, the QNMs determine how fast a thermal state in the boundary theory will reach thermal equilibrium according to the AdS/CFT correspondence [32], where the relaxation time of a thermal state of the boundary thermal theory is proportional to the inverse of the imaginary part of the QNMs of the dual gravity background [33]. In the context of black hole thermodynamics, the QNMs allow the quantum area spectrum of the black hole horizon to be studied, as well as the mass and the entropy spectrum. In this regard, Bekenstein [34] was the first to propose the idea that in quantum gravity the area of black hole horizon is quantized, leading to a discrete spectrum which is evenly spaced. Then, Hod [35] conjectured that the asymptotic QNF is related to the quantized black hole area, by identifying the vibrational frequency with the real part of the QNFs. However, it is not universal for every black hole background. Then, Kunstatter [36] propose that the black hole spectrum can be obtained by imposing the Bohr-Sommerfeld quantization condition to an adiabatic invariant quantity involving the energy and the vibrational frequency. Furthermore, Maggiore [37] argued that in the large damping limit the identification of the vibrational frequency with the imaginary part of the QNF could lead to the Bekenstein universal bound. Then, the consequences of these proposals were studied in several spacetimes, for instance, see [31], where the authors comment on AdS/CFT correspondence and entropy/area spectrum for the new type black holes. Besides, in [38]- [41] the authors discuss a connection between Hawking radiation and black hole quasinormal modes, which is important in the route to quantize gravity, because one can naturally interpret black hole quasinormal modes in terms of quantum levels.
The paper is organized as follows. In Sec. II we give a brief review of the new type black holes in three-dimensional New Massive Gravity. In Sec. III we calculate the exact QNMs of fermionic perturbations for the new type black holes. Finally, conclusions are presented in Sec. IV.

II. NEW TYPE BLACK HOLES
In this work we will consider a matter distribution outside the event horizon of three-dimensional black holes, which are solutions of NMG and conformal gravity in three dimensions [42]. The action of NMG theory is given by where K is given by the following quadratic terms of the curvature where R is the Ricci scalar, and λ is the cosmological constant. At the special case µ 2 = λ the theory admits the following metric as solution where f (r) = Ar 2 + Br + C. In the case that B = 0 there is a curvature singularity due to the Ricci scalar diverges at r = 0 The metric (3) is conformally flat, and hence solve the field equations of conformal gravity in three dimensions too, as we mentioned. The parameter A is proportional to the cosmological constant λ, B is a kind of "gravitational hair" and C is related to the mass of the spacetime. Note that depending of the values of the parameters A, B and C this metric may represents a black hole or not, being the black holes solutions known as new type black holes. The roots of f (r) are for r + > 0 the metric represents new type black holes with an event horizon located at r + . Also, in some cases may exist a cosmological horizon r − , with r − > r + . For A > 0, the new type black holes are asymptotically anti-de Sitter, for A < 0 are asymptotically de Sitter and for A = 0 are asymptotically locally flat. In the case A = 0 the black holes are solutions of the theory that only contains the K term in the Lagrangian. Additionally, for A > 0 and B = 0 this metric reduces to the non-rotating BTZ black hole.

III. DIRAC QUASINORMAL MODES OF NEW TYPE BLACK HOLES
The QNMs of fermionic perturbations in the background of new type black hole are given by the fermionic field solution to the Dirac equation in curved spacetime with suitable boundary conditions. This means that there is only ingoing waves on the event horizon and we will consider that the fermionic field vanishes at spatial infinity for asymptotically AdS black holes, known as Dirichlet boundary conditions, and that there are only outgoing waves at the cosmological horizon or at spatial infinity for asymptotically dS black holes and for asymptotically locally flat black holes, respectively. The Dirac equation is given by where the covariant derivative is defined as and the generators of the Lorentz group J ab are The gamma matrices in curved spacetime γ µ are defined by where γ a are the gamma matrices in flat spacetime. In order to solve the Dirac equation we use the diagonal vielbein and, from the null torsion condition we obtain the spin connection Now, using the following representation of the gamma matrices where σ i are the Pauli matrices, and using the following ansatz for the fermionic field we arrive at the following coupled system of differential equations where m is the mass of the fermionic field ψ, which is minimally coupled to curvature. Decoupling the above system of equations we obtain the following equation for ψ 1 Now, in order to obtain analytical solutions we will consider some special cases.

A. Null angular momentum
In this section we will consider the case κ = 0. So, making the change of variables y = 1 − r+ r , the equation (16) We note that this equation correspond to a Riemann's differential equation, whose general form is [43] where r, s, t are the singular points, and the exponents α, α ′ , β, β ′ , γ, γ ′ are subject to the condition The complete solution of (18) is denoted by the symbol where the P symbol denotes the Riemann's P function, which can be reduced to the hypergeometric function through where the P function is now the Gauss' hypergeometric function.

Asymptotically AdS New Type Black Holes
In this case r + > r − (A > 0). So, considering equations (17) and (18) we can identify the regular singular points r, s and t as Therefore, the exponents are given by and the solution to equation (17) can be written as where we have defined the constants a, b and c as In the near horizon limit, the above expression behaves as Now, we impose as boundary condition that classically nothing can escape from the event horizon. Note that Q/A (Q − 1) > 0, therefore we must take C 2 = 0 in order to have only ingoing waves at the horizon. So, the solution simplifies to Now, we implement boundary conditions at spatial infinity y → 1. In order to do this, we employ the Kummer's relations [43], which allow us to write the solution as In the limit y → 1, the above expression becomes So, imposing that the scalar field be null at spatial infinity, we can determine the quasinormal frequencies. For m > 0 the second term of equation (31) blows up when y → 1 unless we impose the condition a = −n or b = −n, these conditions give the following QNFs Due to the imaginary part of the QNFs is negative the asymptotically AdS new type black holes are stable under fermionic perturbations, at least for the mode with the lowest angular momentum. Also, in a similar way the quasinormal frequencies associated to ψ 2 can be obtained, note that ψ 2 satisfies a similar equation that ψ 1 but making the changes κ → −κ and ω → −ω.

Asymptotically dS New Type Black Holes
In this case, besides the event horizon r + , we have a cosmological horizon r − , where r − > r + (A < 0). So, considering equations (17) and (18) we can identify the regular singular points r, s and t as Therefore, the exponents are given by and the solution to equation (17) can be written as where we have defined the constants a, b and c as In the near horizon limit, the above expression behaves as Now, we impose as boundary condition that classically nothing can escape from the event horizon. Note that Q/A (Q − 1) > 0 as in the case 1, therefore we must take C 2 = 0 in order to have only ingoing waves at the horizon. So, the solution simplifies to Now, we implement boundary conditions at the cosmological horizon y → 1 − Q. In order to do this, we employ the Kummer's relations [43], which allow us to write the solution as In the limit y → 1 − Q, the above expression becomes So, in order to have only outgoing waves at the cosmological horizon, we must impose c − a = −n or c − b = −n. These conditions yield the following set of quasinormal frequencies Due to the imaginary part of the QNFs is negative the asymptotically dS new type black holes are stable under fermionic perturbations, at least for the mode with the lowest angular momentum. As in the previous case, in a similar way the quasinormal frequencies associated to ψ 2 can be obtained, note that ψ 2 satisfies a similar equation that ψ 1 but making the changes κ → −κ and ω → −ω.

B. Massless Dirac QNMs of Asymptotically dS and Asymptotically Locally Flat New Type Black Holes
In this section we compute the QNMs for asymptotically dS (A < 0) and asymptotically locally flat (A = 0) new type black holes for massless fermionic perturbations (m = 0).

Asymptotically dS New Type Black Holes
Performing the change of variables y = 1 − r+ r , the equation (16) reduces to where we have defined Q = r + /r − < 1. It is worth mentioning that in these coordinates the event horizon is located at y = 0, and the cosmological horizon at y = 1 − Q. Now, making the change of variables we arrive at the following equation where α and β are given by In the next, we will use the values α = iQω Ar+(1−Q) and β = − iQω Ar+(1−Q) . Then, performing another change of variable (1 − Q) z = y we obtain and, in these new coordinates the event horizon is located at z = 0, and the cosmological horizon at z = 1. We recognize equation (48) as the hypergeometric equation where The general solution of the hypergeometric equation is which has three regular singular points at z = 0, z = 1 and z = ∞. Here, 2 F 1 (a, b, c; z) is a hypergeometric function and C 1 , C 2 are constants. Then, the solution for the radial function ψ 1 (z) is So, in the vicinity of the event horizon z = 0, and using the property 2 F 1 (a, b, c; 0) = 1, the function ψ 1 (z) behaves as ψ 1 (z) = C 1 e α ln z + C 2 e (1/2−α) ln z .
Now, imposing boundary conditions at the event horizon, that there is only ingoing modes, implies that C 2 = 0. Thus, the solution can be written as On the other hand, using Kummer's formula for hypergeometric functions, the radial function ψ 1 (z) can be written as So, in order to have only outgoing waves at the cosmological horizon z = 1, we must impose c − a = −n or c − b = −n. These conditions yield the following set of quasinormal frequencies In a similar way the quasinormal frequencies associated to ψ 2 can be obtained, note that ψ 2 satisfies a similar equation that ψ 1 but making the changes κ → −κ and ω → −ω.

Asymptotically Locally Flat New Type Black Holes
Under the change of variables y = 1 − r+ r , the equation (16) becomes Now, making the change of variables we arrive at the following equation where α and β are given by In the next, we will use the values α = − iω B and β = iω B . We recognize equation (60) as the hypergeometric equation where As in the previous case, the general solution of the hypergeometric equation is which has three regular singular points at y = 0, y = 1 and y = ∞. Here, 2 F 1 (a, b, c; y) is a hypergeometric function and C 1 , C 2 are constants. Then, the solution for the radial function ψ 1 (y) is So, in the vicinity of the event horizon, y = 0 and using the property 2 F 1 (a, b, c; 0) = 1, the function ψ 1 (y) behaves as Now, imposing boundary conditions at the horizon, that there is only ingoing modes, implies that C 2 = 0. Thus, the solution can be written as On the other hand, using Kummer's formula for hypergeometric functions (55), the radial function ψ 1 (y) can be written as So, in order to have only outgoing waves at the spatial infinity y = 1, we must impose c − a = −n or c − b = −n. These conditions yield the following set of quasinormal frequencies

C. Massless Dirac QNMs of Extremal New Type Black Holes
In this section we consider the extremal case r + = r − (A > 0). First, we take κ = 0. So, using the change of variable y = 1 − r+ r as in the previous sections, equation (16) reduces to The solution is given in terms of Whittaker functions In order to have a regular scalar field at spatial infinity y → 1, we must set C 2 = 0. Therefore, the solution reduces to However, for y → 1 this expression becomes null, which shows the absence of QNMs. Furthermore, for m = 0 we obtain the equation whose solution is given in terms of confluent Heun functions Then, at infinity y → 1, the radial function vanishes. Therefore, there are not quasinormal modes for extremal new type black holes.

IV. CONCLUSIONS
In this work we have calculated analytically the QNMs of fermionic perturbations for some special cases for new type black holes, which are solutions of three-dimensional NMG, and also are solutions of conformal gravity in three dimensions. The first case that we have analyzed are massive fermionic fields perturbations without angular momentum (κ = 0) in the black hole background, and we have found the QNFs for asymptotically AdS (dS) new type black holes. For asymptotically AdS new type black holes the QNFs are purely imaginary and negative, which ensures the stability of the black hole under fermionic perturbations. For asymptotically dS new type black holes the QNFs have a real and imaginary part, and the imaginary part is negative, which ensures the stability of the black hole under fermionic perturbations. It is worth mentioning that in these cases the Dirac equation can be written as a Riemann differential equation, as in [44]. Other cases, where is possible to find the QNFs analytically is for asymptotically dS and asymptotically locally flat new type black holes for massless fermionic field perturbations, in these cases we have found that the QNFs have a real and imaginary part, and the imaginary part is negative, which ensures the stability of the black hole under fermionic perturbations.
Finally, we have analyzed fermionic field perturbations in the extremal new type black holes for some special cases, and we found that there are not quasinormal modes as occurs, for instance, in [44,45], where the authors have showed the absence of QNMs in the extremal BTZ black hole and the extremal four-dimensional Lifshitz Black Hole in Conformal Gravity. However, it was shown that it is possible to construct the QNMs of three-dimensional extremal black holes in an algebraic way as the descendents of the highest weight modes [46], with hidden conformal symmetry being an intrinsic property of the extremal black hole. Also, it is worth to mention that the absence of QNMs for extremal black holes does not always occur, for instance see [47], where the authors have showed the presence of QNMs for the extremal BTZ black holes in TMG.