Starting from a recently constructed stealth Kerr solution of higher order scalar tensor theory involving scalar hair, we analytically construct disformal versions of the Kerr spacetime with a constant degree of disformality and a regular scalar field. While the disformed metric has only a ring singularity and asymptotically is quite similar to Kerr, it is found to be neither Ricci flat nor circular. Non-circularity has far reaching consequences on the structure of the solution. As we approach the rotating compact object from asymptotic infinity we find a static limit ergosurface similar to the Kerr spacetime with an enclosed ergoregion. However, the stationary limit of infalling observers is found to be a timelike hypersurface. A candidate event horizon is found in the interior of this stationary limit surface. It is a null hypersurface generated by a null congruence of light rays which are no longer Killing vectors. Under a mild regularity assumption, we find that the candidate surface is indeed an event horizon and the disformed Kerr metric is therefore a black hole quite distinct from the Kerr solution.

In a subclass of scalar-tensor theories, it has been shown that standard general relativity solutions of neutron stars and black holes with trivial scalar field profiles are unstable. Such an instability leads to solutions which are different from those of general relativity and have non-trivial scalar field profiles, in a process called scalarization. In the present work we focus on scalarization due to a non-minimal coupling of the scalar field to the Gauss-Bonnet curvature invariant. The coupling acts as a tachyonic mass for the scalar mode, thus leading to the instability of general relativity solutions. We point out that a similar effect may occur for the scalar modes in a cosmological background, resulting in the instability of cosmological solutions. In particular, we show that a catastrophic instability develops during inflation within a period of time much shorter than the minimum required duration of inflation. As a result, the standard cosmological dynamics is not recovered. This raises the question of the viability of scalar-Gauss-Bonnet theories exhibiting scalarization.

We develop a new technique for finding black hole solutions in modified gravity that have "stealth" hair, i.e., hair whose only gravitational effect is to tune the cosmological constant. We consider scalar-tensor theories in which gravitational waves propagate at the speed of light, and show that Einstein metrics can be painted with stealth hair provided there exists a family of geodesics always normal to spacelike surfaces. We also present a novel scalar-dressed rotating black hole that has finite scalar field at both the black hole and cosmological event horizons.

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 gravitational wave event GW170817 together with its electromagnetic counterparts constrains the speed of gravity to be extremely close to that of light. We first show, on the example of an exact Schwarzschild-de Sitter solution of a specific beyond-Horndeski theory, that imposing the strict equality of these speeds in the asymptotic homogeneous Universe suffices to guarantee so even in the vicinity of the black hole, where large curvature and scalar-field gradients are present. We also find that the solution is stable in a range of the model parameters. We finally show that an infinite class of beyond-Horndeski models satisfying the equality of gravity and light speeds still provide an elegant self-tuning: The very large bare cosmological constant entering the Lagrangian is almost perfectly counterbalanced by the energy-momentum tensor of the scalar field, yielding a tiny observable effective cosmological constant.

We consider theories describing the dynamics of a four-dimensional metric, whose Lagrangian is diffeomorphism invariant and depends at most on second derivatives of the metric. Imposing degeneracy conditions we find a set of Lagrangians that, apart form the Einstein-Hilbert one, are either trivial or contain more than two degrees of freedom. Among the partially degenerate theories, we recover Chern-Simons gravity, endowed with constraints whose structure suggests the presence of instabilities. Then, we enlarge the class of parity violating theories of gravity by introducing new "chiral scalar-tensor theories". Although they all raise the same concern as Chern-Simons gravity, they can nevertheless make sense as low energy effective field theories or, by restricting them to the unitary gauge (where the scalar field is uniform), as Lorentz breaking theories with a parity violating sector.

We show that the two additional Lagrangians that appear in theories beyond Horndeski can be reexpressed in terms of simple generalizations of the "John" and "Paul" terms of the Fab Four theories. We find that these extended Fab Four satisfy the same properties of self-tuning as the original Fab Four.

We have recently proposed a special class of scalar tensor theories known as the Fab Four. These arose from attempts to analyse the cosmological constant problem within the context of Horndeski's most general scalar tensor theory. The Fab Four together give rise to a model of self-tuning, with the relevant solutions evading Weinberg's no-go theorem by relaxing the condition of Poincare invariance in the scalar sector. The Fab Four are made up of four geometric terms in the action with each term containing a free potential function of the scalar field. In this paper we rigorously derive this model from the general model of Horndeski, proving that the Fab Four represents the only classical scalar tensor theory of this type that has any hope of tackling the cosmological constant problem. We present the full equations of motion for this theory, and give an heuristic argument to suggest that one might be able to keep radiative corrections under control. We also give the Fab Four in terms of the potentials presented in Deffayet et al's version of Horndeski.

Starting from the most general scalar-tensor theory with second order field equations in four dimensions, we establish the unique action that will allow for the existence of a consistent self-tuning mechanism on FLRW backgrounds, and show how it can be understood as a combination of just four base Lagrangians with an intriguing geometric structure dependent on the Ricci scalar, the Einstein tensor, the double dual of the Riemann tensor and the Gauss-Bonnet combination. Spacetime curvature can be screened from the net cosmological constant at any given moment because we allow the scalar field to break Poincar\'e invariance on the self-tuning vacua, thereby evading the Weinberg no-go theorem. We show how the four arbitrary functions of the scalar field combine in an elegant way opening up the possibility of obtaining non-trivial cosmological solutions.

In order for a modified gravity model to be a candidate for cosmological dark energy it has to pass stringent local gravity experiments. We find that a Brans-Dicke (BD) theory with well-defined second order corrections that include the Gauss-Bonnet term possess this feature. We construct the generic second order theory that gives, to linear order, a BD metric solution for a point-like mass source. We find that these theories interpolate between general relativity (GR) and BD gravity. In particular it is found that the relevant Eddington parameter, that is commonly heavily constrained by time delay experiments, can be arbitrarily close to the GR value of 1, with an arbitrary BD parameter. We find the region where the solution is stable to small timelike perturbations.

We discuss a particular higher order gravity theory, Lovelock theory, that generalises in higher dimensions, general relativity. After briefly motivating modifications of gravity, we will introduce the theory in question and we will argue that it is a unique, mathematically sensible, and physically interesting extension of general relativity. We will see, by using the formalism of differential forms, the relation of Lovelock gravity to differential geometry and topology of even dimensional manifolds. We will then discuss a generic staticity theorem, which will give us the charged static black hole solutions. We will examine their asymptotic behavior, analyse their horizon structure and briefly their thermodynamics. We will then examine the distributional matching conditions for Lovelock theory. We will see how induced 4 dimensional Einstein-Hilbert terms result on the brane geometry from the higher order Lovelock terms. With the junction conditions at hand, we will go back to the black hole solutions and give applications for braneworlds: perturbations of codimension 1 braneworlds and the exact solution for braneworld cosmology as well as the determination of maximally symmetric codimension 2 braneworlds. In both cases, the staticity theorem evoked beforehand will give us the general solution for braneworld cosmology in codimension 1 and maximal symmetry warped branes of codimension 2. We will then end with a discussion of the simplest Kaluza-Klein reduction of Lovelock theory to a 4 dimensional vector-scalar-tensor theory which has the unique property of retaining second order field equations. We will conclude by listing some open problems and common difficulties.

Although the Gauss-Bonnet term is a topological invariant for general relativity, it couples naturally to a quintessence scalar field, modifying gravity at solar system scales. We determine the solar system constraints due to this term by evaluating the post-Newtonian metric for a distributional source. We find a mass dependent, 1/r^7 correction to the Newtonian potential, and also deviations from the Einstein gravity prediction for light-bending. We constrain the parameters of the theory using planetary orbits, the Cassini spacecraft data, and a laboratory test of Newton's law, always finding extremely tight bounds on the energy associated to the Gauss-Bonnet term. We discuss the relevance of these constraints to late-time cosmological acceleration.

We show how to construct consistent braneworld models which exhibit late time acceleration. Unlike self-acceleration, which has a de Sitter vacuum state, our models have the standard Minkowski vacuum and accelerate only in the presence of matter, which we dub ``stealth-acceleration''. We use an effective action for the brane which includes an induced gravity term, and allow for an asymmetric set-up. We study the linear stability of flat brane vacua and find the regions of parameter space where the set-up is stable. The 4-dimensional graviton is only quasi-localised in this set-up and as a result gravity is modified at late times. One of the two regions is strongly coupled and the scalar mode is eaten up by an extra symmetry that arises in this limit. Having filtered the well-defined theories we then focus on their cosmology. When the graviton is quasi-localised we find two main examples of acceleration. In each case, we provide an illustrative model and compare it to LambdaCDM.

Models with a scalar field coupled to the Gauss-Bonnet Lagrangian appear naturally from Kaluza-Klein compactifications of pure higher-dimensional gravity. We study linear, cosmological perturbations in the limits of weak coupling and slow-roll, and derive simple expressions for the main observable sub-horizon quantities: the anisotropic stress factor, the time-dependent gravitational constant, and the matter perturbation growth factor. Using present observational data, and assuming slow-roll for the dark energy field, we find that the fraction of energy density associated with the coupled Gauss-Bonnet term cannot exceed 15%. The bound should be treated with caution, as there are significant uncertainies in the data used to obtain it. Even so, it indicates that the future prospects for constraining the coupled Gauss-Bonnet term with cosmological observations are encouraging.

We systematically explore the spectrum of gravitational perturbations in codimension-1 DGP braneworlds, and find a 4D ghost on the self-accelerating branch of solutions. The ghost appears for any value of the brane tension, although depending on the sign of the tension it is either the helicity-0 component of the lightest localized massive tensor of mass $0<m^2 < 2H^2$ for positive tension, the scalar `radion' for negative tension, or their admixture for vanishing tension. Because the ghost is gravitationally coupled to the brane-localized matter, the self-accelerating solutions are not a reliable benchmark for cosmic acceleration driven by gravity modified in the IR. In contrast, the normal branch of solutions is ghost-free, and so these solutions are perturbatively safe at large distance scales. We further find that when the $\mathbb{Z}_2$ orbifold symmetry is broken, new tachyonic instabilities, which are much milder than the ghosts, appear on the self-accelerating branch. Finally, using exact gravitational shock waves we analyze what happens if we relax boundary conditions at infinity. We find that non-normalizable bulk modes, if interpreted as 4D phenomena, may open the door to new ghost-like excitations.

Dark matter caustics have specific density profiles and, therefore, precisely calculable gravitational lensing properties. We present a formalism which simplifies the relevant calculations, and apply it to four specific cases. In the first three, the line of sight is tangent to a smooth caustic surface. The curvature of the surface at the tangent point is positive, negative or zero. In the fourth case the line of sight passes near a cusp. For each we derive the map between the image and source planes. In some cases, a point source has multiple images and experiences infinite magnification when the images merge. Unfortunately, for the dark matter caustics expected in realistic galactic halo models, the angular resolution required to resolve the multiple images is not presently achievable. A more promising approach aims to observe the distortions caused by dark matter caustics in the images of extended sources such as radio jets.

We investigate the spacetime of a thick gravitating domain wall for a general potential $V(\Phi)$. Using general analytical arguments we show that all nontrivial solutions fall into two categories: those interpretable as an isolated domain wall with a cosmological event horizon, and those which are pure false vacuum de Sitter solutions. Although this latter solution is always unstable to the field rolling coherently to its true vacuum, we show that there is an additional instability to wall formation if the scalar field does not couple too strongly to gravity. Using the $\lambda \Phi^4$ and sine-Gordon models as illustrative examples, we investigate the phase space of the gravitating domain wall in detail numerically, following the solutions from weak to strong gravity. We find excellent agreement with the analytic work. Then, we analyse the domain wall in the presence of a cosmological constant finding again the two kinds of solutions, wall and de Sitter, even in the presence of a negative cosmological constant.