Perturbations of black holes, initially considered in the context of possible observations of astrophysical effects, have been studied for the past ten years in string theory, brane-world models and quantum gravity. Through the famous gauge/gravity duality, proper oscillations of perturbed black holes, called quasinormal modes (QNMs), allow for the description of the hydrodynamic regime in the dual finite temperature field theory at strong coupling, which can be used to predict the behavior of quark-gluon plasmas in the nonperturbative regime. On the other hand, the brane-world scenarios assume the existence of extra dimensions in nature, so that multidimensional black holes can be formed in a laboratory experiment. All this stimulated active research in the field of perturbations of higher-dimensional black holes and branes during recent years. In this review recent achievements on various aspects of black hole perturbations are discussed such as decoupling of variables in the perturbation equations, quasinormal modes (with special emphasis on various numerical and analytical methods of calculations), late-time tails, gravitational stability, AdS/CFT interpretation of quasinormal modes, and holographic superconductors. We also touch on state-of-the-art observational possibilities for detecting quasinormal modes of black holes.

Following previous work of ours in spherical symmetry, we here propose a new parametric framework to describe the spacetime of axisymmetric black holes in generic metric theories of gravity. In this case, the metric components are functions of both the radial and the polar angular coordinates, forcing a double expansion to obtain a generic axisymmetric metric expression. In particular, we use a continued-fraction expansion in terms of a compactified radial coordinate to express the radial dependence, while we exploit a Taylor expansion in terms of the cosine of the polar angle for the polar dependence. These choices lead to a superior convergence in the radial direction and to an exact limit on the equatorial plane. As a validation of our approach, we build parametrized representations of Kerr, rotating dilaton, and Einstein-dilaton-Gauss-Bonnet black holes. The match is already very good at lowest order in the expansion and improves as new orders are added. We expect a similar behavior for any stationary and axisymmetric black-hole metric.

The WKB approach for finding quasinormal modes of black holes, suggested in [1] by Schutz and Will at the first order and later developed to higher orders [2][3][4], became popular during the past decades, because, unlike more sophisticated numerical approaches, it is automatic for different effective potentials and mostly provides sufficient accuracy. At the same time, the seeming simplicity of the WKB approach resulted in appearance of a big number of partially misleading papers, where the WKB formula was used beyond its scope of applicability. Here we review various situations in which the WKB formula can or cannot bring us to reliable conclusions. As the WKB series converges only asymptotically, there is no mathematically strict criterium for evaluation of an error. Therefore, here we are trying to introduce a number of practical recipes instead and summarize cases in which higher WKB orders improve accuracy. We show that averaging of the Padé approximations, suggested first by J. Matyjasek and M. Opala [4], leads to much higher accuracy of the WKB approach, estimate the error and present the automatic code [5] which computes quasinormal modes and grey-body factors. WKB orderPublication 1 B. Mashhoon [8], B. Schutz and C. Will [1] 2-3 S. Iyer and C. Will [2] 4-6 R. Konoplya [3] 7-13 J. Matyjasek and M. Opala [4] 3 0.