Scalar field quasinormal modes on asymptotically locally flat rotating black holes in three dimensions

TL;DR

This paper analytically computes the quasinormal modes of a massless scalar field around rotating black holes in three-dimensional New Massive Gravity, revealing unique damping behaviors as angular momentum varies.

Contribution

It provides the first closed-form solutions for scalar quasinormal frequencies on these asymptotically locally flat rotating black holes.

Findings

Quasinormal frequencies are obtained analytically.

Imaginary parts tend to minus infinity as angular momentum approaches zero.

Static black holes in this family do not admit quasinormal modes.

Abstract

The pure quadratic term of New Massive Gravity in three dimensions admits asymptotically locally flat, rotating black holes. These black holes are characterized by their mass and angular momentum, as well as by a hair of gravitational origin. As in the Myers-Perry solution in dimensions greater than five, there is no upper bound on the angular momentum. We show that, remarkably, the equation for a massless scalar field on this background can be solved in an analytic manner and that the quasinormal frequencies can be found in a closed form. The spectrum is obtained requiring ingoing boundary conditions at the horizon and an asymptotic behavior at spatial infinity that provides a well-defined action principle for the scalar probe. As the angular momentum of the black hole approaches zero, the imaginary part of the quasinormal frequencies tends to minus infinity, migrating to the north pole of the Riemann Sphere and providing infinitely damped modes of high frequency. We show that this is consistent with the fact that the static black hole within this family does not admit quasinormal modes for a massless scalar probe.

Equations

$$I\left[g\right]=\frac{1}{{\kappa }^2}\int d^{3}x\sqrt{-g}\left[-R+\frac{1}{m^2}\left(R_{\mu \nu }{R}^{\mu \nu }-\frac{3}{8}{R}^2\right)\right].$$

$$K_{\mu \nu }=2\square R_{\mu \nu }-\frac{1}{2}\left({\nabla }_{\mu }{\nabla }_{\nu }R+g_{\mu \nu }\square R\right)-8R_{\mu \rho }{R}_{\nu }^{\rho }+\frac{9}{2}R{R}_{\mu \nu }+g_{\mu \nu }\left(3R^{\alpha \sigma }{R}_{\alpha \sigma }-\frac{13}{8}{R}^2\right)=0.$$

$$d{s}^2=-\left(br-\mu \right)d{t}^2+\frac{d{r}^2}{br-\mu }-a\left(br-\mu \right)dtd\varphi +\left(r^2+r_0^2\right)d{\varphi }^2,$$

$$r_0^2=\frac{a^4{b}^2+16a^2\mu }{64}.$$

$$M=\frac{b^2}{32G},J=Ma.$$

$$T=\frac{b}{4\pi },S=\frac{\pi b}{4G}.$$

$$d{s}^2$$

$$+\left(r^2+r_0^2\right){\left(d\psi -\frac{32a\left(br-\mu \right)}{64r^2+16a^2\mu +a^4{b}^2}du\right)}^2,$$

$$du$$

$$d\psi$$

$$\underset{r\to +\infty }{\lim }{R}_{\alpha \beta }^{\mu \nu }=0.$$

$$\square \Phi =0.$$

$$\Phi \left(t,\varphi ,r\right)=\Re \left(\sum _{n=-\infty }^{+\infty }\int d\omega e^{-i\omega t+in\varphi }{R}_{\omega ,n}\left(r\right)\right).$$

$$r=\frac{x+\mu }{b},$$

$$A\left(x\right)\frac{d^{2}R\left(x\right)}{d{x}^2}+B\left(x\right)\frac{dR\left(x\right)}{dx}+C\left(x\right)R\left(x\right)=0,$$

$$A\left(x\right)$$

$$B\left(x\right)$$

$$C\left(x\right)$$

$$R\left(x\right)=C_1{x}^{-i\frac{\omega }{b}}\left(1+\mathcal{O}\left(x\right)\right)+C_2{x}^{i\frac{\omega }{b}}\left(1+\mathcal{O}\left(x\right)\right),$$

$$R\left(x\right)=\frac{{\tilde{C}}_1}{x^{{\delta }_{+}}}\left(1+\mathcal{O}\left(\frac{1}{x}\right)\right)+\frac{{\tilde{C}}_2}{x^{{\delta }_{-}}}\left(1+\mathcal{O}\left(\frac{1}{x}\right)\right),$$

$${\delta }_{\pm }=\frac{1}{2}\left(1\pm \sqrt{1-\frac{4{\omega }^2}{b^2}}\right).$$

$$R\left(x\right)$$

$$+C_2{\left(-\frac{8x}{a^2{b}^2+8\mu }\right)}^{1-\gamma }F\left(\alpha +1-\gamma ,\beta +1-\gamma ,2-\gamma ,-\frac{8x}{a^2{b}^2+8\mu }\right)],$$

$$\alpha$$

$$\beta$$

$$\gamma$$

$$R\left(x\right)=C_1{\left(8x+8\mu +a^2{b}^2\right)}^{-i\frac{\sqrt{2}\left(a\omega +2n\right)}{\sqrt{a^2{b}^2+8\mu }}}{x}^{-i\frac{\omega }{b}}F\left(\alpha ,\beta ,\gamma ,-\frac{8x}{a^2{b}^2+8\mu }\right).$$

$$R\left(x\right)\sim \frac{\Gamma \left(\gamma \right)\Gamma \left(\beta -\alpha \right)}{\Gamma \left(\gamma -\alpha \right)\Gamma \left(\beta \right)}\frac{1}{x^{{\delta }_{+}}}\left(1+\mathcal{O}\left(\frac{1}{x}\right)\right)+\xi \frac{\Gamma \left(\gamma \right)\Gamma \left(\alpha -\beta \right)}{\Gamma \left(\alpha \right)\Gamma \left(\gamma -\beta \right)}\frac{1}{x^{{\delta }_{-}}}\left(1+\mathcal{O}\left(\frac{1}{x}\right)\right),$$

$$I\left[\varphi \right]=-\frac{1}{2}\int \sqrt{-g}{d}^{3}x{\nabla }_{\nu }\Phi {\nabla }^{\nu }\Phi ,$$

$$\delta I\left[\varphi \right]=-\int \sqrt{-g}{d}^{3}x{\partial }_{\nu }\Phi {\partial }^{\nu }\delta \Phi =-\int \sqrt{-g}{d}^{3}x\left({\nabla }_{\nu }\left({\nabla }^{\nu }\Phi \delta \Phi \right)-\square \Phi \delta \Phi \right),$$

$${\delta }_{\text{on-shell}}I\left[\varphi \right]=-\underset{r_0\to +\infty }{\lim }{\int d^{2}x\sqrt{-g}{g}^{rr}{\partial }_r\Phi \delta \Phi |}_{r=r_0}.$$

$${\left(\frac{\Gamma \left(\beta -\alpha \right)}{\Gamma \left(\gamma -\alpha \right)\Gamma \left(\beta \right)}\right)}^2\frac{{\delta }_{+}}{x^{2{\delta }_{+}-1}}+\xi \left(\frac{\Gamma \left(\alpha -\beta \right)\Gamma \left(\beta -\alpha \right)}{\Gamma \left(\alpha \right)\Gamma \left(\gamma -\beta \right)\Gamma \left(\gamma -\alpha \right)\Gamma \left(\beta \right)}\right)\frac{\left({\delta }_{+}+{\delta }_{-}\right)}{x^{{\delta }_{+}+{\delta }_{-}-1}}+$$

$$\displaystyle\xi^{2}\left(\frac{\Gamma\left(\alpha-\beta\right)}{\Gamma\left(\alpha\right)\Gamma\left(\gamma-\beta\right)}\right)^{2}\frac{\delta_{-}}{x^{2\delta_{-}-1}}=0\$$

$$\alpha =-p\text{or }\gamma -\beta =-q,$$

Full text

I Introduction

Three dimensional gravity has been a useful arena to explore gravitational models with simpler properties than their counterpart in four dimensions. Since the Weyl tensor identically vanishes in three dimensions, Einstein equations imply that the spacetime is locally flat or (A)dS depending on whether the theory has a null, negative or positive cosmological term, respectively. This implies that General Relavitity has no local degrees of freedom in vacuum, but nevertheless in the case with a negative cosmological constant, its spectrum contains black holes (the BTZ black hole) (BTZ, ). These black holes have been of fundamental importance on the tests of the holographic correspondence between physics in Anti de Sitter spacetime and that of a dual Conformal Field Theory living at the boundary (Maldacena:1997re, ), (Witten:1998qj, ), (Gubser:1998bc, ). Just to mention two examples of this relation, the entropy of the BTZ black holes can be obtained by a microscopic counting of microstates in the dual theory (Strominger, ) and their quasinormal ringing correlates precisely with relaxation time in the dual field theory at finite temperature (Horowitz:1999jd, ), (Birmingham:2001pj, ). By the end of the last decade it was also realized that the lack of local degrees of freedom of General Relavitiy in three dimensions allows to construct a parity invariant, self-interacting theory for a massive particle of spin 2, i.e., a massive gravity theory. The Einstein-Hilbert Lagrangian is supplemented by quadratic terms in the curvature, and despite the fact that the field equations are of fourth order, their linearization around flat space correctly reproduces the Fierz-Pauli equation for a massive spin 2 excitation (NMG, )-(Deser:2009hb, ). The addition of a cosmological term allows to further explore the ideas of the holographic correspondence in the presence of a massive graviton in the bulk. As it generically occurs in theories containing quadratic powers in the curvature, such theory possesses two maximally symmetric (and therefore of constant curvature) solutions. For a precise relation of the couplings both vacua coincide and actually in this case one can find asymptotically locally (A)dS black holes that are characterized by the mass, angular momentum and an extra parameter that stands for a gravitational hair (OTT, )-(Bergshoeff:2009aq, ). As shown in (Giribet:2009qz, ) the entropy of such black holes can be reproduced by counting microstates in the boundary theory, providing a new test for such relation on a theory with local degrees of freedom in the bulk. Recently, there has been a revived interest in the study of asymptotic symmetries to null surfaces, where such surface could be null infinity or the event horizon of a black hole. It is expected that these studies might shed some light on the information paradox (see e.g.(Arcioni:2003td, )-(Donnay:2015abr, ) and references therein). New Massive Gravity provides for a simple setup to carry on such studies, since the spectrum of the purely quadratic theory (that is healthy and intrinsically of fourth-order (Deser:2009hb, )) contains asymptotically locally flat, rotating black holes in 2+1 dimensions (OTT, ). Such black holes do not exist in General Relativity in vacuum. Within the realm of NMG, these black holes can be generalized to construct non-circular black objects, dubbed black flowers (Barnich:2015dvt, ), a family of metrics whose simplest representative is the rotating black hole constructed in (OTT, ) in the massless limit of NMG111It is also worth mentioning that these solutions can be embedded in the Born-Infeld extensions of New Massive Gravity BINMG by removing the Einstein-Hilbert term BFINBINMG .. The aim of this paper is to show that the equation for the massless scalar, remarkably, can be solved analytically on these backgrounds, and that the quasinormal frequencies can be found in a closed manner. As usual, we require ingoing boundary condition at the horizon. We show that at infinity the natural boundary condition that makes the action principle for the scalar to be well-defined is that the field must vanish sufficiently fast. This is one of the few known, rotating black holes, that admit such integration.

This paper is organized as follows: Section I introduces the theory of New Massive gravity and presents the asymptotically locally flat, rotating black hole solution. In Section II we consider a massless scalar field perturbation and solve it in an exact and analytic manner. We find the quasinormal modes in a closed form in the rotating case and show that as the angular momentum of the black hole decreases the imaginary part of the frequencies approach minus infinity. This is consistent with the fact that in the static case, the algebraic equations that determine the spectrum cannot be fulfilled. We therefore explore the behavior of the quasinormal frequencies in terms of the global charges of the background solution. Even though our interests is in the propagation of a scalar probe on the black hole background, let us mentions that our results could shed some light on the stability of the rotating black hole in NMG. Since the full non-linear field equations of NMG are of fourth order, their linearized version around a generic background will be in general of the same order. Nevertheless, there could be an affective quantity, constructed with second derivatives of the perturbation, which may have a second order dynamics, as it occurs in the Teukolsky equation, where the unknowns are linearized expressions for the components of the Weyl tensor in a null tetrad Chandrasekhar . It is also worth mentioning that in some particular cases in GR it is known that the dynamics of some of the modes of the full gravitational perturbation coincide with the dynamics of a scalar probe with a given mass. This happens for example in the massless topological black hole Birmingham:2006zx . A less striking connection between a scalar probe and a gravitational perturbation occurs with the Regge-Wheeler equation, which has the same functional form than the equation for a massless scalar probe in Schwarzschild background, but in this case there is a single numeric coefficient of difference in the effective Schroedinger potential (see Chandrasekhar ). Section III contains conclusions and further comments.

II Massless limit of NMG and its rotating black hole

The theory of New Massive Gravity (NMG, ), has the following action

[TABLE]

Note that the Einstein term has the wrong sign, which is necessary to obtain a ghost free theory for a massive graviton, at the linearized level around flat space. Since General Relativity does not propagate local degrees of freedom in three dimensions, we can choose the sign of the Einstein term at convenience. The massless limit is described by the purely quadratic theory which acquires a linearized conformal invariance (Bergshoeff:2009aq, ), and despite the fact of having fourth order field equations, it defines a healthy theory (Deser:2009hb, ). We take the limit in the following manner $m\rightarrow 0$ and $\kappa\rightarrow+\infty$ with $16\pi G:=\kappa^{2}m^{2}$ fixed. The field equations in vacuum are therefore given by

[TABLE]

Diffeomorphism invariance of the original action ensures the identity $\nabla^{\mu}K_{\mu\nu}\equiv 0$. This is the unique quadratic combination in three dimensions that leads to field equations of fourth order, whose trace reduces to a second order constraint (Nakasone:2009bn, ). This feature is of fundamental importance when proving the healthiness of the theory. The theory defined by the field equations (2) has the following rotating black hole solution

[TABLE]

with

[TABLE]

It can be checked that this spacetime has a vanishing Cotton tensor, and it is therefore a solution of Conformal Gravity in $2+1$ dimensions (the massless limit of Topologically Massive Gravity (TMG1, ), (TMG2, )). In such context, the static solution was originally reported in (Oliva:2009hz, ). The spacetime (3), is characterized by three parameters $b$, $a$ and $\mu$ which determine the global charges associated to the asymptotic Killing vectors $\partial_{t}$ and $\partial_{\phi}$, i.e. the mass and the angular momentum, respectively, which are given by (Barnich:2015dvt, )

[TABLE]

Note that the parameter $\mu$ does not appear in the global charges, and therefore can be interpreted as a hair of gravitational origin. Further thermodynamical properties are the temperature and the entropy of the black hole that respectively read

[TABLE]

The remaining chemical potential, the angular velocity of the horizon, vanishes. The metric (3) can be obtained from the rotating solution of Cosmological New Massive Gravity at the point of coinciding vacua (OTT, ), in the limit $l\rightarrow+\infty$, by a suitable rescalling of the integration constants. This solution can also be written in null coordinates, adapted to future null infinity as

[TABLE]

where the change of coordiantes that relates (3) and (7) is

[TABLE]

Hereafter we will use the metric in the coordinate system defined in equation (3). The event horizon is located at $r=r_{+}=\mu/b$, and the metric is asymptotically locally flat, i.e.

[TABLE]

III Massless scalar on the rotating black hole

Consider a massless scalar $\Phi$, on the rotating black hole metric (3)

[TABLE]

As usual, the symmetries of the background allow to consider the following mode decomposition

[TABLE]

Hereafter we denote $R_{\omega,n}\left(r\right)=R\left(r\right)$. One therefore obtains an ODE for the radial dependence, which acquires a simple fashion in terms of the coordinate $x$, such that

[TABLE]

that maps the domain $r\in[r_{+},+\infty[$ to the range $x\in[0,+\infty[$. The equation for $R\left(x\right)$ therefore reads

[TABLE]

with

[TABLE]

Near the event horizon, located at $x=0$, the solutions of equation (15) allow for the following asymptotic behaviors

[TABLE]

while at infinity the scalar behaves as

[TABLE]

where

[TABLE]

Here $C_{1,2}$ an $\tilde{C}_{1,2}$ are integration constants. Note that since the angular velocity of the horizon vanishes, the leading dependence of the scalar on the near horizon region does not depend on the rotation parameter $a$.

A remarkable fact is that even though (15) describes the radial dependence of a scalar on a rotating black hole, it can be solved in an analytic manner in terms of Hypergeometric functions. The solution of (15) is given by

[TABLE]

where $F$ stands for a Hypergeometric Function ${}_{2}F_{1}$, and

[TABLE]

For quasinormal modes we must require that the energy flux at the horizon must be ingoing, therefore using the fact that ${}_{2}F_{1}\left(\mu,\nu,\sigma,0\right)=1$ one sees that $C_{2}$ in (22) must vanish. Consequently, one obtains

[TABLE]

Using Kummer relations for the Hypergeometric Functions AbramowitzStegun , leads to the following asymptotic behavior at infinity ($x\rightarrow+\infty$) for (26):

[TABLE]

where $\xi$ is a non-vanishing constant. In general $\delta_{\pm}=\frac{1}{2}\left(1\pm\sqrt{1-\frac{4\omega^{2}}{b^{2}}}\right)$ will be complex numbers, related by $\delta_{+}+\delta_{-}=1$.

The dynamics of the scalar field on the rotating black hole background comes from the action principle

[TABLE]

whose variation reads

[TABLE]

which on-shell reduces to a pure boundary term

[TABLE]

Requiring the action principle to be well-defined implies that $\delta_{\text{on-shell}}I\left[\phi\right]=0$, which using the asymptotic expansion (27) implies

[TABLE]

In the limit $x\rightarrow\infty$ the first term goes to zero since $\Re(2\delta_{+}-1)>0$, the second term is a constant and the third term is divergent since $\Re(2\delta_{-}-1)<0$. The last two terms must vanish independently, and therefore we need to impose that their common factor vanishes by imposing $\Gamma(\alpha)\rightarrow\infty$ or $\Gamma(\gamma-\beta)\rightarrow\infty$. These condition lead to the equations for the spectrum

[TABLE]

where $p$ and $q$ take integer values 0,1,2,3,...\. Conditions (32) lead to the determination of the spectrum of quasinormal frequencies for ingoing boundary condition at the horizon and on-shell vanishing boundary term of the variation of the action functional at infinity (30). We can introduce the reduced quantities

[TABLE]

in terms of which the constraints (32) that determine the spectrum read,

[TABLE]

The union of the loci in the complex $\omega$-plane defined by the equations (34) and (35), defines the quasinormal spectrum of the rotating black hole. The integers $p$ and $q$, determine the overtones. As mentioned above, note that $\mathcal{C}_{1}(\hat{J},n)=\mathcal{C}_{2}(-\hat{J},-n)$ with $p$ interchanged by $q$. We observe that the constraint $\mathcal{C}_{1}=0$($\mathcal{C}_{2}=0$) can be solved only for positive(negative) values of the black hole angular momentum $\hat{J}$, but due to the symmetry mentioned above, as expected, both constraints will lie on the same locus.

Consequently, in terms of the reduced angular momentum $\hat{J}=J/\sqrt{M}$, the reduced frequencies $\hat{\omega}=\omega/\sqrt{M}$ can be conveniently summarized as

[TABLE]

where we have defined $\nu:=\sqrt{2\hat{J}^{2}+\mu}$ and set Newton’s constant $G=1/2$. $\sigma_{x}$ stands for the sign of $x$, defined such that $\sigma_{0}=1$.

Some comments on the boundary conditions are in order. For asymptotically flat black holes in four dimensions, as for example in the Kerr family, the scalar field has two possible behaviors at infinity which are plane waves that describe the out-going and the in-coming modes, and therefore the leading contribution goes as $e^{-i\omega\left(t\pm r\right)}$ with the plus(minus) sign standing for the in-coming(out-going) mode. The background of the rotating asymptotically, locally flat, black hole of New Massive Gravity allows different behaviors at infinity, which are given in (27). Even in terms of the proper radial distance at infinity $\rho=\frac{2}{\sqrt{b}}\sqrt{r}$, one would have $x\sim\rho^{2}$, and the behavior (27) would still have a power law fashion and would not be compatible with the interpretation of an out-going or in-coming plane wave. The situation at hand is actually more similar to the analysis of quasinormal modes in asymptotically AdS black holes where the leading asymptotic behaviors are $r^{-\Delta_{\pm}}$ and $\Delta_{\pm}$ are given in terms of the mass of the field and the dimension of the spacetime, defining the conformal weights of the operator dual to the scalar on the boundary CFT. Therefore, with the asymptotic behaviors given by equation (27) a natural physical condition that allows to compute the quasinormal frequencies is to require that the field at infinity decays fast enough as to lead to a well defined variational principle (30) and this is in fact the strategy we have followed.

Below we present plots of the QNM spectrum of the scalar probe in terms of the physical variables, as well as closed analytic expressions for the asymptotic behavior of the frequencies.

III.1 Static case and asymptotic frequencies

Let’s analyze first the static case. We set $J=0$ in (34) and (35). Then from (34) one obtains two equations for the real and imaginary part of the constraint, whose real part reads

[TABLE]

This equation does not have solution for any value of $\omega=\omega_{r}+i\omega_{i}$ (with $\omega_{r,i}$ finite), and therefore (34) cannot be fulfilled. A similar analysis of the real part of (35) leads to the same equation (with $p$ changed by $q$). This is a peculiar feature of the static solution: there are no massless scalar field quasinormal modes for the asymptotically flat, static black hole with gravitational hair in NMG. The constraints (34) and (35) can be mapped to the Riemann sphere, where it can be seen that at the leading order when $\omega_{i}\rightarrow-\infty$ (for constant $\omega_{r}$) the constraint is indeed solved since (37) reads

[TABLE]

which vanishes when $\omega_{i}\rightarrow-\infty$. This result is recovered from the quasinormal modes of the rotating case. We will see that as $J\rightarrow 0$, the imaginary part of the quasinormal frequencies tend to $-\infty$, and therefore the quasinormal modes acquire an infinity damping in this limit, and migrate to the north pole of the Riemann Sphere222Same conclusion is obtained if one starts with the asymptotically locally flat, static black hole from the scratch. Quasinormal modes of scalar and Dirac fields on the asymptotically AdS, hairy black hole of NMG have been explored numerically and perturbatively in Kwon:2011ey and Gonzalez:2014voa ..

III.2 The quasinormal spectrum

In this section we present different plots that clarify the behavior of the spectra as a function of the parameters of the problem, with descriptive captions.

III.3 Asymptotic frequencies

Finally let’s analyze two informative asymptotic behaviors for the spectrum. First, for large angular momentum of the field the asymptotic expansion $n\rightarrow+\infty$ reads

[TABLE]

It is interesting to notice that this behavior is reminiscent of that of the quasinormal modes of a scalar field with large angular momentum $j$ on Schwarzschild black hole. In that case one has Ferrari:1984zz

[TABLE]

It is also interesting to put special attention to the case with arbitrarily large $J$, since this cannot be reached within the Kerr family. (It can be reached in Myers-Perry in $d>5$ Myers:1986un ). Quasinormal frequencies $\hat{\omega}$ in the limit $\hat{J}\rightarrow\infty$ becomes purely imaginary (damped) and depend only on the overtone integer

[TABLE]

The fundamental mode in this regime is $\hat{\omega}_{\hat{J}\rightarrow\infty}^{0}=-2\,i$. The function $\hat{\omega}_{\hat{J}\rightarrow\infty}(p)$ at large $p$ approaches a straight line with slope $2\sqrt{2}$.

IV Conclusions

In this paper we have shown that a massless scalar probe on the asymptotically locally flat, rotating black hole of pure NMG can be analytically solved leading to an exact expression for the quasinormal mode frequencies. We have imposed ingoing boundary conditions at the horizon and a decay at infinity which is fast enough leading to a well defined variational principle for the scalar probe. For a given black hole, parametrized by the mass, angular momentum and gravitational hair, the quasinormal spectrum is defined by the angular momentum of the field as well as an integer counting overtones. We showed that under these boundary conditions the imaginary part of the frequencies is always negative, not leading to any superradiance, as expected for a massless scalar. Even more, when the angular momentum of the black hole approaches zero, the imaginary part of the frequencies tends to minus infinity leading to infinitely damped modes. Remarkably, this is one of the few examples of rotating black holes in which the quasinormal frequencies for a scalar probe can be obtained in an analytic, closed form333The other examples being BTZ (ExactBTZ, ), the so-called subtracted geometries (subtracted, ) and near horizon geometries of rotating black holes (nh, ). For a thorough review on QNMs see (Berti:2009kk, )..

Since the rotating black hole, even in the case with cosmological constant, has vanishing Cotton tensor (OTT, ), (Bergshoeff:2009aq, ), it is a conformally flat spacetime, therefore the solution for a conformal scalar could also be obtained in an analytic manner by applying a conformal transformation to the solution of the free, massless scalar of 3D Minkowski or AdS space (since the latter is also conformally flat). Once the mapping is done on the general solution, one would have to impose the boundary conditions at the horizon and infinity taking care of the potential divergences induced by the conformal mapping that leads to the black hole.

The stable propagation of scalar probes on the rotating background may lead to the existence of rotating black holes with both scalar and gravitation hair. In NMG static black holes with self-interacting scalar hair already exist, with Lifshitz and AdS asymptotics (Correa:2014ika, ), Ayon-Beato:2015jga , and it would be interesting to construct rotating black holes with these hairs.

Aknowledgements

The authors want to thank Gaston Giribet and Nicolas Grandi for enlightening comments. This work was partially supported by FONDECYT grants 1170279 and 1181047. O.F. is supported by the Project DINREG 19/2018 of the Dirección de Investigación of the Universidad Católica de la Santísima Concepción.