Quasinormal modes of a charged spherical black hole with scalar hair for scalar and Dirac perturbations

The quasinormal modes of charged and uncharged massive scalar fields and also of charged Dirac fields against the background of a charged spherical black hole endowed with a scalar hair have been investigated. Special emphasis has been given to the case where negative scalar charge dominates over the electric charge of the black hole which mimics an Einstein-Rosen bridge. Except for the complete monotonic behaviour of the damping (imaginary part of the quasinormal mode) against the charge of the black hole as opposed to the existence of a peak for the pure RN case, the qualitative behaviour does not appreciably change due to the presence of scalar hair.


Introduction
After the dream of detecting gravity waves came true, and that too from a merger of two black holes [1], the importance of a thorough investigation of the quasinormal modes in connection with the black hole perturbations cannot perhaps be exaggerated. These investigations started a long way back, through the work of Regge and Wheeler [2] and Vishveshwara [3,4]. The response of a black hole to a perturbation of an external field or the perturbation of the metric is manifested in the form of a damped wave emitted by the black hole, characterized by a complex frequency, called the quasinormal frequency. The real part of the frequency corresponds to the actual frequency of the wave motion while the imaginary part takes care of the damping factor. For excellent reviews, we refer to the works of Nollert [5], Kokkotas and Schmidt [6] and Konoplya and Zhidenko [7]. a e-mail: ac13ip001@iiserkol.ac.in b e-mail: narayan@iiserkol.ac.in Quasinormal modes (QNM) for a Schwarzchild black hole has been studied by Vishveshwara [4] and also by Davis, Ruffini, Press and Price [8]. QNMs for a Reissner-Nordström black hole was first investigated by Gunter [9]. Investigations regarding QNMs for various kind of black holes are already there in the literature. Dreyer discussed the QNMs, area spectrum and entropy of a black hole and also fixed the value of the Immirizi parameter which arises in Loop quantum gravity [10]. Cardoso and Lemos discussed the QNMs of a BTZ black hole [11] and also Schwarzchild-AdS black holes [12]. The latter had been discussed by Horowitz and Hubeny [13] also. QNMs for a near extremal black hole has been investgated by Starinets [14] and by Cardoso and Lemos [15]. QNMs for a Gauss-Bonnet black hole has been discussed by Chakrabarti [16].
The purpose of the present work is to investigate the QNMs of a black hole endowed with a scalar hair. We pick up the example which has been very recently given by Astorino [17]. The black hole has both electric charge and scalar charge. The scalar field part is basically the same one as that given by Bekenstein [18,19], but in a much more useful form. The metric is qualitatively same as the Reissner-Nordström (RN) metric. The "scalar charge" comes only as an additive correction to the electric charge, so nothing new comes out of it as such to start with. But the scalar charge comes with a power unity in the metric as opposed to the quadratic appearence of the electric charge. Thus one can set the scalar charge s < −e 2 so as to get the 1 r 2 term with a negative coefficient. This is a very simple realization of a "mutated Reissner-Nordtsröm" metric leading to the Einstein-Rosen bridge [20] or a so called wormhole. The present work deals with such a metric, with primarily negative values of s. The perturbation of massless and massive uncharged/charged scalar particles and massless charged Dirac particles and the QNMs generated by the perturbations in the vicinity of a mutated Reissner-Nordström black hole are dis-cussed in the present work. The continued fraction method (see Refs. [21][22][23]) has been adopted. The fundamental mode is the dominating one in the signal and only that mode is dealt with.
In almost all the cases both the frequency and the damping rate decrease with the magnitude of the negative scalar charge. For massive scalar field, the damping rate falls off sharply compared to the massless case whereas the real frequency falls off at a much slower rate. For charged fields, the oscillation frequency and the damping rate is more for higher values of the field charge.
The paper is organized in the following way. We start with a brief description of the background spacetime in section 2. In section 3, we briefly describe the continued fraction technique and discuss the QNMs for both uncharged and charged scalar fields close to a mutated RN black hole. Section 4 includes a discussion on the the QNMs of massless charged Dirac field around an RN black hole endowed with a scalar hair. The fifth and final section contains a summary and discussion on the results obtained. As it is already mentioned, the work is done using the continued fraction method. It has also been worked out using a 3rd order WKB method (see Refs. [24,25]), but not mentioned in the text. To facilitate a comparison, we include a table showing the results of the two methods for one example as an appendix.

Background Spacetime
Starting from the action of general relativity coupled to a Maxwell field F µν and conformally coupled to a scalar field ψ, Astorino [17] arrived at the Reissner-Nordström black hole of mass M and charge e endowed with a scalar hair s, where The net stress-energy tensor looks like It is interesting to note that the scalar field ψ has a constant value. Still it gives a non-trivial contribution to the metric because of its nonminimal coupling with geometry in the form R 6 ψ 2 in the action (1). The scalar hair s is a primary hair since the scalar field ψ survives even in the absence of the electromagnetic field. It is easy to note from relation (5) that the trace of the energy momentum tensor due to the scalar field alone is also zero. Thus the existence of this hair is completely consistent with the theorem given in Ref. [26]. In the range 0 > s > −e 2 , the scalar field has an imaginary value. In such a case, the kinetic part in the action should have been written as ∇ µ ψ * ∇ µ ψ. However, it hardly matters in the present case as the kinetic part becomes trivial as ψ is a constant. As discussed earlier, one of the principal motivation of the work is to look at the QNMs for a mutated Reissner-Nordström black hole which requires s < −e 2 , the question of a complex scalar field will not arise.
The present black hole given by the solution (2) and (3), henceforth referred to as the Reissner-Nordström-scalar hair (shRN) black hole, is also characterised by an inner Cauchy horizon (r − ) and an outer event horizon (r + ). The horizons of the shRN black hole are located at The maximum value of the scalar (or electric) charge is determined by the extremality condition, For s = −e 2 the shRN spacetime reduces to a Schwarzschild black hole with the event horizon at r + = 2M. The mutated RN spacetime is also characterised by a single event horizon as r − is negative and of no physical significance.
3 Massive scalar field around a charged black hole with scalar hair In this section we discuss the dynamics of a massive charged scalar field in the background of an shRN black hole and study the fundamental (n = 0) mode of the quasinormal spectrum of the field around the black hole.

Field dynamics
The dynamics of a massive charged test scalar field Φ of mass µ and electric charge q in the background (2) is governed by the Klein-Gordon equation, where A ν = −δ 0 ν e/r is the electromagnetic vector potential of the black hole. We can decompose the field Φ as where ω is the conserved frequency, l is the spherical harmonic index and m (−l ≤ m ≤ l) is the azimuthal harmonic index. Hereafter we will drop the subscripts l and m for brevity.
With the decomposition (10) one can separate the Klein-Gordon equation (9) into a radial and an angular equation with the separation constant K l = l (l + 1). The radial Klein-Gordon equation is given by where ∆ = r 2 f (r) and If we define a new radial function ζ = rR and adopt the tortoise coordinate r * ( defined by, dr * = dr/ f (r) ), mapping the semi infinite region [r + , ∞) into (−∞, ∞), then the radial Klein-Gordon equation (11) becomes where In the asymptotic limit equation (13) can be solved analytically with the quasinormal mode (QNM) boundary conditions of purely ingoing waves at the horizon (r * → −∞) and purely outgoing waves at spatial infinity (r * → ∞), where Ω = ω 2 − µ 2 . Equation (13) together with the boundary conditions (15) becomes an eigenvalue problem with complex eigenvalues ω representing the quasinormal frequencies.

Continued Fraction technique
In 1985, Leaver [21][22][23] inspired by a seminal work of Jaffé [27] on the calculation of the electronic spectra of hydrogen molecular ion, proposed a very accurate method for finding out the QNM frequencies of black holes.
To implement Leaver's method we start with equation (11) and observe that it has two regular singularities at r + and r − and an irregular singularity as r → ∞ .
We can write a solution to equation (11) with the desired behaviour at the boundaries as where Substituting the ansatz (16) into equation (11) we arrive at the following three term recurrence relations, satisfied by the coefficient a n α 0 a 1 + β 0 a 0 = 0, (17) α n a n+1 + β n a n + γ n a n−1 = 0, where α n , β n and γ n are given by, 4 The convergence of the series (16) requires the recursion coefficients to satisfy an infinite continued fraction relation The solution to this infinite continued fraction equation gives the QNM frequencies. The continued fraction relation (22) can be inverted any number of times. Numerically, the n th QNM frequency is defined to be the most stable root of the n th inversion of the continued fraction relation, In practice the infinite continued fractions in equations (22, 23) are truncated at some large truncation index, N. Nollert [28] has shown that the "error" due to truncation can be minimised and the convergence of the method can be improved by a wise choice of the "remaining" part of the infinite continued fraction , a N , which in turn satisfies the recurrence equation, Assuming that R N can be expanded in a power series of we obtain the first three coefficients C k as,

Numerical Results
We first study the fundamental QNMs due to uncharged massive scalar field then we add electric charge to the perturbing field and study the effect of the scalar hair on the QNMs. For the sake of numerical simplicity we scale the mass of the black hole to unity.

Uncharged massive scalar field
For an uncharged scalar field (q = 0) in the shRN background we assume, without any loss of generality, the constant electric charge of the black hole to be zero,(e = 0). The function W appearing in equation (13) can now be written as The potential of the of the shRN black hole for different values of the scalar charge, field mass and multipole index are shown in Fig. 1.
In Figs. 2(a) and (b), we show the behaviour of the real and imaginary parts of the fundamental QN frequency with the scalar charge for a particular multipole index (l = 1) and different field masses. We observe that for s < 0, the magnitude of both the real and imaginary parts of the QN frequency decrease with the absolute value of the scalar charge. This implies that the real oscillation frequency as well as the damping rate decrease with increasing magnitude of the negative scalar charge. For s > 0, the spacetime (2) effectively behaves as an RN black hole of unit mass and electric charge, e = √ s, showing a distinct peak in the magnitude of the imaginary part of the fundamental quasinormal frequency (see Refs. [29,30]). We also observe that for a particular value of the scalar charge, the real part of the QN frequency increases with the field mass whereas the magnitude of the imaginary part decreases. This behaviour is manifested more clearly in Figs. 2(c) and (d), where we observe that for sufficiently large field masses, the imaginary part of the QN frequency becomes vanishingly small. This results in long lived, purely real modes in the quasinormal spectrum, called quasi-resonance modes [31]. We also note that as the scalar charge changes from positive to negative,quasiresonance occurs at lower field masses with smaller real frequencies. Fig. 3(a) shows a compact view of the behaviour of the real and imaginary parts of the QN frequency with the scalar charge for different values of the multipole number and field mass. We note that, as the multipole number increases the real part of the fundamental QN frequency increases and so does the imaginary part, but only for higher field masses. This behaviour of the imaginary part of the QN frequency can be seen more clearly in Fig. 3(b) where we note that for lower field masses, the damping rate decreases with the multipole index whereas for higher field masses, it increases with the multipole index. For large values of the multipole number, the damping rate is almost insensitive to the field mass.
Following Ref. [32] we define the Quality Factor as Q.F. ∼ | ω Re ω Im |. In Fig. 4(a), we observe that for a given multipole index and for large positive values of the scalar charge the quality factor decreases sharply, however for smaller values of the scalar charge it decreases very gradually. For massless field, the gradual decrease of the quality factor continues to persist for negative values of the scalar charge as well. However, beyond a certain value of the field mass µ, it tends to increase for high negative values of s, the plot corresponding to µ = 0.2 in Fig. 4(a) reveals this feature. The quality factor has higher values for higher multipole indices (see Fig. 4(b)). The quality factor is a measure of the product of the frequency and the ring down time of a black hole radi-  ation, and is an important tool to figure out the black hole parameters [33].

Charged scalar field
The presence of the scalar hair changes the frequency and damping rate of the QN spectrum of charged scalar fields as well. In Figs. 5(a) and (b), we observe that compared to the RN black hole, for fixed non-zero values of e and q, the magnitude of both the real and imaginary parts of the fundamental QN frequency are higher for positive values of the scalar charge and lower for negative values. Konoplya [29] observed that for an RN black hole, the imaginary part of the QN frequency, for any given value of the field charge, approaches that for the uncharged field as the extremal limit is approached. Apart from a similar observation in the presence of a scalar hair (see Figs. 5(c) and (d)) as well, we note from Fig. 5(f) that such convergence of the imaginary part of the fundamental QN frequency occurs for any given value of the black hole electric charge, as the maximal value of the scalar charge is approached. This maximal value is determined by the extremality condition (8). We further observe that, for s > −e 2 , the magnitude of the imaginary part of the fundamental QN frequency shows a distinct peak whereas for s < −e 2 , it decreases monotonically with the magnitude of s. The corresponding behaviour of the real part of the QN frequency with the scalar charge is shown in Fig. 5(e).  The symmetry of the QNMs with respect to the transformation (eq → −eq, ω → −ω * ) is depicted in Fig. 6. Fig. 6(a) also highlights the existence of a critical value of |eq| at which the real part of the QN frequency vanishes. However, such a behaviour of the QN frequency is not new and has been previously observed for the RN black hole for charged scalar and Dirac fields (see Refs. [34,35]). We note in particular, that the critical value of |eq| is almost unaffected by the presence of the scalar hair and does not change with the the black hole electric charge. For an RN black hole with unit multipole index, the critical value is |eq| ≈ 1.3.

Charged Dirac field around a charged black hole with scalar hair
The dynamics of a massless charged Dirac field propagating in the shRN spacetime is given by the Dirac equation, where Ψ is the Dirac four-spinor, γ µ are the coordinate dependent Dirac four-matrices and D µ is spinor covariant derivative defined by, Here q is the charge of the Dirac field and Γ µ are the spinor connection matrices. Following Refs. [36,37] we decompose the Dirac four-spinor as and use the canonical orthonormal (symmetric) tetrad, proposed by Carter [38], to yield two pairs of coupled first order differential equation, and where K = ωr 2 − eqr, − j ≤ m ≤ j and λ = j + 1/2 ( with j = 1/2, 3/2...) is the separation constant. The radial equations (30,31) can then be combined to yield, If we define a new radial function, ξ = ∆ −1/4 rP, then equation (34) can be written in a Schrödinger like form in terms of the tortoise coordinate as wherẽ In the asymptotic limits of the tortoise coordinate equation (36) can be solved analytically with the QNM boundary conditions yielding Equation (34), similar to equation (11), also has two regular singularities at r + and r − and an irregular singularity as r → ∞. So proceeding as before we introduce an ansatz, consistent with the boundary conditions (37), where u = r−r + r−r − ,ρ = −ieq + iω (r + + r − ) and (38) back into equation (34) we again arrive at the three term recurrence relations, α n b n+1 +β n b n +γ n b n−1 = 0, wherẽ α n = (n + 1) The convergence of the series (38) demands the recurrence coefficients to satisfy an infinite continued fraction relation , respectively, show the real and imaginary parts of the fundamental (Dirac) QN frequency as a function of eq for j = 1/2 and e = 0.5. Each curve in each panel corresponds to a particular value of s as indicated. The dashed lines represent the symmetric curves due to the simultaneous transformation eq → −eq and ω → −ω * . similar to equations (22,23). Applying Nollert's improvement, we now get the first three coefficients of the series (25) as, C 0 = −1, C 1 = 2iω (r − − r + ) and C 2 = 5 4 − ieq + 2iωr + .

Numerical Results
The behaviour of the real and imaginary parts of the fundamental QN frequency with the scalar charge is shown in Figs. 7(a) and (b). As before, we observe that for a fixed value of the black hole electric charge, the magnitude of the imaginary part of the QN frequency for a given value of the field charge increases as the scalar charge changes from neg-ative to positive and ultimately approaches the neutral one in the extremal limit. Thus, in the extremal limit the damping rate is independent of s. Here also we observe a peak in the magnitude of the imaginary part of the fundamental QN frequency for s > −e 2 whereas for s < −e 2 , it decreases monotonically with the magnitude of the scalar charge. The real QN frequency on the other hand continues to grow with the scalar charge. Away from the extremal value of the scalar charge this causes the quality factor to grow steadily (see Fig. 7(c)), however as the extremal value of the scalar charge is reached the growth of the quality factor becomes quite rapid.
Similar to the scalar case, we observe in Fig. 8 that the QN frequency is symmetric with respect to the transformation (eq → −eq, ω → −ω * ) and note that the critical value of electromagnetic interaction (|eq| = 0.7) at which the real part of the QN frequency vanishes, is almost unaffected by the presence of the scalar hair.

Summary and Discussion
In the present work we discussed the QN spectrum of massless and massive uncharged as well as charged scalar fields and massless charged Dirac fields in the vicinity of a charged spherically symmetric black hole with a scalar hair dubbed as the "shRN" black hole. We mainly focussed on negative values of the scalar charge with s < −e 2 , for which the metric (2) represents a "mutated RN" spacetime mimicking the Einstein-Rosen bridge.
Unlike the appearance of a distinct peak in the magnitude of the imaginary part of the fundamental QN frequency of an shRN black hole for scalar and Dirac fields with s > −e 2 , the mutated RN spacetime (s < −e 2 ) is characterised by monotonically decreasing |Im(ω)| (see Figs. 2(b), 5(f) and 7(b)). For uncharged fields, the shRN black hole effectively behaves as an RN black hole with effective electric charge, e e f f = √ e 2 + s, provided s lies in M 2 −e 2 ≥ s > −e 2 . For massive scalar field, the phenomenon of quasiresonance, characterised by vanishingly small |Im(ω)| is observed. We also showed the behaviour of the quality factor with the scalar charge for both the scalar and Dirac fields. As the extremal limit is approached either by increasing the electric charge for a fixed s or vice-versa, we find that the imaginary part of ω for neutral and charged scalar or Dirac perturbations to be coincident.
In the presence of electric charge of the perturbing fields, we observe the existence of a critical value of |eq|, above which the real part of the QN frequency vanishes, for both the scalar and Dirac fields. This value is completely unaffected by the presence of scalar hair.
Following the method of Cho [39], we start with the asymptotic form of the Dirac QN frequency and calculated the area spectrum of the shRN black hole based on the proposals of Kunstatter [40] and Maggiore [41]. We obtain the area quantum as ∆ A = 8πh. This being the same as that of an RN black hole [42], we refrain from including a detailed calculation of the same.
Very recently Saleh, Thomas and Kofane [43,44] discussed the QN spectrum of massless uncharged scalar and Dirac fields in the vicinity of a "quantum-corrected" Schwarzschild black hole [45] using 3rd order WKB approximation. The metric used by them is effectively similar to that of the mutated RN spacetime discussed in the present work. The results obtained by us for the massless uncharged scalar and Dirac fields in the shRN background (with s < −e 2 ) using the more accurate continued fraction method, is qualitatively similar to them. The present work is, however, much more general as it includes charge for both the scalar and the Dirac fields and mass for the scalar field.
The QN spectrum analysis was also carried out with the 3rd order WKB approximation which is generally believed to be less accurate. We add a table in the appendix comparing the results for one example, namely that for an uncharged massless scalar perturbation of the shRN black hole. The WKB approximation is known to yield more and more accurate results for higher and higher values of multipole (l). The table contains the values given by the Leaver method and that by the WKB approximation for the real and the imaginary parts of the quasinormal mode frequencies for l = 1, 2, 3. It is apparent from the table that the difference in the results given by the two methods reduces for higher values of l.
We point out that in the present work we have not observed any QNM with positive imaginary part indicating the stability of the shRN black hole under massive (and massless) charged (and uncharged) scalar perturbations as well as under massless charged (and uncharged) Dirac perturbations for both s in [−e 2 , M 2 − e 2 ) and in the "mutated" regime, s < −e 2 . This implies that the mutated RN spacetime is also stable under all the above mentioned types of perturbation.
It has already been pointed out in section 2 that the scalar field in this case is a constant (ψ = ± 6 8πG s s+e 2 ). From equation (1), it is clear that the action can be transformed into an Einstein-Maxwell system, with a different value for the effective Newtonian constant of gravity G e f f = However, the spacetime is qualitatively different if s is negative, and |s| > e 2 . The resulting metric component is in the standard RN spacetime. The distinctive features of this metric has already been discussed at length. Finally, we mention that the qualitative features for the quasinormal modes for the charged or uncharged massive scalar fields and also that for the charged Dirac field for a mutated RN (shRN with s < −e 2 ) background are qualitatively same as that for a usual RN black hole except for the complete monotonic behaviour of the damping (imaginary part of the QNM) in the case of the former as opposed to the existence of a peak in the latter case, an RN black hole. Table 1: Fundamental QN frequencies of massless uncharged scalar fields in the background of an shRN black hole of mass, M = 1 and electric charge, e = 0 for different values of the multipole index and scalar charge. For each value of s, the first line is obtained using the continued fraction method with 175 terms and the second line is obtained using the 3rd order WKB approximation.