The general parametrization of a black-hole spacetime in arbitrary metric theories of gravity includes an infinite set of parameters. It is natural to suppose that essential astrophysically observable quantities, such as quasinormal modes, parameters of shadow, electromagnetic radiation and accreting matter in the vicinity of a black hole, must depend mostly on a few of these parameters. Starting from the parametrization for spherically symmetric configurations in the form of infinite continued fraction, we suggest a compact representation of the asymptotically flat spherically symmetric and slowly rotating black holes in terms of only three and four parameters respectively. A subclass of arbitrarily rotating black holes belonging to the Carter family can also be parametrized by only four parameters. This approximate representation of a black-hole metric should allow one to describe physical observables in the region of strong gravity.

A $(3+1)$-dimensional Einstein-Gauss-Bonnet theory of gravity has been recently formulated in [D. Glavan and C. Lin, Phys. Rev. Lett. 124, 081301 (2020)] which is different from the pure Einstein theory, i.e. bypasses the Lovelock's theorem, and avoids Ostrogradsky instability. The theory was formulated in $D > 4$ dimensions and its action consists of the Einstein-Hilbert term with a cosmological constant, and the Gauss-Bonnet term multiplied by a factor $1/(D-4)$. Then, the four-dimensional theory is defined as the limit $D \to 4$. Here we generalize this approach to the four-dimensional Einstein-Lovelock theory and formulate the most general static $4D$ black-hole solution allowing for a $\Lambda$-term (either positive or negative) and electric charge $Q$. As the metric functions cannot be found in a closed form in the general case, we develop and share publicly the code which constructs the metric functions for every given set of parameters.

A (3+1)-dimensional Einstein-Gauss-Bonnet theory of gravity has been recently formulated as the $D \to 4$ limit of the higher dimensional field equations after the rescaling of the coupling constant. This theory is a nontrivial generalization of General Relativity, it bypasses the Lovelock's theorem, and avoids the Ostrogradsky instability. This approach has been extended to the four-dimensional Einstein-Lovelock gravity. Here we study the eikonal gravitational instability of asymptotically flat, de Sitter and anti-de Sitter black holes in the four dimensional Einstein-Gauss-Bonnet and Einstein-Lovelock theories. We find parametric regions of the eikonal instability for various orders of the Lovelock gravity, values of coupling and cosmological constants, and share the code which allows one to construct the instability region for an arbitrary set of parameters.

Recently, numerical solutions to the field equations of Einstein-scalar-Gauss-Bonnet gravity that correspond to black-holes with non-trivial scalar hair have been reported. Here, we employ the method of the continued-fraction expansion in terms of a compact coordinate in order to obtain an analytical approximation for the aforementioned solutions. For a wide variety of coupling functionals to the Gauss-Bonnet term we were able to obtain analytical expressions for the metric functions and the scalar field. In addition we estimated the accuracy of these approximations by calculating the black-hole shadows for such black holes. Excellent agreement between the numerical solutions and analytical approximations has been found.

The Schwarzschild star is an ultracompact object beyond the Buchdahl limit, which has Schwarzschild geometry outside its surface and positive pressure in the external layer which vanishes at the surface. Recently it has been shown that the Schwarzschild star is stable against spherically symmetric perturbations. Here we study arbitrary axial nonspherical perturbations, and show that the observable quasinormal modes can be as close to the Schwarzschild limit as one wishes, what makes the Schwarzschild star a very good mimicker of a black hole. The decaying time-domain profiles prove that the Schwarzschild star is stable against nonspherical perturbations as well. Another peculiar feature is the absence of echoes at the end of the ringdown. Instead we observe a nonoscillating mode which might belong to the class of algebraically special modes. At asymptotically late times, Schwarzschildian power-law tails dominate in the signal.

In [P. Kanti, B. Kleihaus, J. Kunz, Phys. Rev. Lett. 107, 271101 (2011)] it was shown that the four-dimensional Einstein-dilaton-Gauss-Bonnet theory allows for wormholes without introducing any exotic matter. The numerical solution for the wormhole was obtained there and it was claimed that this solution is gravitationally stable against radial perturbations, what, by now, would mean the only known theoretical possibility for existence of an apparently stable, four-dimensional and asymptotically flat wormhole without exotic matter. Here, more detailed analysis of perturbations shows that the Kanti-Kleihaus-Kunz wormhole is unstable against small perturbations for any values of its parameters. The exponential growth appears in the time domain after a long period of damped oscillations, in the same way as it takes place in the case of unstable higher-dimensional black holes in the Einstein-Gauss-Bonnet theory. The instability is driven by the purely imaginary mode, which is nonperturbative in the Gauss-Bonnet coupling$\alpha$.

Quasinormal modes of a massive Dirac field were calculated for various static black hole backgrounds with the help of the WKB formula. These estimations, however, are rough and valid only for very small values of $\mu M$, where $M$ and $\mu$ are mass of the black hole and field respectively. Thus, no accurate calculations of massive Dirac modes are known even for the Schwarzschild black hole and this is all the more so for the Kerr solution. Here we calculate quasinormal modes of a massive Dirac field in the Kerr background. We have shown that the infinitely long-lived quasinormal modes (quasiresonances), which exist for boson fields, appear also in the fermions' quasinormal spectrum. Two chiralities of massive fermions lead to an additional "fine structure" in the spectrum. We discuss the effect of this fine structure on the behavior of quasiresonances and the stability. The analysis is also extended to a charged massive field in the Kerr-Newman background.

We study quasinormal modes and scattering properties via calculation of the $S$-matrix for scalar and electromagnetic fields propagating in the background of spherically and axially symmetric, traversable Lorentzian wormholes of a generic shape. Such wormholes are described by the Morris-Thorne ansatz and its axially symmetric generalization. The properties of quasinormal ringing and scattering are shown to be determined by the behavior of the wormhole's shape function $b(r)$ and shift factor $\Phi(r)$ near the throat. In particular, wormholes with the shape function $b(r)$, such that $b'(r) \approx 1$, have very long-lived quasinormal modes in the spectrum. We have proved that the axially symmetric traversable Lorentzian wormholes, unlike black holes and other compact rotating objects, do not allow for superradiance. As a by product we have shown that the 6th order WKB formula used for scattering problems of black or wormholes provides high accuracy and thus can be used for quite accurate calculations of the Hawking radiation processes around various black holes.

We construct an analytical approximation for the numerical black hole metric of P. Kanti, et. al. [PRD54, 5049 (1996)] in the four-dimensional Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The continued fraction expansion in terms of a compactified radial coordinate, used here, converges slowly when the dilaton coupling approaches its extremal values, but for a black hole far from the extremal state, the analytical formula has a maximal relative error of a fraction of one percent already within the third order of the continued fraction expansion. The suggested analytical representation of the numerical black hole metric is relatively compact and good approximation in the whole space outside the black hole event horizon. Therefore, it can serve in the same way as an exact solution when analyzing particles' motion, perturbations, quasinormal modes, Hawking radiation, accreting disks and many other problems in the vicinity of a black hole. In addition, we construct the approximate analytical expression for the dilaton field.

Higher derivative extensions of Einstein gravity are important within the string theory approach to gravity and as alternative and effective theories of gravity. H. L\"u, A. Perkins, C. Pope, K. Stelle [Phys.Rev.Lett. 114 (2015), 171601] found a numerical solution describing a spherically symmetric non-Schwarzschild asymptotically flat black hole in the Einstein gravity with added higher derivative terms. Using the general and quickly convergent parametrization in terms of the continued fractions, we represent this numerical solution in the analytical form, which is accurate not only near the event horizon or far from black hole, but in the whole space. Thereby, the obtained analytical form of the metric allows one to study easily all the further properties of the black hole, such as thermodynamics, Hawking radiation, particle motion, accretion, perturbations, stability, quasinormal spectrum, etc. Thus, the found analytical approximate representation can serve in the same way as an exact solution.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

Collaborative international efforts under the name of the Event Horizon Telescope project, using sub- mm very long baseline interferometry, are soon expected to provide the first images of the shadow cast by the candidate supermassive black hole in our Galactic center, Sagittarius A*. Observations of this shadow would provide direct evidence of the existence of astrophysical black holes. Although it is expected that astrophysical black holes are described by the axisymmetric Kerr solution, there also exist many other black hole solutions, both in general relativity and in other theories of gravity, which cannot presently be ruled out. To this end, we present calculations of black hole shadow images from various metric theories of gravity as described by our recent work on a general parameterisation of axisymmetric black holes [R. Konoplya, L. Rezzolla and A. Zhidenko, Phys. Rev. D 93, 064015 (2016)]. An algorithm to perform general ray-tracing calculations for any metric theory of gravity is first outlined and then employed to demonstrate that even for extremal metric deformation parameters of various black hole spacetimes, this parameterisation is both robust and rapidly convergent to the correct solution.

Recently LIGO and VIRGO collaborations reported about observation of gravitational-wave signal corresponding to the inspiral and merger of two black holes, resulting into formation of the final black hole. It was shown that the observations are consistent with the Einstein theory of gravity with high accuracy limited mainly by the statistical error. Angular momentum and mass of the final black hole were determined with rather large allowance of tens of percents. Here we shall show that this indeterminacy in the range of the black-hole parameters allows for some not negligible deformations of the Kerr spacetime leading to the same frequencies of black-hole ringing. This means that at the current precision of the experiment there remain some possibilities for alternative theories of gravity.

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.

We find the quasinormal modes of the charged scalar and Dirac fields in the background of the rotating charged black holes, described by the Kerr-Newman-de Sitter solution. The dependence of the quasinormal spectrum upon the black hole parameters mass M, angular momentum a, charge Q, as well as on values of the \Lambda-term and field charge q is investigated. Special attention is given to the near extremal limit of the black hole charge. In particular, we find that for both scalar and Dirac fields, charged perturbations decay quicker for q>0 and slower for q<0 for values of black holes charge Q less than than some threshold value, which is close to the extremal value of charge and depend on parameters of the black holes.

In our earlier work (PRL 103 (2009) 161101) it was shown that nonextremal highly charged Reissner-Nordstrom-de Sitter black holes are gravitationally unstable in D>6-dimensional space-times. Here, we find accurate threshold values of the $\Lambda$-term at which the instability of the extremally charged black holes starts. The larger $D$ is, the smaller is the threshold value of $\Lambda$. We have shown that the ratio $\rho = r_{h}/r_{cos}$ (where $r_{cos}$ and $r_{h}$ are the cosmological and event horizons) is proportional to $e^{-(D-4)/2}$ at the onset of instability for D=7,8,...11, implying that the same law should fulfill for arbitrary D. This is numerical evidence that extremally charged Reissner-Nordstrom-de Sitter black holes are gravitationally unstable for D>6, while asymptotically flat extremally charged Reissner-Nordstrom black holes are stable for all D. The instability is not connected to the horizon instability discussed recently in the literature, and, unlike the later one, develops also outside the event horizon, that is, it can be seen by an external observer. In addition, for the nonextremal case through fitting of the numerical data we obtained an approximate analytical formula which relates values of charge and the $\Lambda$-term at the onset of instability.

So far analysis of the quasinormal spectrum of a massive charged scalar field in the black hole background has been limited by the regime of small \mu M and qQ, where \mu, q (M, Q) are mass and charge of the field (black hole). Here we shall present a comprehensive picture of quasinormal modes, late-time tails and stability of a massive charged scalar field around Kerr-Newman black holes for any physically meaningful values of the parameters. We shall show that despite presence of the two mechanisms of superradiance (owing to black hole's rotation and charge) and the massive term creating growing bound states, there is no indication of instability under quasinormal modes' boundary conditions. We have shown that for some moderate values of qQ dominant quasinormal modes may have arbitrarily small real oscillation frequencies Re(\omega). An analytic formula for the quasinormal modes has been derived in the regime of large qQ. The larger the field's charge, the sooner asymptotic tails dominate in a signal, making it difficult to extract quasinormal frequencies from a time-domain profile. Analytic expressions for intermediate and asymptotically late-time tails have been found for the Reissner-Nordstr\"om black hole. For the near extremal Kerr-Newman black holes we have obtained a more general picture of the mode branching found recently for massless fields [arXiv:1212.3271] in the Kerr background.

Usually alternative theories of gravity imply deviations from the well-known Kerr space-time, a model of an isolated black hole in General Relativity. In the dominant order, the deformed Kerr metric, free of closed time-like curves outside the event horizon, has been suggested recently by Johannsen and Psaltis. It has a single deformation parameter which is not constrained by the current observations, allowing, thereby, for a kind of unified and simple phenomenological description of black holes in various theories of gravity. Here we consider a number of classical and quantum phenomena of radiation in the vicinity of such deformed Schwarzschild-like and Kerr-like black holes: spiralling of particles into black holes, decay of fields propagating in the black hole's background, Hawking radiation. In particular, we calculate some quantitative characteristics of the above phenomena, such as the binding energy of particles, quasinormal modes, late-time tails of fields of various spin, intensity of Hawking radiation. The binding energy released when a particle goes over from a given stable orbit in the equatorial plane to the innermost stable one is calculated for such non-Kerr black holes. Due to inseparability of the wave equations in the general case, the perturbations and stability of scalar, Dirac, and electromagnetic fields are analyzed for vanishing rotation only. The dependence of the radiation phenomena on the deformation parameter is discussed.

In the present paper we analyze the spectrum of quasinormal modes for massive scalar and Dirac fields, allowing for both tardyonic ($\mu^2 >0$) and tachyonic ($\mu^2 <0$) masses, in the expanding and rotating cosmological background. The spectrum found shows a number of peculiar features, which are absent in the Minkowski space-time. A hypothetical particle that moves faster than light, \emph{a tachyon}, is known to be classically unstable in the Minkowski space-time. This instability has its analog at the quantum level: small vacuum fluctuations of the field lead to the unbounded growth of the amplitude, so that the appearance of the real tachyons in the spectrum means that there is catastrophic instability in the theory. It has been conjectured a long time ago that possibly the lightest particles with nonzero mass, the neutrino, may be a tachyon. Here we shall show that in the rotating and expanding Universe tachyons are stable if their mass is less than some constant, which is related to the Universe's rotation and expansion scales. The current upper bound on the rotation scale gives us a very small upper bound on the tachyon's mass which is many orders less than the mass of the electron. This might explain why only very light particles have the chance of becoming tachyons. It is shown that the spectrum of the "normal" ($\mu^2 >0$) Dirac field has a discontinuity as a function of the z-component of the wave vector $k_z$ at $k_z = 0$. Keywords: OPERA, superluminal neutrino.

Dictated by the string theory and various higher dimensional scenarios, black holes in $D>4$-dimensional space-times must have higher curvature corrections. The first and dominant term is quadratic in curvature, and called the Gauss-Bonnet (GB) term. We shall show that although the Gauss-Bonnet correction changes black hole's geometry only softly, the emission of gravitons is suppressed by many orders even at quite small values of the GB coupling. The huge suppression of the graviton emission is due to the multiplication of the two effects: the quick cooling of the black hole when one turns on the GB coupling and the exponential decreasing of the grey-body factor of the tensor type of gravitons at small and moderate energies. At higher $D$ the tensor gravitons emission is dominant, so that the overall lifetime of black holes with Gauss-Bonnet corrections is many orders larger than was expected. This effect should be relevant for the future experiments at the Large Hadron Collider (LHC). Keywords: Hawking radiation, black hole evaporation.

We study the stability of static, spherically symmetric solutions to the Einstein equations with a scalar field as the source. We describe a general methodology of studying small radial perturbations of scalar-vacuum configurations with arbitrary potentials V(\phi), and in particular space-times with throats (including wormholes), which are possible if the scalar is phantom. At such a throat, the effective potential for perturbations V_eff has a positive pole (a potential wall) that prevents a complete perturbation analysis. We show that, generically, (i) V_eff has precisely the form required for regularization by the known S-deformation method, and (ii) a solution with the regularized potential leads to regular scalar field and metric perturbations of the initial configuration. The well-known conformal mappings make these results also applicable to scalar-tensor and f(R) theories of gravity. As a particular example, we prove the instability of all static solutions with both normal and phantom scalars and V(\phi) = 0 under spherical perturbations. We thus confirm the previous results on the unstable nature of anti-Fisher wormholes and Fisher's singular solution and prove the instability of other branches of these solutions including the anti-Fisher "cold black holes".

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.

We have shown that higher dimensional Reissner-Nordstr\"om-de Sitter black holes are gravitationally unstable for large values of the electric charge and cosmological constant in $D \geq 7$ space-time dimensions. We have found the shape of the slightly perturbed black hole at the threshold point of instability. Why only $D=4, 5$ and 6 dimensional worlds are favorable as to the black stability remains unknown.

It is well-known that the perturbations of Schwarzschild black holes are governed by a wave equation with some effective potential. We consider perturbations of a gas in a tube called Laval nozzle, which is narrow in the middle and has a sonic point in the throat. By equating the wave equation in a Laval nozzle of an arbitrary form with the wave equation of spin-s perturbations of Schwarzschild black holes, we find the exact expression for the form of the Laval nozzle, for which acoustic perturbations of the gas flow corresponds to the general form of perturbations of Schwarzschild black holes. This allows observation, in a laboratory, of the acoustic waves, which are analogue of damping quasinormal oscillations of Schwarzschild black holes. The found exact acoustic analog allows to observe also some other phenomena governed by the wave equation, such as the wave scattering and tunneling.

Using the Frobenius method, we find high overtones of the Dirac quasinormal spectrum for the Schwarzschild black hole. At high overtones, the spacing for imaginary part of $\omega_{n}$ is equidistant and equals to $\Im{\omega_{n+1}}-\Im{\omega_{n}} =i/8M$, ($M$ is the black hole mass), which is twice less than that for fields of integer spin. At high overtones, the real part of $\omega_{n}$ goes to zero. This supports the suggestion that the expected correspondence between quasinormal modes and Barbero-Immirzi parameter in Loop Quantum Gravity is just a numerical coincidence.