Quasinormal modes of the four-dimensional black hole in Einstein-Weyl gravity

In this paper we consider quasi-normal modes of test scalar and electromagnetic fields in the background of a non-Schwarzschild black hole in the Einstein-Weyl gravity. We find the dependence of the modes on the new parameter $p$, which is related to the coupling constant $\alpha$, for different values of the angular parameter $l$ and show that when $l$ is increasing, the frequencies tend to linear dependence on $p$. We have also examined the case of a massive scalar field for which we have shown that the arbitrarily long lived modes (called quasi-resonances) exist in the spectrum.


I. INTRODUCTION
As is known, the general theory of relativity leaves open some fundamental questions, such as: the singularity problem, construction of non-contradictory quantum gravity, the problems of dark matter and dark energy. The search for answers to these questions led to the emergence of new alternative theories. Some of these alternative theories are inspired by string theory, which suggests its ways of resolution of some of the above problems. One of the most popular approaches is related to adding higher than the first order in curvature terms to the Einsteinin action. If no other matter fields are included for four-dimensional space-time D = 4, the most general theory up the second order in curvature is described by the following lagrangian: here g is determinant of the matrix of a metric tensor; α, β and γ are constants; C µνρσ is the Weyl tensor. This theory has been recently considered in a number of papers [1], [2]. On the other hand, in the last time interest in studying perturbations of black holes is growing [3], [4], [5]. This was greatly contributed by observations of LIGO and VIRGO collaborations of the gravitational-wave signal at the merger of two black holes [6]. Although the results of observations are in concordance with General Relativity, a significant uncertainty of tens of percents in our knowledge of angular momentum and mass of the newly formed black hole leaves the window for alternative theories widely open [7], [8], [9]. This enhanced consideration of various alternative theories of gravity, including theories with higher curvature corrections.
Perturbations of black holes in theories with higher curvature corrections have been recently studied in a number of papers when considering higher dimensional black holes in the context of gauge/gravity duality and higher dimensional gravity [10], in the four-dimensional Einstein-dilaton-Gauss-Bonnet theory [11], in the Chern-Simons gravity [12] and in the quadratic [2], [13], [14], [15]. The latter theory will be considered in our paper, because apart from its generality, the latter quadratic theory is interesting also because, unlike many other theories of gravity, it allows for two asymptotically flat solutions representing a black hole: one is the usual Schwarzschild solution and the other is essentially non-Schwazrschild solution. The non-Schwarzschild solution has a remarkable property: when the coupling constant approaches its extremal value, the black hole approaches the massless state, having a non-zero radius at the same time. Thus, near extremal alternatives to the Schwarzschild solution are a kind of "light black holes" which are qualitatively different from its Einsteinian analogues.
The eikonal formula for quasi-normal modes of test fields was obtained for the black hole in this theory in [16]. However, no analysis of the fundamental quasi-normal modes was performed there. At the same time, the analysis of the fundamental mode suggested in [13] only for a massless test field suffers, in our opinion, from numerical inaccuracy. In this paper we investigate the quasi-normal spectrum of a scalar and a electromagnetic test fields and find the analytic expressions for the quasi-normal modes for different values of the multipole number l. We also consider the case of a massive scalar field. The interest of this case is connected with the possibility of a existence of undamped oscillations, i. e. quasi-resonance, for sufficiently large masses of the field. First this phenomenon was observed for the charged field of the Reissner-Nordstrom metric in the paper [17] and further studied for Schwarzschild and Kerr cases in [18], [19]. Moreover, recently in [20] it was shown that the quasi-resonances exist even at non-minimal coupling of the field, that is, they are not artifacts of the minimal coupling. Here we have shown the existence of a quasi-resonance for the massive scalar field for the non-Schwarzschild black hole solution in the quadratic gravity.
The structure of this work is as follows. In Sec. II we consider the basic concepts, introduce a metric and a wave equation for a spherically symmetric black hole. In Sec. III A, the basic principles of the WKB method are briefly considered, quasi-normal modes for test scalar and electromagnetic fields of massless (Sec. III B) and massive (Sec. III C) are found, a comparative analysis is made for the result obtained for the scalar field with the results in other studies. The possibility of the existence of quasi-resonance is also shown.

II. BLACK HOLE METRIC AND ANALYTICS FOR THE WAVE EQUATION
The metric for a spherically symmetric black hole was found in [16] and in general case is: The metric functions A(r) and B(r) were obtained numerically in [14], [15] and an analytical approximation was found for this numerical solution in [16].
We return to equation (1). As mentioned above, in the framework of the quadratic gravity, we should have two analytical solutions: the usual Schwarzschild solution and a solution with the quadratic correction. As shown in [14], [15], we can put β = 0. We also set γ = 1. The remaining undefined constant α we rewrite through the dimensionless parameter p, which parameterizes the solutions up to the rescaling and it can be determined through α as follows: where r 0 is a black hole radius. The Schwarzschild metric is valid for all p, however the second solution for an asymptotically flat black hole appears at some intervals of the values of this parameter, for small values [14], [15]. The approximate minimum and maximum values of p are: We will consider quasi-normal modes of test fields, scalar and electromagnetic, in the background of the black hole. Their perturbations can be represented in terms of the general relativistic Klein-Gordon equation: for a massive scalar field equation (5) has the form: for an electromagnetic field the general covariant Maxwell equations read as follows: Here A µ is a vector potential and The function Φ can be represented in terms of spherical harmonics and radial part as follows: Here R ωl (r) is a radial part, Y l (θ, φ) is a usual spherical harmonics, l is a angular harmonic index. For convenience, we rewrite the radial component R ωl (r) through the new function: Thus, taking into account the replacement (8), equations (6) and (7) in terms of the "tortoise coordinate" r * can be represented in the general form: III. QUASI-NORMAL MODES FOR TEST FIELDS

A. The WKB method
The solution of equation (10) can be found with the help of the WKB approximation [21]. This approach has application in quantum mechanics for finding the first terms in the expansion of the wave function, provided that the amplitude or frequency of the oscillation changes smoothly. When solving the wave equation, strictly fixed boundary conditions are imposed: for a functions Ψ there are only incoming waves at the horizon (r * → −∞) and only the outgoing waves at the infinity (r * → +∞). Finally, we obtain a complex number for the frequency ω, in which the real part is the real oscillation frequency, and imaginary part is the damping rate of oscillations in terms of the parameters of the black hole. It is known that the WKB series converges only asymptotically. Nevertheless, it is observed in many papers [22], [23], [24] that higher WKB orders usually produce more accurate results. In our work we will use the sixth-order wkb-corrections obtained in [24]: Here, the Λ i are the correction term of the i-th order and λ depends on the value of the potentials V and its derivative at the maximum and are explicitly represented in [22], [24]; n = 0, 1, 2, .. is a overtone number. It is sufficient to use the 6th order WKB formula for small n, while n ≤ l. Here we are limited only by the fundamental mode n = 0, because the contribution of higher overtones in the resultant signal is negligible [5]. In addition, since gravitational waves are very weak, only the fundamental mode of the spectrum can be extracted from a signal with sufficient accuracy.

B. Quasi-normal modes for massless fields
We shall consider a test scalar and electromagnetic fields. The effective potentials for the general metric (2) can be written in the forms the following forms [25]: V e (r) = l(l + 1)A(r) r 2 .
We find the quasi-normal modes for given potentials with the help of the sixth order WKB formula [25] with the values of the multipole number l = 0, 1, 2 for the scalar field (s = 0) and l = 1, 2, 3 for the electromagnetic field (s = 1). The results are presented in table III. From the values obtained, we construct the dependencies of the real and, separately, imaginary parts of the frequency ω on the values of the parameter p for various multipole numbers l. As can be seen in figures 1 and 2, the dependence ω(p) tends to be linear. Thus, it is possible to obtain approximate laws of variation of the imaginary and real parts as functions of p for different values of the angular number l = 0, l = 1, l = 2. From the above data one can see that the non-Schwarzschild branch of the Einstein-Weyl theory is always characterized by higher oscillation frequency and quicker damping of the fundamental mode.
Accordingly, for the electromagnetic field (FIG. 2) l = 1, l = 2, l = 3, we have: In formulas (13) and (14), the linear approximation for l = s = 0 for the scalar and l = s = 1 for the electromagnetic field, is carried out for values of p which are larger than 1.04443. It is possible that we do not observe the strict linear behavior at small l and relatively small values of p, because the WKB methods is less accurate in this regime. This is also clearly visible in figures 1a and 2a.  1 have a significant difference from the data provided in the work [13]. The reason of this discrepancy between our results and the data presented in [13] is probably related to potentially insufficient accuracy when integrating the black hole solution: the time domain integration [26] used in [13] requires very high precision for discrete data for A(r) and B(r). In our case, it is examined an analytical metric, which makes it possible to improve the accuracy of the calculations. The accuracy of our calculations is confirmed by the smoothness of ω as a function of p and  The dependence of the real and imaginary parts of the quasinormal frequencies on the dimensionless parameter p for the electromagnetic field (s = 1), l = 3.
reproducing the Schwarzschild limit when p → p min . Indeed, for example, for p min (given by [4a]), l = 2 and r 0 = 2 we have ω = 0.483642 − 0.096766i which is very close to the Schwarzschild limit ω = 0.4836 − 0.0968i [24]. In our case, the function ω(p) tends to be linear, while in the case with the data obtained in [13], the dependence of ω on p is not smooth (see fig. 1), which clearly indicates an unacceptably large numerical error in their calculations. In addition, the data presented in [13] does not allow one to see that the Schwaraschild limit can be reproduced as the mode with the smallest p is rather far from the Schwarzschild mode. The correctness of our results is also indicated by "the convergence" of the obtained values for the frequency at higher WKB-orders. That is demonstrated in table I and in table II.

C. Frequency dependence of quasinormal modes with the mass-term
Let us consider the effective potential for a massive scalar field: We consider how the behavior of quasi-normal modes changes with the increase of the mass-term. Let us fix p = 1 and l = 20, so that the spectrum is in its eikonal regime and can be described by the WKB approximation very well. We find values for the quasi-normal frequencies of this effective potential. Recently it has been shown that strictly speaking the WKB formula cannot describe the regime of quasi-resonances accurately [27] but can clearly show the tendency of the imaginary part to approach zero. Then, the extrapolation of the WKB data will indicate the existence of quasi-resonances. We choose, as an example, p = 1, l = 20. As can be seen in figure 3, the real oscillation frequency exponentially increases with increasing field mass. At the same time, the damping rate of the oscillations tends to zero smoothly. This suggests that mass modes are scattered more slowly than massless modes and, under certain conditions, for massive fields, the oscillations become undamped, which are quasi-resonances.

IV. CONCLUSIONS
We considered quasi-normal modes of the black hole in the quadratic gravity. The dependence of the modes on the new parameter p was determined for different values of the angular parameter l. From figures 1 and 2 it is obvious, that l is increasing, then the real and imaginary parts of ω tend to the linear dependence on p. It is found that the quasinormal modes of the non-Schwarzschild black hole solution are characterized by larger real oscillation frequency and larger damping rate than the Schwarzschild branch. In III C we added a massive term to the test scalar field. For such a field, for sufficiently large masses, the existence of undamped oscillations is shown, which are quasi-resonances discussed in the earlier literature for Schwarzschild [19], Reissner-Nordstrom [17] and Kerr solutions [18]. We have also shown that the previously published results of [13] for quasinormal modes of a massless scalar field are inaccurate.
It would be also interesting to consider perturbations of the gravitational field in the same black hole background. Using the WKB approach Hawking radiation for the above non-Schwarzschild black hole can be considered via solution of the classical scattering problem.   (11) with n = 0, r0 = 1 ; l = 0, 1, 2 for s = 0 and l = 1, 2, 3 for s = 1.