Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

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.

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers.

With a view to developing a formalism that will be applicable at second perturbative order, we devise a new practical scheme for computing the gravitational self-force experienced by a point mass moving in a curved background spacetime. Our method works in the frequency domain and employs the effective-source approach, in which a distributional source for the retarded metric perturbation is replaced with an effective source for a certain regularized self-field. A key ingredient of the calculation is the analytic determination of an appropriate puncture field from which the effective source and regularized residual field can be calculated. In addition to its application in our effective-source method, we also show how this puncture field can be used to derive tensor-harmonic mode-sum regularization parameters that improve the efficiency of the traditional mode-sum procedure. To demonstrate the method, we calculate the first-order-in-the-mass-ratio self-force and redshift invariant for a point mass on a circular orbit in Schwarzschild spacetime.

Building on substantial foundational progress in understanding the effect of a small body's self-field on its own motion, the past 15 years has seen the emergence of several strategies for explicitly computing self-field corrections to the equations of motion of a small, point-like charge. These approaches broadly fall into three categories: (i) mode-sum regularization, (ii) effective source approaches and (iii) worldline convolution methods. This paper reviews the various approaches and gives details of how each one is implemented in practice, highlighting some of the key features in each case.

The equations of motion of a point particle interacting with its own field are defined in terms of a certain regularized self-field. Two of the leading methods for computing this regularized field are the mode-sum and effective-source approaches. In this work we unite these two distinct regularization schemes by generalizing traditional frequency-domain mode-sum calculations to incorporate effective-source techniques. For a toy scalar-field model we analytically compute an appropriate puncture field from which the regularized residual field can be calculated. To demonstrate the method, we compute the self-force for a scalar particle on a circular orbit in Schwarzschild spacetime. We also demonstrate the relation between the worldtube and window function approaches to localizing the puncture field to the neighborhood of the worldline and show how the method reduces to the well-known mode-sum regularization scheme in a certain limit. This new computational scheme can be applied to cases where traditional mode-sum regularization is inadequate, such as in calculations at second perturbative order.

This is the second in a series of papers aimed at developing a practical time-domain method for self-force calculations in Kerr spacetime. The key elements of the method are (i) removal of a singular part of the perturbation field with a suitable analytic "puncture" based on the Detweiler--Whiting decomposition, (ii) decomposition of the perturbation equations in azimuthal ($m$-)modes, taking advantage of the axial symmetry of the Kerr background, (iii) numerical evolution of the individual $m$-modes in 2+1-dimensions with a finite difference scheme, and (iv) reconstruction of the physical self-force from the mode sum. Here we report an implementation of the method to compute the scalar-field self-force along circular equatorial geodesic orbits around a Kerr black hole. This constitutes a first time-domain computation of the self force in Kerr geometry. Our time-domain code reproduces the results of a recent frequency-domain calculation by Warburton and Barack, but has the added advantage of being readily adaptable to include the back-reaction from the self force in a self-consistent manner. In a forthcoming paper---the third in the series---we apply our method to the gravitational self-force (in the Lorenz gauge).

We present the first successful application of the method of Matched Expansions for the calculation of the self-force on a point particle in a curved spacetime. We investigate the case of a scalar charge in the Nariai spacetime, which serves as a toy model for a point mass moving in the Schwarzschild black hole background. We discuss the singularity structure of the Green function beyond the normal neighbourhood and the interesting effect of caustics on null wave propagation.

We present the first application of the Poisson-Wiseman-Anderson method of matched expansions, to compute the self-force acting on a point particle moving in a curved spacetime. The method uses two expansions for the Green function, valid in `quasilocal' and `distant past' regimes, which are matched within the normal neighbourhood. We perform our calculation in a static region of the spherically symmetric Nariai spacetime (dS_2 x S^2), on which scalar perturbations are governed by a radial equation with a P\"oschl-Teller potential. We combine (i) a very high order quasilocal expansion, and (ii) an expansion in quasinormal modes, to determine the Green function globally. We show it is singular everywhere on the null wavefront (even outside the normal neighbourhood), and apply asymptotic methods to determine its singular structure. We find the Green function undergoes a transition every time the null wavefront passes through a caustic: the singular part follows a repeating four-fold sequence $\delta(\sigma)$, $1/\pi \sigma$, $-\delta(\sigma)$, $-1/\pi \sigma$ etc., where $\sigma$ is Synge's world function. The matched expansion method provides new insight into the non-local properties of the self-force; we find the contribution from the segment of the worldline lying outside the normal neighbourhood is significant. We compute the scalar self-force acting on a static particle, and validate against an alternative method. Finally, we discuss wave propagation on black hole spacetimes (where any expansion in quasinormal modes will be augmented by a branch cut integral) and predict that the Green function in Schwarzschild spacetime will inherit the four-fold singular structure found here.