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

This paper provides an extended exploration of the inverse-chirp gravitational-wave signals from stellar collapse in massive scalar-tensor gravity reported in [Phys. Rev. Lett. 119, 201103]. We systematically explore the parameter space that characterizes the progenitor stars, the equation of state, and the scalar-tensor theory of the core collapse events. We identify a remarkably simple and straightforward classification scheme of the resulting collapse events. For any given set of parameters, the collapse leads to one of three end states: a weakly scalarized neutron star, a strongly scalarized neutron star, or a black hole, possibly formed in multiple stages. The latter two end states can lead to strong gravitational-wave signals that may be detectable in present continuous-wave searches with ground-based detectors. We identify a very sharp boundary in the parameter space that separates events with strong gravitational-wave emission from those with negligible radiation.

We study how the gravitational wave signal from stellar collapse in scalar-tensor gravity varies under the influence of scalar self-interaction. To this end, we extract the gravitational radiation from numerical simulations of stellar collapse for a range of potentials with higher-order terms in addition to the quadratic mass term. Our study includes collapse to neutron stars and black holes and we find the strong inverse-chirp signals obtained for the purely quadratic potential to be exceptionally robust under changes in the potential at higher orders; quartic and sextic terms in the potential lead to noticeable differences in the wave signal only if their contribution is amplified, implying a relative fine-tuning to within 5 or more orders of magnitude between the mass and self-interaction parameters. arXiv:1903.09704v1 [gr-qc]

We present new numerical-relativity simulations of eccentric merging black holes with initially anti-parallel spins lying in the orbital plane (the so-called superkick configuration). Binary eccentricity boosts the recoil of the merger remnant by up to 25%. The increase in the energy flux is much more modest, and therefore this kick enhancement is mainly due to asymmetry in the binary dynamics. Our findings might have important consequences for the retention of stellar-mass black holes in star clusters and supermassive black holes in galactic hosts.

We compute families of spherically symmetric neutron-star models in two-derivative scalar-tensor theories of gravity with a massive scalar field. The numerical approach we present allows us to compute the resulting spacetimes out to infinite radius using a relaxation algorithm on a compactified grid. We discuss the structure of the weakly and strongly scalarized branches of neutron-star models thus obtained and their dependence on the linear and quadratic coupling parameters α0, β0 between the scalar and tensor sectors of the theory, as well as the scalar mass μ. For highly negative values of β0, we encounter configurations resembling a “gravitational atom”, consisting of a highly compact baryon star surrounded by a scalar cloud. A stability analysis based on binding-energy calculations suggests that these configurations are unstable and we expect them to migrate to models with radially decreasing baryon density and scalar field strength.