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

TL;DR

This paper calculates the quasinormal modes of a scalar field in a rotating three-dimensional Hořava AdS black hole, revealing stability and insights into the system's thermalization process under Lorentz invariance breaking.

Contribution

It provides the first analysis of quasinormal modes in a rotating Hořava AdS black hole, exploring the effects of Lorentz invariance violation on stability and thermalization.

Findings

Scalar field propagation is stable in this background.

Different solutions reach thermal equilibrium at different times.

Lorentz invariance breaking influences QNM characteristics.

Abstract

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.

Figures

Equations

$$d{s}^2=-{\sinh }^2\mu {\left(r_{+}dt-r_{-}d\varphi \right)}^2+d{\mu }^2+{\cosh }^2\mu {\left(-r_{-}dt+r_{+}d\varphi \right)}^2.$$

$$T_L=(r_{+}-r_{-})/2\pi ,T_R=(r_{+}+r_{-})/2\pi .$$

$$h_R+h_L=\Delta ,h_R-h_L=\pm s,$$

$$\Delta =1+\sqrt{1+m^2}.$$

$${\omega }_L$$

$${\omega }_R$$

$$S_H=\frac{1}{16\pi G_H}\int dT{d}^{2}xN\sqrt{g}\left[L_2+L_4\right],$$

$$L_2=K_{ij}{K}^{ij}-\lambda K^2+\xi \left({}^{(2)}R-2\Lambda \right)+\eta a_i{a}^i,$$

$$d{s}^2=N^{2}d{T}^2-g_{ij}(d{x}^i+N^{i}dT)(d{x}^j+N^{j}dT).$$

$$u_{\alpha }=\frac{{\partial }_{\alpha }T}{\sqrt{g^{\mu \nu }{\partial }_{\mu }T\partial \nu T}},$$

$$\frac{G_H}{G_{ae}}=\xi =\frac{1}{1-c_{13}},\frac{\lambda }{\xi }=1+c_2,\frac{\eta }{\xi }=c_{14},$$

$$S_{ae}=\frac{1}{16\pi G_{ae}}\int d^{3}x\sqrt{-g}(-R-2\Lambda +L_{ae}),$$

$$L_{ae}=-M^{\alpha \beta \mu \nu }{\nabla }_{\alpha }{u}_{\mu }{\nabla }_{\beta }{u}_{\nu },$$

$$M^{\alpha \beta \mu \nu }=c_1{g}^{\alpha \beta }{g}^{\mu \nu }+c_2{g}^{\alpha \mu }{g}^{\beta \nu }+c_3{g}^{\alpha \nu }{g}^{\beta \mu }+c_4{u}^{\alpha }{u}^{\beta }{g}^{\mu \nu }.$$

$$d{s}^2=Z(r{)}^{2}d{t}^2-\frac{1}{F(r{)}^2}d{r}^2-r^2(d\varphi +\Omega (r)dt{)}^2,$$

$$F(r{)}^2=Z(r{)}^2=-M+\frac{{\bar{J}}^2}{4r^2}-\bar{\Lambda }{r}^2,$$

$${\bar{J}}^2=\frac{J^2+4a^2(1-\xi )}{\xi },\Omega (r)=-\frac{J}{2r^2},\bar{\Lambda }=\Lambda -\frac{b^2(2\lambda -\xi -1)}{\xi },$$

$${\xi }_e$$

$${\xi }_c=\frac{4a^2+J^2}{4a^2},$$

$$R=-6\bar{\Lambda }+\frac{1}{2r^2}\left({\bar{J}}^2-J^2\right)$$

$$r_{\pm }^2=-\frac{M}{2\bar{\Lambda }}\left(1\pm \sqrt{1+\frac{{\bar{J}}^2\bar{\Lambda }}{M^2}}\right).$$

$$T_H=\frac{-\bar{\Lambda }(r_{+}^2-r_{-}^2)}{2\pi r_{+}}.$$

$$\square \psi =\frac{1}{\sqrt{-g}}{\partial }_{\mu }\left(\sqrt{-g}{g}^{\mu \nu }{\partial }_{\nu }\right)\psi =-m^2\psi ,$$

$$\left(-\frac{1}{F(r{)}^2}{\partial }_t^2+F(r{)}^2{\partial }_r^2+\frac{1}{r}{\partial }_r(rF(r{)}^2){\partial }_r-\frac{J}{r^{2}F(r{)}^2}{\partial }_t{\partial }_{\varphi }+\frac{1}{r^{2}F(r{)}^2}(F(r{)}^2-\frac{J^2}{4r^2}){\partial }_{\varphi }^2-m^2\right)\Psi =0.$$

$$-M+\frac{{\bar{J}}^2}{4r^2}-\bar{\Lambda }{r}^2-\frac{J^2}{4r^2}$$

$$z(1-z){\partial }_z^{2}R(z)+\left(1-z\right){\partial }_{z}R(z)+$$

$$z(1-z){\partial }_z^{2}R(z)+\left(1-z\right){\partial }_{z}R(z)+\left[-\frac{{\omega }^2(z{r}_{-}^2-r_{+}^2)}{4{\bar{\Lambda }}^2(r_{+}^2-r_{-}^2{)}^{2}z}+\frac{m^2}{4(1-z)\bar{\Lambda }}\right]R(z)=0,$$

$$z(1-z){\partial }_z^{2}R(z)+\left(1-z\right){\partial }_{z}R(z)+\left(A+\frac{B}{z}+\frac{C}{1-z}\right)R(z)=0,$$

$$A=-\frac{{\omega }^2{r}_{-}^2}{4{\bar{\Lambda }}^2(r_{+}^2-r_{-}^2{)}^2},B=\frac{{\omega }^2{r}_{+}^2}{4{\bar{\Lambda }}^2(r_{+}^2-r_{-}^2{)}^2},C=\frac{m^2}{4\bar{\Lambda }}.$$

$$z(1-z)K^{\prime \prime }(z)+\left[c_1-(1+a_1+b_1)z\right]K^{\prime }(z)-a_1{b}_{1}K(z)=0,$$

$$a_1=\alpha +\beta \mp \sqrt{A},b_1=\alpha +\beta \pm \sqrt{A},c_1=1+2\alpha ,$$

Full text

I Introduction

If a dynamical system is perturbed, it will return to equilibrium, and this process is completely determined by the poles of the retarded correlation function of the perturbation. In gravity theories, black holes are thermodynamical systems and perturbations of them at equilibrium are described by the quasi-normal modes (QNMs) Regge:1957td ; Zerilli:1971wd ; Zerilli:1970se ; Kokkotas:1999bd ; Nollert:1999ji ; Berti:2009kk ; Konoplya:2011qq . The QNMs are determined by solving the wave equation of an incident wave with the right boundary conditions. Then the solution of the wave equation determines the complex frequencies, the real part of which gives the rate of oscillations of the wave while their complex part gives the required decay time for the system to reach thermal equilibrium.

The QNMs and quasi-normal frequencies (QNFs) have been the subjects of study for a long time and have recently acquired great interest due to the detection of gravitational waves Abbott:2016blz . Despite the detected signal being consistent with the Einstein gravity TheLIGOScientific:2016src , there are great uncertainties in mass and angular momenta of the ringing black hole, which leaves open possibilities for alternative theories of gravity Konoplya:2016pmh like $f(R)$ gravity Starobinsky:1980te ; DeFelice:2010aj ; Nojiri:2005jg and Galileon gravity theories Nicolis:2008in ; Deffayet:2009wt ; Kolyvaris:2011fk . Also, the QNMs and the QNFs were extensively studied in connection with the stability of black holes in Einstein gravity Wang:2000gsa ; Wang:2004bv and in modified gravity theories Konoplya:2018qov ; Abdalla:2018ggo ; Abdalla:2019irr .

The gauge/gravity duality which results from the AdS/CFT correspondence Maldacena:1997re ; Aharony:1999ti stimulated the interest in calculating the QNMs and QNFs of black holes in AdS spacetime. It was shown in Birmingham:2001pj that this holographic principle leads to the existence of a correspondence between the QNMs in AdS black holes and linear response theory in scale invariant finite temperature field theory. This correspondence of the decay of perturbations in the dual conformal field theory and the QNMs in the gravity bulk was first discussed in Horowitz:1999jd . Considering the $(2+1)$-dimensional AdS black hole BTZ , it was shown analytically Birmingham:2001pj that there is an agreement between its QNFs and the location of the 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 Sotiriou:2014gna , 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 Siopsis:2008xz ). Exact solutions for QNMs of black holes in three spacetime dimensions have been obtained in exactQNM .

The motivation for considering the Hořava gravity theory is twofold. Considering a condensed matter dynamical system, it was argued in Nicolis:2015sra 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 Papantonopoulos:2011zz .

In this context, the QNMs for four-dimensional non-reduced Einstein-aether theory was studied, and it was found in Konoplya:2006rv 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 Ding:2017gfw 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

[TABLE]

The angular coordinate $\phi$ has period $2\pi$, 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 Birmingham:2001pj . Consequently, the conformal field theory splits into two independent sectors at thermal equilibrium with temperatures

[TABLE]

According to the AdS3/CFT2 correspondence, to each field of spin $s$ propagating in AdS3 there corresponds an operator $\cal O$ in the dual conformal field theory characterised by conformal weights $(h_{L},h_{R})$ with Aharony:1999ti

[TABLE]

and $\Delta$ is determined in terms of the mass $m$ of the scalar field,

[TABLE]

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 ${\cal 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 Birmingham:2001pj

[TABLE]

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 Birmingham:2001pj .

We now discuss the three-dimensional Hořava gravity, the action of which is given by Sotiriou:2011dr

[TABLE]

where $G_{H}$ is a coupling constant with dimensions of a length squared and the Lagrangian $L_{2}$ has the following form

[TABLE]

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}=\partial_{i}\ln{N}$, being the line element in the preferred foliation

[TABLE]

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 $\lambda=\xi=1$ and $\eta=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

[TABLE]

and

[TABLE]

where $c_{ij}=c_{i}+c_{j}$. Being the action of Einstein-aether:

[TABLE]

where $G_{ae}$ is a coupling constant with dimensions of a length square, $g$ is the determinant of $g_{\mu\nu}$, $\Lambda$ is the cosmological constant, $R$ is the $3D$ Ricci scalar,

[TABLE]

and

[TABLE]

Another important characteristic of this theory is that only in the sector $\eta=0$, Horava gravity admits Asymptotically AdS solution Sotiriou:2014gna . 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

[TABLE]

where

[TABLE]

with

[TABLE]

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\neq 0$, the aether does not align with the timelike Killing vector asymptotically. Note that when $\xi=1$ and $\lambda=1$, the solution becomes BTZ black holes, and for $\xi=1$, the solution becomes BTZ black holes with a shifted cosmological constant, $\bar{\Lambda}=\Lambda-2b^{2}(\lambda-1)$. Also, $\bar{J}^{2}$ can be negative, when either $\xi<0$ or $\xi>1$, $a^{2}>J^{2}/(4(\xi-1))$. The sign of $\bar{\Lambda}$ determine the asymptotic behavior (flat, dS, or AdS) of the metric Sotiriou:2014gna .

In Sotiriou:2014gna was argued that if $\xi>0$, and $\lambda>0$ the aether represents a well defined folation at large $r$ for any value of $b$. Moreover, if $\lambda\geq(1+\xi)/2$, then $\bar{\Lambda}$ is always negative for any $b$. Also, if the coupling constants are such that $\xi>0$, $\lambda>1/2$ and $\lambda<(1+\xi)/2$, then $\bar{\Lambda}$ÃÑ 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 $\xi$ increases (left panel), we observe a region where there is no a horizon until ${\bf{\xi=\xi_{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 $\xi=\xi_{c}$, and finally there is a region $\xi>\xi_{c}$ where there is only one horizon. Also, we can observe the same behavior when $\lambda$ decreases (right panel). Also, in Figs. 2 and 3, we plot the behavior of $\bar{\Lambda}$ 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 $\bar{\Lambda}$ is negative, as mentioned, if the coupling constants are such that $\xi>0$, $\lambda>1/2$ and $\lambda<(1+\xi)/2$, then $\bar{\Lambda}$ will switch sign at some value of b.

The value of $\xi$ for which the black hole is extremal is given by

[TABLE]

the value of $\xi$ for which the black hole passes from having two horizons to having one horizon, an is given by

[TABLE]

and the value of $\xi$ for which the effective cosmological constant $\bar{\Lambda}$ changes sign is given by $\xi=\frac{(2\lambda-1)b^{2}}{b^{2}+\lambda}$. In Fig. (4) we plot the different regions defined by $\xi_{e}$, $\xi_{c}$, and $\bar{\Lambda}$ for a choice of parameters.

In the case $\bar{J}\neq J$ ($\xi\neq 1$), there is a curvature singularity due to the Ricci scalar

[TABLE]

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_{\pm}$, are given by

[TABLE]

Also, $M$ and $\bar{J}$ can be written as $M=-\bar{\Lambda}(r_{+}^{2}+r_{-}^{2})$ and $\bar{J}=2r_{+}r_{-}\sqrt{-\bar{\Lambda}}$, respectively. The Hawking temperature $T_{H}$ is given by

[TABLE]

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

[TABLE]

where $m$ is the mass of the scalar field $\psi$. Which can be written as

[TABLE]

The term $F(r)^{2}-\frac{J^{2}}{4r^{2}}$ is given by

[TABLE]

It is worth mentioning that the second term in the above expression vanishes for $\xi=1$. Performing the change of variables $z=\frac{r^{2}-r_{+}^{2}}{r^{2}-r_{-}^{2}}$ along with the ansatz $\psi=R(z)e^{-i\omega t}e^{i\kappa\phi}$, Eq. (23) yields

[TABLE]

III.1 Massive scalar field with circular symmetry

For a massive scalar field with circular symmetry ($\kappa=0$) Eq. (III) is

[TABLE]

which can be written as

[TABLE]

where

[TABLE]

Under the decomposition $R(z)=z^{\alpha}(1-z)^{\beta}K(z)$, Eq. (27) can be written as a hypergeometric equation for $K$

[TABLE]

where the coefficients $a_{1}$, $b_{1}$, and $c_{1}$ are given by

[TABLE]

and the exponents $\alpha$ and $\beta$ are

[TABLE]

The general solution of Eq. (29) takes the form

[TABLE]

which has three regular singular points at $z=0$, $z=1$, and $z=\infty$. 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 $\alpha$, and the solution for the radial function $R(z)$ is

[TABLE]

According to our change of variables at the vicinity of the horizon $r\rightarrow r_{+}$, $z\rightarrow 0$, and at infinity $r\rightarrow\infty$, $z\rightarrow 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

[TABLE]

and the scalar field $\psi$ can be written in the following way

[TABLE]

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

[TABLE]

In order to implement boundary conditions at infinity ($z=1$), we shall apply in Eq. (36), the Kummer’s formula, for the hypergeometric function M. Abramowitz ,

[TABLE]

With this expression, the radial function results in

[TABLE]

Therefore, by imposing that the scalar field at infinity is null, for $m^{2}/\bar{\Lambda}<0$ ($\beta_{-}<0$ and $c_{1}-a_{1}-b_{1}+\beta_{-}=\beta_{+}>0$), then the term proportional to $(1-z)^{\beta}$ in Eq. (III.1) diverges. So, we obtain that the scalar field is null only upon the following additional restriction $(c_{1}-a_{1})|_{\beta_{-}}=-n$ or $(c_{1}-b_{1})|_{\beta_{-}}=-n$. Then, the QNFs yields

[TABLE]

[TABLE]

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 $\xi$. Note that, as mentioned, for $\xi<\xi_{e}$, there is no horizon. So, for $\xi>\xi_{e}$, we observe that for $\omega_{1}$ (continuous line) there is a range where $Re(\omega_{1})$ is null and then takes positive values, when the coupling constant $\xi$ increases, while $\lvert Im(\omega_{1})\rvert$ decreases when $\xi$ 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 $\omega_{2}$ (dashed line) and $\xi>\xi_{e}$, $Re(\omega_{2})$ is null and then takes negative values, while its imaginary part increases and then decreases when the coupling constant $\xi$ increases, showing that the relaxation time can decrease or increase depending on the value of $\xi$. It is interesting to note that when $Im(\omega_{2})$ decreases $Im(\omega_{1})=Im(\omega_{2})$. If we consider the BTZ black hole, $\xi=1$ and $\lambda=1$, the real part is null and $Im(\omega_{1})\neq Im(\omega_{2})$. In the following, we will analyze the two branches of QNFs for different values of $\xi$. 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 $\xi=1$ the QNFs coincide and correspond to the QNFs of the BTZ black hole. On the other hand, for $\xi=\xi_{c}$, both purely imaginary branches converge to the same value. For values of $\xi$ near 1, $|Im(\omega_{1})|$ decreases when $\xi$ increases, decreasing faster for large values of $b$, while $|Im(\omega_{2})|$ increases when $\xi$ 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 $\xi$ near $\xi_{c}$ with $\xi>\xi_{c}$, only one branch exists, and $|Im(\omega)|$ decreases for small values of $b$ while it increases for large values of $b$. Notice that for $b=3$ the effective cosmological constant $\bar{\Lambda}$ becomes positive before reaching the value $xi_{c}$.

III.1.1 $\xi_{e}<\xi<\xi_{c}$

Now, in order to observe the behavior of the QNFs (39) and (40), in the range $\xi_{e}<\xi<\xi_{c}$, that is, where $r_{\pm}$ are positive, first we plot $r_{\pm}$ versus $J$ in Fig. 7, for different values of the parameters $\xi$ and $\lambda$, in order to see for which values of the parameter $J$ the horizons $r_{\pm}$ are positive. Then, for this range, we plot the imaginary part $Im(\omega)$ of the QNFs in Fig. 8, and we observe that for $\omega_{1}$, $\lvert Im(\omega_{1})\rvert$ increases when the parameter $J$ (or equivalently $\bar{J}$, see Fig. 9) increases, see left panel of Fig. 8, so the relaxation time decreases. However, for $\omega_{2}$, the behavior is the opposite. $\lvert Im(\omega_{2})\rvert$ decreases when the parameter $J$ increases, see right panel of Fig. 8, so the relaxation time increases. Note that in this range $Re(\omega)$ is null. Also, the sectors $T_{R}$ and $T_{L}$ of the conformal field theory are well defined. Furthermore, $\lvert Im(\omega_{1})\rvert$ decreases and $\lvert Im(\omega_{2})\rvert$ increases, when the coupling constant $\xi$ 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 $\xi=1$, the imaginary part of the fundamental QNF is constant; however, for a asymptotical misalignment of the aether with the timelike Killing vector, $b\neq 0$, the fundamental QNFs depend on $b$. $\lvert Im(\omega_{1})\rvert$ increases when $\xi$ decreases. For $\lvert Im(\omega_{2})\rvert$ the behavior is the opposite since it decreases when $\xi$ decreases.

III.1.2 $\xi>\xi_{c}$

As mentioned, for $\xi>\xi_{c}$ there is only one horizon. $r_{+}>0$ and $r_{-}$ become imaginary. This occurs for $\xi>(J^{2}+4a^{2})/(4a^{2})$, and consequently there is a gap in $J$, see Fig. 7, left panel, and for which $\bar{J}^{2}$ is negative, see Fig. 9, which occurs when $J^{2}<4a^{2}(\xi-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(\omega_{1})=-Re(\omega_{2})$, $\lvert Re(\omega)\rvert$ decreases when $J$ increases, and $Im(\omega_{1})=Im(\omega_{2})$, and it is negative. In this case, the two branches have converge to one branch and when the coupling constant $\xi$ increases $\lvert Im(\omega)\rvert$ decreases, see Fig. 5.

III.1.3 $\xi=\xi_{c}$

Finally, for $\xi=\xi_{c}$, that is $r_{-}=0$, the two sectors converge, this occurs for $\bar{J}=0$, that is $J^{2}+4a^{2}(1-\xi)=0$. In this case the QNFs are given by

[TABLE]

where $r_{+}^{2}=-\frac{M}{\bar{\Lambda}}$. In Fig. 12 we show the behavior of the fundamental QNFs as a function of $\xi=\xi_{c}$. We observe that $\lvert Im(\omega)\rvert$ decreases when $\xi_{c}$ increases whereas $Re(\omega)$ 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(\omega)$ for different values of the constants $a$ and $b$. We observe that for the range $\xi>\xi_{e}$, $\lvert Im(\omega)\rvert$ 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 $\Delta_{+}$ and $\Delta_{-}$, where the second set of QNFs matches exactly the dual operators with $\Delta_{-}$ Birmingham:2001pj . So, by using the condition that the flux, which is given by

[TABLE]

vanishes at asymptotic infinity, we obtain for $\beta=\beta_{-}$ and $0>m^{2}>\bar{\Lambda}$, that the flux vanishes if $(a)|_{\beta_{-}}=-n$ or $(b)|_{\beta_{+}}=-n$, which leads to

[TABLE]

[TABLE]

respectively.

III.2 Massive scalar field

For massive scalar field, the Klein-Gordon equation can be rewritten as

[TABLE]

where

[TABLE]

and

[TABLE]

Under the decomposition $R(z)=z^{\alpha}(1-z)^{\beta}(\bar{a}-z)K(z)$, with

[TABLE]

Eq.(45) can be written as

[TABLE]

where

[TABLE]

[TABLE]

that corresponds to Heun’s differential equation. The condition $2=\epsilon_{1}+\epsilon_{2}-(1+2\alpha)-2\beta+1$, ensures regularity of the point at $\infty$, and $q$ corresponds to the accessory parameter. Heun’s equation has four regular singular points: $0,1,a$, and $\infty$ with exponents $(0,-2\alpha)$, $(0,1-2\beta)$, $(0,-1)$, and $(\epsilon_{1},\epsilon_{2})$. Now, in order to obtain the QNFs, we proceed to perform a numerical analysis by using the pseudospectral Chebyshev method Boyd , which has been applied to compute the QNFs in other geometries, for instance see Gonzalez:2017shu ; Gonzalez:2018xrq ; Finazzo:2016psx . In Fig. 14 we plot the numerical results obtained for the real and imaginary parts of the fundamental QNF of the branch $\omega_{2}$ as a function of $\xi$ for different values of $\kappa$. We observe that for $\xi>1$, the absolute value of the imaginary part decreases as $\kappa$ increases; however, for $\xi<1$, the behavior is the opposite; the absolute value of the imaginary part increases as $\kappa$ 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.

III.3 Case: $(-1+\bar{a})^{2}J^{2}+4\bar{a}\bar{\Lambda}\left(r_{-}^{2}-r_{+}^{2}\right)^{2}=0$ or $J=\bar{J}$

In this case, Eq. (45) becomes

[TABLE]

where $A$, $B$, and $C$ are given by Eq. (46). Under the decomposition $R(z)=z^{\alpha}(1-z)^{\beta}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

[TABLE]

and the exponents $\alpha$ and $\beta$ are

[TABLE]

Following the same procedure detailed in the case of massive radial scalar field we obtain for $m^{2}/\bar{\Lambda}<0$ ($\beta_{-}<0$ and $\beta_{+}>0$) that the field at infinity is null if the gamma function $\Gamma(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})|_{\beta_{-}}=-n$ or $(c_{1}-b_{1})|_{\beta_{-}}=-n$. These conditions determine the form of the QNFs as

[TABLE]

[TABLE]

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 $\beta=\beta_{+}$ and $0>m^{2}>\bar{\Lambda}$, with $(b_{1})|_{\beta_{+}}=-n$ or $(a_{1})|_{\beta_{-}}=-n$, the second set of QNFs, which yields

[TABLE]

[TABLE]

respectively. In this case, the QNFs correspond to $\omega=\omega_{mr}\mp\kappa\sqrt{-\bar{\Lambda}}$, where $\omega_{mr}$ are the QNFs for massive radial scalar fields. So, the condition ($(-1+\bar{a})^{2}J^{2}+4\bar{a}\bar{\Lambda}\left(r_{-}^{2}-r_{+}^{2}\right)^{2}=0$ or $J=\bar{J}$) only has effect on the real part of the QNFs. As mentioned, when $\xi=1$ ($J$=$\bar{J}$), the solution results in the BTZ black holes with a shifted cosmological constant, $\bar{\Lambda}=\Lambda-2b^{2}(\lambda-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 $\xi$ the parameter that differentiates Horava gravity from General Relativity which is obtained for $\xi=1$. For various values of $J$ and $\xi$, we obtained various branches of solutions with different properties of QNMs and QNFs. In particular:

For positive inner and outer horizons $r_{\pm}$, the range $\xi_{e}<\xi<\xi_{c}$ gives $Re(\omega)$ null for the fundamental QNFs, and the two sectors $T_{R}$ and $T_{L}$ are well defined. $\lvert Im(\omega_{1})\rvert$ 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 $\xi$ increases, $\lvert Im(\omega_{1})\rvert$ decreases and therefore the relaxation time increases. However, for $\omega_{2}$, the behavior is opposite. $\lvert Im(\omega_{2})\rvert$ decreases when the parameter $J$ increases and therefore the relaxation time increases. Furthermore, $\lvert Im(\omega_{2})\rvert$ increases when $\xi$ increases; so, the relaxation time decreases.

In the range $\xi>\xi_{c}$ the black hole only has one horizon $r_{+}>0$. For this case $J^{2}<4a^{2}(\xi-1)$ and the QNFs acquire a real part, with $Re(\omega_{1})=-Re(\omega_{2})$, $\lvert Re(\omega)\rvert$ decreases when $J$ increases, and $Im(\omega_{1})=Im(\omega_{2})$, and it is negative. Also, when $\xi$ increases, $\lvert Im(\omega)\rvert$ decreases. Finally, for $\xi=\xi_{c}$, that is, $r_{-}=0$, the two sectors converge. This occurs when $J^{2}+4a^{2}(1-\xi)=0$. In this case $\lvert Im(\omega)\rvert$ decreases when $\xi_{c}$ increases whereas $Re(\omega)$ is null. Therefore, the relaxation time increases. Also, we have considered different values of the constants $a$ and $b$, and $\lvert Im(\omega)\rvert$ increases when the constant $a$ or $b$ 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 $\omega_{2}$ and $\xi>1$, the absolute value of the imaginary part decreases as $\kappa$ increases; however, for $\xi<1$ the behavior is the opposite. The absolute value of the imaginary part increases as $\kappa$ 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 $\xi$ 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 Koutsoumbas:2010pt 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.

Acknowledgments

This work was funded by the Dirección de Investigación y Desarrollo de la Universidad de La Serena (Y.V.). E.P., Y.V. and R.B would like to thank the Facultad de Ingeniería y Ciencias, Universidad Diego Portales for its hospitality, where part of this work was carried out.