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’.

The idea of stable, localized bundles of energy has strong appeal as a model for particles. In the 1950s, John Wheeler envisioned such bundles as smooth configurations of electromagnetic energy that he called geons, but none were found. Instead, particle-like solutions were found in the late 1960s with the addition of a scalar field, and these were given the name boson stars. Since then, boson stars find use in a wide variety of models as sources of dark matter, as black hole mimickers, in simple models of binary systems, and as a tool in finding black holes in higher dimensions with only a single Killing vector. We discuss important varieties of boson stars, their dynamic properties, and some of their uses, concentrating on recent efforts.

A general covariant extension of Einstein´s field equations is considered with a view to Numerical Relativity applications. The basic variables are taken to be the metric tensor and an additional four-vector Zµ. Einstein's solutions are recovered when the additional four-vector vanishes, so that the energy and momentum constraints amount to the covariant algebraic condition Zµ = 0. The extended field equations can be supplemented by suitable coordinate conditions in order to provide symmetric hyperbolic evolution systems: this is actually the case for either harmonic coordinates or normal coordinates with harmonic slicing.

Solitonic boson stars are stable objects made of a complex scalar field with a compactness that can reach values comparable to that of neutron stars. A recent study of the collision of identical boson stars produced only non-rotating boson stars or black holes, suggesting that rotating boson stars may not form from binary mergers. Here we extend this study to include an analysis of the gravitational waves radiated during the coalescence of such a binary, which is crucial to distinguish these events from other binaries with LIGO and Virgo observations. Our studies reveal that the remnant's gravitational wave signature is mainly governed by its fundamental frequency as it settles down to a non-rotating boson star, emitting significant gravitational radiation during this postmerger state. We calculate how the waveforms and their post-merger frequencies depend on the compactness of the initial boson stars and estimate analytically the amount of energy radiated after the merger.