Quasinormal modes of Gauss-Bonnet-AdS black holes: towards holographic description of finite coupling

Here we shall show that there is no other instability for the Einstein-Gauss-Bonnet-anti-de Sitter (AdS) black holes, than the eikonal one and consider the features of the quasinormal spectrum in the stability sector in detail. The obtained quasinormal spectrum consists from the two essentially different types of modes: perturbative and non-perturbative in the Gauss-Bonnet coupling $\alpha$. The sound and hydrodynamic modes of the perturbative branch can be expressed through their Schwazrschild-AdS limits by adding a linear in $\alpha$ correction to the damping rates: $\omega \approx Re(\omega_{SAdS}) - Im(\omega_{SAdS}) (1 - \alpha \cdot ((D+1) (D-4) /2 R^2)) i$, where $R$ is the AdS radius. The non-perturbative branch of modes consists of purely imaginary modes, whose damping rates unboundedly increase when $\alpha$ goes to zero. When the black hole radius is much larger than the anti-de Sitter radius $R$, the regime of the black hole with planar horizon (black brane) is reproduced. If the Gauss-Bonnet coupling $\alpha$ (or used in holography $\lambda_{GB}$) is not small enough, then the black holes and branes suffer from the instability, so that the holographic interpretation of perturbation of such black holes becomes questionable, as, for example, the claimed viscosity bound violation in the higher derivative gravity. For example, $D=5$ black brane is unstable at $|\lambda_{GB}|>1/8$ and has anomalously large relaxation time when approaching the threshold of instability.


Introduction
Quasinormal modes of asymptotically anti-de Sitter black holes play an essential role in the description of strongly coupled processes in the dual conformal field theory via the AdS/CFT correspondence [1]. The poles of the retarded Green functions in the strongly coupled Conformal Field Theory in D-dimensions coincide with the proper oscillation frequencies (called quasinormal modes) of large (in comparison with the AdS radius) (D + 1)dimensional asymptotically anti-de Sitter black holes [2]. The damping rate of the dominant quasinormal mode determines the timescale at which a black hole or quantum system in the dual field theory approaches equilibrium [3].
In 2005 Kovtun, Son, and Starinets showed [4] that the AdS/CFT correspondence predicts the following universal behavior of quark-gluon plasmas when modeled in the Conformal Field Theory through various gravitational backgrounds: where η is the shear viscosity, s is volume density of entropy. They also suggested that (1.1) determines the lower bound on the shear viscocity/entropy density ratio. Soon the theoretical prediction (1.1) was confirmed through comparisons with the data obtained when observing quark-gluon plasmas at RHIC [5].
However, the above results were obtained at the assumption that the dual field theory has large 't Hooft coupling λ [3]. The regime of weak coupling can be described by the kinetic theory, while the transition from weak to strong coupling, i. e., the regime of intermediate coupling is not yet understood. Recently, an approach allowing one to obtain corrections to the regime of infinite coupling has been suggested in [6,7]. There it is shown that higher curvature corrections to the gravitational action, such as the Gauss-Bonnet, Lovelock, R 4 or other terms, may give us a hint on what happens in the regime of intermediate coupling.
In [7,8,9] quasinormal modes of black branes in various theories with higher curvature corrections were investigated and several new interesting findings were reported. The main striking difference of the black-hole spectrum in presence of higher curvature corrections is appearance of a new, purely imaginary (i. e., non-oscillating), non-perturbative (in terms of the Gauss-Bonnet coupling constant α) modes. There was also observed a kind of breakdown of the hydrodynamic regime. It consists in the duplication of the shear, purely imaginary modes, which acquire non-vanishing real parts at some critical value of α [7,8]. These modes may interfere with the hydrodynamic modes in some range of parameters. However, in [10] it was shown that the effect of acquiring non-zero real part occurs already in the region of gravitational instability, so that it is not clear whether this effect may actually have some holographic meaning.
It is natural to learn whether the above new features are artifacts of the particular gravitational backgrounds or appropriate to a broader class of models or, possibly, even universal for a broad class of theories with higher curvature corrections. The straightforward task is to extend the above analysis from planar to spherically symmetric black holes. Although by now there are a lot of papers on quasinormal modes of black holes and branes in Gauss-Bonnet and Lovelock theories, quasinormal modes of gravitational perturbations of Gauss-Bonnet-AdS black holes have not been studied so far. Therefore, here we shall analyze quasinormal modes of spherically symmetric Einstein-Gauss-Bonnet-anti-de Sitter black holes in detail. In our previous paper [10] it was shown that Gauss-Bonnet-AdS black holes can be unstable and the exact parametric region of instability was found. This instability, counter-intuitively, develops at high multipole numbers ℓ and was called therefore eikonal instability. It is intrinsically related to the breakdown of the well-posedness of the initial value problem [11]. A similar instability was found for asymptotically flat [12,13] and de Sitter black holes [14] in Gauss-Bonnet theories, as well as for planar Einstein-Lovelock [15], spherical pure Lovelock (i.e. without the Einstein term) black holes [16], and for small charged Lovelock black holes [17]. The full analysis of the eikonal instability of black holes and branes in the generic Lovelock theory has been recently performed in [18]. However, test fields in the background of the Einstein-Gauss-Bonnet black hole do not show such an instability [19,20]. At the same time, when the non-zero positive Λ-term is turned on, there is another, non-eikonal gravitational instability [14] of Gauss-Bonnet black holes, which may be similar in nature with the gravitational instability of higher dimensional Ressiner-Nordström-de Sitter black holes in Einstein gravity [21]. The search for eikonal instability can be performed by analysis of the dominant (at high multipole numbers ℓ) term of the effective potential only, so that other types of instability at lower ℓ cannot be excluded. Thus, further detailed investigation of the quasinormal spectrum of Gauss-Bonnet-AdS black holes is also motivated, as it allows one to exclude or detect potential non-eikonal instabilities owing to non-zero (negative) Λ-term.
In the present paper we show that the spectrum consists of two essentially different, relatively the Gauss-Bonnet coupling α, types of modes. First type of modes approaches Schwarzschild-AdS values when α → 0. In particular, the most interesting for the holography hydrodynamic and sound modes, can be expressed through their Einsteinian limits by adding a linear (in α) correction. The analytical form of the correction is found by fitting the numerical data. The second type of modes are purely imaginary and their damping rates increase, approaching infinity when α → 0, i.e. these modes are non-perturbative in α and do not exist in the Einsteinian (Schwarzschild-AdS) limit. In the regime of infinitely large, relatively the AdS radius, black holes we reproduce the black brane regime and discuss parametric region of instability of black branes and properties of quasinormal modes/thermalization spectra for stable region and near the threshold of instability. This way we also complement recent study of the black brane's spectrum in [7,8].
The paper is organized as follows. Sec. 2 gives the basic information on the background black-hole metric, perturbation equations and diversity in terminology and designations between AdS/CFT and gravitational communities. Sec. 3 is devoted to analysis of the form of the effective potentials and to the search of a potential non-eikonal instability. Sec. 4 relates the details of the shooting method, which we used for finding the quasinormal modes. Sec. 5 and 6 discuss the found quasinormal modes for large and intermediate/small black holes respectively. Sec. 7 discusses the timescale for relaxation of perturbations. In sec. 8 we summarize the obtained results and mention prospects for future investigation.

The background metric and perturbation equations
The Lagrangian of the D-dimensional Einstein-Gauss-Bonnet theory has the form: The Gauss-Bonnet term corresponds to the full divergence in D = 4, so that it contributes only to the higher dimensional space-times. For D > 6 there are higher than the secondorder curvature corrections given by the Lovelock theory [22]. An exact solution for a static spherically symmetric black hole in the D-dimensional Einstein-Gauss-Bonnet theory (2.1) was found by Boulware and Deser in [23]. The metric has the form: where dΩ 2 n is a line element of the (n = D − 2)-dimensional sphere, and such that it satisfies the following relation: Here Λ is a cosmological constant, µ is a constant, proportional to mass, and the Gauss-Bonnet coupling constant is The black hole solution of (2.4), which goes over into the known Tangherlini solutions [24] allowing, in general, for a non-zero Λ-term, is First of all, we are interested in this branch of solutions, because it has the known asymptotically flat, de Sitter, and anti-de Sitter analogues 1 . Let us describe the range of parameters corresponding to such a black hole with the required AdS asymptotic at a negative Λ-term. As we shall have to investigate the spacetime behavior outside the black hole only, it is useful to express the black hole mass in terms of its size by introducing the radius of the event horizon r H > 0. This choice comes at a price that the parametric space must be taken carefully: one must ensure that the black-hole mass is positive and the event horizon exists. For a given negative value of α, there is a lower bound on the mass parameter µ: for which there is an event horizon r H > 0. From (2.6) it follows that for α < 0 we have The righthand side of (2.6) is positive for any Λ, implying that the black hole exists only for a positive asymptotic mass.
In order to measure all the quantities in the units of the same dimension we express µ as a function of the event horizon r H , as µ = n r n−1 Here we shall measure Λ in units of the AdS radius R, defined by relation Then, we have implying that α < R 2 .
(2.11) It turns out that, when R 2 /2 < α < R 2 , the solution (2.5) does not satisfy eq. (2.9) and describes a black hole, which is identical to the one with α < R 2 /2 after some re-scaling of parameters. Therefore, here we shall consider black holes with α ≤ R 2 /2 only. Other basic properties of Einstein-Gauss-Bonnet-AdS black holes were considered in [25]. After decoupling of the angular variables the perturbation equations can be reduced to the second-order master differential equations [26] where r * is the tortoise coordinate, , (2.13) and i stands for t (tensor ), v (vector ), and s (scalar ) perturbations. The explicit forms of the effective potentials V s (r), V v (r), and V t (r) [14] are given by where ℓ = 2, 3, 4, . . . is the multipole number and functions T (r) and R(r) can be written as follows designations used in gravity designations used in holography scalar sound vector shear tensor scalar Table 1: Terms used for types of gravitational perturbations in AdS/CFT-and gravity-oriented papers.
In the fields of gravity and holography different terms for the types of gravitational perturbations are used. Gravitationists work mostly with spherically symmetric black holes and are used to distinguish the three channels of perturbations according to the irreducible representations of the rotation group on (D − 2)-sphere. Thus, harmonics are called scalar, vector, and tensor. Holographic community analyzes mostly black holes with planar horizons and, following the hydrodynamic analogies, calls the corresponding channels: sound, shear and scalar (see Table 1).
The Gauss-Bonnet coupling λ GB used usually in holography, and the coupling α used in gravity are related as follows: In the regime of large, relatively the AdS-radius R, spherical black holes, the results for planar black holes (black branes) must be reproduced. For instance, the quasinormal modes of large AdS black holes found in [27] coincide with those of black branes [28] after the proper rescaling. Therefore, in order to compare the numerical data, it is necessary to know the relation between the multipole number ℓ in spherical symmetry and momentum k in the planar one: One can easily check that, in this limit formulas (44), (51), and (53) of [20], obtained for black holes, coincide, respectively, with (2.79), (2.80), and (2.81) of [8] for black branes.

The effective potentials and stability
The exact regions of the eikonal instability of the Einstein-Gauss-Bonnet-AdS black holes were found in [10] (see fig. 1 and 2 therein). These regions were found by analysis of regime of dominance of a negative gap in the effective potential at high ℓ. Thus, strictly speaking, it is not guaranteed that no instability exists for lower ℓ, especially, taking into consideration that such non-eikonal instability at the lowest multipole ℓ = 2 was found for the Gauss-Bonnet-de Sitter black holes [14]. Therefore, here we shall start from the analysis of positiveness of the effective potentials at ℓ = 2, concentrating on those regions which are free from the eikonal instability. If the effective potentials are positive definite everywhere outside the event horizon, then the black hole is stable, because the differential operator is a positive-definite self-adjoint operator in the Hilbert space of square integrable functions L(r * , dr * ), i. e. there are no negative-mode solutions for well-behaved initial data (smooth data of compact support). In other words, all solutions are bounded all of the time. In the parametric regions of the eikonal instability, no such well-behaved initial data can be provided, because of absence of convergence in ℓ. Nevertheless, the well-posed initial value problem is expected when there is no eikonal instability. We construct effective potentials and investigate regions of their positiveness for these eikonal-instability-free cases.
Thus, for example, the region for the eikonal instability (see eqs. (4.3) and (4.5) in [10]) in the scalar sector for D = 5 is and in the tensor channel is For D = 6 the region of instability obeys the following inequality (eq. 4.6 in [10]): Construction of effective potentials for the lowest multipole numbers ℓ shows that the scalar and tensor D = 5 potentials are positive-definite, while the vector potential has a negative gap near the event horizon for some values of the parameters. This gap does not lead to the eikonal instability, because here higher ℓ stabilize the system. At the same time, the vector potential can be easily deformed in such a way that the negative gap disappears, while the minimal allowed value of the damping rate of quasinormal frequencies does not become lower. This procedure is called S-deformation and was used in a number of works, devoted to analysis of stability of various black holes and fields [29]. The S-deformation for Gauss-Bonnet black holes was done in [15] implying, first of all, asymptotically flat black holes, but the whole procedure is the same for f (r) with any asymptotic behavior. Thus, according to [15], vector perturbations of Guass-Bonnet-(A)dS black holes are stable. Summarizing the D = 5 case: D = 5 Gauss-Bonnet-AdS black holes do not have other instability than the eikonal one, whose parametric region was found in [10].
The D = 6 case is different, because the scalar potential has a big negative gap, which is the deeper, the larger r H is. In this case we need to make an extensive search for quasinormal modes in the region of parameters, which provides negativeness of the potential, and make sure that no growing modes appear in the spectrum. Unstable modes of a spherically symmetric background do not oscillate, as it was shown in [30], so that, we need to search for unstable modes only on the imaginary axis. A restriction upon the possible values of unstable modes comes from the depth of the negative potential gap: since V − ω 2 > 0 guarantees stability, where V min is the minimal value of the effective potential at r H ≤ r < ∞. With the help of the shooting method, which will be described in the next section, we found no growing purely imaginary modes in the above region. This situation is not new for asymptotically anti-de Sitter space-times, as, for example, the effective potential of the Reissner-Nordström-AdS black hole has large and even infinite region of negativeness of the effective potential, but, owing to the asymptotic non-flatness of spacetime and specific behavior of r * coordinate, no instability develops in this case [32].

The shooting method
In order to analyse gravitational stability of a black hole whose perturbations are governed by cumbersome effective potentials, numerical analysis of black holes' quasinormal spectrum is the only feasible way. For asymptotically AdS black holes, the quasinormal mode boundary conditions for gravitational perturbations require purely incoming waves at the event horizon and the Dirichlet boundary conditions at infinity. If a growing mode is found, then the considered system is unstable. Although, usually, damped quasinormal modes have both real and imaginary parts, i. e. are oscillating, the growing modes for spherically symmetric background are known to be non-oscillating, that is, pure imaginary [30]. Thus, it is sufficient to look for unstable modes only along the imaginary axis.
The shooting method is useful for finding of quasinormal modes of asymptotically Gauss-Bonnet-AdS black holes. Its general features are related in [31], while here, in a similar fashion with [32], we shall adopt it for asymptotically anti-de Sitter space-times. Equation (2.12) has at least two regular singular points: at spacial infinity and at the event horizon. At the event horizon The quasinormal boundary conditions at the event horizon imply that At spatial infinity the two linear independent solutions of Ψ(r) are Let us consider the new function If Ψ(r) satisfies the quasinormal boundary conditions, then y(r H ) = Z 0 . Since y(r) satisfies the linear equation, one can choose Z 0 = 1 to fix the scale. Then y ′ (r H ) can be found from equation (2.12) in the following way Using these initial conditions at the event horizon, we solve equation (2.12) numerically for each ω with the help of the Wolfram Mathematica R built-in function N DSolve for r H ≤ r ≤ r f , where r f ≫ r H (the typical value of r f is ∼ 10 4 r H ).
In the general case the behavior of Ψ(r) at infinity is a superposition of the two solutions (4.1), Ψ C (r) and Ψ D (r): where Ψ C (r) satisfies the quasinormal boundary condition (4.2). If ω is the quasinormal frequency, then the corresponding solution must satisfy the boundary conditions (4.2) at spatial infinity and, thereby, Z D = 0. In order to calculate Z C and Z D we obtain expansions of Ψ C (r) and Ψ D (r) far from the black hole in analytic form. The expansion for Ψ C (r) contains only negative powers of r, . . . Since the difference between the roots of the indicial equation is integer, the expansion for Ψ D (r) contains not only powers of r, but also contributions, proportional to ln(r). In particular, for D = 5, one has while for D = 6, As the series (4.6, 4.7) are convergent, we have used only the first three terms of the expansions. The coefficients A 1 , A 2 , A 3 , B 1 , B 2 , B 3 , C 1 , C 2 , and C 3 are determined by substituting the expansions (4.6, 4.7) into (2.12). Thus, our numerical procedure is the following. For any given value of ω we integrate the equation (2.12) numerically, by imposing quasinormal boundary condition at the event horizon (4.4). At a large distance from the black hole, near r f , we fit the numerical values of Ψ(r) by the analytical expression (4.5) and find Z C and Z D by solving the least squares problem. The quasinormal modes correspond to the roots of the equation Z D (ω) = 0. (4.8) One should be careful when applying the shooting method to finding of higher overtones or situations when the purely imaginary mode acquires a non-zero real part, because in these cases the shooting method requires extra precision in order to reproduce the known analytical results accurately [20].

Large (r H /R ≫ 1) black holes
Once the regions of instability are clearly determined, we are in position to analyze quasinormal modes in the stability sector. There are three channels of perturbations (scalar, vector and tensor) and two qualitatively different regimes of black holes: large and small (relatively the anti-de Sitter radius R). We shall discuss the features of the quasinormal spectrum for all these cases, concentrating on the most interesting from the point of view of the holography issues: the hydrodynamic and sound modes of large black holes.  Table 2: Quasinormal modes of the Gauss-Bonnet-AdS black hole: D = 5, scalar channel, ℓ = 2, r H /R = 6. The frequencies are: ω 0 is the fundamental sound mode, ω 1 is the branch of modes with spacing linear in r H /R, ω h is the lowest mode on the imaginary axis.

Scalar (sound ) channel
There are three essentially different kinds of modes in the scalar channel (see table 2): 1. perturbative in α, sound mode, which usually has the lowest damping rate and whose real part represents the speed of sound in a medium (ω 0 in Table 2); 2. perturbative (in α) branch of modes which scale as r H /R and go over into their Schwarzschild-AdS limits [33] when α = 0. The dominant mode of this series is shown as ω 1 in Table 2. At higher n this branch of the spectrum is equidistant, as it takes place for Schwazrschild-AdS case. Though, detecting of high overtones n is difficult by the shooting method.
3. branch of purely imaginary modes, which do not exist in the α = 0 limit, because they "come from infinity" at small α along the negative imaginary axis, approaching and crossing the origin at the value of α corresponding to the threshold of eikonal instability. Thus, these modes are non-perturbative in α. Notice, that in Table (2) data for some of the modes in the non-perturbative branch are absent, because at sufficiently small α the frequencies become too large and are difficult to detect by the shooting method.
The values of the sound mode given in fig. 1 for ℓ = 2, r H /R = 6 can be very well fitted by the linear law: From fig. 2 one can see that, at a fixed α, Re(ω) of the sound mode is growing and at sufficiently large ℓ this growth looks linear. The damping rate, given by imaginary part of the frequency, approaches some constant. Notice, that when ℓ in fig. 2 is larger than r H /R, such regime is irrelevant for holography, because it corresponds to a very large momentum of matter propagating on a relatively small sphere. In the other limit, of large Im Ω R Figure  Re Ω R Im Ω R mode equals c s k (where c s is speed of sound and k is momentum), this may indicate that c s is unaffected by the Gauss-Bonnet coupling α. Indeed, this is also confirmed by the calculations for black branes (see eq. 3.25 in [8]). By fitting a lot of numerical data, we can guess that when α is nonzero, the imaginary part of ω obeys the following relation: where A(D) ≈ 3 for D = 5, A(D) ≈ 2.3 for D = 6, A(D) ≈ 2 for D = 7, etc. Thus, we can suppose that A(D) ≈ (D+1)/(D−3). This relation is satisfied with good accuracy not only for very large, but also for moderately large (e. g. r H /R = 6) black holes. For sufficiently large black holes, approaching the black brane regime, the quasinormal frequencies are where ω SAdS is the sound mode of the Schwarzschild-AdS black hole (α = 0).

Vector (shear ) channel
The vector channel has also three different types of modes: • purely imaginary "hydrodynamic" mode; • purely imaginary non-perturbative in α modes (similar and numerically close to those found for the scalar channel); • modes with non-zero real part, which scale as r H /R and have the Einsteinian analogues in the limit α = 0 [33]. These frequencies are also numerically close to those of the scalar and tensor channels.
The vector modes allow one to compute the η/s ratio [3,34]. The purely imaginary hydrodynamic modes for large Schwarzschild-anti-de Sitter black hole have the form (see formula (16) in [33]) In order to see how expression (5.5) is modified when adding the non-zero α, we plotted the numerically found, purely imaginary mode ω as a function of α for various values of ℓ = 2, 3, 4, 5, r H /R = 20 (see fig. 4). The data fits very well the following formulas for D = 5, and for D = 6, The factors in these fits approximate the SAdS expression (5.5). For example, for r H /R = 20, ℓ = 5 and D = 5, the exact formula (5.5) gives ωR = −0.4i, while the numerical data is ωR = −0.403i. Thus, we conclude that when α is nonzero, formula (5.5) can be generalized as follows: Comparing the formulas for the sound and hydrodynamic modes (5.4) and (5.6), we observe that the damping rates of these modes are related as follows: Im(ω sound ) = 2 3 · Im(ω hydro ).  The relation (5.7) also occurs for the Einstein-Gauss-Bonnet black branes [34].
The purely imaginary non-perturbative modes for all three types of gravitational perturbations have damping rate which increases as α is deceasing (fig. 5). Thus, for sufficiently small α they reach arbitrarily large damping rates and their contribution to the spectrum is negligible, while at α = 0 they disappear from the spectrum.
At some sufficiently large α the purely imaginary, non-perturbative modes acquire the non-zero real part at some lower ℓ. This occurs in the region of eikonal instability, so that higher ℓ modes are unboundedly growing [10]. One can see how the purely imaginary "mode" becomes oscillating though the exact solutions of the wave equations found in [20] (see, for example eq. (28) therein) for the unstable case α = R 2 /2. Tensor channel has only two kinds of modes: the same non-perturbative branch as for scalar and vector channels and the one scaling as r H /R and having the Einsteinian limit when α → 0. Tensor perturbative modes do not lay on the linear fit, as can be seen from fig. 6, while the non-perturbative branch behaves similar with those for the other two channels and numerical values of the frequencies are close for all three types of perturbations. The non-linear behavior depicted on fig. 6 is different from the linear in α sound and hydrodynamic modes. As we cannot completely exclude that the shooting method simply becomes inaccurate in this case 2 , calculations with the help of an alternative method would be desirable.

Moderate (r H /R ∼ 1) and small black holes
Small black holes are not interesting from the point of view of possible applicability within AdS/CFT correspondence. Nevertheless, in order to have the full picture, we shall consider this regime as well.
Quasinormal modes of Schwarzschild-AdS black holes approach the normal modes of the empty AdS space-time when r H → 0 [35]. This is evidently not applicable to the Gauss-Bonnet case, because the limit r H → 0 is not achievable at any non-zero value of α, since all sufficiently small black holes are unstable (see eqs. 3.2, 3.3, 3.4).
The eikonal instability develops due to the purely imaginary nonperturbative frequencies of the scalar or tensor channel. When the absolute value of α grows, the damping rates decrease for sufficiently large ℓ. Then, at a fixed ℓ the non-perturbative frequency with the slowest decay rate becomes dominant and, at some critical value of α, it turns into the unstable (growing) mode. When increasing ℓ, the critical value of α (at a given ℓ) approaches the threshold of the eikonal instability [10], so that for any given α inside the region of instability there is always a finite value of the multipole number ℓ inst for which a growing mode appears in the spectrum. The lowest unstable multipole number can be arbitrarily large, if the value of α is sufficiently close to the threshold of the eikonal instability. While perturbations of higher than ℓ inst multipoles are always unstable, perturbations of lower multipoles do not necessary contain a growing mode. Thus, in [10] it was observed that within the region of eikonal instability the lower-multipole non-perturbative frequencies are damped and may "duplicate" by acquiring a real part (with positive and negative signs).

Thermalization timescale
According to the AdS/CFT correspondence, perturbation of a spherical, D-dimensional asymptotically AdS black hole is related to the perturbation of the thermal state propagating on the sphere in the (D − 1)-dimensional world. When the black hole is sufficiently large, we approach the results of the black brane perturbations and the quantum system can be considered as the one not limited by a sphere. Damping rate of the dominant quasinormal mode characterizes the time, which is necessary for the dual quantum system to return to the thermal equilibrium. Thus, the thermalization timescales can be expressed as τ = 1 |Im(ω)| (7.1) When the Gauss-Bonnet coupling is turned on, the Hawking temperature of large black holes can be approximated as follows: 2) being, thereby, linear in α and r H . Therefore, using (5.4) and (5.6) for the fundamental scalar (sound) and vector (hydrodynamic) modes, one has Taking into account that | α| 0.1R 2 , we can neglect α in the denominator of the above expression, so that The above relation connects the damping rates to the temperature ratios for large black holes in the Einstein and Einstein-Gauss-Bonnet theories. The sound mode has the longest lifetime, so that it determines the characteristic time for relaxation of perturbations. An essential constrain on possible applicability of the Einstein-Gauss-Bonnet-AdS black hole and brane backgrounds must be imposed by the observed eikonal instability. If the GB coupling is large enough, black holes (and branes) are unstable [10], so that the black hole does not approach an equilibrium state, which could be characterized by a given set of black-hole parameters, Hawking temperature of the event horizon, etc. Thus, no relaxation of perturbations and thermalization should occur in such an unstable system.
In our opinion, when considering holographic applications one should not turn the blind eye to the eikonal instability, justifying it by the fact that the phenomena occurs at large momentum k, while one can be constrained by small momentum in the field theory side. First, when slightly increasing the Gauss-Bonnet coupling constant, the instability occurs at lower and lower momenta, reaching, for the large black hole case, the minimal ℓ = 2 multipole. For example, taking limit of r H → ∞ in (3.2) and (3.3) and substituting into (2.16), we find that the five-dimensional black brane is gravitationally unstable outside the following region Once an instability takes place, then the following issue arises. The perturbation contains various values of momenta and in order to ignore the non-equilibration at high momenta, one should develop mathematically consistent cut-off for this case, if that is possible. Thus, when discussing the possible holographic model for quantum liquids possessing superfluidity one can formally derive the α-corrected formula for the viscosity to entropy density ratio [34], η s ≈ 4πk (1 − 4λ GB ). (7.4) In this context the regime λ GB 1/4 in which the viscosity could be made small was considered in the literature [34,9]. Here from (2.16) (n = 3) we can see that the above regime corresponds to α R 2 /2, which is unstable. Thus, the holographic description of zero viscosity regime and the conclusions made within the Einstein-Gauss-Bonnet-AdS model look questionable and, should not be applied to quantum liquids with high momentum.
Another interesting question is related to the causality violation in theories with higher curvature corrections. For the black brane regime, the causality is violated for an evidently larger region (λ GB > 9/100) [36] than the one suffering from instability (α > 0.146R 2 , λ GB > 0.125). At the same time the analysis of causality suggested in [36] implies a well defined eikonal regime (see e.g. eq. (26) in [36]), valid Fourier transformations (eq. (7) in [36]), etc., in order to work with quasi-particles and show that they can travel faster than the speed of light. All these constituents are not well defined when the eikonal instability occurs.
Recently, an analytical deduction of the viscosity to the entropy density ratio has been found for the system with translational symmetry by representing the black brane within a classical membrane paradigm [37]. It would be interesting to deduce the η/µ ratio from the quasinormal spectrum obtained here and compare it with that found in [37].

Discussion
Here we have studied quasinormal modes of gravitational perturbations of the Einstein-Gauss-Bonnet-AdS black holes in detail. The numerical analysis of the quasinormal spectrum has shown that there is no other than eikonal instability for such black holes. The numerical data for the sound and hydrodynamic modes are shown to be very well described by an analytical formula with a linear in α correction to the corresponding modes of the Schwarzschild-AdS black hole. At the same time, when λ GB is not small enough, the black holes and branes suffer from the instability, so that the holographic interpretation of perturbation of such black holes become questionable, as, for example, the claimed viscosity bound violation [34]. The eikonal instability phenomenon should be taken into account when considering perturbations of various black holes and branes in the Gauss-Bonnet theory [38].
The non-perturbative in α quasinormal modes exist not only for considered here gravitational perturbations of the Einstein-Gauss-Bonnet-AdS black holes, but also in a number of other cases, which we would to like summarize here: • Gravitational perturbations of asymptotically flat Einstein-Gauss-Bonnet black holes in various D 3 [13,39].
• Test scalar field perturbations in the background of Einstein-Gauss-Bonnet-AdS spherical black holes [20].
All the above cases may indicate possible independence of existing of purely imaginary non-perturbative in α quasinormal modes on such peculiarities as: spin of a perturbed field, number of spacetime dimensions, asymptotic behavior (flat, dS, AdS), and, possibly, even character of higher curvature corrections (Gauss-Bonnet, Lovelock, R 4 , etc.). The broadness of this phenomena may have implications for D = 4 large astrophysical black holes. Thus, it is reasonable to check the presence of the non-perturbative branch of modes for Einstein-dilaton-Gauss-Bonnet black holes, dynamical Chern-Simons black holes, various black hole solutions in f(R) gravity etc. The gravitational spectra of the latter cases were investigated recently only from the point of view of small deviations from their Schwarzschild and Kerr values [41], so that it seems reasonable to reconsider these works taking into account possible existence of non-perturbative modes. As nowadays, gravitational wave experiments do not suggest a strict constrain on possible deviations from Kerr geometry [42], the coupling constants in higher curvature corrections can be sufficiently large, so that the non-perturbative modes have a chance to be dominating in the signal at some values of the parameters. These cases will be studied in our future publications [43].
Our work can be extended in a number of other ways: by adding an axion parameter [44], considering higher curvature corrections and higher D.