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.

Black hole spacetimes, like the Kerr spacetime, admit both stable and plunging orbits, separated in parameter space by the separatrix. Determining the location of the separatrix is of fundamental interest in understanding black holes, and is of crucial importance for modeling extreme mass-ratio inspirals. Previous numerical approaches to locating the Kerr separatrix were not always efficient or stable across all of parameter space. In this paper we show that the Kerr separatrix is the zero set of a single polynomial in parameter space. This gives two main results. First, we thoroughly analyze special cases (extreme Kerr, polar orbits, etc.), finding strict bounds on the limits of roots, and unifying a number of results in the literature. Second, we pose a stable numerical method which is guaranteed to quickly and robustly converge to the separatrix. This new approach is implemented in the Black Hole Perturbation Toolkit, and results in a ~45x speedup over the prior robust approach.

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.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

We present the first fully relativistic study of gravitational radiation from bodies in circular equatorial orbits around the massive black hole at the Galactic Center, Sgr A* and we assess the detectability of various kinds of objects by the gravitational wave detector LISA. Our computations are based on the theory of perturbations of the Kerr spacetime and take into account the Roche limit induced by tidal forces in the Kerr metric. The signal-to-noise ratio in the LISA detector, as well as the time spent in LISA band, are evaluated. We have implemented all the computational tools in an open-source SageMath package, within the Black Hole Perturbation Toolkit framework. We find that white dwarfs, neutrons stars, stellar black holes, primordial black holes of mass larger than $10^{-4} M_\odot$, main-sequence stars of mass lower than $\sim 2.5\, M_\odot$ and brown dwarfs orbiting Sgr A* are all detectable in one year of LISA data with a signal-to-noise ratio above 10 during at least $10^5$ years in the slow inspiral towards either the innermost stable circular orbit (compact objects) or the Roche limit (main-sequence stars and brown dwarfs). The longest times in-band, of the order of $10^6$ years, are achieved for primordial black holes of mass $\sim 10^{-3} M_\odot$ down to $10^{-5} M_\odot$, depending on the spin of Sgr A*, as well as for brown dwarfs, just followed by white dwarfs and low mass main-sequence stars. The long time in-band of these objects makes Sgr A* a valuable target for LISA. We also consider bodies on close circular orbits around the massive black hole in the nucleus of the nearby galaxy M32 and find that, among them, compact objects and brown dwarfs stay for $10^3$ to $10^4$ years in LISA band with a 1-year signal-to-noise ratio above 10.

We calculate the evolution and gravitational-wave emission of a spinning compact object inspiraling into a substantially more massive (non-rotating) black hole. We extend our previous model for a non-spinning binary [Phys. Rev. D 93, 064024] to include the Mathisson-Papapetrou-Dixon spin-curvature force. For spin-aligned binaries we calculate the dephasing of the inspiral and associated waveforms relative to models that do not include spin-curvature effects. We find this dephasing can be either positive or negative depending on the initial separation of the binary. For binaries in which the spin and orbital angular momentum are not parallel, the orbital plane precesses and we use a more general osculating element prescription to compute inspirals.

We model the inspiral of a compact object into a more massive black hole rotating very near the theoretical maximum. We find that once the body enters the near-horizon regime the gravitational radiation is characterized by a constant frequency, equal to (twice) the horizon frequency, with an exponentially damped profile. This contrasts with the usual "chirping" behavior and, if detected, would constitute a "smoking gun" for a near-extremal black hole in nature.

We compute the radiation emitted by a particle on the innermost stable circular orbit of a rapidly spinning black hole both (a) analytically, working to leading order in the deviation from extremality and (b) numerically, with a new high-precision Teukolsky code. We find excellent agreement between the two methods. We confirm previous estimates of the overall scaling of the power radiated, but show that there are also small oscillations all the way to extremality. Furthermore, we reveal an intricate mode-by-mode structure in the flux to infinity, with only certain modes having the dominant scaling. The scaling of each mode is controlled by its conformal weight, a quantity that arises naturally in the representation theory of the enhanced near-horizon symmetry group. We find relationships to previous work on particles orbiting in precisely extreme Kerr, including detailed agreement of quantities computed here with conformal field theory calculations performed in the context of the Kerr/CFT correspondence.

We present a calculation of the scalar field self-force (SSF) acting on a scalar-charge particle in a strong-field orbit around a Kerr black hole. Our calculation specializes to circular and equatorial geodesic orbits. The analysis is an implementation of the standard mode-sum regularization scheme: We first calculate the multipole modes of the scalar-field perturbation using numerical integration in the frequency domain, and then apply a certain regularization procedure to each of the modes. The dissipative piece of the SSF is found to be consistent with the flux of energy and angular momentum carried by the scalar waves through the event horizon and out to infinity. The conservative (radial) component of the SSF is found to be attractive (inward pointing) for $r_0>r_{\rm c}(a)$ and repulsive (outward pointing) for $r_0<r_{\rm c}(a)$, where $a$ is the Kerr spin parameter, $r_0$ is the Boyer-Lindquist orbital radius, and $r_{\rm c}$ is a critical $a$-dependent radius at which the conservative SSF vanishes. When the motion is retrograde the conservative SSF is repulsive for all $r_0$ (as in the Schwarzschild case). The dominant conservative effect of the SSF in Schwarzschild spacetime is known to be of 3rd post-Newtonian (PN) order (with a logarithmic running). Our numerical results suggest that the leading-order PN correction due to the black hole's spin arises from spin-orbit coupling at 3PN, which dominates the overall SSF effect at large $r_0$. In PN language, the change-of-sign of the radial SSF is attributed to an interplay between the spin-orbit term ($\propto -ar_0^{-4.5}$) and the "Schwarzschild" term ($\propto r_0^{-5}\log r_0$).

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.

We present a numerical code for calculating the self force on a scalar charge moving in a bound (eccentric) geodesic in the equatorial plane of a Kerr black hole. We work in the frequency domain and make use of the method of extended homogeneous solutions [Phys.\ Rev.\ D {\bf 78}, 084021 (2008)], in conjunction with mode-sum regularization. Our work is part of a program to develop a computational architecture for fast and efficient self-force calculations, alternative to time-domain methods. We find that our frequency-domain method outperforms existing time-domain schemes for small eccentricities, and, remarkably, remains competitive up to eccentricities as high as $\sim 0.7$. As an application of our code we (i) compute the conservative scalar-field self-force correction to the innermost stable circular equatorial orbit, as a function of the Kerr spin parameter; and (ii) calculate the variation in the rest mass of the scalar particle along the orbit, caused by the component of the self force tangent to the four-velocity.

Accurately modeling astrophysical extreme-mass-ratio-insprials requires calculating the gravitational self-force for orbits in Kerr spacetime. The necessary calculation techniques are typically very complex and, consequently, toy scalar-field models are often developed in order to establish a particular calculational approach. To that end, I present a calculation of the scalar-field self-force for a particle moving on a (fixed) inclined circular geodesic of a background Kerr black hole. I make the calculation in the frequency-domain and demonstrate how to apply the mode-sum regularization procedure to all four components of the self-force. I present results for a number of strong-field orbits which can be used as benchmarks for emerging self-force calculation techniques in Kerr spacetime.

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.

For a self-gravitating particle of mass \mu in orbit around a Kerr black hole of mass M >> \mu, we compute the O(\mu/M) shift in the frequency of the innermost stable circular equatorial orbit (ISCEO) due to the conservative piece of the gravitational self-force acting on the particle. Our treatment is based on a Hamiltonian formulation of the dynamics in terms of geodesic motion in a certain locally-defined effective smooth spacetime. We recover the same result using the so-called first law of binary black-hole mechanics. We give numerical results for the ISCEO frequency shift as a function of the black hole's spin amplitude, and compare with predictions based on the post-Newtonian approximation and the effective one-body model. Our results provide an accurate strong-field benchmark for spin effects in the general relativistic two-body problem.

We present results from calculations of the orbital evolution in eccentric binaries of nonrotating black holes with extreme mass-ratios. Our inspiral model is based on the method of osculating geodesics, and is the first to incorporate the full gravitational self-force (GSF) effect, including conservative corrections. The GSF information is encapsulated in an analytic interpolation formula based on numerical GSF data for over a thousand sample geodesic orbits. We assess the importance of including conservative GSF corrections in waveform models for gravitational-wave searches.

Bound geodesic orbits around a Kerr black hole can be parametrized by three constants of the motion: the (specific) orbital energy, angular momentum and Carter constant. Generically, each orbit also has associated with it three frequencies, related to the radial, longitudinal and (mean) azimuthal motions. Here we note the curious fact that these two ways of characterizing bound geodesics are not in a one-to-one correspondence. While the former uniquely specifies an orbit up to initial conditions, the latter does not: there is a (strong-field) region of the parameter space in which pairs of physically distinct orbits can have the same three frequencies. In each such isofrequency pair the two orbits exhibit the same rate of periastron precession and the same rate of Lense-Thirring precession of the orbital plane, and (in a certain sense) they remain "synchronized" in phase.