Charge superradiance on charged BTZ black holes

We study superradiant scattering for a charged scalar field subject to Robin (mixed) boundary conditions on a charged BTZ black hole background. Scalar field modes having a real frequency do not exhibit superradiant scattering, independent of the boundary conditions applied. For scalar field modes with a complex frequency, no superradiant scattering occurs if the black hole is static. After exploring some regions of the parameter space, we provide evidence for the existence of superradiantly scattered modes with complex frequencies for a charged and rotating BTZ black hole. Most of the superradiantly scattered modes we find satisfy Robin (mixed) boundary conditions, but there are also superradiantly scattered modes with complex frequencies satisfying Dirichlet and Neumann boundary conditions. We explore the effect of the black hole and scalar field charge on the outgoing energy flux of these superradiantly scattered modes, and also investigate their stability.


Introduction
In the scattering of waves incident on a black hole, superradiant scattering occurs if the reflected wave has greater amplitude than the incident wave [1].For example, lowfrequency bosonic waves interacting with a rotating Kerr black hole exhibit superradiant scattering [2][3][4][5].This phenomenon is the wave analogue of the Penrose process for particles [6], with the wave extracting rotational energy from the black hole.A similar effect occurs for charged scalar field waves scattered by a charged Reissner-Nordström black hole [7][8][9][10][11].In this case it is electrostatic rather than rotational energy which is extracted from the black hole by the wave.The amplification of superradiantly scattered scalar field modes on a Kerr black hole is comparatively small, for example for a Kerr black hole with angular momentum parameter a = 0.99M (where M is the black hole mass), the maximum amplification factor for a massless scalar field mode is about 0.4% [1].In contrast, charged scalar field modes undergoing superradiant scattering on a Reissner-Nordström black hole can achieve significantly larger amplification factors, up to 40% [1].For this reason, charged superradiant scattering is a very useful toy model for numerical studies.The fact that the Reissner-Nordström black hole is static and spherically symmetric also simplifies the mathematical analysis of charged superradiant scattering compared with that on rotating black holes.
One consequence of superradiant scattering is the possibility that it can lead to a superradiant instability, if there is some mechanism confining a superradiantly amplified field [3].For example, a massive scalar field on a rotating Kerr black hole has a potential which is confining and exhibits a superradiant instability [12][13][14].However, superradiantly unstable scalar field modes on a Kerr black hole grow only very slowly in time, which presents major challenges for numerical studies of the evolution of the superradiant instability [15,16].
For a charged scalar field on a Reissner-Nordström black hole, the presence of a mass is not sufficient to confine the field and cause an instability [17,18].An alternative confining mechanism is therefore required if there is to be a superradiant instability.In the "black hole bomb" scenario [3], the confining mechanism is provided by a perfectly reflecting mirror.For a charged scalar field on a Reissner-Nordström black hole, there is an instability if the mirror is sufficiently far from the event horizon of the black hole [19,20].This instability can grow many times faster than the corresponding instability of a massive scalar field on a Kerr black hole, which enables the numerical evolution of the instability, leading to the formation of a black hole with scalar hair [21,22].
An alternative confining mechanism can arise when considering asymptotically anti-de Sitter (adS) rather than asymptotically flat black holes.In this situation superradiant scattering persists, and the adS boundary can act like a confining mirror.A massless scalar field on a sufficiently small rotating asymptotically-adS black hole is superradiantly unstable, and the same is true for a charged scalar field on a sufficiently small Reissner-Nordström-adS black hole (see, for example,  for a selection of papers from the extensive literature on this topic).
The three-dimensional, asymptotically adS, BTZ black hole [48][49][50] has been widely studied as a toy model for phenomena on black holes in four or more space-time dimensions.In particular, one can consider whether superradiant scattering, or a superradiant instability, is possible on a BTZ black hole.The analysis of these questions for a neutral scalar field is greatly simplified by the fact that the field modes can be written in terms of hypergeometric functions (see, for example, [51,52]).Since null infinity is a timelike hypersurface in adS, boundary conditions have to be applied to a perturbing field.For a neutral scalar field, applying the simplest boundary conditions, namely Dirichlet boundary conditions, superradiant scattering is absent [53].Superradiant scattering also does not occur for a fermion field vanishing on the boundary [54].However, one can impose more general (Neumann or Robin) boundary conditions on the scalar field [55], and, for at least some Robin boundary conditions, there exist modes which exhibit superradiant scattering [56].The impact of Robin boundary conditions on superradiant scattering has also been considered on four-dimensional Reissner-Nordström-adS [44] and Kerr-adS black holes [57].
In this paper we examine whether superradiant scattering exists for a charged scalar field subject to Robin (mixed) boundary conditions on a charged analogue of the BTZ black hole.We seek to understand whether the analogy between superradiant scattering behaviour on a BTZ black hole and four dimensional black holes persists for a charged scalar field.In particular, we ascertain whether there is charged superradiant scattering on a static charged BTZ black hole.We also examine whether superradiant amplification is enhanced in this set-up compared with superradiant scattering on a rotating BTZ black hole.
We begin, in Section 2, by reviewing the black hole metric and separable solutions of the charged scalar field equation on this spacetime background, paying particular attention to the boundary conditions satisfied by the field far from the black hole.In Section 3 we study the possibility of su-perradiant scattering using two approaches: firstly a Wronskian condition (arising from the conservation of the charge current) which is valid for waves having real frequency, and secondly, considering the fluxes of energy and charge down the horizon for an ingoing wave.In particular, an ingoing wave which is superradiantly scattered will have an outgoing energy flux.Using the Wronskian, we find that there is no superradiant scattering for field modes having real frequency.If the frequency is complex, by considering the energy flux, we show that superradiant scattering is absent if the black hole is nonrotating.This leaves open the possibility of superradiant scattering for charged, rotating black holes, which is studied in Section 4. Using a simple numerical method, valid for frequencies in the superradiant regime, we provide evidence for the existence of charged scalar field modes which are superradiantly scattered.We explore the effect of increasing either the black hole or scalar field charge on the energy fluxes of these superradiantly scattered modes, as well as investigating their stability.Our conclusions are presented in Section 5. Two appendices give further details of our numerical procedure and the complex frequencies of the superradiantly scattered modes.
2 Charged scalar field on a charged BTZ black hole

Charged BTZ black holes
The neutral BTZ black hole [48][49][50] is a solution of the three-dimensional Einstein equations with a negative cosmological constant Λ = −ℓ −2 , having metric where with M the mass and J the angular momentum of the black hole.
In the static case (J = 0), the black hole acquires an electric charge Q by introducing the electromagnetic potential A µ = A 0 δ t µ , where and an arbitrary length scale r 0 has been introduced to render the argument of the logarithm dimensionless.The metric for a static charged BTZ black hole is then [48][49][50]] and the lapse function takes the form The generalization of the rotating BTZ black hole to include an electric charge is far from straightforward [50,58].Various three-dimensional, charged, rotating black holes are presented in [59][60][61][62][63].In this paper we consider the following charged generalization of the BTZ metric [63]: where N(r) is the same as in the static case (2c) and the other functions appearing in the metric are: where M, Q and Ω ∈ [0, 1) are constants.The mass M and angular momentum J of the black hole are given in terms of the parameters M, Q, and Ω as follows [63]: Unlike the neutral BTZ black hole, the metric (3a) cannot be obtained by identifying points in three-dimensional adS space-time.In particular, the scalar curvature R is not constant: When Ω > 0, the electromagnetic potential A µ acquires a nonzero magnetic part and is given by [63]: The nonzero components of the electromagnetic field strength tensor For comparison, the electrostatic potential for a four-dimensional Reissner-Nordström black hole with charge Q RN is, for a suitable choice of gauge, from which we find giving the usual Coulomb law for the electric field due to a charge Q RN located at the origin in three spatial dimensions.
Therefore Q RN is the charge of the Reissner-Nordström black hole.Comparing (7,9), it can be seen that there is a sign difference in F rt .The electric field component F rt in (7) has the expected r −1 behaviour for a point charge at the origin in two spatial dimensions, but is proportional to −Q rather than Q.As a result, the charge Q of the black hole is given by The fact that the black hole charge Q has the opposite sign to Q will have important consequences for the interpretation of the charged scalar field modes which we study later in this paper.
In the limit Ω → 0, the metric (3) and gauge field potential (6) reduce to those in the static case (2).If we set Q = 0, the metric (3a) does not reduce to the original rotating BTZ metric (1a) in (t, r, ϕ) coordinates.However, using R as the radial coordinate, when Q = 0 the metric (3a) becomes where we have defined the function We therefore have a metric of the form (1a) with mass M = M 1 + Ω 2 / 1 − Ω 2 2 and angular momentum J = 2MΩ ℓ/ 1 − Ω 2 , in accordance with (4).The horizons of the black hole are located at those values of the radial coordinate r for which N(r) vanishes.If M < Q 2 [1 − 2 ln (Qℓ/2r 0 ) ] /4 there is a naked singularity at r = 0; we do not consider this possibility further.For M > Q 2 [1 − 2 ln (Qℓ/2r 0 )] /4 there is an event horizon at r = r h , the largest zero of N(r) and an inner horizon at the smaller positive zero of N(r).These two horizons coincide when M = Q 2 [1 − 2 ln (Qℓ/2r 0 )] /4 and in this case we have an extremal black hole [64].In this paper we focus on the case where the black hole is nonextremal.
By making a gauge transformation of the form we may set r 0 = r h without loss of generality.We then find, by considering the zeros of (2c), that At the horizon, we have R(r h ) 2 = r 2 h / 1 − Ω 2 and N ϕ (r h ) = −Ω /ℓ, so that Ω /ℓ is the angular speed with which the event horizon rotates.
In our analysis of superradiant scattering, we will be interested in the flux of energy down the event horizon of the black hole.For this analysis, we will require suitable coordinates which are regular across the horizon.We will employ ingoing Eddington-Finkelstein (EF) coordinates.First we define an ingoing null coordinate v by and a new angular coordinate ϕ v by Then the coordinates (v, r, ϕ v ) are ingoing EF coordinates, in terms of which the metric (3a) becomes The metric ( 15) is regular when r = r h and N(r) = 0, as required.Near the horizon, the ingoing EF coordinates take the form where r * is the usual tortoise coordinate, defined by Similarly, we can define outgoing EF coordinates (u, r, ϕ u ) as follows: Near the horizon, the outgoing EF coordinates take the form

Charged scalar field
We consider a scalar field Φ with charge q and mass m propagating on the rotating charged black hole (3a), and satisfying the charged scalar field equation where D µ = ∇ µ − iqA µ is the covariant derivative.We assume that the scalar field is minimally coupled to the geometry.The stress-energy tensor for the charged scalar field is where ℜ is the real part and a star is used to denote complex conjugation.The scalar field has an associated current density J µ given by In this paper we treat the scalar field as a test field.However, if we were to go beyond this approximation, and consider the effect of the charged scalar field on the electromagnetic field, the current density J µ (22) would act as source term in Maxwell's equations: In Section 3 we shall study in particular the radial component J r of the current density (22).Mode solutions of the scalar field equation ( 20) take the form where ω is the frequency of the wave (which may be complex) and k ∈ Z is the azimuthal quantum number.In terms of the tortoise coordinate r * (17), the radial function X ωk (r) satisfies the equation where the potential V ωk (r) takes the form As r → r h and the event horizon is approached, we have r * → −∞ and where we have defined Therefore, near the horizon, the radial function X ωk (r) takes the form X ωk (r) ∼ A ωk e i ωr * + B ωk e −i ωr * (28) where A ωk and B ωk are complex constants.The frequency of the wave has effectively been shifted due to the rotation of the black hole.The fact that ω does not depend on the charge stems from our choice of gauge, in that the electromagnetic gauge potential (6) vanishes at the horizon since we have taken r 0 = r h .Far from the black hole, as r → ∞, the leading-order behaviour of the potential (25b) is, in general, This leading-order behaviour is the same as for the neutral scalar field, and does not depend on the frequency ω or the azimuthal quantum number k.In this regime the tortoise coordinate has the following form: yielding the equidimensional differential equation where µ 2 is a constant given by Let us assume for the moment that µ 2 ̸ = 0.The solutions of (31) are X ωk ∼ r p * ∼ r −p , where For 4µ 2 > −1, the values of p are real and for complex constants C ωk , D ωk .The second term gives a scalar field mode ( 24) which diverges at infinity when 4µ 2 > 3 and therefore we set D ωk = 0 in this case.In this situation there is no choice of boundary conditions which can be imposed on the scalar field at infinity.For −1 < 4µ 2 < 3, both solutions in (34) give regular scalar field modes (24), resulting in some freedom in the choice of boundary conditions at infinity [65,66].The solution with D ωk = 0 satisfies Dirichlet boundary conditions, while, following [65], we define Neumann boundary conditions to be such that C ωk = 0.If both C ωℓ and D ωℓ are nonzero, then we have Robin (mixed) boundary conditions.In this situation we write the solution (34) in the form [56,65] where E ωk is a complex constant and the real angle ζ (which we term the "Robin parameter") can be taken to lie in the interval 0 ≤ ζ < π (we could equally well take ζ ∈ (− π 2 , π 2 ]).Setting ζ = 0 yields Dirichlet boundary conditions, while ζ = π 2 corresponds to Neumann boundary conditions.When 4µ 2 = −1, we have Both solutions are square-integrable in this case, so again we have a choice of boundary conditions.For 4µ 2 < −1, the exponent p (33) is complex and X ωk (r) is oscillatory.Once again both linearly independent solutions of the radial equation are square-integrable at infinity.However, these values of µ 2 violate the Breitenlohner-Freedman bound [67,68] and therefore we do not consider them further in this work.
The above discussion of the boundary conditions at infinity is valid only when µ 2 ̸ = 0, in which case the behaviour of the charged scalar field at infinity is identical to that for a neutral scalar field.In the special case µ 2 = 0 the leading order behaviour of the potential (25b) is no longer (29), but instead we have, as r → ∞, In this case it is not possible to solve the asymptotic form of the radial equation exactly in terms of elementary functions.However, it is possible to perform an asymptotic expansion for the radial function X ωk (r) in this case.The first couple of terms in this asymptotic expansion are: The second solution gives a mode which is not square integrable at infinity, so we set D ωk = 0 in this case.The behaviour at infinity of a massless and conformally coupled charged scalar field is thus rather different from that seen in the neutral case.

Criterion for superradiant scattering
We now explore whether superradiant scattering occurs for a charged scalar field, examining separately the cases where the frequency ω is real or complex.

Real frequency
We first consider the situation in which the frequency ω is real and follow an approach of [44,53].Near the horizon of an eternal charged BTZ black hole, a scalar field wave will have a radial function of the form (28), giving Using (16,19) this can be rewritten in terms of ingoing and outgoing EF coordinates: The second term corresponds to an ingoing wave with amplitude |B ωk |, and the first term to a wave which is outgoing close to the horizon and has amplitude |A ωk |.Since we are interested in superradiantly scattered modes, we set A ωk = 0 so that the wave is purely ingoing at the horizon.
When the frequency ω is real, the potential V ωk (r) (25b) is also real and therefore the Wronskian is a constant for any solution X ωk of the radial equation (25).Near the horizon, using (28) with A ωk = 0 we find The quantity ( 42) is proportional to the radial component of the current (22) evaluated near the horizon for a scalar field mode having a radial function of the form (28) with real frequency and A ωk = 0, namely The fact that the Wronskian is a constant reflects the conservation of the current J µ .The value of W ωk as r → ∞ depends on the form of the radial function X ωk (r) in this regime.Consider first the solution (34) valid when 4µ 2 > −1.In this case we have where ℑ denotes the imaginary part.Therefore, if we specify Dirichlet boundary conditions (for which D ωk = 0) or Neumann boundary conditions (for which C ωk = 0), equating (42,44) gives that |B ωk | 2 = 0, and it is therefore not possible to have a purely ingoing mode with real frequency at the horizon, generalizing the result of [53] to the charged case.
When −1 < 4µ 2 < 3, µ 2 ̸ = 0, Dirichlet and Neumann boundary conditions are not the only possibility, we can also impose Robin boundary conditions for which both C ωk and D ωk are nonzero.Using the parameterization (35), we have C ωk = E ωk cos ζ and D ωk = E ωk sin ζ and the Wronskian (44) becomes Therefore, when Robin boundary conditions are applied, it is again the case that |B ωk | 2 = 0.There are two special cases which need to be considered separately.First, when 4µ 2 = −1, the radial function X ωk (r) has the form (36) as r → ∞, whence Our conclusions are however unchanged: (46) vanishes for Dirichlet, Neumann and Robin boundary conditions and Finally we have the case 4µ 2 = 0, for which the radial function takes the form (38) as r → ∞, with D ωk = 0 to ensure square integrability.In this situation the Wronskian tends to zero as r → ∞, so that again it must be the case that The inclusion of a scalar field charge has made no difference to the analysis of the Wronskian.In particular, we find that for charged scalar field modes having real frequency, any mode which is purely ingoing at the event horizon vanishes identically, irrespective of the boundary conditions.In particular, this means that there can be no superradiant scattering for modes with real frequency.
Similar results were found in [44] for a charged scalar field on a four-dimensional Reissner-Nordström-anti-de Sitter black hole.As in [44], the Wronskian (and hence the radial component of the current J r (22)) vanishes on the boundary for Robin boundary conditions.As a result, the radial component of the current J r (43) must also vanish on the event horizon.Therefore there can be no net positive flux from the horizon, and no superradiant scattering.
Our results in this section are in accordance with those of [56] for the neutral scalar field.In that case there are superradiantly scattered modes when Robin boundary conditions are applied, but these modes have complex frequencies.Our analysis in this subsection is valid only when the modes have real frequencies.In the next subsection we therefore examine the possibility of superradiant scattering for modes with complex frequencies using a different strategy.

Fluxes of energy and charge down the horizon
If the frequency ω is complex, it is no longer the case that the Wronskian ( 41) is a constant and the approach taken in Subsection 3.1 is not applicable.In this subsection, we therefore take an alternative perspective to investigate whether there are superradiantly scattered modes having complex frequency ω.Following, for example [5,56], we consider the energy flux down the event horizon due to an ingoing mode.For the remainder of this paper the frequency ω will be complex.
Consider an ingoing mode given, in terms of ingoing EF coordinates (14), by where B ωk is a complex constant, and whose radial part is therefore In terms of the Killing vectors where Ω /ℓ is the angular speed of the event horizon, the flux of energy down the black hole is [56] where the stress-energy tensor for the charged scalar field is given by (21).Evaluating the required components of the stress-energy tensor (21) gives where Thus the flux of energy down the horizon due to an ingoing mode will be positive unless This is exactly the same condition for superradiant scattering as in the neutral case [56].Following the analysis in [44], we also consider the radial component of the scalar current (22).For an ingoing mode of the form (47), near the horizon this is given by Assuming that q > 0, this will be negative if ℜ(ω) > kΩ ℓ −1 , corresponding to an ingoing flux of charge.However, if ℜ(ω) < kΩ ℓ −1 and q > 0, we have J r > 0, and there is an outgoing flux of charge.An ingoing mode with q > 0 and satisfying the criterion ( 52) is therefore extracting both charge and energy from the black hole.In this situation the scalar field mode has been superradiantly scattered.
For a nonrotating charged black hole we have and if q > 0 and ℜ(ω) > 0. Therefore an ingoing mode has ingoing fluxes of both energy and charge and there is no superradiant scattering for a charged scalar field on a static, charged BTZ black hole, irrespective of the boundary conditions applied to the field.This result is in contrast to that for a charged scalar field on a four-dimensional, nonrotating, asymptotically flat Reissner-Nordström black hole [44], for which superradiant scattering can occur.The key difference between these two situations is that, for charged scalar field scattering on a Reissner-Nordström black hole, the electrostatic potential at the horizon plays the same role as kΩ ℓ −1 does in (53), resulting in a range of ingoing mode frequencies for which superradiant scattering is possible.In our case, the electromagnetic potential (6) vanishes on the horizon.Therefore the black hole charge Q does not appear in the fluxes of energy (50) or charge (53) at the event horizon.

Superradiant scattering on charged rotating BTZ black holes
We now provide evidence for the existence of superradiantly scattered modes satisfying the condition (52), using a numerical method.We restrict attention to the regime −1 < 4µ 2 < 3, µ 2 ̸ = 0 for which there is a choice of boundary conditions that can be applied to the scalar field at infinity.
Our numerical method involves integrating the radial equation (25) to find complex frequencies ω for which the radial function X ωℓ (r) satisfies ingoing boundary conditions (48) at the horizon and Robin boundary conditions (35) at infinity.In other words, we are seeking quasi-normal modes (QNMs).There are many methods in the literature for the accurate computation of QNM frequencies (see, for example [69][70][71][72] for reviews and [73][74][75][76][77][78][79][80][81][82][83][84][85] for a selection of references concerning QNMs on BTZ black holes).However, the form of the potential (25b) (in particular, the presence of the nonanalytic ln r r h term) hinders implementing these in our situation.Our aims in this section are rather less ambitious than the high-precision computation of QNM frequencies.Instead, we are looking for numerical evidence for the existence of superradiantly scattered modes, and some qualitative information about the energy flux (50) for these modes.With this in mind, we employ a rather naive direct integration method, which is sufficiently accurate for our purposes for modes lying in the region (52) for which superradiant scattering is possible.Our computations are described in Appendix A and are implemented in MATHEMATICA.
According to (52), superradiantly scattered modes can exist only for frequencies having a real part ℜ(ω) satisfying 0 < ℜ(ω) < kΩ /ℓ for k > 0. We therefore consider only frequencies whose real parts lie in this interval.For a general frequency ω, the constants C ωk and D ωk (34) describing the behaviour of the scalar field mode at infinity will be complex.To apply Robin boundary conditions (35), we require the ratio D ωk /C ωk to be real, which will only occur at particular frequencies.When the imaginary part of D ωk /C ωk vanishes, we determine the parameter ζ (35) governing the Robin boundary conditions from and the energy flux F E /F using (50).We mostly take a branch of the arctan function such that ζ ∈ [0, π); however in part of our analysis it will be helpful to consider instead ].We find that our numerical method yields satisfactory results only when either the scalar field charge vanishes (q = 0) or for reasonably large values of at least one of the charges |Q|, q.In order to obtain good results for a wider range of values of the charges, and for nonsuperradiantly scattered modes, a more sophisticated method would be needed to find the QNM frequencies.However, our method is sufficiently accurate to give a selection of superradiantly scattered modes which enables us to qualitatively explore the effect of black hole and/or scalar field charge on superradiant scattering.Below we consider first the energy flux (50) from superradiantly scattered modes and then the real and imaginary parts of the mode frequencies.
For our discussion of superradiantly scattered modes, it will be important to recall, from (10), that the black hole charge Q has the opposite sign to the parameter Q.Therefore, if qQ > 0, we have q Q < 0 and the scalar field charge has the opposite sign to the black hole charge, while if qQ < 0, then q Q > 0 and the scalar field charge has the same sign as the black hole charge.The radial equation ( 25) depends only on the product qQ (explicit in the first term in the potential V ωk (25b)) and Q 2 (which appears in the metric function N(r) (2c) and its derivative).Changing the sign of the black hole charge parameter Q (with the scalar field charge q fixed) will only impact the second term in the square brackets in the first line of the potential (25b).All other terms in the potential are unaffected on changing the sign of Q.In this section we therefore consider positive scalar field charge q > 0 and both positive and negative values of the parameter Q.

Energy flux
In Figures 1-4 we present our numerical results providing evidence for the existence of superradiant scattering for a charged scalar field on a charged, rotating BTZ black hole.To aid comparison with the results for a neutral scalar field in [56], we set m 2 = −0.65,ℓ = 1, M = 16 (which corresponds to Figure 2 in [56]) and consider only modes with azimuthal quantum number k = 1.In all our plots we show the energy flux F E /F (50)   We begin by setting the scalar field charge q = 0, see Figure 1, where we have also fixed the rotation parameter Ω = 0.6 and varied the black hole charge parameter Q.For a neutral scalar field with q = 0, the radial equation (25) does not depend on the sign of the charge parameter Q so it is sufficient to consider Q ≥ 0. When Q = 0, we reproduce the results in [56], which provides verification of our numerical method.QNM for Q ̸ = 0 and q = 0 were studied in [64], although the focus in that work was the mode frequencies rather than the energy flux, as is the case here.It can be seen in Figure 1 that, as the black hole charge parameter Q increases, the maximum magnitude of the energy flux F E /F in superradiantly scattered modes decreases, and that superradiantly scattered modes exist for smaller values of the Robin parameter ζ .The effects of superradiant scattering are small in this situation: there is only a narrow interval of values of ζ for which there are superradiantly scattered Energy Flux (q=1)  Energy Flux (Q=2) q = 0.7 q = 0.9 q = 1.2 q = 1.5 q = 1.7 (c) Energy Flux (Q= -2) q = 0.7 q = 0.9 q = 1.2 q = 1.5 q = 1.7 modes, and the resulting fluxes of energy have small magnitudes.We also note that all values of the Robin parameter ζ for the superradiantly scattered modes when q = 0 shown in Figure 1 are greater than π/2, the value corresponding to Neumann boundary conditions.We now examine the effect of the scalar field charge on the energy flux for superradiantly scattered modes.
In Figure 2(a-b) we set the scalar field charge q = 1 and consider a selection of values of the black hole charge parameter Q, again for fixed rotation parameter Ω = 0.6.We have qQ > 0 in Figure 2(a) and qQ < 0 in Figure 2(b).For fixed Q, superradiantly scattered modes exist only in a narrow interval of values of the Robin parameter ζ , with the width of this interval decreasing as |Q| increases.The widths of these intervals are slightly greater when qQ < 0 compared to qQ > 0. As Q varies, the possible values of ζ for which there are superradiantly scattered modes is much broader than in the case q = 0, and we find superradiantly scattered modes with boundary conditions close to Dirichlet (ζ = 0) when |Q| is large.The magnitude of the energy flux F E /F for the superradiantly scattered modes with q = 1 in Figure 2(a-b) is roughly an order of magnitude greater than those in Figure 1 for q = 0, indicating a significant enhancement in superradiant scattering due to the scalar field charge.The maximum magnitude of the energy flux is slightly larger for superradiantly scattered modes with qQ < 0 than for the same values of q and |Q| but with qQ > 0. Therefore there is a slight additional enhancement in superradiant scattering Energy Flux (Q=-3, q=1) The black hole charge parameter is fixed to be Q = ±3, and the scalar field charge is q = 1.A selection of values of the rotation parameter Ω are considered.
when the scalar field charge q has the same sign as the black hole charge Q (10) compared to the case where q has the opposite sign to Q.The values of the Robin parameter ζ for the superradiantly scattered modes in Figure 2(a-b) mostly lie between the Dirichlet value ζ = 0 and that for Neumann boundary conditions ζ = π 2 .In Figure 2(c-d), with the rotation parameter again set to be Ω = 0.6, we fix the black hole charge parameter Q = ±2 and vary the scalar field charge q.In Figure 2(c-d), the interval of values of the Robin parameter ζ for which superradiantly scattered modes exists shrinks as the scalar field charge q increases, and moves to smaller values of ζ .The widths of these intervals are slightly greater for Q = −2 than for Q = 2 with the same value of q.At the same time, the maximum magnitude of the energy flux for superradiantly scattered modes F E /F in Figure 2(c-d) increases as q increases for both positive and negative Q.The maximum magnitude of the energy flux is slightly larger for Q = −2 than it is for Q = 2 with the same value of q.Combining the results in Figures 2, we provide evidence that increasing the magnitudes of either the black hole or scalar field charges gives a narrower interval of values of ζ yielding superradiantly scattered modes, with that interval being closer to Dirichlet boundary conditions.The maximum magnitude of the energy flux for superradiantly scattered modes shown in Figure 2 generally increases as either q or |Q| increases.
So far we have studied superradiantly scattered modes with the rotation parameter Ω fixed.In Figure 3 we fix the black hole charge parameter Q = ±3 and scalar field charge q = 1, and consider a selection of values of Ω .In Figure 3, increasing the rotation parameter results in large increases in both the width of the interval of values of ζ for which there are superradiantly scattered modes, and the maximum magnitude of the energy flux for superradiantly scattered modes.These effects in Figure 3 are significantly larger than those in Figure 2 resulting from changing either the scalar field or black hole charges.We therefore provide evidence that the most important factor influencing superradiant scattering is the rotation of the black hole.
Comparing the results for Q = ±3 in Figure 3, there is again a slight enhancement in the maximum magnitude of the energy flux when Q < 0 compared to Q > 0. The difference in the values of the Robin parameter ζ for superradiantly scattered modes when Q = 3 compared to Q = −3 is also striking.
We close our discussion of the energy flux by exploring, in Figure 4, some results for a large value of the black hole charge parameter, namely |Q| = 5, again with the rotation parameter Ω = 0.6.In Figure 4 we find behaviour which is qualitatively different from that shown in Figure 2.
Consider first Q = 5 and the superradiantly scattered modes shown in Figure 4(a).For larger values of the scalar field charge q ≥ 2.2, in Figure 4 Energy Flux (Q=5) q = 1.8 q = 1.9 q = 2.0 q = 2.1 q = 2.2 q = 2.5 Energy Flux (Q= -5) q = 1.2 q = 1.3 q = 1.5 q = 1.7 q = 1.8 q = 1.9 ular, for q = 2.0, 2.1 and 2.2 we show, in Figure 4(a), superradiantly scattered modes for which ζ = π 2 , corresponding to Neumann boundary conditions.For 1.8 < q < 2.2, Figure 4(a) shows that the value of the Robin parameter ζ at which the energy flux has its maximum magnitude shifts from a location close to ζ = 0 to a location close to ζ = π.On decreasing q further, for fixed q we find two "branches" of superradiantly scattered modes, one in a neighbourhood of ζ = π and one in a neighbourhood of ζ = 0.These superradiantly scattered modes are depicted in Figure 4(b), where we have chosen a branch of the arctan function in (56) for which − π 2 < ζ < π 2 instead of 0 < ζ < π as used elsewhere.In Figure 4(b), we can see that for q = 1.5 and 1.7, there are superradiantly scattered modes for which ζ = 0, corresponding to Dirichlet boundary conditions.Superradiantly scattered charged scalar field modes for Q = −5 are shown in Figure 4(c-d).Again we find that for large values of the scalar field charge q, superradiantly scattered modes lie in a narrow interval of values of the Robin parameter ζ , and these values are close to the Dirichlet value ζ = 0. Decreasing q again results in superradiantly scattered modes in a much larger interval of values of ζ , but the width of the interval of values of ζ shrinks again at a larger value of q than for the Q = 5 case.In Figure 4(c), we find superradiantly scattered modes for which ζ = π 2 when q = 2.2.As q decreases further, as for Q = 5 we also find two "branches" of superradiantly scattered modes, which are shown in Figure 4(d), where again we use a branch of the arctan function in (56) for which − π 2 < ζ < π 2 .The value of q below which we have two "branches" is higher for Q = −5 than it is for Q = 5.From Figure 4(d), we find superradiantly scattered modes for which ζ = 0 (corresponding to Dirichlet boundary conditions) when q = 1.8 and 1.9.In Figure 4, it can also be seen that the maximum value of the energy flux is slightly larger for superradiantly scattered modes with Q = −5 than for modes with the same value of the scalar charge q but with Q = 5.Comparing Figure 4(b) and (d), it can be seen that the range of values of the Robin parameter We therefore provide evidence that, unlike the situation for a neutral scalar field, for a charged scalar field on a charged BTZ black hole background, at least for a small subset of the (Q, q)-parameter space, there are superradiantly scattered modes with complex frequencies satisfying either Dirichlet or Neumann boundary conditions at infinity.

Mode frequencies
Having provided evidence for superradiant scattering of a charged scalar field on a charged, rotating, BTZ black hole, we now consider the frequencies of the superradiant modes in order to determine whether or not a superradiant instability arises.
We begin, in Figure 5, by considering the frequencies of neutral scalar field modes.The curves for a neutral black hole with Q = 0 match those in [56], and those for Q ̸ = 0 have qualitatively similar properties (recall that for q = 0 the radial equation (25) does not depend on the sign of Q).The real parts of the frequencies decrease as the Robin parameter ζ increases.It can be seen that the values of the frequencies do not change very much as the black hole charge Q varies, but the Robin parameter does vary more significantly.This is in agreement with the energy fluxes shown in Figure 1.The imaginary parts of the frequencies increase as the Robin parameter increases.All the superradiantly scattered neutral scalar field modes we consider here have positive imaginary part, and are therefore unstable, as in [56] for Q = 0.The range of values of the imaginary parts of the frequencies is not significantly changed by including a black hole charge Q, and therefore the timescale for the superradiant instability is also not significantly affected.
We now explore the frequencies of superradiantly scattered charged scalar field modes.In Figures 6-7 we show the real and imaginary parts of the frequencies of the superradiantly scattered modes whose energy fluxes are depicted in Figure 4.The features of the frequencies for the remaining superradiantly scattered modes shown in Figures 2-3 are qualitatively similar so we relegate these to Appendix B.
Consider first qQ > 0, so that the scalar field charge q has the opposite sign to the black hole charge Q (10).The frequencies in this case are shown in Figure 6.In contrast to the neutral scalar case, the real parts of the frequencies (left-hand plots in Figure 6) increase as the Robin parameter ζ increases, for all values of the black hole charge parameter Q > 0, scalar field charge q and rotation parameter Ω considered.
When qQ > 0, the imaginary parts of the frequencies (right-hand plots in Figure 6) are negative for all parameter values examined, indicating that these modes are stable.This is a surprising result, since, for example, superradiantly scattered charged scalar field modes on a four-dimensional Reissner-Nordström-adS black hole are unstable [44].We may understand this result heuristically as follows.
From (50,53), an ingoing charged scalar field mode having q > 0 and which is superradiantly scattered gives outgoing fluxes of both energy and charge at the event horizon.A superradiantly scattered charged scalar field mode is therefore extracting both energy and charge from the black hole.At infinity, as in [44], the fluxes of both energy and charge vanish for all Robin boundary conditions.When qQ > 0, the black hole charge Q (10) has the opposite sign to the scalar field charge q.A scalar field mode with q > 0, which is extracting charge from the black hole, therefore results in an increase in the magnitude of the black hole charge | Q| (since Q < 0 when Q > 0).Since the magnitude of the black hole charge is increasing, the energy in the electromagnetic field is also increasing as a result of the superradiant scattering.In this case the energy extracted from the black hole by the superradiantly scattered charged scalar field mode is absorbed by the electromagnetic field rather than the scalar field.Consequently, the scalar field does not grow with time and the superradiantly scattered mode is stable.
Turning now to qQ < 0, the frequencies for superradiantly scattered modes in this case are shown in Figure 7.For the values of the scalar field charge q > 0, black hole charge parameter Q < 0 and rotation parameter Ω considered, the real parts of the frequencies (left-hand plots in Figure 7) decrease as the Robin parameter ζ increases, as is the case for the neutral scalar field mode frequencies shown in Figure 5.
When qQ < 0, the imaginary parts of the frequencies (right-hand plots in Figure 7) are positive for all parameter values studied, indicating that these superradiantly scattered modes are unstable, in contrast to the superradiantly scattered modes for qQ > 0. A heuristic argument for the instability of the superradiantly scattered charged scalar field modes with qQ < 0 is as follows.
As discussed above, a superradiantly scattered ingoing charged scalar field mode with q > 0 gives outgoing fluxes of both energy and charge, and therefore is extracting both energy and charge from the black hole.When qQ < 0, the black hole charge Q (10) has the same sign as the scalar field Real Part (q=0) Imaginary Part (q=0) charge q.Therefore a scalar field mode with q > 0 which is extracting charge from the black hole results in a decrease in the magnitude of the black hole charge | Q|, and consequently a decrease in the energy in the electromagnetic field.As a result, the energy extracted from the black hole is absorbed by the scalar field.Since the boundary conditions at infinity are reflecting, the scalar field wave is repeatedly superradiantly scattered off the black hole, resulting in an instability.For all values of the scalar field charge q, black hole charge parameter Q and rotation parameter Ω studied in Figures 6-7, we find that the imaginary part of the frequency of superradiantly scattered charged scalar field modes increases as the Robin parameter ζ increases.Therefore the stable superradiantly scattered modes for qQ > 0 decay more slowly in time as ζ increases, while the unstable superradiantly scattered modes for qQ < 0 are more rapidly growing with time as ζ increases.

Conclusions
In this paper we have explored the effect of black hole and scalar field charge on superradiant scattering on three-dimensional BTZ black holes.We sought to compare superradiant scattering of a charged scalar field on a charged BTZ black hole with that on a four-dimensional Reissner-Nordström black hole.In particular, we examined whether charge superradiance is present on a nonrotating charged BTZ black hole.Superradiant amplification due to charge superradiant scattering on a Reissner-Nordström black hole is up to two orders of magnitude greater than that due to superradiant scattering on a rotating Kerr black hole, and we investigated whether similar enhancement also exists on BTZ black holes.
We considered separable mode solutions of the charged scalar field equation on the charged generalization of the rotating BTZ black hole metric [63].Working in the frequency domain, we find, as in the neutral scalar field case, that modes with real frequency do not exhibit superradiant scattering.For modes with complex frequencies, following [5,56], we define superradiant scattering as occurring if the ingoing flux of energy due to an ingoing scalar field mode is negative (in other words, if an ingoing mode results in an outgoing flux of energy it is said to be superradiantly scattered).We find that it is necessary for the black hole to be rotating in order for superradiant scattering to occur.Therefore, there is no superradiant scattering for nonrotating BTZ black holes, unlike the situation for four-dimensional, Reissner-Nordström black holes.Superradiantly scattered modes can only lie in a region of the complex frequency plane satisfying the inequality (52).We use a simple numerical method, applicable to modes in the superradiant scattering regime, to provide evidence for the existence of superradiantly scattered charged scalar field modes when the black hole charge is nonzero.We have not performed an exhaustive search of the parameter space, but instead considered a sample of black holes.For the region of parameter space explored, the presence of black hole and scalar field charges results in a flux of outgoing energy which is about an order of magnitude larger than in the uncharged case.However, the dominant parameter affecting Real Part (Q= -5) q = 2.0 q = 2.1 q = 2.2 q = 2.5 (a) Imaginary Part (Q= -5) q = 2.0 q = 2.1 q = 2.2 q = 2.5 Real Part (Q= -5) q = 1.2 q = 1.3 q = 1.5 q = 1.7 q = 1.8 q = 1.9 Imaginary Part (Q= -5) q = 1.2 q = 1.3 q = 1.5 q = 1.7 q = 1.8 q = 1.9 Q > 0 studied, are stable.For both signs of qQ, the superradiantly scattered modes are ingoing at the event horizon but extract energy from the black hole.When qQ < 0, superradiantly scattered modes decrease the magnitude of the black hole charge and thereby the energy in the electromagnetic field also decreases.This means that the energy extracted from the black hole is absorbed by the scalar field, resulting in an instability.In contrast, for qQ > 0, superradiantly scattered modes increase the magnitude of the black hole charge and thereby the energy in the electromagnetic field also increases.This means that the energy extracted from the black hole is absorbed by the electromagnetic field rather than the scalar field with no consequent scalar field instability.
Our numerical method has limited us to exploring a comparatively small region of the parameter space.In particular, we find reliable numerical results only when at least one of the scalar field charge q or black hole charge parameter |Q| is comparatively large.We have also fixed the black hole mass parameter M and azimuthal quantum number k, as well as the scalar field mass m.Furthermore, we have restricted our attention to a charged scalar field minimally coupled to the spacetime curvature.With a more sophisticated numerical method (such as that employed in [84,85]), it would be interesting to probe the parameter space more widely.In particular, we conjecture that there will be other superradiantly scattered charged scalar field modes that are unstable, as well as other unstable (but nonsuperradiantly scattered) modes.Testing this conjecture requires an alternative numerical method and hence is an avenue for future work.
In this paper we have studied a classical charged scalar field.A natural application of our work would be to consider a quantum charged scalar field.The study of a massless, conformally coupled quantum scalar field on a neutral BTZ black hole is comparatively straightforward due to the construction of the BTZ metric by identifying points in adS space-time [48,49].In particular, when either Dirichlet or Neumann boundary conditions are applied, the maximal symmetry of adS can be exploited to enable the computation of the renormalized expectation value of the stress-energy tensor using the method of images [87,88], see also [89][90][91].This method is not applicable when Robin boundary conditions are applied as these break the maximal symmetry of the underlying adS geometry [92][93][94].The ground state Green's function for a neutral scalar field with Robin boundary conditions applied is constructed in [65] using a mode sum decomposition.
It would be interesting to explore what effect the superradiantly scattered modes we have found in this paper have on the definition of quantum states for a charged scalar field on a charged BTZ black hole.On four-dimensional asymptotically flat black holes, the presence of superradiantly scattered modes introduces subtleties in the construction of quantum states, both in the rotating [95][96][97][98][99] and charged scenarios [100], and one might anticipate similar challenges on a BTZ black hole.Neutral scalar field modes on a neutral BTZ black hole are given by hypergeometric functions which simplifies the analysis [56,65].For a charged scalar field on a charged BTZ background there appears to be no simple closed-form expression for the modes, which will complicate the construction.We therefore postpone further consideration of the quantum charged scalar field to future work.BY) licence.This version of the article has been accepted for publication, after peer review but is not the Version of Record and does not reflect post-acceptance improvements, or any corrections.The Version of Record is available online at: http://dx.doi.org/10.1140/epjc/s10052-024-12910-7.

A Numerical method
Our goal is to find complex frequencies ω such that the corresponding solution X ωℓ (r) of the radial equation ( 25) satisfies ingoing boundary conditions (48) at the horizon and Robin boundary conditions (35) at infinity.There are two key components of our method, integrating the radial equation for a given complex frequency ω, and a root-finding algorithm to locate frequencies for which the boundary conditions at infinity are satisfied.
Given a complex frequency ω, we impose ingoing boundary conditions (48) on the radial function X ωk (r) at r = r h + ε, where ε ≪ 1.For r ≫ r h , the function X ωk (r) takes the form (34) for complex constants C ωk , D ωℓ .We rewrite the radial equation (25) in terms of the radial coordinate r and a new dependent variable Y ωk (r) = r X ωk (r).We numerically integrate this new radial equation from r = r h + ε to r = r max , where r max ≫ r h , using MATHEMATICA's built-in NDSolve command.We have tested the numerical integration by using a sample of both stable and unstable QNM frequencies for a charged scalar field on nonrotating charged BTZ black hole given in [85].The value of D ωk can be found as the limit of Y ωk (r) as r → ∞.The value of C ωk is found from the limit of r 1+ √ 1+4µ 2 Y ′ ωk (r) as r → ∞.We require the numerical integration of the radial equation (25) to very high precision in order to extract the constants C ωk and D ωk to a reasonable accuracy.
For a general frequency ω, the constants C ωk and D ωk thus found will be complex.To apply Robin boundary conditions (35), we require the ratio D ωk /C ωk to be real.For fixed ℜ(ω) in the interval 0 < ℜ(ω) < kΩ /ℓ for which superradiant scattering is possible, we use MATHEMATICA's inbuilt root-finding command FindRoot to find the value of ℑ(ω) for which the imaginary part of D ωk /C ωk vanishes.
From [56], the spectrum of QNM frequencies for the neutral scalar field on a neutral BTZ black hole is very complicated, with many modes located very close together in (ω, ζ ) space (where ζ is the parameter governing the Robin boundary conditions ( 56)).We suspect that the same is true for a charged scalar field on a charged BTZ black hole.This presents challenges for the root-finding algorithm, which can jump from one branch of QNM frequencies to another as ℜ(ω) is varied.
When Q = 0 = q, the root-finding algorithm reproduces the results presented in [56] (see the Q = 0 curves in Figures 1 and 5).This is a highly nontrivial check, since we are directly integrating the radial equation (25) while the authors of [56] exploit the fact that the radial functions are given in terms of hypergeometric functions.
For q ̸ = 0, the root-finding algorithm gives robust results for sufficiently large q and |Q|, as presented in Section 4 and Appendix B. We have verified that changing the parameters ε and r max in the numerical integration procedure yields plots which are indistinguishable from those presented in Section 4 and Appendix B, and estimate the maximum relative error in the QNM frequencies to be of the order of O(10 −5 ).

Fig. 1 :
Fig. 1: Energy flux F E /F (50) for superradiantly scattered neutral scalar field modes as a function of the Robin parameter ζ (56).A negative energy flux corresponds to superradiant scattering.We have fixed m 2 = −0.65,ℓ = 1, M = 16, Ω = 0.6, q = 0 and k = 1.The values of the black hole charge parameter Q are as given in the legends.

Fig. 2 :
Fig. 2: Energy flux F E /F (50) for superradiantly scattered charged scalar field modes as a function of the Robin parameter ζ (56).A negative energy flux corresponds to superradiant scattering.We have fixed m 2 = −0.65,ℓ = 1, M = 16, Ω = 0.6 and k = 1.The values of the black hole charge parameter Q and scalar field charge q are as given in the legends.

Fig. 3 :
Fig. 3: Energy flux F E /F (50) for superradiantly scattered charged scalar field modes as a function of the Robin parameter ζ (56).A negative energy flux corresponds to superradiant scattering.We have set m 2 = −0.65,ℓ = 1, M = 16 and k = 1.The black hole charge parameter is fixed to be Q = ±3, and the scalar field charge is q = 1.A selection of values of the rotation parameter Ω are considered.
(a) we find a narrow interval of values of the Robin parameter ζ which yield superradiantly scattered modes, and furthermore these values of ζ lie close to the Dirichlet value ζ = 0, similarly to the results in Figure2(c) for Q = 2.However, as q decreases (again in Figure4(a)), the interval of values of ζ for which there are superradiantly scattered modes widens significantly, and comprises the majority of the interval 0 < ζ < π.In partic-

Fig. 4 :
Fig. 4: Energy flux F E /F (50) for superradiantly scattered charged scalar field modes as a function of the Robin parameter ζ (56).A negative energy flux corresponds to superradiant scattering.We have set m 2 = −0.65,ℓ = 1, M = 16, Ω = 0.6 and k = 1.The black hole charge parameter is fixed to be Q = ±5, and a selection of values of the scalar field charge q are considered.

Fig. 5 :
Fig. 5: Real (left) and imaginary (right) parts of the frequencies of neutral scalar field modes as a function of the Robin parameter ζ (56).We have fixed m 2 = −0.65,ℓ = 1, M = 16, Ω = 0.6, q = 0 and k = 1.The values of the black hole charge parameter Q are as given in the legends.

Fig. 7 :
Fig. 7: Real (left) and imaginary (right) parts of the frequencies of superradiantly scattered charged scalar field modes as a function of the Robin parameter ζ (56).We have fixed m 2 = −0.65,ℓ = 1, M = 16, Ω = 0.6 and k = 1.The black hole charge parameter is fixed to be Q = −5, and a selection of values of the scalar field charge q are considered.
as a function of the Robin parameter ζ (56).A negative energy flux corresponds to superradiant scattering.