Search

We compute the (displacement) gravitational wave memory due to a quasicircular inspiral of two black holes using a variety of perturbative techniques. Within post-Newtonian theory, we extend previous results for non-spinning binaries to 3.5PN order. Using the gravitational self-force approach, we compute the memory at first order in the mass ratio for inspirals into a Kerr black hole. We do this both numerically and via a double post-Newtonian--self-force expansion which we carry out to 5PN order. At second order in the self-force approach, near-zone calculations encounter an infrared divergence associated with memory, which is resolved through matching the near-zone solution to a post-Minkowskian expansion in the far zone. We describe that matching procedure for the first time and show how it introduces nonlocal-in-time memory effects into the two-body dynamics at second order in the mass ratio, as was also predicted by recent 5PN calculations within the effective field theory approach. We then compute the gravitational-wave memory through second order in the mass ratio (excluding certain possible memory distortion effects) and find that it agrees well with recent results from numerical relativity simulations for near-comparable-mass binaries., Comment: 66 pages plus appendices and bibliography, 10 figures

Worldtube excision is a method of reducing computational burden in Numerical Relativity simulations of binary black holes in situations where there is a good analytical model of the geometry around (one or both of) the objects. Two such scenarios of relevance in gravitational-wave astronomy are (1) the case of mass-disparate systems, and (2) the early inspiral when the separation is still large. Here we illustrate the utility and flexibility of this technique with simulations of the fully self-consistent radiative evolution in the model problem of a scalar charge orbiting a Schwarzschild black hole under the effect of scalar-field radiation reaction. We explore a range of orbital configurations, including inspirals with large eccentricity (which we follow through to the final plunge and ringdown) and hyperbolic scattering., Comment: 8 pages, 4 figures

In the last decade, the first law of binary black hole mechanics played an important unifying role in the gravitational two-body problem. More recently, binary black hole scattering and the application of high-energy physics methods have provided a new avenue into this classical problem. In this Letter, we connect these two themes by extending the first law to the case of scattering orbits. We present derivations based on classical S-matrix, Hamiltonian, and pseudo-Hamiltonian methods, the last of which allows us to include dissipative effects for the first time. Finally, a "boundary to bound" map links this first law to the traditional bound-orbit version. Through this map a little-known observable for scatter orbits, the elapsed proper time, is mapped to the Detweiler redshift for bound orbits., Comment: 5 pages + supplemental material

The past three years have seen two significant advances in models of gravitational waveforms emitted by quasicircular compact binaries in two regimes: the weak-field, post-Newtonian regime, in which the gravitational wave energy flux has now been calculated to fourth-and-a-half post-Newtonian order (4.5PN) [Phys. Rev. Lett. 131, 121402 (2023)]; and the small-mass-ratio, gravitational self-force regime, in which the flux has now been calculated to second perturbative order in the mass ratio (2SF) [Phys. Rev. Lett. 127, 151102 (2021)]. We compare these results and find excellent agreement for the total flux, showing consistency between the two calculations at all available PN and SF orders. However, although the total fluxes agree, we find disagreements in the fluxes due to individual spherical-harmonic modes of the waveform, strongly suggesting the two waveforms might be in different asymptotic frames., Comment: 15 pages, 4 figures, modes of PN flux provided in ancillary file

In this paper we present an approach to compute analytical post-Minkowskian corrections to unbound two-body scattering in the self-force formalism. Our method relies on a further low-velocity (post-Newtonian) expansion of the motion. We present a general strategy valid for gravitational and non-gravitational self-force, and we explicitly demonstrate our approach for a scalar charge scattering off a Schwarzschild black hole. We compare our results with recent calculations in [Barack et al., PRD 108, 024025 (2023)], showing complete agreement where appropriate and fixing undetermined scale factors in their calculation. Our results also extend their results by including in our dissipative sector the contributions from the flux into the black hole horizon., Comment: 24 pages, minor changes

Quasi-normal modes (QNMs) uniquely describe the gravitational-wave ringdown of post-merger black holes. While the linear QNM regime has been extensively studied, recent work has highlighted the importance of second-perturbative-order, quadratic QNMs (QQNMs) arising from the nonlinear coupling of linear QNMs. Previous attempts to quantify the magnitude of these QQNMs have shown discrepant results. Using a new hyperboloidal framework, we resolve the discrepancy by showing that the QQNM/QNM ratio is a function not only of the black hole parameters but also of the ratio between even- and odd-parity linear QNMs: the ratio QQNM/QNM depends on what created the ringing black hole, but only through this ratio of even- to odd-parity linear perturbations., Comment: 8 pages, 1 figure

Compact binaries with asymmetric mass ratios are key expected sources for next-generation gravitational wave detectors. Gravitational self-force theory has been successful in producing post-adiabatic waveforms that describe the quasi-circular inspiral around a non-spinning black hole with sub-radian accuracy, in remarkable agreement with numerical relativity simulations. Current inspiral models, however, break down at the innermost stable circular orbit, missing part of the waveform as the secondary body transitions to a plunge into the black hole. In this work we derive the transition-to-plunge expansion within a multiscale framework and asymptotically match its early-time behaviour with the late inspiral. Our multiscale formulation facilitates rapid generation of waveforms: we build second post-leading transition-to-plunge waveforms, named 2PLT waveforms. Although our numerical results are limited to low perturbative orders, our framework contains the analytic tools for building higher-order waveforms consistent with post-adiabatic inspirals, once all the necessary numerical self-force data becomes available. We validate our framework by comparing against numerical relativity simulations, surrogate models and the effective one-body approach., Comment: Clarifications added regarding the dependence on $\delta M$ and $\delta J$ of metric perturbations and self-force

Second-order self-force calculations will be critical for modelling extreme-mass-ratio inspirals, and they are now known to have high accuracy even for binaries with mass ratios $\sim 1:10$. Many of the challenges facing these calculations are related to slow convergence of spherical-harmonic (or spheroidal harmonic) mode sums in a region containing the small companion. In this paper, we begin to develop a multi-domain framework that can evade those problems. Building on recent work by Osburn and Nishimura, in the problematic region of spacetime we use a puncture scheme and decompose the punctured field equations into a basis of Fourier and azimuthal $m$ modes, avoiding a harmonic decomposition in the $\theta$ direction. Outside the problematic region, we allow for a complete spherical- or spheroidal-harmonic decomposition. As a demonstration, we implement this framework in the simple context of a scalar charge in circular orbit around a Schwarzschild black hole. Our implementation utilizes several recent advances: a spectral method in each region, hyperboloidal compactification, and an extremely high-order puncture., Comment: 33 pages; 21 figures

Most self-force calculations rely, in one way or another, on representations of a particle's Detweiler- Whiting singular field. We present a simple method of calculating the singular field to high order in a local expansion in powers of distance from the particle. As a demonstration, we compute the singular field to 14th order in distance, 10 orders beyond the previous state of the art, in the simple case of a scalar charge in circular orbit around a Schwarzschild black hole. We provide the result in both a 4-dimensional form and a decomposed form suitable for use in an m-mode puncture scheme. Our method should have applications in overcoming bottlenecks in current self-force calculations at both first and second order in perturbation theory., Comment: 24 pages, 5 figures

Post-adiabatic models of extreme- and intermediate-mass-ratio inspirals will require calculations of second-order gravitational self-force effects in the spacetime of a spinning, Kerr black hole. We take a step toward such calculations by implementing the recently formulated Teukolsky puncture scheme with Green-Hollands-Zimmerman metric reconstruction [CQG 39, 015019 (2022)]. This scheme eliminates the critical obstacle of gauge singularities that arise in the standard no-string metric reconstruction. Our first proof-of-principle implementation is limited to the simple case of circular orbits in Schwarzschild spacetime, but the method also applies to generic orbits on a Kerr background. We conclude with a discussion of various approaches to the second-order self-force problem in Kerr., Comment: 31 pages, 9 figures

This is a third installment in a program to develop a method for alleviating the scale disparity in binary black hole simulations with mass ratios in the intermediate astrophysical range, where simulation cost is prohibitive while purely perturbative methods may not be adequate. The method is based on excising a "worldtube" around the smaller object, much larger than the object itself, replacing it with an analytical model that approximates a tidally deformed black hole. Previously (arXiv:2304.05329) we have tested the idea in a toy model of a scalar charge in a fixed circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field in 3+1 dimensions on the SpECTRE platform. Here we take the significant further step of allowing the orbit to evolve radiatively, in a self-consistent manner, under the effect of back-reaction from the scalar field. We compute the inspiral orbit and the emitted scalar-field waveform, showing a good agreement with perturbative calculations in the adiabatic approximation. We also demonstrate how our simulations accurately resolve post-adiabatic effects (for which we do not have perturbative results). In this work we focus on quasi-circular inspirals. Our implementation will shortly be publicly accessible in the SpECTRE numerical relativity code., Comment: 20 pages, 15 figures

We construct retarded and advanced Green's functions for gravitational perturbations in Kerr in an ingoing radiation gauge. Our Green's functions have a frequency domain piece that has previously been obtained by Ori [Phys. Rev. D 67 (2003)] based on the Chrzanowski-Cohen-Kegeles metric reconstruction method. As is well known, this piece by itself is not sufficient to obtain an actual Green's function. We show how to complete it with a piece based on a method by Green et al. [Class. Quant. Grav. 37 (2020)]. The completion piece has a completely explicit form in the time-domain and is supported on pairs of points on the same outgoing principal null geodesic which are in the appropriate causal order. We expect our Green's functions to be useful for gravitational self-force calculations and other perturbation problems on Kerr spacetime., Comment: 43 pages

The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase., Comment: 155 pages, with executive summary and table of contents

Second-order gravitational self-force theory has recently led to the breakthrough calculation of ``first post-adiabatic'' (1PA) compact-binary waveforms [Phys. Rev. Lett. 130, 241402 (2023)]. The computations underlying those waveforms depend on a method of solving the perturbative second-order Einstein equation on a Schwarzschild background in the Fourier domain. In this paper we present that method, which involves dividing the domain into several regions. Different regions utilize different time slicings and allow for the use of ``punctures'' to tame sources and enforce physical boundary conditions. We demonstrate the method for Lorenz-gauge and Teukolsky equations in the relatively simple case of calculating parametric derivatives (``slow time derivatives'') of first-order fields, which are an essential input at second order., Comment: 41 pages, 10 figures. Minor changes in response to referee

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome., Comment: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

We investigate the impact of post-adiabatic (1PA) terms on parameter estimation for extreme and intermediate mass-ratio inspirals using state-of-the-art waveform models. Our analysis is the first to employ Bayesian inference to assess systematic errors for 1PA waveforms. We find that neglecting 1PA terms introduces significant biases for the (small) mass ratio $\epsilon \gtrsim 10^{-5}$ for quasi circular orbits in Schwarzschild spacetime, which can be mitigated with resummed 3PN expressions at 1PA order. Moreover, we show that the secondary spin is strongly correlated with the other intrinsic parameters, and it can not be constrained for $\epsilon \lesssim 10^{-5}$. Finally, we highlight the need for addressing eccentric waveform systematics in the small-mass-ratio regime, as they yield stronger biases than the circular limit in both intrinsic and extrinsic parameters., Comment: 25 + 11 pages, 10 figures, 5 tables, accepted for publication in PRD. Minor corrections to figures and new appendix on 1PA corrections to the EMRI waveform amplitudes

High-accuracy gravitational-wave modeling demands going beyond linear, first-order perturbation theory. Particularly motivated by the need for second-order perturbative models of extreme-mass-ratio inspirals and black hole ringdowns, we present practical spherical-harmonic decompositions of the Einstein equation, Regge-Wheeler-Zerilli equations, and Teukolsky equation at second perturbative order in a Schwarzschild background. Our formulations are covariant on the $t$--$r$ plane and on the two-sphere, and we express the field equations in terms of gauge-invariant metric perturbations. In a companion Mathematica package, PerturbationEquations, we provide these invariant formulas as well as the analogous formulas in terms of raw, gauge-dependent metric perturbations. Our decomposition of the second-order Einstein equation, when specialized to the Lorenz gauge, was a key ingredient in recent second-order self-force calculations [Phys. Rev. Lett. 124, 021101 (2020); ibid. 127, 151102 (2021); ibid. 130, 241402 (2023)]., Comment: 41 pages. v2 contains additional material and corrects minor errors. v3 corrects additional minor errors. v4 makes minor changes in response to referee comments

To fully exploit the capabilities of next-generation gravitational wave detectors, we need to significantly improve the accuracy of our models of gravitational-wave-emitting systems. This paper focuses on one way of doing so: by taking black hole perturbation theory to second perturbative order. Such calculations are critical for the development of nonlinear ringdown models and of gravitational self-force models of extreme-mass-ratio inspirals. In the most astrophysically realistic case of a Kerr background, a second-order Teukolsky equation presents the most viable avenue for calculating second-order perturbations. Motivated by this, we analyse two second-order Teukolsky formalisms and advocate for the one that is well-behaved for gravitational self-force calculations and which meshes naturally with recent metric reconstruction methods due to Green, Hollands, and Zimmerman [CQG 37, 075001 (2020)] and others. Our main result is an expression for the nonlinear source term in the second-order field equation; we make this available, along with other useful tools, in an accompanying Mathematica notebook. Using our expression for the source, we also show that infrared divergences at second order can be evaded by adopting a Bondi--Sachs gauge., Comment: 22 pages, 0 figures. Accompanying Mathematica notebook available at https://github.com/DrAndrewSpiers/NP-and-GHP-Formalisms-for-2nd-order-Teukolsky

Binary black hole simulations become increasingly more computationally expensive with smaller mass ratios, partly because of the longer evolution time, and partly because the lengthscale disparity dictates smaller time steps. The program initiated by Dhesi et al. (arXiv:2109.03531) explores a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysical range ($10^{-4} \lesssim q \lesssim 10^{-2}$), where purely perturbative methods may not be adequate. A region ("worldtube") much larger than the small black hole is excised from the numerical domain, and replaced with an analytical model approximating a tidally deformed black hole. Here we apply this idea to a toy model of a scalar charge in a fixed circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field. This is a first implementation of the worldtube excision method in full 3+1 dimensions. We demonstrate the accuracy and efficiency of the method, and discuss the steps towards applying it for evolving orbits and, ultimately, in the binary black-hole scenario. Our implementation is publicly accessible in the SpECTRE numerical relativity code., Comment: 19 pages, 10 figures; remove typo in Eq. (9c)

We leverage recent breakthrough calculations using second-order gravitational self-force (2GSF) theory to improve both the gravitational-mode amplitudes and radiation-reaction force in effective-one-body~(EOB) waveform models. We achieve this by introducing new calibration parameters in the SEOBNRv5HM mode amplitudes, and matching them to the newly available 2GSF energy-flux multipolar data for quasicircular nonspinning binary black holes. We find that this significantly improves the SEOBNRv5HM energy flux, when compared to numerical-relativity (NR) simulations of binary black holes with mass ratios between 1:1 and 1:20. Moreover, we find that, once the conservative part of the SEOBNRv5 dynamics is calibrated, the SEOBNRv5HM waveform model with 2GSF information reproduces the binding energy of NR simulations more accurately, providing a powerful check of the consistency and naturalness of the EOB approach. While we only include nonspinning 2GSF information, the more accurate binding energy and energy flux carry over to the SEOBNRv5 waveform models for spinning binary black holes. Thus, our results improve the latest generation of SEOBNR waveform models (i.e., SEOBNRv5), which has been recently completed for use in the upcoming fourth observing (O4) run of the LIGO-Virgo-KAGRA Collaboration.

We compare recently computed waveforms from second-order gravitational self-force (GSF) theory to those generated by a new, GSF-informed, effective one body (EOB) waveform model for (spin-aligned, eccentric) inspiralling black hole binaries with large mass ratios. We focus on quasi-circular, nonspinning, configurations and perform detailed GSF/EOB waveform phasing comparisons, either in the time domain or via the gauge-invariant dimensionless function $Q_\omega\equiv \omega^2/\dot{\omega}$, where $\omega$ is the gravitational wave frequency. The inclusion of high-PN test-mass terms within the EOB radiation reaction (notably, up to 22PN) is crucial to achieve an EOB/GSF phasing agreement below 1~rad up to the end of the inspiral for mass ratios up to 500. For larger mass ratios, up to $5\times 10^4$, the contribution of horizon absorption becomes more and more important and needs to be accurately modeled. Our results indicate that our GSF-informed EOB waveform model is a promising tool to describe waveforms generated by either intermediate or extreme mass ratio inspirals for future gravitational wave detectors, Comment: 14 pages, 8 figures, published on Phys. Rev. D

We present the first systematic comparison between gravitational waveforms emitted by inspiralling, quasi-circular and nonspinning black hole binaries computed with three different approaches: second-order gravitational self-force (2GSF) theory, as implemented in the 1PAT1 model; numerical relativity (NR), as implemented by the SXS collaboration; and the effective one body (EOB) formalism, as implemented in the TEOBResumS waveform model. To compare the models we use both a standard, time-domain waveform alignment and a gauge-invariant analysis based on the dimensionless function $Q_\omega(\omega)\equiv \omega^2/\dot{\omega}$, where $\omega$ is the gravitational wave frequency. We analyse the domain of validity of the 1PAT1 model, deriving error estimates and showing that the effects of the final transition to plunge, which the model neglects, extend over a significantly larger frequency interval than one might expect. Restricting to the inspiral regime, we find that, while for mass ratios $q = m_1/m_2\le 10$ TEOBResumS is largely indistinguishable from NR, 1PAT1 has a significant dephasing $\gtrsim 1$rad; conversely, for $q\gtrsim 100$, 1PAT1 is estimated to have phase errors $<0.1$rad on a large frequency interval, while TEOBResumS develops phase differences $\gtrsim1$rad with it. Most crucially, on that same large frequency interval we find good agreement between TEOBResumS and 1PAT1 in the intermediate regime $15\lesssim q\lesssim 64$, with $<0.5$rad dephasing between them. A simple modification to the TEOBResumS flux further improves this agreement for $q\gtrsim 30$, reducing the dephasing to $\approx0.27$rad even at $q=128$. Our results pave the way for the construction of GSF-informed EOB models for both intermediate and extreme mass ratio inspirals for the next generation of gravitational wave detectors., Comment: 31 pages, 19 figures, published on Phys. Rev. D

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

We compute the linear metric perturbation to a Schwarzschild black hole generated by a spinning compact object, specialising to circular equatorial orbits with an (anti-)aligned spin vector. We derive a two-timescale expansion of the field equations, with an attendant waveform-generation framework, that includes all effects through first post-adiabatic order, and we use the Regge-Wheeler-Zerilli formalism in the frequency domain to generate waveforms that include the complete effect of the spin on the waveform phase. We perform the calculations using expansions at fixed orbital frequency, increasing the computational efficiency and simplifying the procedure compared to previous approaches. Finally, we provide the first fully relativistic, first-principles regularisation procedure for gauge invariant self-force quantities to linear order in spin. We use this procedure to produce the first strong-field, conservative self-force calculation including the spin of the secondary -- computing Detweiler's redshift invariant., Comment: 27 pages, 6 figures

We produce gravitational waveforms for nonspinning compact binaries undergoing a quasicircular inspiral. Our approach is based on a two-timescale expansion of the Einstein equations in second-order self-force theory, which allows first-principles waveform production in tens of milliseconds. Although the approach is designed for extreme mass ratios, our waveforms agree remarkably well with those from full numerical relativity, even for comparable-mass systems. Our results will be invaluable in accurately modelling extreme-mass-ratio inspirals for the LISA mission and intermediate-mass-ratio systems currently being observed by the LIGO-Virgo-KAGRA Collaboration.

Scientific analysis for the gravitational-wave detector LISA will require theoretical waveforms from extreme-mass-ratio inspirals (EMRIs) that extensively cover all possible orbital and spin configurations around astrophysical Kerr black holes. However, on-the-fly calculations of these waveforms have not yet overcome the high dimensionality of the parameter space. To confront this challenge, we present a user-ready EMRI waveform model for generic (eccentric and inclined) orbits in Kerr spacetime, using an analytical self-force approach. Our model accurately covers all EMRIs with arbitrary inclination and black hole spin, up to modest eccentricity ($\lesssim 0.3$) and separation ($\gtrsim2$--$10M$ from the last stable orbit). In that regime, our waveforms are accurate at the leading `adiabatic' order, and they approximately capture transient self-force resonances that significantly impact the gravitational-wave phase. The model fills an urgent need for extensive waveforms in ongoing data-analysis studies, and its individual components will continue to be useful in future science-adequate waveforms., Comment: 9+5 pages,4+1 figures,5PN-e10 Teukolsky data are all available at the Black Hole Perturbation Club: https://sites.google.com/view/bhpc1996/data v2:updated references,expanded supplemental discussion of initial phases.v3:contains minor corrections and more references.Matches version published in PRL

In recent years, gravitational-wave astronomy has motivated increasingly accurate perturbative studies of gravitational dynamics in compact binaries. This in turn has enabled more detailed analyses of the dynamical black holes in these systems. For example, Pound et al. [Phys. Rev. Lett. 124, 021101 (2020)] recently computed the surface area of a Schwarzschild black hole's apparent horizon, perturbed by an orbiting body, to second order in the binary's mass ratio. In this paper, we take that as the starting point for a comprehensive study of a perturbed Schwarzschild black hole's apparent and event horizon at second perturbative order, deriving generic formulas for the first- and second-order corrections to the horizons' radial profiles, surface areas, Hawking masses, and intrinsic curvatures. We find that the two horizons are remarkably similar, and that any teleological behavior of the event horizon is suppressed in several ways. Critically, we establish that at all orders, the perturbed event horizon in a small-mass-ratio binary is effectively localized in time. Even more pointedly, the event horizon is identical to the apparent horizon at linear order regardless of the source of perturbation, implying that the seemingly teleological "tidal lead", previously observed in linearly perturbed event horizons, is not genuinely teleological in origin. The two horizons do generically differ at second order, but their Hawking masses remain identical, implying that the event horizon obeys the same energy-flux balance law as the apparent horizon. At least in the case of a binary system, the difference between their surface areas remains extremely small even in the late stages of inspiral. In the course of our analysis, we also numerically illustrate puzzling behaviour in the black hole's motion around the binary's center of mass., Comment: 39 pages, 4 figures. v2 modifies the title, abstract, and Conclusion to better highlight the content. v3 corrects minor typos

The computational cost of inspiral and merger simulations for black-hole binaries increases in inverse proportion to the square of the mass ratio $q:=m_2/m_1\leq 1$. One factor of $q$ comes from the number of orbital cycles, which is proportional to $1/q$, and another is associated with the required number of time steps per orbit, constrained (via the Courant-Friedrich-Lewy condition) by the need to resolve the two disparate length scales. This problematic scaling makes simulations progressively less tractable at smaller $q$. Here we propose and explore a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysical range ($10^{-4} \lesssim q\lesssim 10^{-2}$), where purely perturbative methods may not be adequate. A region of radius much larger than $m_2$ around the smaller object is excised from the numerical domain, and replaced with an analytical model approximating a tidally deformed black hole. The analytical model involves certain a priori unknown parameters, associated with unknown bits of physics together with gauge-adjustment terms; these are dynamically determined by matching to the numerical solution outside the excision region. In this paper we develop the basic idea and apply it to a toy model of a scalar charge in a circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field in a 1+1D framework. Our main goal here is to explore the utility and properties of different matching strategies, and to this end we develop two independent implementations, a finite-difference one and a spectral one. We discuss the extension of our method to a full 3D numerical evolution and to gravity., Comment: 24 pages, 18 figures

Inspirals of stellar-mass objects into massive black holes will be important sources for the space-based gravitational-wave detector LISA. Modelling these systems requires calculating the metric perturbation due to a point particle orbiting a Kerr black hole. Currently, the linear perturbation is obtained with a metric reconstruction procedure that puts it in a "no-string" radiation gauge which is singular on a surface surrounding the central black hole. Calculating dynamical quantities in this gauge involves a subtle procedure of "gauge completion" as well as cancellations of very large numbers. The singularities in the gauge also lead to pathological field equations at second perturbative order. In this paper we re-analyze the point-particle problem in Kerr using the corrector-field reconstruction formalism of Green, Hollands, and Zimmerman (GHZ). We clarify the relationship between the GHZ formalism and previous reconstruction methods, showing that it provides a simple formula for the "gauge completion". We then use it to develop a new method of computing the metric in a more regular gauge: a Teukolsky puncture scheme. This scheme should ameliorate the problem of large cancellations, and by constructing the linear metric perturbation in a sufficiently regular gauge, it should provide a first step toward second-order self-force calculations in Kerr. Our methods are developed in generality in Kerr, but we illustrate some key ideas and demonstrate our puncture scheme in the simple setting of a static particle in Minkowski spacetime., Comment: 73 pages in CQG format, 2 figures. v2 has an expanded introduction

Within the framework of self-force theory, we compute the gravitational-wave energy flux through second order in the mass ratio for compact binaries in quasicircular orbits. Our results are consistent with post-Newtonian calculations in the weak field and they agree remarkably well with numerical-relativity simulations of comparable-mass binaries in the strong field. We also find good agreement for binaries with a spinning secondary or a slowly spinning primary. Our results are key for accurately modelling extreme-mass-ratio inspirals and will be useful in modelling intermediate-mass-ratio systems., Comment: 5 pages + supplemental material, 7 figures; Updated to reflect published version

Extreme-mass-ratio inspirals (EMRIs) will be key sources for LISA. However, accurately extracting system parameters from a detected EMRI waveform will require self-force calculations at second order in perturbation theory, which are still in a nascent stage. One major obstacle in these calculations is the strong divergences that are encountered on the worldline of the small object. Previously, it was shown by one of us [Phys. Rev. D 95, 104056 (2017)] that a class of "highly regular" gauges exist in which the singularities have a qualitatively milder form, promising to enable more efficient numerical calculations. Here we derive expressions for the metric perturbation in this class of gauges, in a local expansion in powers of distance $r$ from the worldline, to sufficient order in $r$ for numerical implementation in a puncture scheme. Additionally, we use the highly regular class to rigorously derive a distributional source for the second-order field and a pointlike second-order stress-energy tensor (the Detweiler stress-energy) for the small object. This makes it possible to calculate the second-order self-force using mode-sum regularisation rather than the more cumbersome puncture schemes that have been necessary previously. Although motivated by EMRIs, our calculations are valid in an arbitrary vacuum background, and they may help clarify the interpretation of point masses and skeleton sources in general relativity more broadly., Comment: 35 pages. Version to appear in PRD

Much of the success of gravitational-wave astronomy rests on perturbation theory. Historically, perturbative analysis of gravitational-wave sources has largely focused on post-Newtonian theory. However, strong-field perturbation theory is essential in many cases such as the quasinormal ringdown following the merger of a binary system, tidally perturbed compact objects, and extreme-mass-ratio inspirals. In this review, motivated primarily by small-mass-ratio binaries but not limited to them, we provide an overview of essential methods in (i) black hole perturbation theory, (ii) orbital mechanics in Kerr spacetime, and (iii) gravitational self-force theory. Our treatment of black hole perturbation theory covers most common methods, including the Teukolsky and Regge-Wheeler-Zerilli equations, methods of metric reconstruction, and Lorenz-gauge formulations, presenting them in a new consistent and self-contained form. Our treatment of orbital mechanics covers quasi-Keplerian and action-angle descriptions of bound geodesics and accelerated orbits, osculating geodesics, near-identity averaging transformations, multiscale expansions, and orbital resonances. Our summary of self-force theory's foundations is brief, covering the main ideas and results of matched asymptotic expansions, local expansion methods, puncture schemes, and point particle descriptions. We conclude by combining the above methods in a multiscale expansion of the perturbative Einstein equations, leading to adiabatic and post-adiabatic evolution schemes. Our presentation is intended primarily as a reference for practitioners but includes a variety of new results. In particular, we present the first complete post-adiabatic waveform-generation framework for generic (nonresonant) orbits in Kerr., Comment: 121 pages, 1 figure. Invited chapter for "Handbook of Gravitational Wave Astronomy" (Eds. C. Bambi, S. Katsanevas, and K. Kokkotas; Springer, Singapore, 2021). v4 typos corrected

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., Comment: 41 pages, 3 figures. The third version corrects more typos

Self-force theory is the leading method of modeling extreme-mass-ratio inspirals (EMRIs), key sources for the gravitational-wave detector LISA. It is well known that for an accurate EMRI model, second-order self-force effects are critical, but calculations of these effects have been beset by obstacles. In this letter we present the first implementation of a complete scheme for second-order self-force computations, specialized to the case of quasicircular orbits about a Schwarzschild black hole. As a demonstration, we calculate the gravitational binding energy of these binaries., Comment: 6 pages + supplemental material, 3 figures, revised title and two added references

[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., Comment: 79 pages, 11 figures; invited by Reports on Progress in Physics. v2 contains minor corrections and it is the published version

When a small, uncharged, compact object is immersed in an external background spacetime, at zeroth order in its mass it moves as a test particle in the background. At linear order, its own gravitational field alters the geometry around it, and it moves instead as a test particle in a certain effective metric satisfying the linearized vacuum Einstein equation. In the letter [Phys. Rev. Lett. 109, 051101 (2012)], using a method of matched asymptotic expansions, I showed that the same statement holds true at second order: if the object's leading-order spin and quadrupole moment vanish, then through second order in its mass it moves on a geodesic of a certain smooth, locally causal vacuum metric defined in its local neighbourhood. Here I present the complete details of the derivation of that result. In addition, I extend the result, which had previously been derived in gauges smoothly related to Lorenz, to a class of highly regular gauges that should be optimal for numerical self-force computations., Comment: 31 pages. Corrected typos. To be published in PRD

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., Comment: 27 pages, Minor typos corrected

Second-order self-force computations, which will be essential in modeling extreme-mass-ratio inspirals, involve two major new difficulties that were not present at first order. One is the problem of large scales, discussed in [Phys. Rev. D 92, 104047 (2015)]. Here we discuss the second difficulty, which occurs instead on small scales: if we expand the field equations in spherical harmonics, then because the first-order field contains a singularity, we require an arbitrarily large number of first-order modes to accurately compute even a single second-order mode. This is a generic feature of nonlinear field equations containing singularities, allowing us to study it in the simple context of a scalar toy model in flat space. Using that model, we illustrate the problem and demonstrate a robust strategy for overcoming it., Comment: 18 pages, 6 figures, submitted to PRD

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., Comment: 25 pages. Submitted to PRD

In recent years, asymptotic approximation schemes have been developed to describe the motion of a small compact object through a vacuum background to any order in perturbation theory. The schemes are based on rigorous methods of matched asymptotic expansions, which account for the object's finite size, require no "regularization" of divergent quantities, and are valid for strong fields and relativistic speeds. Up to couplings of the object's multipole moments to the external background curvature, these schemes have established that at least through second order in perturbation theory, the object's motion satisfies a generalized equivalence principle: it moves on a geodesic of a certain smooth metric satisfying the vacuum Einstein equation. I describe the foundations of this result, particularly focusing on the fundamental notion of how a small object's motion is represented in perturbation theory. The three common representations of perturbed motion are (i) the "Gralla-Wald" description in terms of small deviations from a reference geodesic, (ii) the "self-consistent" description in terms of a worldline that obeys a self-accelerated equation of motion, and (iii) the "osculating geodesics" description, which utilizes both (i) and (ii). Because of the coordinate freedom in general relativity, any coordinate description of motion in perturbation theory is intimately related to the theory's gauge freedom. In this review I describe asymptotic solutions of the Einstein equations adapted to each of the three representations of motion, and I discuss the gauge freedom associated with each. I conclude with a discussion of how gauge freedom must be refined in the context of long-term dynamics., Comment: 100 pages, 3 figures, book contribution to Equations of Motion in Relativistic Gravity, edited by D. Puetzfeld et al., Fundamental Theories of Physics 179, Springer, 2015

Through second order in perturbative general relativity, a small compact object in an external vacuum spacetime obeys a generalized equivalence principle: although it is accelerated with respect to the external background geometry, it is in free fall with respect to a certain \emph{effective} vacuum geometry. However, this single principle takes very different mathematical forms, with very different behaviors, depending on how one treats perturbed motion. Furthermore, any description of perturbed motion can be altered by a gauge transformation. In this paper, I clarify the relationship between two treatments of perturbed motion and the gauge freedom in each. I first show explicitly how one common treatment, called the Gralla-Wald approximation, can be derived from a second, called the self-consistent approximation. I next present a general treatment of smooth gauge transformations in both approximations, in which I emphasise that the approximations' governing equations can be formulated in an invariant manner. All of these analyses are carried through second perturbative order, but the methods are general enough to go to any order. Furthermore, the tools I develop, and many of the results, should have broad applicability to any description of perturbed motion, including osculating-geodesic and two-timescale descriptions., Comment: 26 pages, 3 figures. Minor corrections. Equations (120) and (126) are more general than in PRD version

Using black hole perturbation theory and arbitrary-precision computer algebra, we obtain the post-Newtonian (pN) expansions of the linear-in-mass-ratio corrections to the spin-precession angle and tidal invariants for a particle in circular orbit around a Schwarzschild black hole. We extract coefficients up to 20pN order from numerical results that are calculated with an accuracy greater than 1 part in $10^{500}$. These results can be used to calibrate parameters in effective-one-body models of compact binaries, specifically the spin-orbit part of the effective Hamiltonian and the dynamically significant tidal part of the main radial potential of the effective metric. Our calculations are performed in a radiation gauge, which is known to be singular away from the particle. To overcome this irregularity, we define suitable Detweiler-Whiting singular and regular fields in this gauge, and we devise a rigorous mode-sum regularization method to compute the invariants constructed from the regular field.

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., Comment: 17 pages, no figures. Proceedings of the LISA Symposium X, to be published at the Journal of Physics

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., Comment: 20 pages, 1 figure, version to appear in PRD, additional references, minor typos corrected

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., Comment: 23 pages, 1 figure, final version to be published in Phys Rev D

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., Comment: 46 pages, 1 figure, 1 table. Minor corrections. Final version to be published in Phys. Rev. D

In order to extract physical parameters from the waveform of an extreme-mass-ratio binary, one requires a second-order--accurate description of the motion of the smaller of the two objects in the binary. Using a method of matched asymptotic expansions, I derive the second-order equation of motion of a small, nearly spherical and non-rotating compact object in an arbitrary vacuum spacetime. I find that the motion is geodesic in a certain locally defined effective metric satisfying the vacuum Einstein equation through second order, and I outline a method of numerically calculating this effective metric., Comment: 3 pages, contribution to the proceedings of the 13th Marcel Grossman Meeting, summarizes and discusses material in 1201.5089 and 1206.6538