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. This write-up is an initiative taken within the framework of the European Action on ‘Black holes, Gravitational waves and Fundamental Physics’.

Effective field theory methods suggest that some rather general extensions of general relativity include, or are mimicked by, certain higher-order curvature corrections, with coupling constants expected to be small but otherwise arbitrary. Thus, the tantalizing prospect to test the fundamental nature of gravity with gravitational-wave observations, in a systematic way, emerges naturally. Here, we build black hole solutions in such a framework and study their main properties. Once rotation is included, we find the first purely gravitational example of geometries without ${\mathbb{Z}}_2$ symmetry. Despite the higher-order operators of the theory, we show that linearized fluctuations of such geometries obey second-order differential equations. We find nonzero tidal Love numbers. We study and compute the quasinormal modes of such geometries. These results are of interest to gravitational-wave science but also potentially relevant for electromagnetic observations of the galactic center or x-ray binaries.

Black hole solutions in general relativity are simple. The frequency spectrum of linear perturbations around these solutions (i.e., the quasinormal modes) is also simple, and therefore it is a prime target for fundamental tests of black hole spacetimes and of the underlying theory of gravity. The following technical calculations must be performed to understand the imprints of any modified gravity theory on the spectrum: 1. Identify a healthy theory; 2. Find black hole solutions within the theory; 3. Compute the equations governing linearized perturbations around the black hole spacetime; 4. Solve these equations to compute the characteristic quasinormal modes. In this work (the first of a series) we assume that the background spacetime has spherical symmetry, that the relevant physics is always close to general relativity, and that there is no coupling between the perturbation equations. Under these assumptions, we provide the general numerical solution to step 4. We provide publicly available data files such that the quasinormal modes of any spherically symmetric spacetime can be computed (in principle) to arbitrary precision once the linearized perturbation equations are known. We show that the isospectrality between the even- and odd-parity quasinormal mode spectra is fragile, and we identify the necessary conditions to preserve it. Finally, we point out that new modes can appear in the spectrum even in setups that are perturbatively close to general relativity.

Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work [ Phys. Rev. D 99, 104077 (2019)] we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model.

We derive a general formula for the center-of-mass (CM) energy for the near-horizon collision of two particles of the same rest mass on the equatorial plane around a Kerr black hole. We then apply this formula to a particle which plunges from the innermost stable circular orbit (ISCO) and collides with another particle near the horizon. It is found that the maximum value of the CM energy Ecm is given by Ecm/(2m0)$\simeq${}1.40/1-a*24 for a nearly maximally rotating black hole, where m0 is the rest mass of each particle and a* is the nondimensional Kerr parameter. This coincides with the known upper bound for a particle which begins at rest at infinity within a factor of 2. Moreover, we also consider the collision of a particle orbiting the ISCO with another particle on the ISCO and find that the maximum CM energy is then given by Ecm/(2m0)$\simeq${}1.77/1-a*26. In view of the astrophysical significance of the ISCO, this result implies that particles can collide around a rotating black hole with an arbitrarily high CM energy without any artificial fine-tuning in an astrophysical context if we can take the maximal limit of the black hole spin or a*$\rightarrow${}1. On the other hand, even if we take Thorne's bound on the spin parameter into account, highly or moderately relativistic collisions are expected to occur quite naturally, for Ecm/(2m0) takes 6.95 (maximum) and 3.86 (generic) near the horizon and 4.11 (maximum) and 2.43 (generic) on the ISCO for a*=0.998. This implies that high-velocity collisions of compact objects are naturally expected around a rapidly rotating supermassive black hole. Implications to accretion flows onto a rapidly rotating black hole are also discussed.