This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Stationary scalar clouds around a rotating BTZ-like black hole in the Einstein-bumblebee gravity

Fangli Quan1, Fengjiao Li1, Qiyuan Pan1,2111[email protected], Mengjie Wang1222[email protected], and Jiliang Jing1,2333[email protected] 1 Department of Physics, Key Laboratory of Low Dimensional Quantum Structures and Quantum Control of Ministry of Education, Institute of Interdisciplinary Studies, and Synergetic Innovation Center for Quantum Effects and Applications, Hunan Normal University, Changsha, Hunan 410081, China 2 Center for Gravitation and Cosmology, College of Physical Science and Technology, Yangzhou University, Yangzhou 225009, China
Abstract

Abstract

We have studied stationary clouds of massive scalar fields around a rotating BTZ-like black hole in the Einstein-bumblebee gravity, by imposing the Robin type boundary conditions at the AdS boundary. We establish, by scanning the parameter space, the existence of fundamental stationary scalar clouds (i.e.i.e., the overtone number n=0n=0). In particular, we observe that the Lorentz symmetry breaking parameter ss and the quantum number kk play an opposite role in determining scalar clouds, which indicates the existence of degenerate scalar clouds. To illustrate the fact that scalar clouds may only be supported for the n=0n=0 case, we have analyzed the impact of various parameters on scalar quasinormal modes. It is shown that the Lorentz symmetry breaking parameter ss does not change the superradiance condition, and superradiant instabilities only appear for the fundamental modes. Our work shows that the Lorentz symmetry breaking provides richer physics in stationary scalar clouds around black holes.

pacs:
04.70.-s, 04.70.Bw, 97.60.Lf

I Introduction

In 1971, Ruffini and Wheeler proposed the no-hair conjecture [1], which is one of the most important developments in black hole physics, stating that stationary black holes should be characterized by only three parameters: mass, charge and angular momentum [2, 3]. This pioneering work has led to many investigations concerning the new hair of black holes. Based on the superradiance mechanism, Hod first analyzed the dynamics of a test massive scalar field surrounding Kerr black holes [4, 5] and further extended to the case of Kerr-Newman black holes [6, 7], and found out that there exist the stationary and regular scalar field configurations surrounding the realistic rotating black holes which are coined as stationary scalar clouds. Herdeiro et al. presented a family of solutions of Einstein’s gravity minimally coupled to a complex massive scalar field, which describes asymptotically flat, spinning black holes with scalar hair [8], and made a thorough investigation of the scalar clouds, which are bound states at the threshold of the superradiant instability, due to neutral (charged) scalar fields around Kerr(-Newman) black holes [9, 10, 11]. Sampaio et al. considered the Proca fields on the background of Reissner-Nordström black holes and obtained the marginal Proca cloud configurations [12], and the authors of Ref. [13] completed the analysis of bosonic fields on Kerr-anti-de Sitter (Kerr-AdS) black holes by studying Maxwell perturbations and constructed stationary Maxwell clouds with Robin boundary conditions. Extending the investigation to the excited states, Wang et al. found that, in contrast to Kerr black holes with ground state scalar hair, the first excited Kerr black holes with scalar hair have two types of nodes, including radial nr=1n_{r}=1 and angular nθ=1n_{\theta}=1 nodes [14]. For a simple, geometrically elegant, black hole solution of three dimensional general relativity, Ferreira et al. showed that the BTZ black hole can support stationary scalar clouds of a massive scalar field by using appropriate Robin boundary conditions [15, 16]. Along this line, there have been accumulated interest to study various cloud configurations around black holes in the Einstein gravity [17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29] and modified gravity [30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43].

In this work, we focus on the stationary scalar clouds around a black hole in the Einstein-bumblebee gravity, which is one of simple effective theories with the Lorentz violation [44]. In the Einstein-bumblebee model, the Lorentz violation arises from the dynamics of the bumblebee vector field BμB_{\mu} [45, 46]. Casana et al. found the first black hole solution in the Einstein-bumblebee gravity, i.e.i.e., the exact Schwarzschild-like black hole solution, and presented the effects of the spontaneous Lorentz symmetry breaking on some classic tests, including the advance of perihelion, the bending of light, and Shapiro’s time delay [47]. More recently, Ding et al. obtained a three-dimensional rotating BTZ-like black hole solution in the Einstein-bumblebee gravity theory, and found the effects of the Lorentz breaking constant on its thermodynamics and central charges [48]. Considering a scalar perturbation around such a rotating BTZ-like black hole in the Einstein-bumblebee gravity, Chen et al. investigated the quasinormal modes of this rotating BTZ-like black hole and pointed out that the Lorentz symmetry breaking parameter imprints only in the imaginary parts of the quasinormal frequencies for the right-moving and left-moving modes [49]. Moreover, Ge et al. constructed the mass ladder operators for the static BTZ-like black hole in the Einstein-bumblebee gravity, and probed the quasinormal frequencies of the mapped modes by the mass ladder operators for a scalar perturbation under the Dirichlet and Neumann boundary conditions [50]. Since the three-dimensional gravity certainly offers potential insights into quantum gravity, it is of great interest to generalize the investigation to the stationary scalar clouds around a rotating BTZ-like black hole in the Einstein-bumblebee gravity. We will examine the influence of the spontaneous Lorentz symmetry breaking on the stationary scalar clouds, and check whether it is possible to distinguish between the BTZ black hole in the Einstein gravity and BTZ-like black hole in the Einstein-bumblebee gravity based on the scalar field configurations.

This work is organized as follows. In Section II, we briefly review the rotating BTZ-like black holes in the Einstein-bumblebee gravity, and present the massive Klein-Gordon equation and the appropriate boundary conditions for the non-extremal BTZ-like black holes. In Section III, we analyze stationary scalar configurations around the rotating BTZ-like black hole in the Einstein-bumblebee gravity. In Section IV, we investigate quasinormal modes for a massive scalar field which is governed by the Klein-Gordon equation. We conclude in the last section with our main results. For completeness, in Appendix A, we study massive Klein-Gordon equation in the extremal BTZ-like black hole.

II Background geometry, massive scalar fields and boundary conditions

II.1 Rotating BTZ-like black holes

We start with the action of the Einstein-bumblebee gravity with a negative cosmological constant Λ=1/2\Lambda=-1/\ell^{2} in three dimensions [51, 44, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62]

S=d3xg[R2Λ2κ+ξ2κBμBνRμν14BμνBμνV(BμBμ±b2)],\displaystyle S=\int d^{3}x\sqrt{-g}\left[\frac{R-2\Lambda}{2\kappa}+\frac{\xi}{2\kappa}B^{\mu}B^{\nu}R_{\mu\nu}-\frac{1}{4}B^{\mu\nu}B_{\mu\nu}-V\left(B_{\mu}B^{\mu}\pm b^{2}\right)\right], (1)

where RR represents the Ricci scalar, κ=8πG/c3\kappa=8\pi G/c^{3} is a constant related to the three-dimensional Newtonian constant GG, ξ\xi denotes the coupling constant with the dimension of M1M^{-1}, and BμB_{\mu} is the bumblebee vector field with the strength Bμν=μBννBμB_{\mu\nu}=\partial_{\mu}B_{\nu}-\partial_{\nu}B_{\mu}. The potential VV of the bumblebee field has a minimum at BμBμ±b2=0B_{\mu}B^{\mu}\pm b^{2}=0, yielding a nonzero vacuum expectation value Bμ=bμ\langle B_{\mu}\rangle=b_{\mu} with bμbμ=b2b_{\mu}b^{\mu}=\mp b^{2} and leading to the violation of the U(1)U(1) symmetry, where bb is a real positive constant and the signs “±\pm” mean that bμb_{\mu} is timelike or spacelike.

By varying the action (1) with respect to the metric, the authors of Ref. [48] obtained a rotating BTZ-like black hole solution

ds2=f(r)dt2+1+sf(r)dr2+r2(dφj2r2dt)2,ds^{2}=-f(r)dt^{2}+\frac{1+s}{f(r)}dr^{2}+r^{2}\left(d\varphi-\frac{j}{2r^{2}}dt\right)^{2}, (2)

with

f(r)=r22M+j24r2,f(r)=\frac{r^{2}}{\ell^{2}}-M+\frac{j^{2}}{4r^{2}}, (3)

where MM and jj are the mass and spin parameters of the black hole, ss is the spontaneous Lorentz symmetry breaking parameter, and \ell is the AdS radius. The determinant of the metric and the Kretschmann scalar may be obtained as

g=(s+1)r2,RμνρλRμνρλ=124(1+s)2,g=-(s+1)r^{2}\;,\;\;\;R^{\mu\nu\rho\lambda}R_{\mu\nu\rho\lambda}=\frac{12}{\ell^{4}(1+s)^{2}}, (4)

which means that the spacetime (2) has a singularity when s=1s=-1. To keep the spacetime which is Lorentzian and is free of curvature singularity we, therefore, require the constraint s>1s>-1. The black hole solution (2) has an event horizon and an inner horizon, i.e.i.e.,

r±2=22(M±M2j22),r^{2}_{\pm}=\frac{\ell^{2}}{2}\left(M\pm\sqrt{M^{2}-\frac{j^{2}}{\ell^{2}}}\right), (5)

which are exactly the same with the BTZ black hole solution and, therefore, are independent of the spontaneous Lorentz symmetry breaking parameter ss. For the extremal BTZ-like black hole, we have M2=j2/2M^{2}=j^{2}/\ell^{2} which leads to the coincidence of the event and inner horizons, i.e.i.e., r+=r=M/2r_{+}=r_{-}=\ell\sqrt{M/2}. As shown in Appendix A, we observe that extremal BTZ-like black holes do not support stationary scalar clouds. Thus, we will focus on the non-extremal BTZ-like black holes in the following.

In addition, the rotating BTZ-like black hole has an ergoregion rerg=Mr_{erg}=\ell\sqrt{M} and the angular velocity of the horizon

ΩH=rr+,\Omega_{H}=\frac{r_{-}}{\ell r_{+}}, (6)

and the black hole mass parameter may be rewritten as

M=r+2+r22=r+22(1+2ΩH2),M=\frac{r^{2}_{+}+r^{2}_{-}}{\ell^{2}}=\frac{r^{2}_{+}}{\ell^{2}}\left(1+\ell^{2}\Omega^{2}_{H}\right), (7)

which are, again, exactly the same with the BTZ case.

II.2 Massive Klein-Gordon equation

In the rotating BTZ-like black hole background, a massive scalar field Φ\Phi, with the mass μ/l\mu/l, evolves according to the Klein-Gordon equation

1gxμ(ggμνxν)Φμ22Φ=0,\frac{1}{\sqrt{-g}}\frac{\partial}{\partial x^{\mu}}\left(\sqrt{-g}g^{\mu\nu}\frac{\partial}{\partial x^{\nu}}\right)\Phi-\frac{\mu^{2}}{\ell^{2}}\Phi=0, (8)

where μ\mu is dimensionless. The spacetime isometries allow full seperation of variables, and by taking the ansatz

Φ(t,r,φ)=eiωt+ikφϕ(r),\Phi(t,r,\varphi)=e^{-i\omega t+ik\varphi}\phi(r), (9)

the radial equation of massive scalar fields turns into

d2ϕ(r)dr2+[1r+1f(r)df(r)dr]dϕ(r)dr+1+sf(r)[k2r2μ22+(jk2r2ω)24r4f(r)]ϕ(r)=0,\frac{d^{2}\phi(r)}{dr^{2}}+\left[\frac{1}{r}+\frac{1}{f(r)}\frac{df(r)}{dr}\right]\frac{d\phi(r)}{dr}+\frac{1+s}{f(r)}\left[-\frac{k^{2}}{r^{2}}-\frac{\mu^{2}}{\ell^{2}}+\frac{(jk-2r^{2}\omega)^{2}}{4r^{4}f(r)}\right]\phi(r)=0, (10)

where ω\omega and kk are the frequency and the angular quantum number of scalar fields.

To solve the radial equation (10), one may introduce a new radial coordinate zz for non-extremal BTZ-like black holes as

zr2r+2r2r2,r+r,z\equiv\frac{r^{2}-r^{2}_{+}}{r^{2}-r^{2}_{-}},\qquad r_{+}\neq r_{-}, (11)

which transforms the exterior region from r(r+,+)r\in(r_{+},+\infty) to z(0,1)z\in(0,1), and Eq. (10) becomes

z(1z)d2ϕ(z)dz2+(1z)dϕ(z)dz+(Az+B+C1z)ϕ(z)=0,z(1-z)\frac{d^{2}\phi(z)}{dz^{2}}+(1-z)\frac{d\phi(z)}{dz}+\left(\frac{A}{z}+B+\frac{C}{1-z}\right)\phi(z)=0, (12)

with

A=4(ωr+kr/)24(r+2r2)2(1+s),B=4(ωrkr+/)24(r+2r2)2(1+s),C=μ2(1+s)4.A=\frac{\ell^{4}(\omega r_{+}-kr_{-}/\ell)^{2}}{4(r^{2}_{+}-r^{2}_{-})^{2}}(1+s),\;\;B=-\frac{\ell^{4}(\omega r_{-}-kr_{+}/\ell)^{2}}{4(r_{+}^{2}-r_{-}^{2})^{2}}(1+s),\;\;C=-\frac{\mu^{2}(1+s)}{4}. (13)

Then, letting ϕ(z)=zα(1z)βF(z)\phi(z)=z^{\alpha}(1-z)^{\beta}F(z), one obtains the standard hypergeometric equation [63]

z(1z)d2F(z)dz2+[c(1+a+b)z]dF(z)dzabF(z)=0,z(1-z)\frac{d^{2}F(z)}{dz^{2}}+\left[c-(1+a+b)z\right]\frac{dF(z)}{dz}-abF(z)=0, (14)

with

aβi(ω+k)1+s2(r++r),bβi(ωk)1+s2(r+r),c1+2α,αi2r+2(r+2r2)(ωkΩH)1+s,β12[1+1+μ2(1+s)],\begin{split}a\equiv&\beta-i\ell\frac{(\omega\ell+k)\sqrt{1+s}}{2(r_{+}+r_{-})},\qquad b\equiv\beta-i\ell\frac{(\omega\ell-k)\sqrt{1+s}}{2(r_{+}-r_{-})},\qquad c\equiv 1+2\alpha,\\ \alpha\equiv&-i\frac{\ell^{2}r_{+}}{2(r^{2}_{+}-r^{2}_{-})}(\omega-k\Omega_{H})\sqrt{1+s},\qquad\beta\equiv\frac{1}{2}\left[1+\sqrt{1+\mu^{2}(1+s)}\right],\end{split} (15)

where F(z)F(z) is the Gaussian hypergeometric function.

II.3 Boundary conditions

To look for scalar clouds and scalar quasinormal modes by solving Klein-Gordon equation, we have to employ an ingoing wave boundary condition at the event horizon and the Robin type vanishing energy flux boundary condition at infinity, following the Ref. [64, 65, 66, 15, 16].

At the horizon z=0z=0, one obtains the general solution for ϕ(z)\phi(z)

ϕ(z)=Azα(1z)βF(a,b;c;z)+Bzα(1z)βF(ac+1,bc+1;2c;z),\phi(z)=Az^{\alpha}(1-z)^{\beta}F(a,b;c;z)+Bz^{-\alpha}(1-z)^{\beta}F(a-c+1,b-c+1;2-c;z), (16)

where AA and BB are two integration constants, and an ingoing wave boundary condition implies B=0B=0.

At infinity z1z\to 1, the general solution of ϕ(z)\phi(z) is given by

ϕ(z)=C(D)ϕ(D)(z)+C(N)ϕ(N)(z)=A[Γ(c)Γ(cab)Γ(ca)Γ(cb)ϕ(D)(z)+Γ(c)Γ(a+bc)Γ(a)Γ(b)ϕ(N)(z)],\begin{split}\phi(z)=&C^{(D)}\phi^{(D)}(z)+C^{(N)}\phi^{(N)}(z)\\ =&A\left[\frac{\Gamma(c)\Gamma(c-a-b)}{\Gamma(c-a)\Gamma(c-b)}\phi^{(D)}(z)+\frac{\Gamma(c)\Gamma(a+b-c)}{\Gamma(a)\Gamma(b)}\phi^{(N)}(z)\right],\end{split} (17)

with two linearly independent solutions

ϕ(D)(z)=zα(1z)βF(a,b;a+bc+1;1z),ϕ(N)(z)=zα(1z)1βF(ca,cb;cab+1;1z),\begin{split}\phi^{(D)}(z)=&z^{\alpha}(1-z)^{\beta}F(a,b;a+b-c+1;1-z),\\ \phi^{(N)}(z)=&z^{\alpha}(1-z)^{1-\beta}F(c-a,c-b;c-a-b+1;1-z),\end{split} (18)

where C(D)C^{(D)} and C(N)C^{(N)} are two complex constants, and the parameters aa, bb, cc, α\alpha and β\beta are given in Eq. (15). Here ϕ(D)(z)\phi^{(D)}(z) and ϕ(N)(z)\phi^{(N)}(z) are introduced to denote the Dirichlet and Neumann solutions [15], and to keep these two solutions square-integrable near infinity, we will concentrate on the parameter space

1<s<(1+1μ2),μ2<0.-1<s<-\left(1+\frac{1}{\mu^{2}}\right)\;,\;\;\;\mu^{2}<0. (19)

By calculating the flux of energy at rr\to\infty (z1z\to 1)

=limrr𝑑ϕggrrTrt,\mathcal{F}=\lim\limits_{r\to\infty}\int_{\sum_{r}}d\phi\sqrt{-g}g^{rr}T_{rt}\;, (20)

where r\sum_{r} and TμνT_{\mu\nu} are the hypersurface of constant rr and the energy momentum tensor of scalar fields, the vanishing of energy flux [67, 68, 15] leads to

ϕ=cos(ζ)ϕ(D)+sin(ζ)ϕ(N),\phi=\cos(\zeta)\phi^{(D)}+\sin(\zeta)\phi^{(N)}, (21)

where the boundary parameter ζ\zeta satisfies ζ[0,π)\zeta\in\left[0,\pi\right). It is shown clearly that ζ=0\zeta=0 corresponds to the Dirichlet boundary condition and ζ=π/2\zeta=\pi/2 corresponds to the Neumann boundary condition.

By comparing Eq. (17) with Eq. (21), we have the boundary condition

tanζ=Γ(a+bc)Γ(ca)Γ(cb)Γ(cab)Γ(a)Γ(b).\tan\zeta=\dfrac{\Gamma(a+b-c)\Gamma(c-a)\Gamma(c-b)}{\Gamma(c-a-b)\Gamma(a)\Gamma(b)}\;. (22)

III Stationary scalar clouds

In this section, we investigate stationary scalar configurations around the rotating BTZ-like black hole in the Einstein-bumblebee gravity. For this case, the resonance condition should be satisfied, i.e.i.e., the field’s frequency equals to the critical frequency

ω=ωc=kΩH,\omega=\omega_{c}=k\Omega_{H}, (23)

which leads to α=0\alpha=0 and c=1c=1, a=βi1+sk/(2r+)a=\beta-i\ell\sqrt{1+s}k/(2r_{+}) in Eq. (15). Then the solutions in Eq. (18) reduce to

ϕ(D)(z)\displaystyle\phi^{(D)}(z) =\displaystyle= (1z)βF(a,a;2β;1z),\displaystyle(1-z)^{\beta}F(a,a^{*};2\beta;1-z), (24)
ϕ(N)(z)\displaystyle\phi^{(N)}(z) =\displaystyle= (1z)1βF(1a,1a;22β;1z),\displaystyle(1-z)^{1-\beta}F(1-a,1-a^{*};2-2\beta;1-z),

where β\beta is given in Eq. (15), and the boundary condition (22) for scalar clouds becomes

tan(ζ)=Γ(2β1)Γ(1a)Γ(1a)Γ(12β)Γ(a)Γ(a)=Γ(1+(1+s)μ2)|Γ(12121+(1+s)μ2+ik1+s2r+)|2Γ(1+(1+s)μ2)|Γ(12+121+(1+s)μ2+ik1+s2r+)|2.\begin{split}\tan(\zeta)=\frac{\Gamma(2\beta-1)\Gamma(1-a)\Gamma(1-a^{*})}{\Gamma(1-2\beta)\Gamma(a)\Gamma(a^{*})}=\frac{\Gamma(\sqrt{1+(1+s)\mu^{2}})|\Gamma(\frac{1}{2}-\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}+i\frac{k\ell\sqrt{1+s}}{2r_{+}})|^{2}}{\Gamma(-\sqrt{1+(1+s)\mu^{2}})|\Gamma(\frac{1}{2}+\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}+i\frac{k\ell\sqrt{1+s}}{2r_{+}})|^{2}}.\end{split} (25)

By scanning the parameter space of the system for given values of ss, μ2\mu^{2}, kk, r+r_{+} and ζ\zeta, one may obtain scalar clouds. Here, for concreteness we take μ2=0.65\mu^{2}=-0.65, ζ=0.9π\zeta=0.9\pi and, without loss of generality, we take the AdS radius =1\ell=1 to measure all other physical quantities. One should further note that the boundary parameter ζ\zeta is bounded in the range [ζ,π)[\zeta_{*},\pi), where

ζ=arctan(Γ(1+(1+s)μ2)|Γ(12121+(1+s)μ2)|2Γ(1+(1+s)μ2)|Γ(12+121+(1+s)μ2)|2),\displaystyle\zeta_{*}=\arctan\left(\frac{\Gamma(\sqrt{1+(1+s)\mu^{2}})|\Gamma(\frac{1}{2}-\frac{1}{2}\sqrt{1+(1+s)\mu^{2}})|^{2}}{\Gamma(-\sqrt{1+(1+s)\mu^{2}})|\Gamma(\frac{1}{2}+\frac{1}{2}\sqrt{1+(1+s)\mu^{2}})|^{2}}\right), (26)

satisfying ζ(π2,π)\zeta_{*}\in(\frac{\pi}{2},\pi). Since ζ\zeta_{*} increases as the Lorentz symmetry breaking parameter ss increases, the allowed parameter space decreases as the increase of ss.

To explore the impact of the Lorentz symmetry breaking parameter ss on scalar clouds, in Fig. 1 we present scalar clouds for the case of n=0n=0, ζ=0.9π\zeta=0.9\pi and for various Lorentz symmetry breaking parameters ss (s=0.5,0.1,0,0.1,0.5s=-0.5,-0.1,0,0.1,0.5) with fixed k=1k=1 (left panel) and k=2k=2 (right panel). It is shown that the clouds exist for smaller background mass with fixed background angular velocity by increasing ss with fixed kk. In Fig. 2 we exhibit scalar clouds for the case of n=0n=0, ζ=0.9π\zeta=0.9\pi and for various quantum number kk (k=1,2,3,4k=1,2,3,4) with fixed s=0.5s=-0.5 (left panel) and s=0.5s=0.5 (right panel). It is shown that the clouds exist for larger background mass with fixed background angular velocity by increasing kk with fixed ss. The results shown in the above figures indicate the Lorentz symmetry breaking parameter ss and the quantum number kk have entirely opposite effects on the existence lines of nodeless scalar clouds. One should further notice that the Lorentz symmetry breaking parameter ss has no effect on the extremal lines (the dashed black curve in the figure), which is consistent with Eq. (6). This is completely different from the effect of the MOG parameter α\alpha on stationary scalar clouds in Kerr-MOG black holes, where the increase of parameter α\alpha makes the extremal lines lower [37].

Refer to caption
Refer to caption
Figure 1: Existence lines of nodeless scalar clouds (n=0n=0) for various Lorentz symmetry breaking parameter ss with fixed quantum numbers k=1k=1 (left), k=2k=2 (right), in the mass vs horizon angular velocity parameter space of BTZ-like black holes. The dashed black curve corresponds to extremal BTZ-like black holes with ΩH=1/\Omega_{H}=1/\ell and non-extremal BTZ-like black holes exist in the shaded region.
Refer to caption
Refer to caption
Figure 2: Existence lines of nodeless scalar clouds (n=0n=0) for various quantum number kk with fixed s=0.5s=-0.5 (left), s=0.5s=0.5 (right), in the mass vs horizon angular velocity parameter space of BTZ-like black holes. The dashed black curve corresponds to extremal BTZ-like black holes with ΩH=1/\Omega_{H}=1/\ell and non-extremal BTZ-like black holes exist in the shaded region.

Based on the above observation, in Fig. 3 we present scalar clouds in the parameter space by taking a few pairs of (k,sk,s), including (k=1,s=0)(k=1,s=0), (k=2,s=0.36794)(k=2,s=0.36794), (k=3,s=0.42830)(k=3,s=0.42830) and (k=4,s=0.45342)(k=4,s=0.45342). It is interesting to note that the existence lines for all the cases we considered here are the same, and we term such configurations as degenerate clouds. This indicates that infinite degenerate clouds may be constructed by increasing kk and the corresponding ss is bounded by s=(1+1/μ2)s=-(1+1/\mu^{2}). The existence of degenerate clouds raises a question, i.e.i.e. how to distinguish BTZ black holes and BTZ-like black holes. To solve this issue, we further explore the radial profile of degenerate clouds in Fig. 4, by taking the exactly the same parameters given in Fig. 3. It is shown clearly that the radial profiles of degenerate clouds ϕ\phi are completely different, indicating that degenerate clouds may be discerned by the radial profile.

Refer to caption
Figure 3: Degenerate existence lines of nodeless scalar clouds (n=0n=0) for several pairs of the quantum number kk and Lorentz symmetry breaking parameters ss, in the mass vs horizon angular velocity parameter space of BTZ-like black holes, with μ2=0.65\mu^{2}=-0.65, r+=0.40182r_{+}=0.40182, =1\ell=1 and with the Robin boundary condition ζ=0.9π\zeta=0.9\pi at infinity.
Refer to caption
Figure 4: The radial profile of degenerate clouds on BTZ-like black holes is presented in terms of zz, with μ2=0.65\mu^{2}=-0.65, r+=0.40182r_{+}=0.40182, =1\ell=1 and with the Robin boundary condition ζ=0.9π\zeta=0.9\pi at infinity.

Before proceeding, we would like to stress that in the above we only find the nodeless scalar clouds (n=0n=0), i.e.i.e., the ground state stationary clouds, which is different from the Kerr case [8, 14]. Considering that the clouds only exist at the threshold of superradiance [69], we will study the quasinormal modes in the BTZ-like black hole and find out the corresponding reasons.

IV Quasinormal modes

For a massive scalar field Φ\Phi satisfying the Klein-Gordon equation (8), from Eqs. (17) and (21) we can obtain quasinormal frequencies by solving the following equation

tan(ζ)=Γ(1+(1+s)μ2)Γ(12121+(1+s)μ2i1+sωk2(r+r))Γ(1+(1+s)μ2)Γ(12+121+(1+s)μ2i1+sω+k2(r++r))×Γ(12121+(1+s)μ2i1+sω+k2(r++r))Γ((12+121+(1+s)μ2i1+sωk2(r+r)),\begin{split}\tan(\zeta)=&\frac{\Gamma\left(\sqrt{1+(1+s)\mu^{2}}\right)\Gamma\left(\frac{1}{2}-\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}-i\ell\sqrt{1+s}\frac{\omega\ell-k}{2(r_{+}-r_{-})}\right)}{\Gamma\left(-\sqrt{1+(1+s)\mu^{2}}\right)\Gamma\left(\frac{1}{2}+\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}-i\ell\sqrt{1+s}\frac{\omega\ell+k}{2(r_{+}+r_{-})}\right)}\\ &\times\frac{\Gamma\left(\frac{1}{2}-\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}-i\ell\sqrt{1+s}\frac{\omega\ell+k}{2(r_{+}+r_{-})}\right)}{\Gamma\left((\frac{1}{2}+\frac{1}{2}\sqrt{1+(1+s)\mu^{2}}-i\ell\sqrt{1+s}\frac{\omega\ell-k}{2(r_{+}-r_{-})}\right)},\end{split} (27)

where the Robin parameter satisfies ζ[0,π)\zeta\in[0,\pi). Obviously, if ζ=0\zeta=0, i.e.i.e., the so-called Dirichlet boundary condition, we can solve Eq. (27) exactly, and get the left-moving frequencies

ω(D),L=kir++r1+s2[2n+1+1+μ2(1+s)],\omega^{(D),L}=-\frac{k}{\ell}-i\frac{r_{+}+r_{-}}{\sqrt{1+s}~{}\ell^{2}}\left[2n+1+\sqrt{1+\mu^{2}(1+s)}\right], (28)

and the right-moving frequencies

ω(D),R=kir+r1+s2[2n+1+1+μ2(1+s)],\omega^{(D),R}=\frac{k}{\ell}-i\frac{r_{+}-r_{-}}{\sqrt{1+s}~{}\ell^{2}}\left[2n+1+\sqrt{1+\mu^{2}(1+s)}\right], (29)

which reduce to the quasinormal frequencies given in Ref. [16] for the BTZ black hole if the Lorentz symmetry breaking parameter s=0s=0. Similarly, if ζ=π/2\zeta=\pi/2, i.e.i.e., the so-called Neumann boundary condition, we have

ω(N),L\displaystyle\omega^{(N),L} =\displaystyle= kir++r1+s2[2n+11+μ2(1+s)],\displaystyle-\frac{k}{\ell}-i\frac{r_{+}+r_{-}}{\sqrt{1+s}~{}\ell^{2}}\left[2n+1-\sqrt{1+\mu^{2}(1+s)}\right], (30)
ω(N),R\displaystyle\omega^{(N),R} =\displaystyle= kir+r1+s2[2n+11+μ2(1+s)].\displaystyle\frac{k}{\ell}-i\frac{r_{+}-r_{-}}{\sqrt{1+s}~{}\ell^{2}}\left[2n+1-\sqrt{1+\mu^{2}(1+s)}\right]. (31)

For other values of Robin parameter ζ\zeta, we have to solve Eq. (27) numerically and obtain the numerical results for the quasinormal frequencies.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: Imaginary parts of some quasinormal frequencies as a function of ζ/π\zeta/\pi with different Lorentz symmetry breaking parameters ss for BTZ-like black holes and scalar field with μ2=0.65\mu^{2}=-0.65, r+=5r_{+}=5, r=3r_{-}=3, =1\ell=1 and k=1k=1.

In Fig. 5, we present the imaginary parts of some selected quasinormal frequencies as a function of ζ/π\zeta/\pi with different Lorentz symmetry breaking parameters for a given black hole configuration (the fixed \ell, r+r_{+} and rr_{-}) and scalar field configuration (the fixed μ2\mu^{2} and kk). The red dots denote quasinormal frequencies corresponding to the Dirichlet boundary condition (ζ=0)(\zeta=0) and the red triangles correspond to those for the Neumann boundary condition (ζ=π2)(\zeta=\frac{\pi}{2}), which can be served as a check of the numerical computation. We can see clearly that, as a function of ζ\zeta, the mode continuously connects the case of Dirichlet boundary condition with that of Neumann boundary condition for quasinormal frequencies. However, all the imaginary parts of quasinormal frequencies are negative except for the fundamental frequencies of type L, corresponding to the nodeless modes (n=0n=0). For the nodeless modes, the imaginary parts of frequencies always change from negative to positive as the parameter ζ\zeta increases, which means that there exist the so-called superradiance zero-modes, occurring in between decaying modes and superradiantly amplified modes [15]. It should be noted that, regardless of the Lorentz symmetry breaking parameter ss, the superradiance zero-modes only appear in the nodeless modes (n=0n=0) of type L, which indicates that the stationary scalar clouds only exist in this ground state since the resonance condition corresponds precisely to the threshold of the superradiant instability [69]. This supports the observation shown in Figs. 1, 2 and 3 where we can only obtain the existence lines of nodeless scalar clouds (n=0n=0).

Refer to caption
Figure 6: Imaginary parts of some quasinormal frequencies of nodeless modes (n=0n=0) as a function of ζ/π\zeta/\pi with different Lorentz symmetry breaking parameters ss for BTZ-like black holes and scalar field with μ2=0.65\mu^{2}=-0.65, r+=5r_{+}=5, r=3r_{-}=3, =1\ell=1 and k=1k=1. The dashed purple curve corresponds to the case of sapp=0.53846s_{app}=0.53846.

In order to further analyze the superradiance zero-modes, we display the imaginary parts of some quasinormal frequencies of nodeless modes (n=0n=0) as a function of ζ/π\zeta/\pi with different Lorentz symmetry breaking parameters for a given black hole configuration (the fixed \ell, r+r_{+} and rr_{-}) and scalar field configuration (the fixed μ2\mu^{2} and kk) in Fig. 6. We observe that the threshold of the superradiant instability, i.e.i.e., the critical value of ζ\zeta with Im[ω]=0[\omega]=0, increases as the Lorentz symmetry breaking parameter ss increases. In particular, we find this threshold of the superradiant instability is identically equal to the value of ζ\zeta given by Eq. (25) for the existence of the stationary scalar clouds. From Fig. 6, we note that the critical value of ζ0.75π\zeta\rightarrow 0.75\pi when s=sapps=s_{app}. Thus, we argue that, in the case of μ2=0.65\mu^{2}=-0.65, r+=5r_{+}=5, r=3r_{-}=3, =1\ell=1 and k=1k=1, the threshold of the superradiant instability has a range 0.5π<ζ<0.75π0.5\pi<\zeta<0.75\pi, which suggests that the stationary scalar clouds only exist in this range. As a matter of fact, we only take μ2=0.65\mu^{2}=-0.65 and =1\ell=1 in Figs. 1-4 with different valves of r+r_{+}, and so the superradiant instability can appear in the case of ζ=0.9π\zeta=0.9\pi.

V Conclusions

We have investigated the Klein-Gordon equation for a massive scalar field on a rotating BTZ-like black hole in the Einstein-bumblebee gravity. By employing appropriate Robin boundary conditions at the AdS boundary, we have established the existence of stationary fundamental clouds. We found that, by increasing the Lorentz symmetry breaking parameter ss (the quantum number kk), the clouds exist for smaller (larger) background mass with the same background angular velocity, which means that the parameters ss and kk have entirely opposite effects on the existence lines of scalar clouds. In particular, we obtained the degenerate clouds where different parameters ss and kk have the same existence line of scalar clouds, which has not ever been observed as far as we know. Interestingly, we pointed out that the radial profiles of degenerate clouds ϕ\phi may be used to discriminate degenerate clouds. Moreover, we noticed that, below the limiting value s=(1+1/μ2)s=-(1+1/\mu^{2}), there exist infinite degenerate clouds for any initial values of ss and kk.

By analyzing quasinormal modes in BTZ-like black holes, we observed that, regardless of the Lorentz symmetry breaking parameter ss, superradiant instabilities are only appeared in the nodeless mode (n=0n=0) of type L, which illustrates the observation, shown in Figs. 1, 2 and 3, that only fundamental clouds are allowed. Moreover, we observed that the threshold of the superradiant instability, which determines scalar clouds, depends on the parameter ss. Our results indicate that the Lorentz symmetry breaking parameter ss yields richer physics in the stationary clouds of massive scalar fields around black holes.

Acknowledgements.
This work was supported by the National Key Research and Development Program of China (Grant No. 2020YFC2201400), National Natural Science Foundation of China (Grant Nos. 12275079, 12475050 and 12035005), the Scientific Research Fund of Hunan Provincial Education Department (Grant No. 22A0039) and innovative research group of Hunan Province (2024JJ1006).

Appendix A Massive Klein-Gordon equation in the extremal BTZ-like black hole

In this appendix, we consider the massive Klein-Gordon equation in the extremal BTZ-like black hole of the Einstein-bumblebee gravity. For this case, we have r+=r=M/2r_{+}=r_{-}=\ell\sqrt{M/2} since M2=j2/2M^{2}=j^{2}/\ell^{2}. Thus, the angular velocity of the event horizon is completely determined by the AdS radius, ΩH=1/\Omega_{H}=1/\ell, and the black hole mass can be expressed as M=2r+2ΩH2M=2r_{+}^{2}\Omega_{H}^{2}. We still take the ansatz (9) but change the radial coordinate (11) into

zr+2r2r+2,z\equiv\frac{r^{2}_{+}}{r^{2}-r^{2}_{+}}, (32)

which maps the exterior region from r(r+,+)r\in(r_{+},+\infty) to z(0,+)z\in(0,+\infty) with the AdS boundary at z0z\to 0. Then Eq. (10) turns into

d2ϕ(z)dz2+(Az+B+Cz2)ϕ(z)=0,\frac{d^{2}\phi(z)}{dz^{2}}+\left(\frac{A}{z}+B+\frac{C}{z^{2}}\right)\phi(z)=0, (33)

with

A=2(ω22k2)4r+2(1+s),B=2(ωk)24r+2(1+s),C=μ2(1+s)4.\displaystyle A=\frac{\ell^{2}(\omega^{2}\ell^{2}-k^{2})}{4r_{+}^{2}}(1+s),\;\;B=\frac{\ell^{2}(\omega\ell-k)^{2}}{4r_{+}^{2}}(1+s),\;\;C=-\frac{\mu^{2}(1+s)}{4}. (34)

Thererfore, with the Kummer’s confluent hypergeometric function M(a,b,z)M(a,b,z), we get two linearly independent solutions

ϕ(D)(z)=zβeiαzM(a,b,2iαz),ϕ(N)(z)=z1βeiαzM(ab+1,2b,2iαz),\begin{split}\phi^{(D)}(z)=&z^{\beta}e^{i\alpha z}M(a,b,-2i\alpha z),\\ \phi^{(N)}(z)=&z^{1-\beta}e^{i\alpha z}M(a-b+1,2-b,-2i\alpha z),\end{split} (35)

where β\beta is still as in Eq. (15) and

αωk2r+1+s,aβi1+sω+k4r+,b=2β.\alpha\equiv\ell\frac{\omega\ell-k}{2r_{+}}\sqrt{1+s},\qquad a\equiv\beta-i\ell\sqrt{1+s}\frac{\omega\ell+k}{4r_{+}},\qquad b=2\beta. (36)

In order to obtain stationary scalar modes, we impose the resonance condition ω=kΩH=k/\omega=k\Omega_{H}=k/\ell and arrive at

ϕ(D)(z)=zβ,ϕ(N)(z)=z1β.\phi^{(D)}(z)=z^{\beta},\qquad\phi^{(N)}(z)=z^{1-\beta}. (37)

Obviously, in the extremal case, there is no linear combination of the solutions which is regular at the horizon z+z\to+\infty. Thus, just as in the extremal BTZ case [15], one can not find stationary scalar cloud configurations around the extremal BTZ-like black hole.

References

  • [1] R. Ruffini and J.A. Wheeler, Phys. Today 24, 30 (1971).
  • [2] B. Carter, in Black Holes, Proceedings of 1972 Session of Ecole d’ete de Physique Theorique, edited by C. De Witt and B.S. De Witt (Gordon and Breach, New York, 1973).
  • [3] C.W. Misner, K.S. Thorne, and J.A. Wheeler, Gravitation (Freeman, San Francisco, 1973).
  • [4] S. Hod, Stationary scalar clouds around rotating black holes, Phys. Rev. D 86, 104026 (2012); 86, 129902(E) (2012).
  • [5] S. Hod, Stationary resonances of rapidly-rotating Kerr black holes, Eur. Phys. J. C 73, 2378 (2013); arXiv:1311.5298 [gr-qc].
  • [6] S. Hod, Kerr-Newman black holes with stationary charged scalar clouds, Phys. Rev. D 90, 024051 (2014).
  • [7] S. Hod, Analytic treatment of the system of a Kerr-Newman black hole and a charged massive scalar field, Phys. Rev. D 94, 044036 (2016).
  • [8] C. Herdeiro and E. Radu, Kerr black holes with scalar hair, Phys. Rev. Lett. 112, 221101 (2014); arXiv:1403.2757 [gr-qc].
  • [9] C.L. Benone, L.C.B. Crispino, C. Herdeiro, and E. Radu, Kerr-Newman scalar clouds, Phys. Rev. D 90, 104024 (2014).
  • [10] C. Herdeiro and E. Radu, Construction and physical properties of Kerr black holes with scalar hair, Class. Quant. Grav. 32, 144001 (2015); arXiv:1501.04319 [gr-qc].
  • [11] C. Herdeiro, E. Radu, and H. Rúnarsson, Kerr black holes with self-interacting scalar hair: Hairier but not heavier, Phys. Rev. D 92, 084059 (2015).
  • [12] M. Sampaio, C. Herdeiro, and M. Wang, Marginal scalar and Proca clouds around Reissner-Nordström black holes, Phys. Rev. D 90, 064004 (2014).
  • [13] M. Wang and C. Herdeiro, Maxwell perturbations on Kerr-anti-de Sitter black holes: Quasinormal modes, superradiant instabilities, and vector clouds, Phys. Rev. D 93, 064066 (2016).
  • [14] Y.Q. Wang, Y.X. Liu, and S.W. Wei, Excited Kerr black holes with scalar hair, Phys. Rev. D 99, 064036 (2019); arXiv:1811.08795 [gr-qc].
  • [15] H.R.C. Ferreira and C.A.R. Herdeiro, Stationary scalar clouds around a BTZ black hole, Phys. Lett. B 773, 129 (2017); arXiv:1707.08133 [gr-qc].
  • [16] C. Dappiaggi, H.R.C. Ferreira, and C.A.R. Herdeiro, Superradiance in the BTZ black hole with Robin boundary conditions, Phys. Lett. B 778, 146 (2018); arXiv:1710.08039 [gr-qc].
  • [17] C. Herdeiro, E. Radu, and H. Rúnarsson, Non-linear Q-clouds around Kerr black holes, Phys. Lett. B 739, 302 (2014).
  • [18] H. Okawa, Nonlinear evolutions of bosonic clouds around black holes, Class. Quant. Grav. 32, 214003 (2015).
  • [19] Y. Huang and D.J. Liu, Scalar clouds and the superradiant instability regime of Kerr-Newman black hole, Phys. Rev. D 94, 064030 (2016).
  • [20] S. Hod, Spinning Kerr black holes with stationary massive scalar clouds: the large-coupling regime, J. High Energ. Phys. 01, 030 (2017); arXiv:1612.00014 [hep-th].
  • [21] G. García and M. Salgado, Obstructions towards a generalization of no-hair theorems: Scalar clouds around Kerr black holes, Phys. Rev. D 99, 044036 (2019).
  • [22] S. Hod, Analytic treatment of near-extremal charged black holes supporting non-minimally coupled massless scalar clouds, Eur. Phys. J. C 80, 1150 (2020).
  • [23] S.P. Liu, G.H. Liu, and Y. Peng, Stationary scalar clouds outside charged reflecting compact stars, Mod. Phys. Lett. A 35, 2050175 (2020).
  • [24] G. García and M. Salgado, Existence or absence of superregular boson clouds around extremal Kerr black holes and its connection with number theory, Phys. Rev. D 101, 044040 (2020).
  • [25] G. García and M. Salgado, Regular scalar charged clouds around a Reissner-Nordström black hole and no-hair theorems, Phys. Rev. D 104, 064054 (2021).
  • [26] S. Hod, Nonequatorial scalar clouds supported by maximally spinning Kerr black holes, Phys. Rev. D 108, 124028 (2023).
  • [27] G. García and M. Salgado, Regular scalar clouds around a Kerr-Newman black hole: Subextremal and extremal scenarios, Phys. Rev. D 108, 104012 (2023).
  • [28] S. Hod, Spatially regular charged black holes supporting charged massive scalar clouds, Phys. Rev. D 109, 064074 (2024).
  • [29] G.Z. Guo, P. Wang, T.S. Wu, and H.T. Yang, Stationary Scalar Clouds around Kerr-Newman Black Holes, arXiv:2408.09243 [gr-qc].
  • [30] R. Li, J.K. Zhao, X.H. Wu, and Y.M. Zhang, Scalar clouds in charged stringy black hole-mirror system, Eur. Phys. J. C 75, 142 (2015)
  • [31] C.L. Benone, L.C.B. Crispino, C. Herdeiro, and E. Radu, Acoustic clouds: Standing sound waves around a black hole analogue, Phys. Rev. D 91, 104038 (2015).
  • [32] E. Radu, D.H. Tchrakian, and Y.S. Yang, Non-Abelian clouds around Reissner-Nordström black holes: The existence line, Phys. Rev. D 93, 124069 (2016).
  • [33] C. Bernard, Stationary charged scalar clouds around black holes in string theory, Phys. Rev. D 94, 085007 (2016).
  • [34] Y. Huang, D.J. Liu, X.H. Zhai, and X.Z. Li, Scalar clouds around Kerr-Sen black holes, Class. Quant. Grav. 34, 155002 (2017).
  • [35] G. Tokgoz and I. Sakalli, Stationary scalar clouds around maximally rotating linear dilaton black holes, Class. Quant. Grav. 34, 125007 (2017).
  • [36] N. Grandi and I.S. Landea, Scalar hair around charged black holes in Einstein-Gauss-Bonnet gravity, Phys. Rev. D 97, 044042 (2018).
  • [37] X. Qiao, M. Wang, Q. Pan, and J. Jing, Kerr-MOG black holes with stationary scalar clouds, Eur. Phys. J. C 80, 509 (2020).
  • [38] S. Hod, Stationary scalar clouds supported by rapidly-rotating acoustic black holes in a photon-fluid model, Phys. Rev. D 103, 084003 (2021).
  • [39] M. Ciszak and F. Marino, Acoustic black-hole bombs and scalar clouds in a photon-fluid model, Phys. Rev. D 103, 045004 (2021).
  • [40] Y. Huang and H.S. Zhang, True gravitational atoms: Spherical cloud of dilatonic black holes, Phys. Rev. D 105, 124056 (2022).
  • [41] P.H.C. Siqueira and M. Richartz, Quasinormal modes, quasibound states, scalar clouds, and superradiant instabilities of a Kerr-like black hole, Phys. Rev. D 106, 024046 (2022).
  • [42] S.J. Zhang, Nonlinear instability and scalar clouds of spherical exotic compact objects in scalar-Gauss-Bonnet theory, Eur. Phys. J. C 83, 950 (2023).
  • [43] S. Hod, Charged Gauss-Bonnet black holes supporting non-minimally coupled scalar clouds: analytic treatment in the near-critical regime, Eur. Phys. J. C 83, 214 (2023).
  • [44] V.A. Kostelecky and S. Samuel, Gravitational Phenomenology in Higher Dimensional Theories and Strings, Phys. Rev. D 40, 1886 (1989).
  • [45] Z.F. Mai, R. Xu, D.C. Liang, and L.J. Shao, Extended thermodynamics of the bumblebee black holes, Phys. Rev. D 108, 024004 (2023); arXiv:2304.08030 [gr-qc].
  • [46] X. Zhang, M. Wang, and J. Jing, Quasinormal modes and late time tails of perturbation fields on a Schwarzschild-like black hole with a global monopole in the Einstein-bumblebee theory, Sci. China Phys. Mech. Astron. 66, 100411 (2023); arXiv:2307.10856 [gr-qc].
  • [47] R. Casana, A. Cavalcante, F.P. Poulis, and E.B. Santos, Exact Schwarzschild-like solution in a bumblebee gravity model, Phys. Rev. D 97, 104001 (2018); arXiv:1711.02273 [gr-qc].
  • [48] C. Ding, Y. Shi, J. Chen, Y. Zhou, C. Liu, and Y. Xiao, Rotating BTZ-like black hole and central charges in Einstein-bumblebee gravity, Eur. Phys. J. C 83, 573 (2023); arXiv:2302.01580 [gr-qc].
  • [49] C. Chen, Q. Pan, and J. Jing, Quasinormal modes of a scalar perturbation around a rotating BTZ-like black hole in Einstein-bumblebee gravity, Phys. Lett. B 846, 138186 (2023), arXiv:2302.05861 [gr-qc].
  • [50] F. Ge, Q. Pan, S. Chen, and J. Jing, Mass ladder operators and quasinormal modes of the static BTZ-like black hole in the Einstein-bumblebee gravity, accepted by Chin. Phys. C.
  • [51] V.A. Kostelecky and S. Samuel, Spontaneous breaking of Lorentz symmetry in string theory, Phys. Rev. D 39, 683 (1989).
  • [52] R. Bluhm and V.A. Kostelecky, Spontaneous Lorentz violation, Nambu-Goldstone modes, and gravity, Phys. Rev. D 71, 065008 (2005); arXiv:hep-th/0412320.
  • [53] O. Bertolami and J. Paramos, The Flight of the bumblebee: Vacuum solutions of a gravity model with vector-induced spontaneous Lorentz symmetry breaking, Phys. Rev. D 72, 044001 (2005); arXiv:hep-th/0504215.
  • [54] Q.G. Bailey and V.A. Kostelecky, Signals for Lorentz violation in post-Newtonian gravity, Phys. Rev. D 74, 045001 (2006); arXiv:gr-qc/0603030.
  • [55] R. Bluhm, N.L. Gagne, R. Potting, and A. Vrublevskis, Constraints and Stability in Vector Theories with Spontaneous Lorentz Violation, Phys. Rev. D 77, 125007 (2008); Erratum ibid. 79 029902 (2009); arXiv:0802.4071 [hep-th].
  • [56] V.A. Kostelecky and J. Tasson, Prospects for Large Relativity Violations in Matter-Gravity Couplings, Phys. Rev. Lett. 102, 010402 (2009); arXiv:0810.1459 [gr-qc].
  • [57] M.D. Seifert, Generalized bumblebee models and Lorentz-violating electrodynamics, Phys. Rev. D 81, 065010 (2010); arXiv:0909.3118 [hep-ph].
  • [58] R.V. Maluf, C.A.S. Almeida, R. Casana, and M. Ferreira, Einstein-Hilbert graviton modes modified by the Lorentz-violating bumblebee Field, Phys. Rev. D 90, 025007 (2014); arXiv:1402.3554 [hep-th].
  • [59] J. Páramos and G. Guiomar, Astrophysical Constraints on the Bumblebee Model, Phys. Rev. D 90, 082002 (2014); arXiv:1409.2022 [astro-ph].
  • [60] C.A. Escobar and A. Martín-Ruiz, Equivalence between bumblebee models and electrodynamics in a nonlinear gauge, Phys. Rev. D 95, 095006 (2017); arXiv:1703.01171 [hep-th].
  • [61] J.F. Assunão, T. Mariz, J.R. Nascimento, and A.Y. Petrov, Dynamical Lorentz symmetry breaking in a tensor bumblebee model, Phys. Rev. D 100, 085009 (2019); arXiv:1902.10592 [hep-th].
  • [62] D. Capelo and J. Páramos, Cosmological implications of Bumblebee vector models, Phys. Rev. D 91, 104007 (2015); arXiv:1501.07685 [gr-qc].
  • [63] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Academic, New York, 1996).
  • [64] M. Wang, Quantum and classical aspects of scalar and vector fields around black holes, PhD thesis, Aveiro U. (2016); arXiv:1606.00811 [gr-qc].
  • [65] M. Wang, Boundary conditions for Maxwell fields in Kerr-AdS spacetimes, Int. J. Mod. Phys. D 25, 1641011 (2016).
  • [66] M. Wang, C. Herdeiro and M. O. P. Sampaio, Maxwell perturbations on asymptotically anti-de Sitter spacetimes: Generic boundary conditions and a new branch of quasinormal modes, Phys. Rev. D 92, 124006 (2015), arXiv:1510.04713 [gr-qc].
  • [67] M. Wang, Z. Chen, X. Tong, Q.Y. Pan, and J.L. Jing, Bifurcation of the Maxwell quasinormal spectrum on asymptotically anti-de Sitter black holes, Phys. Rev. D 103, 064079 (2021); arXiv:2104.04970 [gr-qc].
  • [68] M. Wang, Z. Chen, Q.Y. Pan, and J.L. Jing, Maxwell quasinormal modes on a global monopole Schwarzschild-anti-de Sitter black hole with Robin boundary conditions, Eur. Phys. J. C 81, 469 (2021); arXiv:2105.10951 [gr-qc].
  • [69] R. Brito, V. Cardoso, and P. Pani, Superradiance: New Frontiers in Black Hole Physics, Lect. Notes Phys. 906, 1 (2015); Lect. Notes Phys. 971, 1 (2020); arXiv:1501.06570 [gr-qc].