Maxwell quasinormal modes on a global monopole Schwarzschild-anti-de Sitter black hole with Robin boundary conditions

We generalize our previous studies on the Maxwell quasinormal modes around Schwarzschild-anti-de-Sitter black holes with Robin type vanishing energy flux boundary conditions, by adding a global monopole on the background. We first formulate the Maxwell equations both in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms and derive, based on the vanishing energy flux principle, two boundary conditions in each formalism. The Maxwell equations are then solved analytically in pure anti-de Sitter spacetimes with a global monopole, and two different normal modes are obtained due to the existence of the monopole parameter. In the small black hole and low frequency approximations, the Maxwell quasinormal modes are solved perturbatively on top of normal modes by using an asymptotic matching method, while beyond the aforementioned approximation, the Maxwell quasinormal modes are obtained numerically. We analyze the Maxwell quasinormal spectrum by varying the angular momentum quantum number $\ell$, the overtone number $N$, and in particular, the monopole parameter $8\pi\eta^2$. We show explicitly, through calculating quasinormal frequencies with both boundary conditions, that the global monopole produces the repulsive force.


I. INTRODUCTION
Black hole quasinormal modes (QNMs), describing the characteristic oscillations of black holes, have attracted a lot of attention recently, see for example reviews [1][2][3][4] and references therein. Due to the existence of an event horizon, the black hole spacetimes are intrinsically dissipative so that quasinormal frequencies are complex in general and the imaginary part is associated with the timescale of the perturbation. QNMs play vital roles on various aspects, ranging from gravitational wave astronomy [5,6] to the application in the context of the anti-de Sitter/conformal field theory (AdS/CFT) correspondence [7][8][9].
The AdS/CFT correspondence states that QNMs of a (D + 1)-dimensional asymptotically AdS black hole or brane are poles of the retarded Green's function in the dual conformal field theory in D dimensions at strong coupling. Horowitz and Hubeny first studied scalar QNMs on Schwarzschild-AdS black holes [10] (see also [11,12]), and numerous works were then followed to explore QNMs of various spin fields on asymptotically AdS black holes, see for example .
Mathematically QNMs are defined as eigenvalues of perturbation equations with physically relevant boundary conditions. Considering a lot of studies already performed in literatures, however, a generic boundary condition is still lacking. Recently, we have proposed the vanishing energy flux principle [38,39], which may be applied both to the Regge-Wheeler-Zerilli and to the Teukolsky formalisms, and leads to two sets of Robin * mjwang@hunnu.edu.cn. † panqiyuan@hunnu.edu.cn ‡ jljing@hunnu.edu.cn type boundary conditions and has been successfully employed to explore QNMs of the Maxwell [38,40] and Dirac fields [41,42]. In this paper, we follow the same rationale and generalize our previous studies of the Maxwell QNMs on Schwarzschild-AdS black holes, by adding a global monopole on the backgrounds. The global monopoles, as a special class of topological defects, may be formed in the early universe through the spontaneous symmetry breaking of the global O(3) symmetry to U(1) [43,44], according to the Grand Unified Theories. The gravitational properties of monopoles have been extensively studied, and an unusual property induced by global monopoles is that it possesses a solid deficit angle. This property makes black holes with a global monopole and without a global monopole topologically different, and thus leads to interesting physical consequences [45][46][47][48][49][50][51][52].
The purpose of this study is twofold. On one hand, we explore the impact of the global monopole on the Maxwell quasinormal spectrum on Schwarzschild-AdS black holes, by imposing vanishing energy flux boundary conditions. On the other hand, it is well known that, on spherically symmetric backgrounds, the Maxwell equations may be written either in the Regge-Wheeler-Zerilli or in the Teukolsky formalisms. As we argued before [38], by imposing vanishing energy flux boundary conditions, the Maxwell equations in both formalisms lead to the same quasinormal spectrum. Here we show explicitly, through calculating normal modes in both formalisms with vanishing energy flux boundary conditions, that it is indeed the case, even if a global monopole is included.
The structure of this paper is organized as follows. In Section II we introduce the Schwarzschild-AdS black holes with a global monopole, and show the Maxwell equations both in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms. In Section III we present the ex-plicit boundary conditions, based on the vanishing energy flux principle, for both the Regge-Wheeler-Zerilli variable and the Teukolsky variable of the Maxwell field. We then perform an analytic matching calculation for small AdS black holes in Section IV, and a numeric calculation in Section V. Final remarks and conclusions are presented in the last section.

II. BACKGROUND GEOMETRY AND THE FIELD EQUATIONS
In this section, we first briefly review the background geometry we shall study, i.e. Schwarzschild-AdS black holes with a global monopole, and then present equations of motion for the Maxwell fields on the aforementioned backgrounds both in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms.

A. The line element
We start by considering the following line element of a Schwarzschild-AdS black hole with a global monopole with the metric function where L is the AdS radius, M is the mass parameter.
Here the dimensionless parameterη 2 is defined bỹ where η is the global monopole parameter, and the Schwarzschild-AdS spacetimes may be recovered when 8πη 2 = 0. The Hawking temperature may be calculated, and one obtains where r + is the event horizon determined by the non-zero real root of ∆ r (r + ) = 0, and where the mass parameter has been expressed in terms of r + as By introducing the following coordinates transformationt and a new mass paramter Eq. (1) becomes Now it becomes clear that the global monopole introduces a solid deficit angle, so that the solid angle of the above spacetime is 4πη 2 .

B. Equations of motion in the Regge-Wheeler-Zerilli formalism
In a spherically symmetric background, one may obtain variable separated and decoupled Maxwell equations by using the Regge-Wheeler-Zerilli method [53,54]. For that purpose, we start from the Maxwell equations where the field strength tensor is defined as We then expand the vector potential A µ in terms of the scalar and vector spherical harmonics [55] A with the definition of the vector spherical harmonics where Y m are the scalar spherical harmonics, m is the azimuthal number, and is the angular momentum quantum number. Note that the first term in the right hand side of Eq. (8) has parity (−1) +1 while the second term has parity (−1) , and we shall call the former (latter) the axial (polar) modes. By substituting Eq. (8) into Eq. (7) with the assumption one obtains the Schrodinger-like radial wave equation where the tortoise coordinate is defined as with Ψ(r) = a m (r) for axial modes, and for polar modes.

C. Equations of motion in the Teukolsky formalism
Equations of motion of the Maxwell fields may be also derived within the Teukolsky formalism [56]. This approach is based on the Newmann-Penrose algorithm [57], and is particularly relevant to study linear perturbations of the massless spin fields on rotating black hole backgrounds. In this subsection we outline the radial equations, which may be obtained following the procedures presented in [58].

III. BOUNDARY CONDITIONS
In order to solve the radial equations, given by Eqs. (9) and (11), one has to impose physically relevant boundary conditions, both at the horizon and at infinity. At the horizon, we impose the commonly used ingoing wave boundary conditions. At infinity, we impose the vanishing energy flux principle, proposed in [38] (see also [39,59]), which have already been employed to study the Maxwell [38,40] and the Dirac [41,42] QNMs on asymptotically AdS spacetimes. Based on this principle, in the following we derive explicit boundary conditions for Eqs. (9) and (11), which are obtained in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms respectively, and we will show both equations with the corresponding boundary conditions lead to the same spectrum in the next section.

A. Boundary conditions in the Regge-Wheeler-Zerilli formalism
We start from the energy-momentum tensor of the Maxwell field, which is given by Then the spatial part of the radial energy flux may be calculated as where denotes the derivative with respect to r. By expanding Eq. (9) asymptotically as Eq. (13) becomes Then the vanishing energy flux principle, i.e. F| r,∞ = 0, leads to B. Boundary conditions in the Teukolsky formalism The explicit boundary conditions for the Teukolsky variables of the Maxwell fields on a global monopole Schwarzschild-AdS black hole can be derived directly, following the similar prescriptions described in [38,59]. Since the monopole parameter does not alter the asymptotic structure of AdS spacetimes, one may get exactly the same boundary conditions as to the Schwarzschild-AdS case, and the results are listed in the following.
To be specific, we focus on the boundary conditions for R −1 . From Eq. (11) one obtains the asymptotic behavior of R −1 as and the vanishing energy flux principle leads to [38,59] IV. ANALYTICS

A. Normal modes
The normal modes of the Maxwell fields on an empty AdS background with a global monopole are calculated analytically in this subsection, both in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms, by solving Eq. (9) with boundary conditions (15) and (16), and Eq. (11) with boundary conditions (18) and (19). These calculations provide a concrete example to show explicitly that, vanishing energy flux is a generic principle, which can be applied to both formalisms and leads to the same spectrum.

Normal modes in the Regge-Wheeler-Zerilli formalism
In a pure AdS spacetime with a global monopole (M = 0), the metric function becomes then the radial equation (9) can be solved, and one obtains Here c 1 , c 2 are two integration constants with dimension of inverse length, F (a, b, c, z) is the hypergeometric function, and where = 1, 2, 3, · · ·. By expanding Eq. (20) at large r, we get relations between c 1 and c 2 , i.e.
which corresponds to the first boundary condition given by Eq. (15), and which corresponds to the second boundary condition given by Eq. (16). Then by expanding Eq. (20) at small r we shall set c 2 = 0 to get a regular solution at the origin. This condition leads to two sets of normal modes from Eq. (22), and from Eq. (23), where N = 0, 1, 2, · · ·. The above two normal modes, by noticing that˜ is not an integer anymore, are different. This is an interesting observation, since for the case without a global monopole, the two sets of the Maxwell normal modes are isospectral up to one mode [38].

Normal modes in the Teukolsky formalism
In this case the radial Teukolsky equation (11) becomes The general solution for Eq. (27) is where F (a, b, c; z) is again the hypergeometric function, c 3 and c 4 are two integration constants with dimension of inverse length. These two constants are related to each other by the boundary conditions through expanding Eq. (28) at large r: • By imposing the first boundary condition given in Eq. (18), one gets a first relation between c 3 and c 4 • By imposing the second boundary condition given in Eq. (19), on the other hand, one gets a second relation between c 3 and c 4 where A 1 =(1 +˜ )F (˜ , 1 +˜ +ωL, 2 + 2˜ ; 2) +ωLF (1 +˜ , 1 +˜ +ωL, 2 + 2˜ ; 2) , Then from the small r behavior of Eq. (28) we have to set c 4 = 0 in order to get a regular solution of R −1 at the origin. This regularity condition picks the normal modes, from Eqs. (29) and (30): F (1 +˜ , 1 +˜ +ωL, 2 + 2˜ ; 2) = 0 , where again N = 0, 1, 2, ···, and two sets of normal modes are different. One may observe that normal modes obtained in the Teukolsky formalism, given in Eqs. (33) and (34), are exactly the same with the counterpart obtained in the Regge-Wheeler-Zerilli formalism, given in Eqs. (25) and (26), which indicates the equivalence of the two formalisms and the universality of the vanishing energy flux boundary conditions.

B. Quasinormal modes for small black holes
In this subsection, we perform an analytic calculation of quasinormal frequencies for the Maxwell fields on a small Schwarzschild-AdS black hole with a global monopole, by using an asymptotic matching method. Note that for this case the analytic calculation is only applicable to the Teukolsky formalism.

Near region
In the near region, and with small black hole approximation r + L , Eq. (11) becomes whereω is defined in Eq. (21). By defining a new dimensionless variable where R −1 ≡ R −1 (z), and˜ is defined in Eq. (21). The above equation can be solved in terms of the hypergeometric function where an ingoing boundary condition has been imposed. In order to match the far region solution, here we shall further expand the near region solution, given in Eq. (38), at large r. To do so, we take the z → 1 limit and use the property of the hypergeometric function [60], then obtain

Far region
In the far region, the black hole effects may be neglected, and the solution is given by Eq. (28). In order to match this solution with the near region solution, we shall expand Eq. (28) at small r, then obtain with where the constants c 3 and c 4 are related with each other by Eqs. (29) and (30), corresponding to the first and second boundary conditions.

Overlap region
In the overlap region the solutions, obtained in the near region given by Eq. (39) and in the far region given by Eq. (41), are the same up to a constant. Then one may impose the matching condition, R near −1,r R far −1,1/r = R far −1,r R near −1,1/r , which gives By imposing the first boundary condition and using the corresponding relation between c 3 and c 4 given by Eq. (29), one obtains while by imposing the second boundary condition and using the corresponding relation between c 3 and c 4 given by Eq. (30), one obtains where A 1 and A 2 are given by Eq. (31). For a small black hole (r + L ), at the leading order of r + /L, the left terms in Eqs. (44) and (45) vanish, and then we shall require the right terms in both equations to vanish as well. These conditions lead to two sets of normal modes, given by Eqs. (33) and (34). Then QNMs of small black holes may be obtained perturbatively by solving Eqs. (44) and (45), on top of normal modes. To achieve this goal, we expand the frequencỹ where j = 1, 2, andω j,N refer to normal modes. Here δ j is complex in general, and its real part, i.e. (δ j ), reflects damping rate of a black hole. The general expression of δ j , which is usually messy and lengthy but can be derived straightforwardly by substituting Eq. (46) into Eqs. (44) and (45).

V. NUMERIC REUSLTS
Beyond the regime where the asymptotic matching method is valid, one has to look for black hole quasinormal spectrum by resorting to numerics. In this part, we utilize a numeric pseudospectral method, adapted from our previous works [42,61], to solve the Maxwell equations given in the Regge-Wheeler-Zerilli formalism (9) with the corresponding boundary conditions given by Eqs. (15) and (16). 1 Before we introduce the pseudospectral method, here goes a few comments on the dimensionless form of Eq. (9) which is essential for numeric calculations. For the case we considered in this paper, one may either take the unit of L or take the unit ofL (with the definition given in Eq. (21)). For the former choice, Eq. (9) may be written as where r is an abbreviation of r L so that it is dimensionless, and where r + is a dimensionless event horizon. By taking the unit ofL, Eq. (9) becomes with where r is an abbreviation of dimensionless radial coordinate r L ,ω,L and˜ are given in Eq. (21).
By noticing that g(r) has the same form with the metric of Schwarzschild-AdS, so Eq. (49) is exactly the same with the Maxwell equation on Schwarzschild-AdS, by replacing with˜ . This is also the Maxwell equation one may obtain by starting from the metric given by Eq. (6). Therefore, in our numeric calculations, we take the unit ofL and setL = 1, and calculate the frequenciesω. As we have checked, in the unit ofL, the quasinormal frequencies have the uniform behaviors for various values of r + , and N , and which is consistent with the physical picture that the global monopole produces the repulsive force.
In order to employ a pseudospectral method conveniently, we first transform Eq. (9), which is a quadratic eigenvalue problem, into a linear eigenvalue problem, by where the tortoise coordinate r * is still defined in Eq. (10). Then changing the coordinate from r to z through which brings the integration domain from r ∈ [r + , ∞] to z ∈ [−1, +1], and discretizing the z coordinate according to the Chebyshev points where n denotes the number of grid points, Eq. (9) turns into an algebraic equation Here M 0 and M 1 are matrices, which may be constructed straightforwardly by discretizing Eq. (9) in terms of the Chebyshev points and Chebyshev differential matrices [62]. Boundary conditions associated to φ(z), may be derived from the transformation given by Eq. (51). At the horizon, since an ingoing wave boundary condition is satisfied automatically, we simply impose a regular boundary condition for φ(z). At infinity, from Eqs. (51) and (14), one obtains corresponding to the condition given in Eq. (15), and corresponding to the condition given in Eq. (16).
One should note that we useω 1 (ω 2 ) to represent the quasinormal frequency corresponding to the first (second) boundary conditions. A few selected data are presented below to demonstrate, in particular, the impact  of global monopole on the spectrum. Also note that in our numeric calculations we focus on black holes with size r + ≤ 1 since in this regime the monopole effects are more relevant. 2 In Fig. 1, we compare the analytic calculations with numeric data, by taking the angular momentum quantum number = 1, the overtone number N = 0 and the monopole parameter 8πη 2 = 0.05, and find a good agreement for small black holes.
A few numeric data are tabulated in Table. I. As one may observe, by taking 8πη 2 = 0.1 and N = 0, the real part of the Maxwell QNMs decreases while the magnitude of the imaginary part increases as the black hole size r + increases, similarly to the Schwarzschild-AdS case. In particular, the isospectrality of the modes for = 1 with the first boundary condition and = 2 with the second The impact of the angular momentum quantum number on the real and imaginary parts of quasinormal modes for the first (red) and second (blue) boundary conditions.
boundary condition is broken, due to the presence of the global monopole.
The effect of the angular momentum quantum number on the Maxwell quasinormal spectrum is presented In Fig. 2, for a black hole with size r + = 1, the global monopole 8πη 2 = 0.1 and with the overtone number N = 0. We observe, similarly to the Schwarzschild-AdS case (i.e. 8πη 2 = 0) reported in [38], that for both boundary conditions the real part of the Maxwell quasinormal frequencies increases while the magnitude of imaginary part decreases, as increases.
As the main goal of this paper, we explore the impact of the global monopole on the Maxwell spectrum in Fig. 3. As an illustrative example, here we take r + = 0.5, = 1, N = 0, and we observe that, for both boundary conditions, the real (the magnitude of imaginary) part of the Maxwell QNMs increases (decreases) as the global monopole 8πη 2 increases. As we have checked for various values of r + , and N , the above mentioned behaviors are held. This may be understood as follows. From Eq. (49), it shows clearly that the monopole parameter only appears in˜ and˜ plays the same role as . From Fig. 4, we observe˜ , by fixing , increases as the monopole parameter 8πη 2 increases, indicating the global monopole produces the repulsive force. This implies that, for larger monopole parameter, the perturbation fields (the Maxwell fields here) live longer around black holes, i.e. decay slower, exactly as shown in Fig. 3.
We have also studied the dependence of the Maxwell quasinormal frequencies on the overtone number N in Fig. 5. For this case, we take r + = 0.5 and = 1. It is shown, from the left and middle panels, that for two boundary conditions, both the real part and the magnitude of imaginary part of the Maxwell frequencies increase as N increases, and the excited modes for both (r + , N , ℓ) = (0.5, 0, 1) branches are approximately evenly spaced in N. In the right panel, we display the imaginary part in terms of the real part of the Maxwell QNMs. It shows interestingly that two branches of QNMs (for excited states) lie on the same line for different N. This phenomenon has also been observed for the Dirac case [42], and indicates that, although two branches of QNMs are different, they are similar in the sense that the excited modes of one branch may be interpolated from the other branch.

VI. DISCUSSION AND FINAL REMARKS
In this paper we have studied the Maxwell quasinormal spectrum on a global monopole Schwarzschild-AdS black hole, by imposing a generic Robin type boundary condition. To this end, we first presented the Maxwell equations both in the Regge-Wheeler-Zerilli and in the Teukolsky formalisms and derived the explicit boundary conditions for the Regge-Wheeler-Zerilli and the Teukol-sky variables, based on the vanishing energy flux principle. Then the Maxwell equations were solved in each formalism, both analytically and numerically.
In a pure AdS space with a global monopole, we have solved the Maxwell equations analytically in the aforementioned two formalisms. We found that two boundary conditions in each formalism lead to two different normal modes, due to the presence of the global monopole. This is very different with the Schwarzschild-AdS case where normal modes obtained from two boundary conditions are the same, up to one mode. In the small black hole and low frequency approximations, we also solved the Maxwell equations in the Teukolsky formalism by using an analytic matching method and we verified that the analytic calculations coincide with the numeric data well.
We then varied the black hole size r + , the angular momentum quantum number , and the overtone number N , in the presence of a global monopole; and analyzed their effects on the two sets of the Maxwell quasinormal spectrum in the numeric calculations. We observed that, the impact of r + , and N on the Maxwell QNMs are very similar to the Schwarzschild-AdS case. In partic-ular, we explored the monopole effects on the Maxwell spectrum, and we found that for both boundary conditions, the real part of the Maxwell spectrum increases while the magnitude of imaginary part decreases as the monopole parameter 8πη 2 increases. These trends are direct consequences of the fact that the global monopole produces the repulsive force.
Finally, we would like to stress that the above mentioned QNMs behaviors were obtained in the unit of L. One may alternatively use the unit of L, and as we have checked for this case the monopole effects on the Maxwell spectrum are more involved. In the former choice, the Maxwell equations on Schwarzschild-AdS black holes with a global monopole may be reformulated to the Maxwell equations without a global monopole but with the modified angular momentum quantum number , so that the repulsive nature of the global monopole becomes more clear.