Quasinormal modes of three-dimensional rotating Ho\v{r}ava AdS black hole and the approach to thermal equilibrium

We compute the quasinormal modes (QNMs) of a massive scalar field in the background of a rotating three-dimensional Ho\v{r}ava AdS black hole, and we analyze the effect of the breaking of the Lorentz invariance on the QNMs. Imposing on the horizon that there are only ingoing waves and at infinity the Dirichlet boundary conditions and the Neumann boundary condition, we calculate the oscillatory and the decay modes of the QNMs. We find that the propagation of the scalar field is stable in this background and employing the holographic principle we find the different times of the perturbed system to reach thermal equilibrium for the various branches of solutions.

poles of the retarded correlation function describing the linear response on the conformal field theory side.
In this work we will consider a matter distribution in the background of three-dimensional rotating Hořava AdS black holes [27], parameterized by a scalar field. We will perturb the scalar field assuming that there is no back reaction on the metric. This will result in the calculation of the QNMs, which are characterized by a spectrum that is independent of the initial conditions of the perturbation and depends only on the black hole and probe field parameters, and on the fundamental constants of the system (for a review see [28]). Exact solutions for QNMs of black holes in three spacetime dimensions have been obtained in [29].
The motivation for considering the Hořava gravity theory is twofold. Considering a condensed matter dynamical system, it was argued in [30] that this condensed matter system breaks Lorentz invariance spontaneously and its excitations, the superfluid's phonons, have to non-linearly realize the spontaneously broken Lorentz boosts forcing their interactions to have a very constrained structure. Then, to holographically describe such a system on the boundary, we need to have a Lorentz braking gravity theory in the bulk. The other motivation is, by calculating the QNMs, to see what is the effect of Lorentz breaking symmetry on the relaxation time of the dynamical system to reach thermal equilibrium on the boundary [31].
In this context, the QNMs for four-dimensional non-reduced Einstein-aether theory was studied, and it was found in [32] that the oscillation and damping rate of QNMs are larger than those of Schwarszschild black holes of the Einstein theory, for an effective potential that is known only numerically. More recently, the QNMs for two kinds of aether black holes were analyzed, and it was shown that quasinormal ringing of the first kind of aether black hole is similar to that of another Lorentz violation model-the QED-extension limit of standard model. Also, it was found in [33] that both the first and the second kind of aether black holes have larger damping rate and smaller real oscillation frequency of QNMs compared to Schwarzschild black hole.
The work is organized as follows. In Sec. II after giving a brief review of the BTZ black hole, we discuss the three-dimensional Hořava gravity and its connection with three-dimensional Einstein-aether theory. In Sec. III we find the QNMs analytically for massive scalar fields with circular symmetry and for a specify value of J. Also, for massive scalar field we show that the Klein-Gordon equation can be written as the Heun's equation, and we find the QNFs numerically by applying the pseudospectral method. Finally, our conclusions are in Sec. IV.

II. THREE-DIMENSIONAL ROTATING HOŘAVA BLACK HOLES
In this Section, after reviewing in brief the BTZ black hole, we discuss the Hořava gravity and its connection with three-dimensional Einstein-aether theory. The metric of the BTZ black hole is given by The angular coordinate φ has period 2π, and the radii of the inner and outer horizons are denoted by r − and r + , respectively. The dual conformal field theory on the boundary is (1 + 1)-dimensional, the conformal symmetry being generated by two copies of the Virasoro algebra acting separately on left-and right-moving sectors [24]. Consequently, the conformal field theory splits into two independent sectors at thermal equilibrium with temperatures According to the AdS 3 /CFT 2 correspondence, to each field of spin s propagating in AdS 3 there corresponds an operator O in the dual conformal field theory characterised by conformal weights (h L , h R ) with [23] and ∆ is determined in terms of the mass m of the scalar field, For a small perturbation, one expects that at late times the perturbed system will approach equilibrium exponentially with a characteristic time-scale. This time-scale is inversely proportional to the imaginary part of the poles, in momentum space, of the correlation function of the perturbation operator O. For a conformal field theory at zero temperature, the 2-point correlation functions can be determined, up to a normalisation, from conformal invariance. Then two sets of poles were found [24] where n takes the integer values (n = 0, 1, 2, ...). This set of poles characterises the decay of the perturbation on the CFT side, and coincides precisely with the quasi-normal frequencies of the BTZ black hole [24]. We now discuss the three-dimensional Hořava gravity, the action of which is given by [35] where G H is a coupling constant with dimensions of a length squared and the Lagrangian L 2 has the following form where K ij , K, and (2) R correspond to extrinsic, mean, and scalar curvature, respectively, and a i is a parameter related to the lapse function N via a i = ∂ i ln N , being the line element in the preferred foliation Also, g is the determinant of the induced metric g ij on the constant-T hypersurfaces. L 4 corresponds to a set of all the terms with four spatial derivatives that are invariant under diffeomorphisms. For λ = ξ = 1 and η = 0, the action reduces to that of General Relativity. In the infrared limit of the theory the higher order terms L 4 (UV regime) can be neglected, and the theory is equivalent to a restricted version of Einstein-aether theory, through and where c ij = c i + c j . Being the action of Einstein-aether: where G ae is a coupling constant with dimensions of a length square, g is the determinant of g µν , Λ is the cosmological constant, R is the 3D Ricci scalar, and M αβµν = c 1 g αβ g µν + c 2 g αµ g βν + c 3 g αν g βµ + c 4 u α u β g µν .
Another important characteristic of this theory is that only in the sector η = 0, Horava gravity admits Asymptotically AdS solution [27]. Therefore, assuming stationary and circular symmetry, the theory will admit the BTZ analogue to the three-dimensional rotating Hořava black holes described by metric where where a and b are constants that can be regarded as measures of aether misalignment, with b as a measure of asymptotically misalignment, for b = 0, the aether does not align with the timelike Killing vector asymptotically. Note that when ξ = 1 and λ = 1, the solution becomes BTZ black holes, and for ξ = 1, the solution becomes BTZ black holes with a shifted cosmological constant,Λ = Λ − 2b 2 (λ − 1). Also,J 2 can be negative, when either ξ < 0 or ξ > 1, a 2 > J 2 /(4(ξ − 1)). The sign ofΛ determine the asymptotic behavior (flat, dS, or AdS) of the metric [27]. In [27] was argued that if ξ > 0, and λ > 0 the aether represents a well defined folation at large r for any value of b. Moreover, if λ ≥ (1 + ξ)/2, thenΛ is always negative for any b. Also, if the coupling constants are such that ξ > 0, λ > 1/2 and λ < (1 + ξ)/2, thenΛ will switch sign at some value of b.
In Fig. 1 we show the behavior of F (r) 2 as a function of r. When ξ increases (left panel), we observe a region where there is no a horizon until ξ = ξ e for which the black hole becomes extremal, this value can be obtained from r − = r + . Then there is a region where r + increases and r − decreases until becomes null when ξ = ξ c , and finally there is a region ξ > ξ c where there is only one horizon. Also, we can observe the same behavior when λ decreases (right panel). Also, in Figs. 2 and 3, we plot the behavior ofΛ as a function of b, and its sign determines the asymptotic behavior (flat, dS, or AdS) of the metric. Note that for −1.4 < b < 1.4, the sign ofΛ is negative, as mentioned, if the coupling constants are such that ξ > 0, λ > 1/2 and λ < (1 + ξ)/2, thenΛ will switch sign at some value of b. The   value of ξ for which the black hole is extremal is given by the value of ξ for which the black hole passes from having two horizons to having one horizon, an is given by and the value of ξ for which the effective cosmological constantΛ changes sign is given by ξ = (2λ−1)b 2 b 2 +λ . In Fig curvature singularity due to the Ricci scalar is divergent at r = 0. This is in contrast to BTZ black holes where the Ricci and Kretschmann scalars are finite and smooth at r = 0. The locations of the inner and outer horizons r = r ± , are given by Also, M andJ can be written as M = −Λ(r 2 + + r 2 − ) andJ = 2r + r − √ −Λ, respectively. The Hawking temperature T H is given by

III. QNMS
The quasi-normal modes of scalar perturbations for a minimally coupled massive scalar field to curvature on the background of three-dimensional Hořava AdS black holes are described by the solution of the Klein-Gordon equation where m is the mass of the scalar field ψ. Which can be written as The term F (r) 2 − J 2 4r 2 is given by It is worth mentioning that the second term in the above expression vanishes for ξ = 1. Performing the change of variables z = A. Massive scalar field with circular symmetry For a massive scalar field with circular symmetry (κ = 0) Eq. (25) is which can be written as where Under the decomposition R(z) = z α (1 − z) β K(z), Eq. (27) can be written as a hypergeometric equation for K where the coefficients a 1 , b 1 , and c 1 are given by and the exponents α and β are The general solution of Eq. (29) takes the form which has three regular singular points at z = 0, z = 1, and z = ∞. Here, 2 F 1 (a 1 , b 1 , c 1 ; z) is a hypergeometric function and C 1 , C 2 are constants. Then, without loss of generality, we choose the negative sign for α, and the solution for the radial function R(z) is According to our change of variables at the vicinity of the horizon r → r + , z → 0, and at infinity r → ∞, z → 1. In the vicinity of the horizon, z = 0 and using the property F (a 1 , b 1 , c 1 , 0) = 1, the function R(z) behaves as and the scalar field ψ can be written in the following way in which the first term represents an ingoing wave and the second term an outgoing wave in the black hole. To compute the QNMs, we have to impose the boundary conditions on the horizon that there exist only ingoing waves. This fixes C 2 = 0. So, the radial solution becomes In order to implement boundary conditions at infinity (z = 1), we shall apply in Eq. (36), the Kummer's formula, for the hypergeometric function [36], With this expression, the radial function results in Therefore, by imposing that the scalar field at infinity is null, for m 2 /Λ < 0 (β − < 0 and c 1 − a 1 − b 1 + β − = β + > 0), then the term proportional to (1 − z) β in Eq. (38) diverges. So, we obtain that the scalar field is null only upon the following additional restriction ( Then, the QNFs yields respectively. Note that the imaginary part of the QNFs is negative, which ensures that the propagation of scalar fields is stable in this background. Now, in order to observe the behavior of the QNFs (39) and (40), we plot in Fig. 5, the behavior of the real (left panel) and imaginary parts (right panel) of the fundamental QNFs as a function of ξ. Note that, as mentioned, for ξ < ξ e , there is no horizon. So, for ξ > ξ e , we observe that for ω 1 (continuous line) there is a range where Re(ω 1 ) is null and then takes positive values, when the coupling constant ξ increases, while |Im(ω 1 )| decreases when ξ increases. So, according to the gauge/gravity duality, the relaxation time in order to reach the thermal equilibrium increases for the right sector. However, for ω 2 (dashed line) and ξ > ξ e , Re(ω 2 ) is null and then takes negative values, while its imaginary part increases and then decreases when the coupling constant ξ increases, showing that the relaxation time can decrease or increase depending on the value of ξ. It is interesting to note that when Im(ω 2 ) decreases Im(ω 1 ) = Im(ω 2 ). If we consider the BTZ black hole, ξ = 1 and λ = 1, the real part is null and Im(ω 1 ) = Im(ω 2 ). In the following, we will analyze the two branches of QNFs for different values of ξ. Fig. (6) is similar to Fig. (5), but in order to see the effect of b on the behavior of the QNFs, we have plotted several curves corresponding to different values of parameter b. For ξ = 1 the QNFs coincide and correspond to the QNFs of the BTZ black hole. On the other hand, for ξ = ξ c , both purely imaginary branches converge to the same value. For values of ξ near 1, |Im(ω 1 )| decreases when ξ increases, decreasing faster for large values of b, while |Im(ω 2 )| increases when ξ increases, increasing faster for small values of b, which implies that the relaxation time of the right sector increases and the relaxation time of the left sector decreases. On the other hand, we observe that for ξ near ξ c with ξ > ξ c , only one branch exists, and |Im(ω)| decreases for small values of b while it increases for large values of b. Notice that for b = 3 the effective cosmological constantΛ becomes positive before reaching the value xi c .

ξe < ξ < ξc
Now, in order to observe the behavior of the QNFs (39) and (40), in the range ξ e < ξ < ξ c , that is, where r ± are positive, first we plot r ± versus J in Fig. 7, for different values of the parameters ξ and λ, in order to see for which values of the parameter J the horizons r ± are positive. Then, for this range, we plot the imaginary part Im(ω) of the QNFs in Fig. 8, and we observe that for ω 1 , |Im(ω 1 )| increases when the parameter J (or equivalentlyJ, see Fig. 9) increases, see left panel of Fig. 8, so the relaxation time decreases. However, for ω 2 , the behavior is the opposite. |Im(ω 2 )| decreases when the parameter J increases, see right panel of Fig. 8, so the relaxation time increases. Note that in this range Re(ω) is null. Also, the sectors T R and T L of the conformal field theory are well defined. Furthermore, |Im(ω 1 )| decreases and |Im(ω 2 )| increases, when the coupling constant ξ increases; so, the relaxation time increases and decreases, respectively. In Fig. 10 we plot the fundamental QNFs as a function of b. We observe that for the BTZ black hole, the central panel with ξ = 1, the imaginary part of the fundamental QNF is constant; however, for a asymptotical misalignment of the aether with the timelike Killing vector, b = 0, the fundamental QNFs depend on b. |Im(ω 1 )| increases when ξ decreases. For |Im(ω 2 )| the behavior is the opposite since it decreases when ξ decreases.

ξ > ξc
As mentioned, for ξ > ξ c there is only one horizon. r + > 0 and r − become imaginary. This occurs for ξ > (J 2 + 4a 2 )/(4a 2 ), and consequently there is a gap in J, see Fig. 7, left panel, and for whichJ 2 is negative, see Fig. 9, which occurs when J 2 < 4a 2 (ξ − 1). In Fig. 11, we plot the fundamental QNFs for the range of values of J in which it is positive, and there is only one horizon. So, we observe that the fundamental QNFs acquire a real part, with Re(ω 1 ) = −Re(ω 2 ), |Re(ω)| decreases when J increases, and Im(ω 1 ) = Im(ω 2 ), and it is negative. In this case, the two branches have converge to one branch and when the coupling constant ξ increases |Im(ω)| decreases, see Fig. 5.

ξ = ξc
Finally, for ξ = ξ c , that is r − = 0, the two sectors converge, this occurs forJ = 0, that is J 2 + 4a 2 (1 − ξ) = 0. In this case the QNFs are given by where r 2 + = − M Λ . In Fig. 12 we show the behavior of the fundamental QNFs as a function of ξ = ξ c . We observe that |Im(ω)| decreases when ξ c increases whereas Re(ω) is null. So, according to the gauge/gravity duality, the relaxation time in order to reach the thermal equilibrium increases. Finally, we plot in Fig. 13 Im(ω) for different values of the constants a and b. We observe that for the range ξ > ξ e , |Im(ω)| increases when the constants a or b increase, so the relaxation time decreases.
Neumann boundary conditions. The frequencies found above for the scalar perturbation have been obtained by imposing the vanishing Dirichlet boundary condition at infinity. It is known that the Dirichlet boundary condition does not lead to any quasinormal modes for m 2 < 0. However, it is also possible to find a second set of QNFs, for negative mass squared, if we consider that the flux of the scalar field vanishes at infinity or vanishing Neumann boundary condition at infinity, which allows us to describe tachyons. Furthermore, it was shown that for negative mass squared there are two sets of dual operators ∆ + and ∆ − , where the second set of QNFs matches exactly the dual operators with ∆ − [24]. So, by using the condition that the flux, which is given by vanishes at asymptotic infinity, we obtain for β = β − and 0 > m 2 >Λ, that the flux vanishes if (a)| β− = −n or (b)| β+ = −n, which leads to respectively.

B. Massive scalar field
For massive scalar field, the Klein-Gordon equation can be rewritten as where Under the decomposition R(z) = z α (1 − z) β (ā − z)K(z), with Eq.(45) can be written as where that corresponds to Heun's differential equation. The condition 2 = 1 + 2 −(1+2α)−2β +1, ensures regularity of the point at ∞, and q corresponds to the accessory parameter. Heun's equation has four regular singular points: 0, 1, a, and ∞ with exponents (0, −2α), (0, 1 − 2β), (0, −1), and ( 1 , 2 ). Now, in order to obtain the QNFs, we proceed to perform a numerical analysis by using the pseudospectral Chebyshev method [37], which has been applied to compute the QNFs in other geometries, for instance see [38][39][40]. In Fig. 14 we plot the numerical results obtained for the real and imaginary parts of the fundamental QNF of the branch ω 2 as a function of ξ for different values of κ. We observe that for ξ > 1, the absolute value of the imaginary part decreases as κ increases; however, for ξ < 1, the behavior is the opposite; the absolute value of the imaginary part increases as κ increases. Note that, for the cases analyzed, the QNFs have a negative imaginary part, which ensures that the propagation of scalar fields is stable in this background.
where A, B, and C are given by Eq. (46). Under the decomposition R(z) = z α (1 − z) β K(z), Eq. (52) can be written as a hypergeometric equation for K, Eq. (29), where the coefficients a 1 , b 1 and c 1 are given by and the exponents α and β are Following the same procedure detailed in the case of massive radial scalar field we obtain for m 2 /Λ < 0 (β − < 0 and β + > 0) that the field at infinity is null if the gamma function Γ(x) has the poles at x = −n for n = 0, 1, 2, .... Then, the wave function satisfies the considered boundary condition only upon the following additional restriction (c 1 − a 1 )| β− = −n or (c 1 − b 1 )| β− = −n. These conditions determine the form of the QNFs as respectively. Now, in a similar manner to the first case, that is, using the condition that the flux vanishes at asymptotic infinity, we obtain for β = β + and 0 > m 2 >Λ, with (b 1 )| β+ = −n or (a 1 )| β− = −n, the second set of QNFs, which yields respectively. In this case, the QNFs correspond to ω = ω mr ∓ κ √ −Λ, where ω mr are the QNFs for massive radial scalar fields. So, the condition ((−1 +ā) 2 J 2 + 4āΛ r 2 − − r 2 + 2 = 0 or J =J) only has effect on the real part of the QNFs. As mentioned, when ξ = 1 (J=J), the solution results in the BTZ black holes with a shifted cosmological constant,Λ = Λ − 2b 2 (λ − 1). Note that the QNFs have real and imaginary parts, with an imaginary part that is negative, which ensures that the propagation of scalar fields is stable in this background.

IV. REMARKS AND CONCLUSIONS
In this work, we computed the QNMs of rotating three-dimensional Hořava AdS black holes and we analyzed the effect of the breaking of the Lorentz invariance on the QNMs. We showed that depending on the parameters, the lapsus function can represent a spacetime without an event horizon, a black hole geometry with one event horizon, an extremal black hole and finally a black hole with two horizons. The QNMs have been obtained by imposing on the horizon that there are only ingoing waves, while at infinity Dirichlet boundary conditions and Neumann boundary conditions were imposed. We found that the propagation of the scalar field is stable in this background, since the imaginary part of the QNFs is negative. Also, we made a systematic study of the QNMs and QNFs for various ranges of parameters J, the angular momentum and ξ the parameter that differentiates Horava gravity from General Relativity which is obtained for ξ = 1. For various values of J and ξ, we obtained various branches of solutions with different properties of QNMs and QNFs. In particular: For positive inner and outer horizons r ± , the range ξ e < ξ < ξ c gives Re(ω) null for the fundamental QNFs, and the two sectors T R and T L are well defined. |Im(ω 1 )| increases when the parameter J increases, so according to the gauge/gravity duality, the relaxation time in order to reach the thermal equilibrium decreases. Furthermore, as the coupling constant ξ increases, |Im(ω 1 )| decreases and therefore the relaxation time increases. However, for ω 2 , the behavior is opposite. |Im(ω 2 )| decreases when the parameter J increases and therefore the relaxation time increases. Furthermore, |Im(ω 2 )| increases when ξ increases; so, the relaxation time decreases.
Moreover, for the general case, that is, a massive scalar field, it was shown that the Klein-Gordon equation can be written as the Heun's differential equation, and we have studied the behavior of the QNMs numerically via the pseudospectral Chebyshev method, and mainly it was found that for ω 2 and ξ > 1, the absolute value of the imaginary part decreases as κ increases; however, for ξ < 1 the behavior is the opposite. The absolute value of the imaginary part increases as κ increases.
As can be seen from the above discussion, the oscillatory and the decay modes of the QNMs have quite different behavior for the various branches, indicating that the time required for a system to reach thermal equilibrium on the boundary is different for the various values of the parameters. An intriguing result is that if the parameter ξ lies between two different critical values, then the time required for the two sectors T R and T L to reach thermal equilibrium is competing in the sense that in one sector the time is increasing while in the other sector it is decreasing. This behavior deserves further study in connection of trying to find a system on the boundary that exhibits such a behavior.
It would be interesting to extent this work to higher dimensional Horava black holes [41] and calculating the QNMs and QNFs of a massive wave to study how a gravity theory in the bulk with broken Lorentz invariance affects the boundary field theory to reach thermal equilibrium.