Black holes have turned into cosmic laboratories to search for ultra-light scalars by virtue of the superradiant instability. In this paper we present a detailed study of the impact of the superradiant evolution on the black hole shadow and investigate the exciting possibility to explore it with future observations of Very Long Baseline Interferometry. We simulated the superradiant evolution numerically, in the adiabatic regime, and derived analytic approximations modelling the process. Driven by superradiance, we evolve the black hole shadow diameter and (i) find that it can change by a few $\mu$as, just below the current resolution of the Event Horizon Telescope, albeit on timescales that are longer than realistic observation times; (ii) show that the shadow diameter can either shrink or grow; and (iii) explore in detail how the shadow's end state is determined by the initial parameters and coupling.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

Ultralight bosonic fields in the mass range $\sim (10^{-20}-10^{-11})\,{\rm eV}$ can trigger a superradiant instability that extracts energy and angular momentum from an astrophysical black hole with mass $M\sim(5,10^{10})M_\odot$, forming a nonspherical, rotating condensate around it. So far, most studies of the evolution and end-state of the instability have been limited to initial data containing only the fastest growing superradiant mode. By studying the evolution of multimode data in a quasi-adiabatic approximation, we show that the dynamics is much richer and depend strongly on the energy of the seed, on the relative amplitude between modes, and on the gravitational coupling. If the seed energy is a few percent of the black-hole mass, a black hole surrounded by a mixture of superradiant and nonsuperradiant modes with comparable amplitudes might not undergo a superradiant unstable phase, depending on the value of the boson mass. If the seed energy is smaller, as in the case of an instability triggered by quantum fluctuations, the effect of nonsuperradiant modes is negligible. We discuss the implications of these findings for current constraints on ultralight fields with electromagnetic and gravitational-wave observations.

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 coupling between scalar and vector fields has a long and interesting history. Axions are one key possibility to solve the strong CP problem and axion-like particles could be one solution to the dark matter puzzle. Given the nature of the coupling, and the universality of free fall, nontrivial important effects are expected in regions where gravity is strong. Here, we show that i. A background EM field induces an axionic instability in flat space, for large enough electric fields. Conversely, a homogeneous harmonic axion field induces an instability in the Maxwell sector. When carried over to curved spacetime, this phenomena translates into generic instabilities of charged black holes (BHs). ii. In the presence of charge, BH uniqueness results are lost. We find solutions which are small deformations of the Kerr-Newman geometry and hairy stationary solutions without angular momentum, which are `dragged' by the axion. Axion fields must exist around spinning BHs if these are immersed in external magnetic fields. The axion profile can be obtained perturbatively from the electro-vacuum solution derived by Wald. iii. Ultralight axions trigger superradiant instabilities of spinning BHs and form an axionic cloud in the exterior geometry. The superradiant growth can be interrupted or suppressed through axionic or scalar couplings to EM. These couplings lead to periodic bursts of light, which occur throughout the history of energy extraction from the BH. We provide numerical and simple analytical estimates for the rates of these processes. iv. Finally, we discuss how plasma effects can affect the evolution of superradiant instabilities.

We study the nonlinear dynamics of black holes that carry scalar hair and binaries composed of such black holes. The scalar hair is due to a linear or exponential coupling between the scalar and the Gauss--Bonnet invariant. We work perturbatively in the coupling constant of that interaction but nonperturbatively in the fields. We first consider the dynamical formation of hair for isolated black holes of arbitrary spin and determine the final state. This also allows us to compute for the first time the scalar quasinormal modes of rotating black holes in the presence of this coupling. We then study the evolution of nonspinning black-hole binaries with various mass ratios and produce the first scalar waveform for a coalescence. An estimate of the energy loss in scalar radiation and the effect this has on orbital dynamics and the phase of the GWs (entering at quadratic order in the coupling) show that GW detections can set the most stringent constraint to date on theories that exhibit a coupling between a scalar field and the Gauss--Bonnet invariant.

Scalar fields coupled to the Gauss-Bonnet invariant evade the known no-hair theorems and have nontrivial configurations around black holes. We focus on a scalar field that couples linearly to the Gauss-Bonnet invariant and hence exhibits shift symmetry. We study its dynamical evolution and the formation of scalar hair in a Schwarzschild background. We show that the evolution eventually settles to the known static hairy solutions in the appropriate limit.

A linear coupling between a scalar field and the Gauss-Bonnet invariant is the only known interaction term between a scalar and the metric that: respects shift symmetry; does not lead to higher order equations; inevitably introduces black hole hair in asymptotically flat, 4-dimensional spacetimes. Here we focus on the simplest theory that includes such a term and we explore the dynamical formation of scalar hair. In particular, we work in the decoupling limit that neglects the backreaction of the scalar onto the metric and evolve the scalar configuration numerically in the background of a Schwarzschild black hole or a collapsing dust star described by the Oppenheimer-Snyder solution. For all types of initial data that we consider, the scalar relaxes at late times to the known, static, analytic configuration that is associated with a hairy, spherically symmetric black hole. This suggests that the corresponding black hole solutions are indeed endpoints of collapse.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

Fundamental fields are a natural outcome in cosmology and particle physics and might therefore serve as a proxy for more complex interactions. The equivalence principle implies that all forms of matter gravitate, and one therefore expects relevant, universal imprints of new physics in strong field gravity, such as that encountered close to black holes. Fundamental fields in the vicinities of supermassive black holes give rise to extremely long-lived, or even unstable, configurations which slowly extract angular momentum from the black hole or simply evolve non-linearly over long timescales, with important implications for particle physics and gravitational-wave physics. Here, we perform a fully non-linear study of scalar-field condensates around rotating black holes. We provide novel ways to specify initial data for the Einstein-Klein-Gordon system, with potential applications in a variety of scenarios. Our numerical results confirm the existence of long-lived bar-modes which act as lighthouses for gravitational wave emission: the scalar field condenses outside the black hole geometry and acts as a constant frequency gravitational-wave source for very long timescales. This effect could turn out to be a potential signature of beyond standard model physics and also a promising source of gravitational waves for future gravitational wave detectors.

Black holes are among the most exciting phenomena predicted by General Relativity and play a key role in fundamental physics. Many interesting phenomena involve dynamical black hole configurations in the high curvature regime of gravity. In these lecture notes I will summarise the main numerical relativity techniques to explore highly dynamical phenomena, such as black hole collisions, in generic $D$-dimensional spacetimes.

The dynamics of black hole spacetimes play a crucial role in astrophysics, high energy physics and fundamental physics. In this thesis I have investigated the dynamics of black holes in generic spacetimes by extending established numerical relativity methods to higher dimensional or non-asymptotically flat spacetimes. Additionally, I have explored BH spacetimes perturbed my massive scalar fields. By developing and further improving NR techniques I have been able to push our knowledge to new grounds.

Light bosonic degrees of freedom have become a serious candidate for dark matter. The evolution of these fields around curved spacetimes is poorly understood but is expected to display interesting effects. In particular, the interaction of light bosonic fields with supermassive black holes, key players in most galaxies, could provide colourful examples of superradiance and nonlinear bosenova-like collapse. In turn, the observation of spinning black holes is expected to impose stringent bounds on the mass of putative massive bosonic fields in our universe. Our purpose here is to present a comprehensive study of the evolution of linearized massive scalar and vector fields in the vicinities of rotating black holes. For a certain boson field mass range, the field can become trapped in a potential barrier outside the horizon and transition to a bound state. Because there are a number of such quasi-bound states, the generic outcome is an amplitude modulated sinusoidal, or beating, signal. We believe that the appearance of such beatings has gone unnoticed in the past, and in fact mistaken for exponential growth. The amplitude modulation of the signal depends strongly on the relative excitation of the overtones, which in turn is strongly tied to the bound-state geography. For the first time we explore massive vector fields in generic BH background which are hard, if not impossible, to separate in the Kerr background. Our results show that spinning BHs are generically strongly unstable against massive vector fields.

Nonlinear dynamics in cosmological backgrounds has the potential to teach us immensely about our universe, and also to serve as prototype for nonlinear processes in generic curved spacetimes. Here we report on dynamical evolutions of black holes in asymptotically de Sitter spacetimes. We focus on the head-on collision of equal mass binaries and for the first time compare analytical and perturbative methods with full blown nonlinear simulations. Our results include an accurate determination of the merger/scatter transition (consequence of an expanding background) for small mass binaries and a test of the Cosmic Censorship conjecture, for large mass binaries. We observe that, even starting from small separations, black holes in large mass binaries eventually lose causal contact, in agreement with the conjecture.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

The puncture method for dealing with black holes in the numerical simulation of vacuum spacetimes is remarkably successful when combined with the BSSN formulation of the Einstein equations. We examine a generalized class of formulations modeled along the lines of the Laguna-Shoemaker system, including BSSN as a special case. The formulation is a two parameter generalization of the choice of variables used in standard BSSN evolutions. Numerical stability of the standard finite difference methods is proven for the formulation in the linear regime around flat space, a special case of which is the numerical stability of BSSN. Numerical evolutions are presented and compared with a standard BSSN implementation. We find that a significant portion of the parameter space leads to stable evolutions and that standard BSSN is located near the edge of the stability region. Non-standard parameter choices typically result in smoother behaviour of the evolution variables close to the puncture and thus hold promise for improved accuracy in, e.g., long-term BH binary inspirals, and for overcoming (numerical) stability problems still encountered in some types of black-hole simulations, e.g., in $D \ge 6$ dimensions.

Numerical relativity has seen incredible progress in the last years, and is being applied with success to a variety of physical phenomena, from gravitational-wave research and relativistic astrophysics to cosmology and high-energy physics. Here we probe the limits of current numerical setups, by studying collisions of unequal mass, non-rotating black holes of mass-ratios up to 1:100 and making contact with a classical calculation in General Relativity: the infall of a point-like particle into a massive black hole. Our results agree well with the predictions coming from linearized calculations of the infall of point-like particles into non-rotating black holes. In particular, in the limit that one hole is much smaller than the other, and the infall starts from an infinite initial separation, we recover the point-particle limit. Thus, numerical relativity is able to bridge the gap between fully non-linear dynamics and linearized approximations, which may have important applications. Finally, we also comment on the "spurious" radiation content in the initial data and the linearized predictions.

Black objects in higher dimensional space-times have a remarkably richer structure than their four dimensional counterparts. They appear in a variety of configurations (e.g. black holes, black branes, black rings, black Saturns), and display complex stability phase diagrams. They might also play a key role in high energy physics: for energies above the fundamental Planck scale, gravity is the dominant interaction which, together with the hoop-conjecture, implies that the trans-Planckian scattering of point particles should be well described by black hole scattering. Higher dimensional scenarios with a fundamental Planck scale of the order of TeV predict, therefore, black hole production at the LHC, as well as in future colliders with yet higher energies. In this setting, accurate predictions for the production cross-section and energy loss (through gravitational radiation) in the formation of black holes in parton-parton collisions is crucial for accurate phenomenological modelling in Monte Carlo event generators. In this paper, we use the formalism and numerical code reported in arXiv:1001.2302 to study the head-on collision of two black holes. For this purpose we provide a detailed treatment of gravitational wave extraction in generic D-dimensional space-times, which uses the Kodama-Ishibashi formalism. For the first time, we present the results of numerical simulations of the head-on collision in five space-time dimensions, together with the relevant physical quantities. We show that the total radiated energy, when two black holes collide from rest at infinity, is approximately (0.089\pm 0.006)% of the centre of mass energy, slightly larger than the 0.055% obtained in the four dimensional case, and that the ringdown signal at late time is in very good agreement with perturbative calculations.

The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modelling black hole production in TeV gravity scenarios, analysis of the stability of exact solutions and tests of Cosmic Censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D-2) isometry group for D\ge 5, or SO(D-3) for D\ge 6. Performing a dimensional reduction on a (D-4)-sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata and Nakamura (BSSN) formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the LEAN code and perform a variety of simulations of non-spinning black hole space-times. Specifically, we present a modified moving puncture gauge which facilitates long term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5,6.

Geodesic motion determines important features of spacetimes. Null unstable geodesics are closely related to the appearance of compact objects to external observers and have been associated with the characteristic modes of black holes. By computing the Lyapunov exponent, which is the inverse of the instability timescale associated with this geodesic motion, we show that, in the eikonal limit, quasinormal modes of black holes in any dimensions are determined by the parameters of the circular null geodesics. This result is independent of the field equations and only assumes a stationary, spherically symmetric and asymptotically flat line element, but it does not seem to be easily extendable to anti-de Sitter spacetimes. We further show that (i) in spacetime dimensions greater than four, equatorial circular timelike geodesics in a Myers-Perry black hole background are unstable, and (ii) the instability timescale of equatorial null geodesics in Myers-Perry spacetimes has a local minimum for spacetimes of dimension d > 5.