Extreme-mass-ratio inspirals, in which a stellar-mass compact object spirals into a supermassive black hole in a galactic core, are expected to be key sources for LISA. Modelling these systems with sufficient accuracy for LISA science requires going to second (or {\em post-adiabatic}) order in gravitational self-force theory. Here we present a practical two-timescale framework for achieving this and generating post-adiabatic waveforms. The framework comprises a set of frequency-domain field equations that apply on the fast, orbital timescale, together with a set of ordinary differential equations that determine the evolution on the slow, inspiral timescale. Our analysis is restricted to the special case of quasicircular orbits around a Schwarzschild black hole, but its general structure carries over to the realistic case of generic (inclined and eccentric) orbits in Kerr spacetime. In our restricted context, we also develop a tool that will be useful in all cases: a formulation of the frequency-domain field equations using hyperboloidal slicing, which significantly improves the behavior of the sources near the boundaries. We give special attention to the slow evolution of the central black hole, examining its impact on both the two-timescale evolution and the earlier self-consistent evolution scheme.

[Abridged] This review surveys the theory of gravitational self-force in curved spacetime and its application to the gravitational two-body problem in the extreme-mass-ratio regime. We first lay the relevant formal foundation, describing the rigorous derivation of the equation of self-forced motion using matched asymptotic expansions and other ideas. We then review the progress that has been achieved in numerically calculating the self-force and its physical effects in the astrophysical scenario of a compact object inspiralling into a (rotating) massive black hole. We highlight the way in which, nowadays, self-force calculations make a fruitful contact with other approaches to the two-body problem and help inform an accurate universal model of binary black hole inspirals, valid across all mass ratios. We conclude with a summary of the state of the art, open problems and prospects. Our review is aimed at non-specialist readers and is for the most part self-contained and non-technical; only elementary-level acquaintance with General Relativity is assumed. Where useful, we draw on analogies with familiar concepts from Newtonian gravity or classical electrodynamics.

Vacuum perturbations of the Kerr metric can be reconstructed from the corresponding perturbation in either of the two Weyl scalars $\psi_0$ or $\psi_4$, using a procedure described by Chrzanowski and others in the 1970s. More recent work, motivated within the context of self-force physics, extends the procedure to metric perturbations sourced by a particle in a bound geodesic orbit. However, the existing procedure leaves undetermined a certain stationary, axially-symmetric piece of the metric perturbation. In the vacuum region away from the particle, this "completion" piece corresponds simply to mass and angular-momentum perturbations of the Kerr background, with amplitudes that are, however, a priori unknown. Here we present and implement a rigorous method for finding the completion piece. The key idea is to impose continuity, off the particle, of certain gauge-invariant fields constructed from the full (completed) perturbation, in order to determine the unknown amplitude parameters of the completion piece. We implement this method in full for bound (eccentric) geodesic orbits in the equatorial plane of the Kerr black hole. Our results provide a rigorous underpinning of recent results by Friedman {\it et al.}\ for circular orbits, and extend them to non-circular orbits.

In general-relativistic perturbation theory, a point mass accelerates away from geodesic motion due to its gravitational self-force. Because the self-force is small, one can often approximate the motion as geodesic. However, it is well known that self-force effects accumulate over time, making the geodesic approximation fail on long timescales. It is less well known that this failure at large times translates to a failure at large distances as well. At second perturbative order, two large-distance pathologies arise: spurious secular growth and infrared-divergent retarded integrals. Both stand in the way of practical computations of second-order self-force effects. Utilizing a simple flat-space scalar toy model, I develop methods to overcome these obstacles. The secular growth is tamed with a multiscale expansion that captures the system's slow evolution. The divergent integrals are eliminated by matching to the correct retarded solution at large distances. I also show how to extract conservative self-force effects by taking local-in-time "snapshots" of the global solution. These methods are readily adaptable to the physically relevant case of a point mass orbiting a black hole.

A compact object moving on a quasicircular orbit about a Schwarzschild black hole gradually spirals inward due to the dissipative action of its gravitational self-force. But in addition to driving the inspiral, the self-force has a conservative piece. Within a second-order self-force formalism, I derive a second-order generalization of Detweiler's redshift variable, which provides a gauge-invariant measure of conservative effects on quasicircular orbits. I sketch a frequency-domain numerical scheme for calculating this quantity. Once this scheme has been implemented, its results may be used to determine high-order terms in post-Newtonian theory and parameters in effective-one-body theory.

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

Accurately modeling an extreme-mass-ratio inspiral requires knowledge of the second-order gravitational self-force on the inspiraling small object. Recently, numerical puncture schemes have been formulated to calculate this force, and their essential analytical ingredients have been derived from first principles. However, the \emph{puncture}, a local representation of the small object's self-field, in each of these schemes has been presented only in a local coordinate system centered on the small object, while a numerical implementation will require the puncture in coordinates covering the entire numerical domain. In this paper we provide an explicit covariant self-field as a local expansion in terms of Synge's world function. The self-field is written in the Lorenz gauge, in an arbitrary vacuum background, and in forms suitable for both self-consistent and Gralla-Wald-type representations of the object's trajectory. We illustrate the local expansion's utility by sketching the procedure of constructing from it a numerically practical puncture in any chosen coordinate system.

Calculations of the gravitational self-force (GSF) in curved spacetime require as input the metric perturbation in a sufficiently regular gauge. A basic challenge in the program to compute the GSF for orbits around a Kerr black hole is that the standard procedure for reconstructing the perturbation is formulated in a class of radiation gauges, in which the particle singularity is non-isotropic and extends away from the particle's location. Here we present two practical schemes for calculating the GSF using a radiation-gauge reconstructed metric as input. The schemes are based on a detailed analysis of the local structure of the particle singularity in the radiation gauges. We identify 3 types of radiation gauges: two containing a radial string-like singularity emanating from the particle, either in one direction ("half-string" gauges) or both directions ("full-string" gauges); and a third type containing no strings but with a jump discontinuity across a surface intersecting the particle. Based on a flat-space example, we argue that the standard mode-by-mode reconstruction procedure yields the "regular half" of a half-string solution, or (equivalently) either of the regular halves of a no-string solution. For the half-string case, we formulate the GSF in a locally deformed radiation gauge that removes the string singularity near the particle. We derive a mode-sum formula for the GSF in this gauge, analogous to the standard Lorenz-gauge formula but with modified regularization parameters. For the no-string case, we formulate the GSF directly, without a local deformation, and we derive a mode-sum formula that requires no correction to the parameters but involves a certain averaging procedure. We explain the consistency of our results with Gralla's invariance theorem, and discuss the correspondence between our method and a related approach by Friedman et al.