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.

The necessary information to distinguish a local inhomogeneous mass density field from its spatial average on a compact domain of the universe can be measured by relative information entropy. The Kullback-Leibler (KL) formula arises very naturally in this context, however, it provides a very complicated way to compute the mutual information between spatially separated but causally connected regions of the universe in a realistic, inhomogeneous model. To circumvent this issue, by considering a parametric extension of the KL measure, we develop a simple model to describe the mutual information which is entangled via the gravitational field equations. We show that the Tsallis relative entropy can be a good approximation in the case of small inhomogeneities, and for measuring the independent relative information inside the domain, we propose the R\'enyi relative entropy formula.

We present a new approach to cosmological perturbations based on the theory of Lie groups and their representations. After re-deriving the standard covariant formalism from SO(3) considerations, we provide a new expansion of the perturbed Friedmann-Lemaitre-Robertson-Walker (FLRW) metric in terms of irreducible representations of the Lorentz group. The resulting decomposition splits into (scalar, scalar), (scalar, vector) and (vector, vector) terms. These equations directly correspond to the standard Lifshitz classification of cosmological perturbations using scalar, vector and tensor modes which arise from the irreducible SO(3) representation of the spatial part of the metric. While the Lorentz group basis matches the underlying local symmetries of the FLRW spacetime better than the SO(3), the new equations do not provide further simplification compared to the standard cosmological perturbation theory. We conjecture that this is due to the fact that the so(3,1) ~ su(2) x su(2) Lorentz algebra has no pair of commuting generators commuting with any of the translation group generators.

We consider general-relativistic rotational perturbations in homogeneous and isotropic Friedman - Robertson - Walker (FRW) cosmologies. Taking linear perturbations of FRW models, the general solution of the field equations contains tensorial, vectorial and scalar functions. The vectorial terms are in connection with rotations in the given model and due to the Sachs - Wolfe effect they produce contributions to the temperature fluctuations of the cosmic microwave background radiation (CMBR). In present paper we obtain the analytic time dependence of these contributions in a spatially flat, FRW model with pressureless ideal fluid, in the presence and the absence of a cosmological constant. We find that the solution can be separated into an integrable and a non-integrable part as is the situation in the case of scalar perturbations. Analyzing the solutions and using the results of present observations we estimate the order of magnitude of the angular velocity corresponding to the rotation tensor at the time of decoupling and today.

The integrated Sachs-Wolfe (ISW) effect in a Lambda dominated universe can be an important factor in the evolution of cosmic microwave background fluctuations. With the inclusion of cosmological constant we present the complete analytic solution of the covariant linear perturbations of the Einstein equations in the Newtonian gauge, in the case of a spatially flat (K=0), Friedman-Robertson-Walker (FRW) universe, filled with pressureless ideal fluid. We use the analytic time dependence of the perturbation potentials to derive the anisotropy power spectrum of the late ISW effect. We choose the scale invariant Harrison-Zeldovich spectrum for obtaining the C_l momenta and show the agreement of our results with earlier, numeric calculations.

We present the complete solution of the first order metric and density perturbation equations in a spatially flat (K=0), Friedmann-Robertson-Walker (FRW) universe filled with pressureless ideal fluid, in the presence of cosmological constant. We use covariant linear perturbation formalism and the comoving gauge condition to obtain the field and conservation equations. The solution contains all modes of the perturbations, i.e. scalar, vector and tensor modes, and we show that our results are in agreement with the Sachs & Wolfe metric perturbation formalism.

Spatially homogeneous and isotropic cosmological models, with a perfect fluid matter source and non-vanishing cosmological constant, are studied. The equations governing linear perturbations of the space-time and the variation of energy density are given. The complete solution of the problem is obtained for $C^{\infty}$ perturbations, using a comoving time. The Sachs-Wolfe fluctuations of the temperature of the cosmic background radiation are obtained for the relatively growing density perturbations. It is found that the observable celestial microwave fluctuation pattern underwent a reversal approximately two billion years ago. What is observed today is a negative image of the last scattering surface with an attenuation of the fluctuations, due to the presence of the cosmological constant.

We point out that, due to the nonlinearity of the Einstein equations, a homogeneous approximation in cosmology leads to the appearance of an additional term in the Friedmann equation. This new term is associated with the spatial inhomogeneities of the metric and can be expressed in terms of density fluctuations. Although it is not constant, it decays much slower (as t^(-2/3)) than the other terms (like density) which decrease as t^(-2). The presence of the new term leads to a correction in the scale factor that is proportional to t^2 and may give account of the recently observed accelerating expansion of the universe without introducing a cosmological constant.