Scalar field perturbation around a rotating hairy black hole: quasinormal modes, quasibound states and superradiant instability

We consider the quasinormal modes, quasibound states and superradiant instability of a rotating hairy black hole, which possesses a Horndeski hair as deviation from Kerr black hole, under the perturbation of massive scalar field. With the use of the matrix method, we mainly calculate the eigenfrequencies related to those modes of the perturbation. Under the perturbation of the massless scalar field, the Horndeski hair and spin parameter have significant influences on the quasinormal frequency, but its imaginary part is always finite negative and no unstable mode is found. Under the perturbation of the massive scalar field, we focus on the eigenfrequencies of quasibound states and find the modes of which the imaginary part of eigenfrequencies is positive, indicating that the black hole undergoes superradiant instability. Then we scan the parameters and figure out a diagram in the space of Horndeski hair and spin parameters to distinguish the rotating hairy black hole with superradiant instability from the stable one.

The stability of black hole solution under certain perturbations is an over-lasting crucial problem in black hole physics.Quasinormal modes (QNMs) with specific eigenfrequencies (Quasinormal frequencies, or QNFs) play a remarkable role in this scenario.QNFs of a field perturbation on the black hole are infinite discrete spectrum of complex frequencies, of which the real part determines the oscillation timescale of the QNMs, while the imaginary part determines their exponential decaying timescale.It is one of the most important characteristics of black hole geometry, and the interest of QNMs in more fundamental physics can be referred into the reviews [1][2][3].Early in the 1950s, Regge and Wheeler studied the stability of black holes in general relativity (GR) by analyzing QNFs of gravitational perturbation in Schwarzschild black hole, and they found that the perturbation field oscillates and decays over time such that the Schwarzschild black hole is stable [4].Later, various external field perturbations in the Schwarzschild black hole were explored and further verified the stability of the Schwarzschild black hole in GR [5][6][7][8] .
However, rotating neutral black holes(i.e.Kerr black holes) in GR may develop instabilities under certain conditions.In fact, if an ingoing bosonic wave is scattered off a Kerr black hole with the horizon angular velocity Ω H , the scattered wave will be amplified once the superradiance condition 0 < ω < mΩ H is satisfied, where ω is the frequency of the wave and m is the azimuthal number [9].The extra energy of the scattered wave is obtained from the rotational energy of the black hole, so this is a wave analogue of the Penrose process [10,11].In this process, if there exists a hypothetic reflective "mirror" that makes the amplified wave be scattered back and forth between the "mirror" and the black hole, the amplifying proceeds continuously and will finally cause the black hole background to suffer from superradiant instability [12].Such a system was also dubbed the "black hole bomb" [13].The Kerr black hole was found to be stable under massless scalar, electromagnetic or gravitational perturbations according to the analysis of those QNMs [14,15].However, once the perturbing scalar field is massive, Kerr black hole was found to undergo superradiant instability because the mass term of the perturbing scalar field acts as the natural mirror [16][17][18][19][20][21][22][23].Moreover, it was addressed in [12,20] that such superradiatively unstable mode may be induced from the quasi-bound state (QBS) of the massive scalar field, which decays exponentially far away from the black hole with certain eigenfrequency.
In GR, the eigenfrequencies for QNM and QBS around a Kerr black hole are only characterized by three parameters: the mass and angular momentum of the black hole and the mass of external field.Though recent observations of gravitational wave generated from binary compact objects [24,25] and black hole shadow [26,27] provide us with remarkable chances to test GR in the strong field regime, the uncertainties in the data leave some space for alternative theories of gravity and non-Kerr metric of black holes.In addition, despite the great success, there are still some challenging problems for GR, such as the explanation of the Universe expansion history, the large scale structure and the understanding of quantum gravity [28][29][30].Therefore, a more general theory of gravity is required from both observational and theoretical aspects.Plenty of modified gravitational theories have been proposed [31][32][33], among which the usual way is to modify the action of GR.The scalar tensor theories, which contains a scalar field in addition to the metric field are considered as the simplest nontrivial extensions of GR.The most general scalar tensor theory with second order equations of motion for both the metric and the scalar field in four dimensions is the Horndeski theory, which is free of the Ostrogradski instabilities [34,35] and is found to be equivalent to the generalized Galileon theory [36,37] derived in a different way.The observational constraints or bounds on Horndeski gravity have been extensively explored in [38][39][40][41][42][43].Moreover, Horndeski gravity attracts lots of attention in both cosmological and astrophysical communities because it has significant consequences in describing the accelerated expansion and other interesting features (please see [35] for a review), as well as holographic applications [44][45][46][47][48][49][50][51][52] and references therein.
In particular, the Horndeski framework has also been used to test the no hair theorem of black holes, which states that the isolated black hole in GR is merely characterized by three parameters: mass, electric charge and angular momentum.Due to the extra scalar field in the action, it is natural to figure out if there are black hole solutions with scalar hair, and this makes sense because similar to GR, the Horndeski theory has diffeomorphism invariance and second order field equations.The no hair theorem for the shift symmetric Horndeski theory [53] has proved that a static and spherically symmetric black hole cannot sustain a nontrivial scalar field with vanishing boundary condition at infinity.But later it turns out that there is a loophole in the demonstration and the counter hairy black hole solution was constructed in [54,55].The no hair theorem has also been considered in [56] for the shift symmetric Horndeski and beyond Horndeski theories, which also admit the hairy static and spherically symmetric black hole solutions with nontrivial scalar profile.Up to date, many hairy black holes have been constructed and analyzed in Horndeski theory, including the radially dependent [44,57,58] and time-dependent scalar hair [59][60][61], among which the black holes with linear time-dependent scalar hair is unstable under perturbation [62], hence ruling out the possibility of this type of hairy solutions.Fortunately, a spherically symmetric hairy black hole solution in shift symmetric Horndeski theory has recently been constructed in [63], which was verified to be stable under various massless external field perturbations by us [64].Some observational investigations have been carried out for this static hairy Horndeski black hole, for example, the strong gravitational lensing [65] and photon rings in the black hole image [66], which shew that the Horndeski hair has a significant influence on the related observable.Especially, with the use of the revised Newman-Janis algorithm [67,68], the authors of [69] soon constructed its axially symmetric counterpart, known as rotating black hole with Horndeski hair.The thermodynamic properties and weak gravitational lensing [69], black hole shadow [70] and superradiant energy extraction [71] have been extensively studied in this rotating hairy black hole.
In this paper, we will study the QNMs and QBSs for a massive perturbing scalar field on the rotating hairy black hole in Horndeski theory.Our study motivates from two aspects.One is that the computation of eigenfrequencies of QNMs for scalar field perturbation is a first step to analyze the dynamical stability of the rotating hairy black hole.The other is that the interesting QBSs could be the potential source of the superradiatively unstable modes as we aforementioned.It is noticed that a brief investigation of superradiant instabilities in the rotating hairy metric was performed in [71] using approximate analytical methods.Our work expands the results of [71] by analyzing in detail the eigenfrequency of these unstable modes as a function of the black hole parameters.We will explore how the hair parameter influences this instability apart from the rotation.To this end, we will use the matrix method to numerically determine the frequency spectrums of QNMs, QBSs and superradiatively unstable modes.
This paper is organized as follows.In Sec.II we present the background geometry of the rotating hairy Horndeski black hole, derive the master equations of the scalar field perturbation and analyze the corresponding boundary conditions.The eigenfrequencies of QNMs, QBSs and superradiant instabilities of the perturbation are explored in Sec.III, in which we mainly analyze the effect of the hair parameter and figure out a 'phase' diagram in black hole parameter space.Sec.IV contains our conclusions.

II. BACKGROUND GEOMETRY AND SETUP OF THE SCALAR FIELD PERTURBATION
The action for the dubbed quartic Horndeski gravity reads as [56] where ) are arbitrary functions of χ and Q i,χ ≡ ∂Q i /∂χ, R is the Ricci scalar and G µν is the Einstein tensor.A static hairy black hole in a specific quartic Horndeski theory, saying that Q 5 in the above action vanishes, has been constructed in [63] with the metric Here, M and h are the parameters related to the black hole mass and scalar hair.The metric reduces to Schwarzschild case as h → 0, and it is asymptotically flat.The hairy black hole described by the above metric was found to be stable under various external field perturbations by computing their QNFs and time evolutions [64].Besides, some observable phenomena affected by the Horndeski hair have also been explored in [65,66].Very recently, the axial symmetric counterpart of metric (2) was constructed, and in the Boyer-Lindquist coordinates the metric of the corresponding rotating hairy black hole reads as [69] with the functions where a is the spin parameter.The axially symmetric spacetime is also asymptotically flat and reduces to the Kerr spacetime as h → 0. The horizons of the rotating hairy black hole are determined by the real positive roots of ∆ = 0, which only exist in a certain parameter region of (a/M, h/M ) shown in FIG. 1.In detail, in the gray parameter region, ∆ = 0 has two real positive roots r ± , in which the larger root r + indicates the event horizon while the smaller one r − denotes the Cauchy horizon.In the white parameter region, ∆ = 0 has no real root meaning there is no horizon such that the metric (3) describes a naked singularity.And the parameters on the black curve that separates these two regions make the rotating hairy black hole become extremal with r − = r + .The horizons r ± as functions of the black hole parameters are explicitly shown in FIG.Here in order to analyze the dynamical (in)stability of the rotating hairy black hole (3), we, as the first step, will introduce the massive scalar field as the probe and study the eigenfrequencies of its master equation.The dynamic of a perturbing scalar field around the rotating hairy black hole is governed by the Klein-Gordon equation where µ is the mass of the scalar field.Then with the decomposition in which ω is the frequency and m is the azimuthal number, we can separate Eq.( 5), yielding the radial master equation and the angular master equation where u = cos θ and λ is the separation constant.Then, we shall analyze the physically relevant boundary conditions of Eq.( 7) and Eq.( 8).The boundary conditions of S(u) can be directly obtained from its asymptotic solutions near the boundaries u → ±1, which are For the radial equation, by introducing the tortoise coordinate and the function transformation we can obtain the equation with It is straightforward to get the asymptotic behavior of V (r) as with Ω H = a r 2 + +a 2 the horizon angular velocity, based on which we can easily obtain the asymptotic solutions as the boundary conditions with ω = ω − mΩ H and k = √ ω 2 − µ 2 .Reminding that the perturbation cannot come out from the event horizon, we only consider the case with the minus sign in the first formula of Eq.( 15), i.e., the ingoing boundary condition near the event horizon.For the second formula of Eq.( 15), we take the plus (minus) sign and ensure Re(k) > 0 (Im(k) < 0) to get the purely outgoing (exponentially decaying) boundary condition at the spatial infinity for the QNMs (QBSs) [20].The eigenfrequencies ω obtained from the master equations are discrete complex numbers, of which the real part describes the oscillation while the imaginary part can be employed to judge the (in)stability of the black hole under the perturbations.QNMs are known to have negative Im(ω), which means that the perturbation decays and the lifetime is evaluated by τ = |Im(ω)| −1 .But the QBSs may become the source of the unstable modes with positive Im(ω), implying an exponentially growing amplitude due to the superradiant instability.Therefore, we shall solve the master equations to obtain the eigenfrequencies, and analyze the possible modes of the perturbation in the rotating hairy black hole.Since it's not possible to solve the master equations analytically, so we shall solve them numerically with the use of the matrix method [72][73][74] and extract the eigenfrequencies of the QNMs and QBSs.For the convenience to employ this method, we shall make some transformations to Eq.( 7) and Eq.( 8) beforehand.Firstly, by rewriting the boundary conditions of the radial equation Eq.( 15) into the form of function R via Eq.( 11), and expanding them at the corresponding asymptotic region, we can obtain where σ = r+(r++h)−a 2 ω, and the upper (lower) signs correspond to the plus (minus) sign in the second formula of Eq. (15).Then, in order to factor out the asymptotic behaviors of the radial and angular functions, we redefine so that R(r) and S(u) will automatically satisfy the aforementioned boundary conditions.Next, we introduce the new coordinates to bring the integration domains into x ∈ [0, 1] and v ∈ [0, 1] to replace r ∈ [r + , ∞) and u ∈ [−1, 1], respectively.Finally, by further performing the transformation we obtain the equations of χ R and χ S , which we can directly solve by using the matrix method, as where the coefficient functions C Ri and C Si (i = 0 ∼ 2) are functions of ω and λ.Moreover, the corresponding boundary conditions are then converted into the homogeneous form The last step is to follow the matrix method algorithm to extract the eigenfrequencies of QNMs and QBSs by numerically solving the algebraic equations where the expressions of the matrix M R or S (ω, λ) are determined by the coefficient functions in the equations ( 20) and ( 21).Here we prefer not to repeat the steps of the matrix method algorithm, and readers can refer to [72][73][74] for more details.

III. RESULTS OF THE PERTURBING MODES AND THEIR EIGENFREQUENCIES
In this section, we shall analyze the results of QNMs, QBSs and superradiant instabilities of the perturbing scalar field on the rotating hairy black hole background (3).It is noted that all the physical quantities in this section are rescaled by the mass parameter M into be dimensionless.

A. Quasinormal modes
We first consider the dynamic behavior of a massless (µM = 0) scalar field perturbing around the rotating hairy Horndeski black hole.We mainly focus on the effects of spin parameter a/M and the Horndeski hair parameter h/M on the fundamental QNMs with overtone number N = 0, so we select samples of the two parameters to show the dependence of the fundamental QNFs on them.The fundamental QNFs with the mode ℓ = m = 0 are shown in FIG. 3.With fixed a/M (see the dashed curves), as the hair parameter h/M decreases, the real part of the QNFs decreases while the imaginary part increases.This means that comparing to the Kerr black hole, the Horndeski hair makes the massless scalar field perturbation oscillate slower and decay slower.With fixed h/M (see the solid curves), the effect of the spin parameter on QNFs could depend on the value of h/M .In detail, for weaker Horndeski hair, i.e., small | h/M |, as a/M increases, the real part of QNFs first increases and then decreases which is similar to that found in Kerr black hole [20], but for large enough | h/M | it decreases monotonically.Meanwhile, the imaginary part first increases with the increasing of a/M and then decreases slightly as the black hole approaches the extreme for all the selected values of h/M , but it tends to a finite negative value for the extremal case.This implies that for the rotating hairy black hole with small | h/M |, the faster spin makes the corresponding perturbation oscillate first faster and then slower as the black hole becomes near extremal, similar to that occurs in Kerr black hole.However, in the rotating hairy black hole with large enough | h/M |, the oscillation always becomes slower as a/M increases, this is because the strong suppression effect of h/M on the real part could balance the enhancement effect of a/M and even dominates when the black hole closes to extreme.(ℓ, m, N ) = (0, 0, 0) FIG. 3: The variation of the real and imaginary parts of eigenfrequencies for fundamental ℓ = m = 0 QNM of the massless scalar field on the rotating hairy Horndeski black hole in terms of the hair (spin) parameter h/M (a/M ) with several selected spin (hair) parameters.For each curve with given a/M (h/M ), the value of h/M (a/M ) starts from 0 (Kerr case) [(static case)] and then extends toward the negative (positive) direction until the black hole becomes near extreme.
To compare the influences of different angular quantum number ℓ and azimuthal number m on QNMs, we compute the fundamental QNFs for ℓ = 1.The results are depicted in FIG. 4 the same selected values of parameters as those in ℓ = 0 case.We can extract the following features.(i) As h/M decreases, the real (imaginary) part of the fundamental QNFs decreases (increases) with any selected a/M for all the ℓ = 1 modes, except for the mode ℓ = m = 1 with large a/M of which the real part first decreases then increases with the decreasing of h/M .(ii) The effect of varying a/M can also be read off.For the modes m = 0 and m = 1, both the real and the imaginary parts of the QNFs increase as the spin becomes faster for the selected h/M .For the mode m = −1, as the rotation increases, the real part decreases monotonically for any selected h/M , but the imaginary part increases for small | h/M | and decreases for large | h/M | on the whole.(iii) According to our study, all the imaginary part of eigenfrequencies are finite negative value and we do not find any unstable mode.

B. Quasibound states and superradiant instabilities
We then turn on the mass of the scalar field and study the existence of QBSs and their eigenfrequencies.As we mentioned in Sec.I, the mass term may trigger the superradiant instability that we will show in this subsection.To be specific, we fix the mass of the scalar field µM = 0.4 as it was done in [75], and compute the fundamental eigenfrequencies of QBSs for the modes ℓ = 0 and ℓ = 1 for the same selected parameters as those in the QNMs.The results are depicted in FIG. 5 and FIG.6, from which we can find some common features of the modes.(i) The real part of the frequencies exhibits a typical feature of the QBSs, i.e., they are always smaller than the the mass of the scalar field.(ii) Similar to the QNFs, the strength of the Horndeski hair also influences the effect of the spin parameter on the QBS frequency, which could be different for different modes.Here we will pay more attention to the effect of Horndeski hair on the eigenfrequencies rather than starchily setting out the effect of rotation.
In FIG. 5 with ℓ = 0, the imaginary part of the eigenfrequency is always negative, indicating that the mode is stable.With a fixed spin parameter, the stronger Horndeski hair suppresses both the real and imaginary parts of the eigenfrequency, which is different from the one that happens in QNFs.Moreover, it is obvious that by comparing FIG. 5 and FIG. 3, the real part of the QBSs of the massive scalar field with µM = 0.4 is larger than those of the QNMs of the massless scalar field, indicating faster oscillation for the former.In FIG. 6 with ℓ = 1, with the selected a/M , the real and the imaginary parts of QBSs from all modes always decrease with the decreasing of h/M , however, the inset of the right panel shows the exception for the mode ℓ = m = 1 with a/M = 0.8, where the imaginary part first decreases slightly then increases and eventually becomes positive.These QBSs with positive Im(ω) means that the amplitude of these states will grow exponentially over time, leading to the superradiant instabilities.Moreover, note that the superradiant instabilities only exist in the m > 0 states, as the result of the superradiance condition 0 < ω < mΩ H .
To clearly show how the parameters of hairy black hole affect the superradiant instabilities, we further plot the eigenfrequencies as a function of the Horndeski hair for several large spin parameters in FIG. 7. In the left plot, the dashed line with Im(M ω) = 0 separates the decaying modes with Im(M ω) < 0 from the superradiatively unstable modes with Im(M ω) > 0. In fact, the modes with Im(M ω) = 0 are bound states which could form as the scalar cloud around the rotating hairy black hole.As addressed in [76], the eigenfrequency of the scalar cloud is purely real ω sc = mΩ H , which means that the scalar cloud is synchronized with the rotation of the hairy black hole such that there is no energy flux towards the black hole [76][77][78].In the right plot, our numerical real part matches well with ω sc = mΩ H marked by the red dots for the scalar cloud.In a word, the Horndeski hair significantly affects the superradiant instability and will make the onset of superradiant instability need a smaller spin parameter.We collect the data for the scalar cloud and superradiatively unstable modes into the (a/M, h/M ) space, depicted in FIG. 8 where we rescale the axes, and the red curve denotes the onset of instability or the black hole with a scalar cloud.The rotating hairy Horndeski black hole with the parameters in the light gray region is stable under the perturbation of the scalar field with mass µM = 0.4, while the black hole with the parameters in the dark gray region could undergo superradiant instability.It is shown that comparing to the Kerr black hole (h = 0), the superradiant instability can occur for slower rotating hairy Horndeski black hole but the parameter region of the existence of superradiant instability is narrower.
Before closing this section, we shall study the effect of the scalar field mass µM on the QBSs and the superradiant instabilities.As addressed in [20], the superradiant instability of Kerr black hole occurs when the mass is in the range of 0 < µM < µ c M , where µ c M depends on the spin of black hole and we name it as the upper bound of the scalar mass corresponding to Im(M ω) = 0.Here we fix the spin parameter a/M = 0.86 and calculate the eigenfrequencies of The fundamental eigenfrequency of the QBSs in the mode ℓ = m = 1 of massive (µM = 0.4) scalar field as a function of h/M with several selected a/M .In the plots, each curve for fixed a/M extends toward the negative h/M direction until the black hole becomes near extreme.In the left plot, there exists positive values of the imaginary part, indicating the superradiant instabilities.In the right plot, the solid curves denote the numerical results of the real part.The red dots denote the eigenfrequencies of the bound state or scalar clouds for which the imaginary part vanishes, so the parameters are extracted from the intersections between each curve and the dashed line in the left plot.The real part of these bound states marked by red dots are directly evaluated by ωsc = mΩ H . FIG. 9: The fundamental eigenfrequency of the QBSs in the mode ℓ = m = 1 as a function of h/M for selected masses of the scalar field.Here we fix a/M = 0.86.In the left plot for imaginary part, we see that the onset of superradiant instability is very sensitive to the scalar mass.Again, the red dots in the right plot denote the eigenfrequencies of bound state or scalar cloud, which we evaluate via ωsc = mΩ H .

IV. CONCLUSIONS
Neutral black holes in general relativity are known to be characterized by two hairs: the mass and spinning parameter.However, theories of gravity beyond GR usually allow additional hairs, which in general modify various properties of the spacetime.In this paper, with the use of matrix method, we investigated the eigenfrequencies associated with the perturbation of a massive scalar field around a rotating hairy black hole in Horndeski theory.
Under the perturbation of massless scalar field, we found that comparing to the Kerr black hole with the same fixed spin parameter, the Horndeski hair could make the perturbation mode ℓ = m = 0 oscillate slower and decay slower, corresponding to smaller real part and larger imaginary part of the QNFs.The effect of the spin parameter on the QNFs could depend on the fixed value of h/M .For small | h/M |, the real part of QNFs first increases and then decreases with the increase of a/M , which is similar to that found in Kerr black hole.While for large enough | h/M |, it decreases monotonically, which is understandable because the strong suppression effect of | h/M | on the real part could balance the enhancement effect of a/M and even dominates when the black hole closes to the extreme.Meanwhile, for all the selected values of h/M , the imaginary part first increases with the increasing of a/M but then decreases slightly as the black hole closes to the extreme.In addition, QNFs with ℓ = 1 and m = −1, 0, 1 were also calculated.The Horndeski hair also has a significant influence on those QNFs, though the rules may be different from those for the mode ℓ = m = 0. Nevertheless, all the QNFs for various modes have finite negative imaginary part, which in some sense suggests that the rotating hairy Horndeski black hole is stable under the perturbation of a massless scalar field.
Then we turned on the mass of the scalar field perturbation, and mainly checked the eigenfrequencies of the quasibound state.For the scalar mass µM = 0.4, we found significant effects of the Horndeski hair and spin parameter on the eigenfrequencies of the quasibound states with ℓ = 0 and ℓ = 1.Especially, the stronger Horndeski hair always suppresses both the real and imaginary parts of the eigenfrequency for all the modes, except for the mode ℓ = m = 1 with a/M = 0.8 of which the imaginary part first decreases slightly and then increases to be positive.This could be associated with the superradiant instability of the black hole under the perturbed mode ℓ = m = 1.Afterward, we scanned the parameter space to fix the onset of superradiant instability, also known as bound state or scalar cloud, which satisfies Im(M ω) = 0 and Re(M ω) = mΩ H .Our results indicate that due to the existence of the Horndeski hair, the rotating hairy black hole with a smaller spin than Kerr black hole could undergo superradiant instability.Finally, we discussed the effect of mass of the scalar field on the quasibound state and superradiant instability.Our study in the selected parameter region, shew that the upper bound of the scalar field mass for the superradiant instability increases when the value h/M is smaller.This implies that comparing to the Kerr black hole, the range of the scalar field mass for superradiant instability of the rotating hairy Horndeski black hole may be extended.
To conclude, this study provided the preliminary step to analyze the (in)stability of the rotating hairy black hole in Horndeski theory, which we believe deserves to extend into the perturbations of fields with higher spin.Our findings on superradiant instability expand the approximate analytical results of [71].Moreover, our results rich the phenomenal properties introduced by the Horndeski hair, which we hope to give some insight in the test of the no-hair theorem of black holes as well as the observational applications of hairy black holes.
geometry and setup of the scalar field perturbation III.Results of the perturbing modes and their eigenfrequencies A. Quasinormal modes B. Quasibound states and superradiant instabilities IV. Conclusions I. INTRODUCTION

2 .FIG. 1 :
FIG.1:The (a/M, h/M ) parameter regions in which the metric (3) describes a black hole (gray region) or a naked singularity (white region), the parameters on the black curve that separates these two regions correspond to a extreme black hole.

6 FIG. 5 :FIG. 6 :
FIG.5:The fundamental eigenfrequencies of the quasibound states in the mode ℓ = m = 0 of perturbing scalar field with the mass µM = 0.4.The setting of black hole parameters are the same as those in QNMs.
FIG. 7:The fundamental eigenfrequency of the QBSs in the mode ℓ = m = 1 of massive (µM = 0.4) scalar field as a function of h/M with several selected a/M .In the plots, each curve for fixed a/M extends toward the negative h/M direction until the black hole becomes near extreme.In the left plot, there exists positive values of the imaginary part, indicating the superradiant instabilities.In the right plot, the solid curves denote the numerical results of the real part.The red dots denote the eigenfrequencies of the bound state or scalar clouds for which the imaginary part vanishes, so the parameters are extracted from the intersections between each curve and the dashed line in the left plot.The real part of these bound states marked by red dots are directly evaluated by ωsc = mΩ H .

FIG. 8 :
FIG.8:The parameter regions for the stable rotating hairy Horndeski black hole and the one with superradiant instability under the perturbation of massive (µM = 0.4) scalar field, of the fundamental ℓ = m = 1 QBSs.The black curve corresponds to the extreme black hole as in FIG.1, and the red curve corresponds to the onset of superradiant instability, or bound state or scalar cloud with the eigenfrequency ωsc = mΩ H which is purely real.