We consider quasinormal modes and Hawking radiation of four dimensional asymptotically flat black holes in the most general up to-cubic-order-in-curvature dimension-independent theory of gravity that shares its graviton spectrum with the Einstein theory on constant curvature backgrounds. We show that damping rate and real oscillation frequencies of quasinormal modes for scalar, electromagnetic and Dirac fields are suppressed once the coupling with the cubic term is on. The intensity of Hawking radiation is suppressed as well, leading to, roughly, one order longer lifetime at a sufficiently large coupling constant.

Recently a $D$-dimensional regularization approach leading to the non-trivial $(3+1)$-dimensional Einstein-Gauss-Bonnet (EGB) effective description of gravity was formulated in [D. Glavan and C. Lin, Phys. Rev. Lett. 124, 081301 (2020)] which was claimed to bypasses the Lovelock's theorem and avoids Ostrogradsky instability. Later Aoki, Gorji and Mukohyama [arXiv:2005.03859] formulated a well-defined EGB theory which is four-dimensional owing to introduction of additional (scalar) degrees of freedom. A number solutions of the first naive approach are also solutions of the well-defined theory. Here we calculate quasinormal modes of scalar, electromagnetic and gravitational perturbations and find the radius of shadow for spherically symmetric and asymptotically flat black holes in this theory. We show that the black hole is gravitationally stable when ($-16 M^2<\alpha \lessapprox 0.6 M^2$). The instability in the outer range is the eikonal one and it develops at high multipole numbers. The radius of the shadow $R_{Sh}$ obeys the linear law with a remarkable accuracy.

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 spherically symmetric deformation of the Schwarzschild solution owing to the quantum corrections to gravity is known as Kazakov-Solodukhin black-hole metric. Neglecting non-spherical deformations of the background the problem was solved non-perturbatively. Here we analyze the basic characteristics of this geometry, such as: quasinormal modes and grey-body factors of fields of various spin and shadow cast by this black hole. The WKB approach and time-domain integration method, which we used for calculation of quasinormal modes, are in a good concordance. The analytical formula for quasinormal modes is deduced in the eikonal regime. The radius of shadow is decreasing when the quantum deformation is turned on.

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.

Classical (quasinormal) and quantum (Hawking) radiations are investigated for test fields in the background of a four dimensional, spherically symmetric and asymptotically flat black hole in the Einstein-dilaton-Gauss-Bonnet (EdGB) theory. The geometry of the EdGB black hole deviates from the Schwarzschild geometry only slightly. Therefore, here we observe that the quasinormal spectrum also deviates from its Schwarzschild limit at most moderately, allowing for a $9\%$ decrease in the damping rate and up to a $6\%$ decrease in the real oscillation frequency. However, the intensity of Hawking radiation of an electromagnetic field turned out to be much more sensitive characteristic than its quasinormal spectrum, allowing for a $54\%$ increase of the energy emission rate. The analytical formula for the eikonal regime of quasinormal modes is derived for test fields and it is shown that the correspondence between the eikonal quasinormal modes and null geodesics is indeed fulfilled for test fields, but is not expected for the gravitational one.

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$.

Here we shall show how to reconstruct the shape function of a spherically symmetric traversable Lorenzian wormhole near its throat if one knows high frequency quasinormal modes of the wormhole. The wormhole spacetime is given by the Morris-Thorne ansatz. The solution to the inverse problem via fitting of the parameters within the WKB approach is unique for arbitrary tideless wormholes and some wormholes with non-zero tidal effects, but this is not so for arbitrary wormholes. As examples, we reproduce the near throat geometries of the Bronnikov-Ellis and tideless Morris-Thorne metrics by their quasinormal modes at high multipole numbers $\ell$.

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.

Einstein-dilaton-Gauss-Bonnet gravity is a well-motivated alternative theory of gravity that emerges naturally from string theory. While black hole solutions have been known in this theory in numerical form for a while, an approximate analytical metric was obtained recently by some of us, which allows for faster and more detailed analysis. Here we test the accuracy of the analytical metric in the context of X-ray reflection spectroscopy. We analyze innermost stable circular orbits (ISCO) and relativistically broadened iron lines and find that both the ISCO and iron lines are determined sufficiently accurately up to the limit of the approximation. We also find that, though the ISCO increases by about 7% as dilaton coupling increases from zero to extremal values, the redshift at ISCO changes by less than 1%. Consequently, the shape of the iron line is much less sensitive to the dilaton charge than expected.

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.

In Phys. Rev. D 79, 064016 (2009) it was claimed that quasinormal modes which any stationary, spherically symmetric and asymptotically flat black hole emits in the eikonal regime are determined by the parameters of the circular null geodesic: the real and imaginary parts of the quasinormal mode are multiples of the frequency and instability timescale of the circular null geodesics respectively. We shall consider asymptotically flat black hole in the Einstein-Lovelock theory, find analytical expressions for gravitational quasinormal modes in the eikonal regime and analyze the null geodesics. Comparison of the both phenomena shows that the expected link between the null geodesics and quasinormal modes is violated in the Einstein-Lovelock theory. Nevertheless, the correspondence exists for a number of other cases and here we formulate its actual limits.

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.

Here we consider two phenomena in the vicinity of a black hole deformed by the tidal gravitational force of surrounding matter and by a strong magnetic field: equatorial motion of charged particles and the decay of a test scalar field. We were able to analyze both phenomena with analytical and simple numerical tools, which was unexpected given the low symmetry of the system. We show that both the tidal gravitational force and the magnetic field strongly enhance the release of the binding energy for the matter spiralling into the black hole. In the presence of the magnetic field, the left and right handed rotations of charged particles are not equivalent and for sufficiently large $|q| B$ there are stable anti-Larmor orbits very close to the event horizon, although Larmor orbits are only stable at some distance from the black hole. The larger the tidal force, the closer the innermost stable orbit to the black hole for both types of rotation. It was also shown that the real oscillation frequencies of the characteristic quasinormal modes are considerably suppressed by the tidal force.

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.

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.

Recently a renormalizable model of gravity has been proposed, which might be a UV completion of General Relativity (GR) or its infra-red modification, probably with a strongly coupled scalar mode. Although the generic vacuum of the theory is anti-de Sitter one, particular limits of the theory allow for the Minkowski vacuum. In this limit (though without consideration of the strongly coupled scalar field) post-Newtonian coefficients of spherically symmetric solutions coincide with those of the General Relativity. Thus the deviations from the convenient GR should be tested beyond the post-Newtonian corrections, that is for a system with strong gravity at astrophysical scales. In this letter we consider potentially observable properties of black holes in the deformed Horava-Lifshitz gravity with Minkowski vacuum: the gravitational lensing and quasinormal modes. We have showed that the bending angle is seemingly smaller in the considered Horava-Lifshitz gravity than in GR. The quasinormal modes of black holes are longer lived and have larger real oscillation frequency in the Horava-Lifshitz gravity than in GR. These corrections should be observable in the near future experiments on lensing and by gravitational antennas, helping to constrain parameters of the Horava-Lifshitz gravity or to discard it.

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.

We consider the bound states of the massive scalar field around a rotating black hole immersed in the asymptotically uniform magnetic field. In the regime of slow black hole rotation, the Klein-Gordon equation allows separation of variables. We show that the growth rate of the instability can be amplified a few times by the magnetic field. The effect occurs because the magnetic field adds the "effective mass" term $B |m|$ to the scalar field potential for a Kerr black hole. In addition, and as a by-product, we discuss the behavior of the quasinormal modes for the magnetized rotating black holes.

We found quasinormal modes, both in time and frequency domains, of the Ernst black holes, that is neutral black holes immersed in an external magnetic field. The Ernst solution reduces to the Schwarzschild solution, when the magnetic field vanishes. It is found that the quasinormal spectrum for massless scalar field in the vicinity of the magnetized black holes acquires an effective "mass" $\mu = 2 B m$, where m is the azimuthal number and B is parameter describing the magnetic field. We shall show that in the presence of a magnetic field quasinormal modes are longer lived and have larger oscillation frequencies. The perturbations of higher dimensional magnetized black holes by Ortaggio and of magnetized dilaton black holes by Radu are considered.

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.

We analyze motion of massless and massive particles around black holes immersed in an asymptotically uniform magnetic field and surrounded by some mechanical structure, which provides the magnetic field. The space-time is described by Preston-Poisson metric, which is the generalization of the well-known Ernst metric with a new parameter, tidal force, characterizing the surrounding structure. The Hamilton-Jacobi equations allow separation of variables in the equatorial plane. The presence of tidal force from surroundings considerably changes parameters of the test particle motion: it increases the radius of circular orbits of particles, increases the binding energy of massive particles going from a given circular orbits to the innermost stable orbit near black hole. In addition, it increases the distance of minimal approach, time delay and bending angle for a ray of light propagating near black hole.

We use the Ernst-Schwarzschild solution for a black hole immersed in a uniform magnetic field to estimate corrections to the bending angle and time delay due-to presence of weak magnetic fields in galaxies and between galaxies, and also due-to influence of strong magnetic field near supermassive black holes. The magnetic field creates a kind of confinement in space, that leads to increasing of the bending angle and time delay for a ray of light propagating in the equatorial plane.

We present here a detailed study of the quasi-normal spectrum of brane-localised Standard Model fields in the vicinity of D-dimensional black-holes. A variety of such backgrounds (Schwarzschild, Reissner-Nordstrom and Schwarzszchild-(Anti) de Sitter) are investigated. The dependence of the quasi-normal spectra on the dimensionality D, spin of the field s, and multipole number l is analyzed. Analytical formulae are obtained for a number of limiting cases: in the limit of large multipole number for Schwarzschild, Schwarzschild-de Sitter and Reissner-Nordstrom black holes, in the extremal limit of the Schwarzschild-de Sitter black hole, and in the limit of small horizon radius in the case of Schwarzschild-Anti de Sitter black holes. We show that an increase in the number of hidden, extra dimensions results in the faster damping of all fields living on the brane, and that the localization of fields on a brane affects the QN spectrum in a number of additional ways, both direct and indirect.

It is presented a thorough analysis of scalar perturbations in the background of Gauss-Bonnet, Gauss-Bonnet-de Sitter and Gauss-Bonnet-anti-de Sitter black hole spacetimes. The perturbations are considered both in frequency and time domain. The dependence of the scalar field evolution on the values of the cosmological constant $\Lambda$ and the Gauss-Bonnet coupling $\alpha$ is investigated. For Gauss-Bonnet and Gauss-Bonnet-de Sitter black holes, at asymptotically late times either power-law or exponential tails dominate, while for Gauss-Bonnet-anti-de Sitter black hole, the quasinormal modes govern the scalar field decay at all times. The power-law tails at asymptotically late times for odd-dimensional Gauss-Bonnet black holes does not depend on $\alpha$, even though the black hole metric contains $\alpha$ as a new parameter. The corrections to quasinormal spectrum due to Gauss-Bonnet coupling is not small and should not be neglected. For the limit of near extremal value of the (positive) cosmological constant and pure de Sitter and anti-de Sitter modes in Gauss-Bonnet gravity we have found analytical expressions.

We investigate the response of the traversable wormholes to the external perturbations through finding their characteristic frequencies and time-domain profiles. The considered solution describes traversable wormholes between the branes in the two brane Randall-Sundrum model and was previously found within Einstein gravity with a conformally coupled scalar field. The evolution of perturbations of a wormhole is similar to that of a black hole and represents damped oscillations (ringing) at intermediately late times, which are suppressed by power law tails (proportional to t^{-2} for monopole perturbations) at asymptotically late times.

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.

We present a thorough analysis for the quasinormal (QN) behavior, associated with the decay of scalar, electromagnetic and gravitational perturbations, of Schwarzschild-anti-de Sitter black holes. As it is known the anti-de Sitter (AdS) QN spectrum crucially depends on the relative size of the black hole to the AdS radius. There are three different types of behavior depending on whether the black hole is large, intermediate, or small. The results of previous works, concerning lower overtones for large black holes, are completed here by obtaining higher overtones for all the three black hole regimes. There are two major conclusions that one can draw from this work: First, asymptotically for high overtones, all the modes are evenly spaced, and this holds for all three types of regime, large, intermediate and small black holes, independently of l, where l is the quantum number characterizing the angular distribution; Second, the spacing between modes is apparently universal, in that it does not depend on the field, i.e., scalar, electromagnetic and gravitational QN modes all have the same spacing for high overtones. We are also able to prove why scalar and gravitational perturbations are isospectral, asymptotically for high overtones, by introducing appropriate superpartner potentials.