The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We disclose remarkable features of the scalar-tensor theory with the derivative coupling of the scalar field to the curvature in the Palatini formalism. Using the disformal transformations, we show that this theory is free from Otrogradski ghosts. For a special relation between two coupling constants, it is disformally dual to the Einstein gravity minimally coupled to the scalar, which opens the way to constructing several exact solutions. The disformal transformation degenerates at the boundary of the physical region, near which the properties of desingularization are revealed, illustrated by exact solutions: non-singular accelerating cosmology and a static spherically symmetric geon. We also construct the exact pp-waves of this theory propagating at the speed of light.

Recently double black hole vacuum and electrovacuum metrics attracted attention as exact solutions suitable for visualization of ultra-compact objects beyond the Kerr paradigm. However, many of the proposed systems are plagued with ring curvature singularities. Here we present a new simple solution of this type which is asymptotically Kerr, has zero electric and magnetic charges, but is endowed with magnetic dipole moment and electric quadrupole moment. It is manifestly free of ring singularities, and contains only a mild string-like singularity on the axis corresponding to a distributional energy-momentum tensor.Its main constituents are two extreme co-rotating black holes carrying equal electric andopposite magnetic and NUT charges.

In many astrophysically relevant situations radiation reaction force acting upon a charge can not be neglected and the question arises about the location and stability of circular orbits in such regime. Motion of point charge with radiation reaction in flat spacetime is described by Lorenz-Dirac (LD) equation, while in curved spacetime -- by DeWitt-Brehme (DWB) equation containing the Ricci term and the tail term. We show that for the motion of elementary particles in vacuum metrics the DWB equation can be reduced to the covariant form of the LD equation which we use here. Generically, the LD equation is plagued by runaway solutions, so we discuss computational ways to avoid this problem in constructing numerical solutions. We also use the first iteration of the covariant LD equation which is the covariant Landau-Lifshitz equation, comparing results of these two approaches and showing smallness of the third-order Schott term in the ultrarelativistic case. We calculate the corresponding energy and angular momentum loss of a particle and study the damping of charged particle oscillations around an equilibrium radius. We find that depending on the orientation of the Lorentz force, the oscillating charged particle either spirals down to the black hole, or stabilizes the circular orbit by decaying its oscillations. The later case leads to an interesting new result of shifting of the particle orbit outwards from the black hole. We also discuss the astrophysical relevance of the presented approach and provide estimations of the main parameters of the model.

We present a one-parameter family of stationary, asymptotically flat solutions of the Einstein-Maxwell equations with only a mild singularity, which are endowed with mass, angular momentum, a dipole magnetic moment and a quadrupole electric moment. We briefly analyze the structure of this solution, which we interpret as a system of two extreme co-rotating black holes with equal masses and electric charges, and opposite magnetic and gravimagnetic charges, held apart by an electrically charged, magnetized string which also acts as a Dirac-Misner string.

We study the massless limit in synchrotron radiation and one-photon pair creation in magnetic field. In this limit Schwinger critical field $H_0=m^2c^3/(e\hbar)$ tends to zero, so two characteristic quantum parameters $\eta=H/H_0,\;\chi=\eta E/mc^2$ are infinite, and the standard approximations used in analytical calculations fail. Applying Schwinger's proper time methods we derive simple expressions for synchrotron radiation spectra emitted by massless charges of spins $s=0,\,1/2$ and the pair creation probability distribution in the quasiclassical (high Landau levels) regime exhibiting simple scaling properties and possessing universal spectral shapes.

We show that supercritically charged black holes with NUT provide a new setting for traversable wormholes. This does not require exotic matter, a price being the Misner string singularities. Without assuming time periodicity to make Misner strings unobservable, we show that, contrary to expectations, geodesics do not stop there. Moreover, since there is no central singularity the space-time turns out to be geodesically complete. Another unpleasant feature of spacetimes with NUTs is the presence of regions where the azimuthal angle $\varphi$ becomes timelike, signalling the appearance of closed timelike curves (CTCs). We show that among them there are no closed timelike or null geodesics, so the freely falling observers should not encounter causality violations. Considering worldlines of charged particles, we find that, although these can become closed in the vicinity of the wormhole throat for large enough charge-to-mass ratio, the non-causal orbits are still disconnected from the distant zones. All these findings support our feeling that wormholes with NUTs deserve to be taken seriously. Integrating the geodesic equations completely, we demonstrate the existence of timelike and null geodesics connecting two asymptotic regions of the wormhole, such that the tidal forces in the throat are reasonably small. We discuss bounds on the NUT charge which follow from the Schwinger pair creation and ionization thresholds and speculate that such NUT wormholes could be present in some galactic centers.

We demonstrate that full description of both electromagnetic and gravitational radiation from massless particles lies outside the scope of classical theory. Synchrotron radiation from the hypothetical massless charge in quantum electrodynamics in external magnetic field has finite total power while the corresponding classical formula diverges in the massless limit. We argue that in both cases classical theory describes correctly only the low-frequency part of the spectra, while the total power diverges because of absence of the UV frequency cutoff. Failure of description of gravitational radiation from massless particles by classical General Relativity may be considered as another appeal for quantization of gravity apart from the problem of singularities.

Classical radiation power from an accelerated massive charge diverges in the zero-mass limit, while some general arguments suggest that strictly massless charge does not not radiate at all. On the other hand, the regularized classical radiation reaction force, though looking odd, is non-zero and finite. To clarify this controversy, we consider radiation problem in massless scalar quantum electrodynamics in the external magnetic field. In this framework, synchrotron radiation is found to be non-zero, finite, and essentially quantum. Its spectral distribution is calculated using Schwinger's proper time technique for {\em ab initio} massless particle of zero spin. Provided $E^2\gg eH$, the maximum in the spectrum is shown to be at $\hbar \omega=E/3$, and the average photon energy is $4E/9$. The normalized spectrum is universal, depending neither on $E$ nor on $H$. Quantum nature of radiation makes classical radiation reaction equation meaningless for massless charge. Our results are consistent with the view (supported by the renormalization group argument) that the correct classical limit of massless quantum electrodynamics is free theory.

We discuss cosmological models involving homogeneous and isotropic Yang-Mills (YM) fields. Such models were proposed recently as an alternative to scalar models of cosmic acceleration. There exists a unique SU(2) YM configuration (generalizable to larger gauge groups) whose energy-momentum tensor is homogeneous and isotropic in space. It is parameterized by a single scalar field with a quatric potential. In the case of the closed universe the coupled YM -- doublet Higgs system admits homogeneous and isotropic configurations too. While pure Einstein-Yang-Mills (EYM) cosmology with the standard conformally invariant YM action gives rise to the hot universe, Einstein-Yang-Mills-Higgs (EYMH) cosmology has a variety of regimes which include inflationary stages, bounces, and exhibits global cycling behavior reminiscent of the Multiverse developed in time. We also discuss other mechanisms of conformal symmetry breaking such as string-inspired Born-Infeld (BI) modification of the YM action or field-theoretical quantum corrections.

We review basic features of cosmological models with large-scale classical non-Abelian Yang-Mills (YM) condensates. There exists a unique SU(2) YM configuration (generalizable to larger gauge groups) compatible with homogeneity and isotropy of the three-space which is parameterized by a single scalar field. In the past various aspects of Einstein-Yang-Mills (EYM) cosmology were discussed in the context of the Early Universe. Due to conformal invariance, solvable EYM FRW models exist both on the classical and quantum levels. To develop the YM model for dark energy one has to find mechanisms of the conformal symmetry breaking. We discuss the Born-Infeld generalization and some phenomenological models motivated by quantum corrections exploring possibility of transient DE and phantom regimes.

Contrary to a recent claim by Anderson ["The Mathematical Theory of Cosmic Strings", I.O.P. Publishing, Bristol 2003], we show that the Letelier-Gal'tsov metric does represent a system of crossed straight infinite cosmic strings moving with arbitrary constant velocities.

We calculate axion radiation emitted in the collision of two straight global strings. The strings are supposed to be in the unexcited ground state, to be inclined with respect to each other, and to move in parallel planes. Radiation arises when the point of minimal separation between the strings moves faster than light. This effect exhibits a typical Cerenkov nature. Surprisingly, it allows an alternative interpretation as bremsstrahlung under a collision of point charges in 2+1 electrodynamics. This can be demonstrated by suitable world-sheet reparameterizations and dimensional reduction. Cosmological estimates show that our mechanism generates axion production comparable with that from the oscillating string loops and may lead to further restrictions on the axion window.

A new mechanism of cosmic axion production is proposed: axion bremsstrahlung from collisions of straight global strings. This effect is of the second order in the axion coupling constant, but the resulting cosmological estimate is likely to be of the same order as that corresponding to radiation from oscillating string loops. This may lead to a further restriction on the axion window.