Within General Relativity, the unique stationary solution of an isolated black hole is the Kerr spacetime, which has a peculiar multipolar structure depending only on its mass and spin. We develop a general method to extract the multipole moments of arbitrary stationary spacetimes and apply it to a large family of horizonless microstate geometries. The latter can break the axial and equatorial symmetry of the Kerr metric and have a much richer multipolar structure, which provides a portal to constrain fuzzball models phenomenologically. We find numerical evidence that all multipole moments are always larger (in absolute value) than those of a Kerr black hole with the same mass and spin. Current measurements of the quadrupole moment of black-hole candidates could place only mild constraints on fuzzballs, while future gravitational-wave detections of extreme mass-ratio inspirals with the space mission LISA will improve these bounds by orders of magnitude.

Einstein's theory of general relativity predicts that the only stationary configuration of an isolated black hole is the Kerr spacetime, which has a unique multipolar structure and a spherical shape when non-spinning. This is in striking contrast to the case of other self-gravitating objects, which instead can in principle have arbitrary deformations even in the static case. Here we develop a general perturbative framework to construct stationary stars with small axisymmetric deformations, and study explicitly compact stars with an intrinsic quadrupole moment. The latter can be sustained, for instance, by crust stresses or strong magnetic fields. While our framework is general, we focus on quadrupolar deformations of neutron stars induced by an anisotropic crust, which continuously connect to spherical neutron stars in the isotropic limit. Deformed neutron stars might provide a more accurate description for stellar remnants formed in supernovae and in binary mergers, and can be used to improve constraints on the neutron-star equation of state through gravitational-wave detections and through the observation of low-mass X-ray binaries. We argue that, if the (dimensionless) intrinsic quadrupole moment is of a few percent or higher, the effect of the deformation is stronger than that of tidal interactions in coalescing neutron-star binaries, and might also significantly affect the electromagnetic signal from accreting neutron stars.

Anisotropic stresses are ubiquitous in nature, but their modeling in General Relativity is poorly understood and frame dependent. We introduce the first study on the dynamical properties of anisotropic self-gravitating fluids in a covariant framework. Our description is particularly useful in the context of tests of the black hole paradigm, wherein ultracompact objects are used as black hole mimickers but otherwise lack a proper theoretical framework. We show that: (i) anisotropic stars can be as compact and as massive as black holes, even for very small anisotropy parameters; (ii) the nonlinear dynamics of the 1+1 system is in good agreement with linearized calculations, and shows that configurations below the maximum mass are nonlinearly stable; (iii) strongly anisotropic stars have vanishing tidal Love numbers in the black-hole limit; (iv) their formation will usually be accompanied by gravitational-wave echoes at late times.

Motivated by the lack of a general parametrization for exotic compact objects, we construct a class of perturbative solutions valid for small (but otherwise generic) multipolar deviations from a Schwarzschild metric in general relativity. We introduce two classes of exotic compact objects, with "soft" and "hard" hair, for which the curvature at the surface is respectively comparable to or much larger than that at the corresponding black-hole horizon. We extend the Hartle-Thorne formalism to relax the assumption of equatorial symmetry and to include deformations induced by multipole moments higher than the spin, thus constructing the most general, axisymmetric quasi-Schwarzschild solution to Einstein's vacuum equations. We explicitly construct several particular solutions of objects with soft hair, which might be useful for tests of quasi-black-hole metrics, and also to study deformed neutron stars. We show that the more compact a soft exotic object is, the less hairy it will be. All its multipole moments can approach their corresponding Kerr values only in two ways as their compactness increases: either logarithmically (or faster) if the moments are spin-induced, or linearly (or faster) otherwise. Our results suggest that it is challenging (but possibly feasible with next-generation gravitational-wave detectors) to distinguish Kerr black holes from a large class of ultracompact exotic objects on the basis of their different multipolar structure.

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 tidal Love numbers (TLNs) encode the deformability of a self-gravitating object immersed in a tidal environment and depend significantly both on the object's internal structure and on the dynamics of the gravitational field. An intriguing result in classical general relativity is the vanishing of the TLNs of black holes. We extend this result in three ways, aiming at testing the nature of compact objects: (i) we compute the TLNs of exotic compact objects, including different families of boson stars, gravastars, wormholes, and other toy models for quantum corrections at the horizon scale. In the black-hole limit, we find a universal logarithmic dependence of the TLNs on the location of the surface; (ii) we compute the TLNs of black holes beyond vacuum general relativity, including Einstein-Maxwell, Brans-Dicke and Chern-Simons gravity; (iii) We assess the ability of present and future gravitational-wave detectors to measure the TLNs of these objects, including the first analysis of TLNs with LISA. Both LIGO, ET and LISA can impose interesting constraints on boson stars, while LISA is able to probe even extremely compact objects. We argue that the TLNs provide a smoking gun of new physics at the horizon scale, and that future gravitational-wave measurements of the TLNs in a binary inspiral provide a novel way to test black holes and general relativity in the strong-field regime.