Greybody factors and quasinormal modes for a nonminimally coupled scalar field in a cloud of strings in (2+1)-dimensional background

We study the propagation of a probe massless nonminimally coupled scalar field in a fixed gravitational background of a cloud of strings in (2+1) dimensions. We obtain exact analytical expressions for the reflection coefficient, the absorption cross section, the decay rate as well as the quasinormal frequencies. The impact of the nonminimal coupling is investigated in detail. Our results show that Universality is not respected in general, and that scalar perturbations are stable.


Introduction
Black holes (BHs hereafter), a generic prediction of Einstein's General Relativity, are objects of paramount importance both for classical and quantum gravity. Greybody factors and quasinormal modes are two topics related to black hole physics of particular interest. On the one hand, Hawking radiation [1,2], since it is as a manifestation of a quantum effect in curved spacetime, has always attracted a lot of interest although it has not been detected in the Universe yet. The emitted particles feel an effective potential that back scatters part of the emitted radiation back into the black hole. The greybody factor is a frequency dependent quantity that measures the deviation from the original black body radiation spectrum, and provides us with valuable information about the black hole horizon structure [3]. The propagation and relativistic scattering of fields has been investigated both in asymptotically flat spacetimes and in background with a non-vanishing cosmological constant. For a partial list see e.g. [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23] and references therein.
On the other hand, LIGO historical direct detection of gravitational waves [24][25][26] from black hole mergers has opened a completely new window to our Universe, and allows us to test gravity and probe strong gravitational fields. Consequently, lately there is an increasing interest in black hole perturbations [27][28][29][30][31][32] and quasinormal modes of black holes, intimately related to the ring down phase after the formation of the distorded object a E-mail: arrincon@uc.cl b E-mail: grigorios.panotopoulos@tecnico.ulisboa.pt during the merging of two black holes. When a black hole is perturbed the geometry of spacetime undergoes dumped oscillations, which are characterized by the quasinormal modes with a non-vanishing imaginary part. Chandrasekhar'e monograph provides us with a comprehensive overview of black hole perturbations [33]. Quasinormal modes of black holes have been extensively studied, and for excellent reviews see e.g. [34][35][36].
The Bañados, Teitelboim and Zanelli (BTZ) black hole solution [37][38][39] in (1+2) dimensions marked the birth of the interest in lower-dimensional gravity. The absence of propagating degrees of freedom as well as its deep connection to the Chern-Simons term only [40][41][42] make threedimensional gravity special, and at the same time a framework which allow us to get insight into realistic black holes in four dimensions by studying a mathematically simpler three-dimensional system. The BTZ black hole is sourced by a negative cosmological constant, but other possibilities, such as scalar or electromagnetic fields [43][44][45][46], also exist. What is more, looking for a complete theory of quantum gravity, black hole solutions that admit scaledependent couplings, have been recently investigated. For an incomplete list see [47][48][49][50][51][52][53][54][55][56] and references therein. a massless canonical scalar field in this particular background were studied.
The relevance of a nonminimal coupling to gravity from a theoretical point of view is well-known. Quantum loop corrections will give rise to a nonminimal coupling, even if it is absent at tree level. Its role, however, in this particular context, to the best of our knowledge, has not been investigated yet. In the present article we extend the work of [64] regarding the propagation of a test massless scalar filed into a fixed gravitational background of a cloud of strings considering also a non-vanishing nonminimal coupling to gravity. Note, however, that here we assume Einstein's General Relativity, although in [64] the authors assumed a modified f (R) gravitational theory. Our work is organized as follows: In the next section we present the model and the equation for scalar perturbations, while in section three we solve the full radial equation exactly in terms of hypergeometric functions. The reflection coefficient, the greybody factor, the decay rate and the quasinormal spectrum are discussed in section four, and finally we summarize our work in the fifth section. We adopt the mostly positive metric signature (−, +, +), and we work in geometrical units where the universal constants are set to unity, = c = k B = 8πG N = 1.

Gravitational background and wave equation
We consider a model described by the action [60,61,64] where the gravitational part S G is given by the Einstein-Hilbert term with R being the Ricci scalar and g being the determinant of the metric tensor g µν , while the second contribution due to a cloud of strings is given by where λ 0 , λ 1 are the string parameters, and h = det(h AB ) with h AB being the string metric. Varying with respect to g µν we obtain Einstein's field equations with R µν being the Ricci tensor. Adopting the coordinate system x µ = (t, r, θ) the energy-momentum tensor corresponding to a cloud of strings takes the form [61] T µ ν = η r diag(1, 1, 0) with η being the coupling constant of the cloud of strings. Assuming a static and circularly symmetric line element of the form the metric function f (r) is found to be [60,64] f (r) = −M + 2ηr (7) with M being the mass of the black hole. In order to avoid a naked singularity the coupling constant η must be positive, and the horizon radius is computed to be r H = M/(2η). The metric function may be written down equivalently as where we trade the mass for the horizon radius. In the rest of the article we shall study the propagation of a probe massless nonminimally coupled scalar field in a fixed gravitational background of a cloud of strings in (2+1) dimensions.

Scalar perturbations: the wave equation
Let us consider a probe scalar field Φ with a non-vanishing nonminimal coupling to gravity ξ described by the action The corresponding wave equation is given by the classical Klein-Gordon equation which reads [16][17][18] where the nonminimal coupling is taken to be positive in order to maintain BH solutions (otherwise we will have naked singularity), and R 3 = −4η/r is the Ricci scalar of the cloud string background. Applying the usual separation of variables where m is the quantum number of angular momentum, we obtain an ordinary differential equation for the radial part or rewritten it explicitly we have In order to read-off the effective potential barrier that the probe scalar field feels, we define new variables as follows x = dr f (r) (15) where x is the so-called tortoise coordinate given by and recast the equation for the radial part into a Schrödingerlike equation of the form Therefore we obtain for the effective potential barrier the expression The effective potential as a function of the radial distance can be seen in Fig. (1) for three different values of the coupling ξ. The effective potential vanishes at r = r H , it reaches a maximum value V max , and tends to a constant V 0 when r → ∞ Since the effective barrier potential vanishes at the horizon, close to the horizon ω 2 V (x), and the solution for the Schrödinger-like equation is given by Requiring purely ingoing solution [3,7,20] we set A + = 0 in the following. Furthermore, when r → ∞ the wave equation admits plane wave solutions of the form ψ ∼ e ±iΩx when ω 2 − V 0 > 0, with Ω = √ ω 2 − V 0 being the modified frequency.

Solution of the radial equation
Given the metric function, we observe that in the far field region, r r H , it takes the simple form f (r) ∼ 2ηr. In addition, the angular momentum term is subdominant compared to the other two terms, and thus the differential equation for the radial part takes the form where the parameter k is defined according to Thus, the algebraic equation (which is Euler's equation) admits power-law solutions of the form R(r) ∼ r ρ . The algebraic equation for ρ satisfies The roots of the above equation are given by which are real when 1 − 4k ≥ 0 and complex when 4k > 1. In the following we shall consider the case of complex roots, and therefore the far field solution is given by where D − , D + are two arbitrary coefficients, while the roots are given by To find the solution of the full radial equation we introduce the dimensionless parameter z which is defined as follow and it takes values between 0 and 1. Then the new differential equation with respect to z becomes where the three constants are given by The last differential equation can be recast in the form of the Gauss' hypergeometric equation by removing the poles in the last term through the ansatz where now F satisfies the following differential equation and the new constants are given bȳ Demanding thatĀ = 0 =B, we determine the parameters α and β as follows and finally we obtain the hypergeometric equation with parameters a, b, c given by Note that the parameters a, b, c satisfy the condition c − a−b = 1−2β. Therefore the general solution for the radial part is given by [20] where C 1 , C 2 are two arbitrary coefficients, and the hypergeometric function can be expanded in a Taylor series as follows [65] F (a, b; c; z) = 1 + ab c z + · · · Setting C 2 = 0 and for the choice for α = −iω/(2η) we recover the purely ingoing solution close to the horizon, R ∼ (r − r H ) α , or ψ ∼ e −iωx . Therefore in the following we consider the first solution only, namely where now we have replaced C 1 by D. The sign in the expression for β does not really matter, and in the following we consider the plus sign.
In order to match with the far field solution obtained earlier (where now z → 1) we use the transformation [65] F and therefore the radial part as z → 1 reads Note that −β = ρ − and β − 1 = ρ + , and since z = 1 − (r H /r) the radial part R(r) for r r H can be written down as follows where we have introduced D − , D + in terms of D as follows 4 Greybody factors and quasinormal modes

Absorption cross section and decay rate
First we compute the reflection coefficient defined by and therefore it is computed by R = |D + /D − | 2 . Using the following identities for the Γ function [7] |Γ (iy)| 2 = π y sinh(πy) (54) we obtain the final expression R = cosh(πy 1 )cosh(πy 2 ) cosh(πy 3 )cosh(πy 4 ) where y i are given by the following expressions We see that one may consider two sectors in the problem separated by the critical value ξ c ≡ 1/8. Therefore, in the rest of the discussion we shall consider two separate cases, namely the case of weak nonminimal coupling, 0 ≤ ξ < ξ c , and strong coupling, ξ ≥ ξ c . The reflection coefficient as a function of the frequency can be seen in Fig. (2) for weak (left) and strong (right) nonminimal coupling, respectively. Notice that in the weak regime there is for ω a minimum allowed value that depend on the nonminimal coupling, The reflection coefficient is a monotonically deceasing function of the frequency, it starts at R ini = 1 and eventually it tends to zero. Furthermore, we see that when ξ increases the curves are shifted downwards.
Then, the absorption cross section is given by the simple formula [7,66] σ abs (ω) = 1 − R(ω) ω (61) which is the three-dimensional version of the optical theorem [3]. In [66] it was shown that for a generic spherically symmetric black hole the absorption cross section of a minimally coupled massless scalar field for vanishing angular momentum in the low energy regime tends to a constant that coincides with the area of the horizon. It is not obvious that this still holds for a massive or a nonminimally coupled scalar field. In fact in [6] it was found that Universality was respected under certain conditions, while in [16] it was shown that for a four-dimensional nonminimally coupled scalar field the absorption cross-section in the low energy regime tends to zero like ω 2 . In this work we find that as a function of ω for m = 0 the absorption cross-section starts from a constant and eventually goes to zero, but this constant does not necessarily coincide with the area of the horizon A H = 2πr H . We found very similar results in a previous work of ours [21]. We define the dimensionless parameter σ abs /A H and we plot it as a function of the frequency in Fig. (3). The constant increases with the coupling and finally acquires a limiting value when the coupling becomes sufficiently large.
In the weak coupling regime as well as in the strong coupling regime, for m > 0, the greybody factor reaches a maximum value, and then it tends to zero decreasing monotonically. What is more, the maximum value increases with ξ, which shifts the curves upwards. In the strong coupling regime for m = 0, however, the greybody factor is a monotonically decreasing function tending to zero starting from its maximum value at the origin. Similarly to the weak coupling regime, when ξ increases the curves are shifted upwards.
Since the spectrum emitted by the black hole (in three spatial dimensions) is given by [16] we define the decay rate of the black hole Γ decay by [7]: where the Hawking temperature of the cloud string black hole is given by T H = η/2π. The decay rate as a function of frequency can be seen in the Fig. (4) for low and large nonminimal coupling (left and right figures), respectively. Our figures are consistent with those of [7] in the high energy regime: the curves asymptotically go to zero. As before, when ξ increases the curves are shifted upwards. In the weak coupling regime the decay rate reaches a maximum value, while in the strong coupling regime the decay rate is a monotonically decreasing function going eventually to zero.

Quasinormal spectrum
First we remark that computing the QNMs of black holes analytically is possible only in some cases, see e.g. [67][68][69][70][71][72][73][74][75][76][77], see also [78] for the influence of the back reaction of the Hawking radiation upon black hole quasinormal modes. To obtain the quasinormal modes we apply the quasinormal boundary condition, according to which at infinity purely outgoing solution is required. Therefore we require that D − = 0 which is satisfied when the Gamma functions in the denominator have a pole where n = 0, 1, 2, ... is the overtone number. Given the time dependence of the scalar field, ∼ e −iωt , the mode is unstable (exponential growth) when ω I > 0 and stable (exponential decay) when ω I < 0. In the latter case the real part determines the frequency of the oscillation, ω R /(2π), while the inverse of |ω I | determines the dumping time, t −1 D = |ω I |. The real part (left panel) and the imaginary part (right panel) of the quasinormal frequencies as a function of the nonminimal coupling are shown in Fig. (5). We see that the slope decreases both with m and with n. The real part vanishes for m = 0, while for m > 0 it is initially positive, but eventually it becomes negative when ξ becomes sufficiently large. Since it is a monotonically decreasing function of the nonminimal coupling, the frequency of the oscillation decreases with ξ. In particular, the period of the oscillation in the case of a non-vanishing ξ is larger than the period corresponding to a canonical scalar field (ξ = 0).
Due to the emission of gravitational waves the spacetime undergoes dumbed oscillations, and this is encoded into the non-vanishing imaginary part. The latter is always negative, and therefore the scalar perturbations studied here are stable. Since the imaginary part is a monotonically decreasing function of the nonminimal coupling, a certain mode for given (m, n) decays faster as ξ increases.

Conclusion
In the present article we have studied the propagation of a test massless nonminimally coupled scalar field into a fixed gravitational background of a cloud of strings, extending a previous work where the probe field was a scalar field with a canonical kinetic term minimally coupled to gravity. We have obtained the expression for the effective potential, and we have solved the full radial equation exactly in terms of hypergeometric functions. We thus have obtained exact analytical expressions for the reflection coefficient, the absorption cross section as well as the decay rate. Finally, applying the quasinormal boundary condition we have obtained the expression for the quasinormal spectrum. The impact of the nonminimal coupling has been investigated in detail. Our results show that Universality is not respected in general, and that scalar perturbations are stable.