Showing 1–50 of 52 results for author: Warburton, N

Comparison of 4.5PN and 2SF gravitational energy fluxes from quasicircular compact binaries

Authors: Niels Warburton, Barry Wardell, David Trestini, Quentin Henry, Adam Pound, Luc Blanchet, Leanne Durkan, Guillaume Faye, Jeremy Miller

Abstract: 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 re… ▽ More 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. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Comments: 15 pages, 4 figures, modes of PN flux provided in ancillary file

Probing fundamental physics with Extreme Mass Ratio Inspirals: a full Bayesian inference for scalar charge

Authors: Lorenzo Speri, Susanna Barsanti, Andrea Maselli, Thomas P. Sotiriou, Niels Warburton, Maarten van de Meent, Alvin J. K. Chua, Ollie Burke, Jonathan Gair

Abstract: Extreme Mass Ratio Inspirals (EMRIs) are key sources for the future space-based gravitational wave detector LISA, and are considered promising probes of fundamental physics. Here, we present the first complete Bayesian analysis of EMRI signals in theories with an additional massless scalar, which could arise in an extension of General Relativity or of the Standard Model of Particle Physics. We dev… ▽ More Extreme Mass Ratio Inspirals (EMRIs) are key sources for the future space-based gravitational wave detector LISA, and are considered promising probes of fundamental physics. Here, we present the first complete Bayesian analysis of EMRI signals in theories with an additional massless scalar, which could arise in an extension of General Relativity or of the Standard Model of Particle Physics. We develop a waveform model accurate at adiabatic order for equatorial eccentric orbits around spinning black holes. Using full Bayesian inference, we forecast LISA's ability to probe the presence of new fundamental fields with EMRI observations. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Fast inspirals and the treatment of orbital resonances

Authors: Philip Lynch, Vojtěch Witzany, Maarten van de Meent, Niels Warburton

Abstract: Extreme mass ratio inspirals (EMRIs), where a compact object orbits a massive black hole, are a key source of gravitational waves for the future Laser Interferometer Space Antenna (LISA). Due to their small mass ratio, ($ε\sim 10^{-4}$--$10^{-7}$), the binary evolves slowly and EMRI signals will be in-band for years. Additionally, astrophysical EMRIs are expected to have complex dynamics featuring… ▽ More Extreme mass ratio inspirals (EMRIs), where a compact object orbits a massive black hole, are a key source of gravitational waves for the future Laser Interferometer Space Antenna (LISA). Due to their small mass ratio, ($ε\sim 10^{-4}$--$10^{-7}$), the binary evolves slowly and EMRI signals will be in-band for years. Additionally, astrophysical EMRIs are expected to have complex dynamics featuring both spin-precession and eccentricity. A standard approach to modelling these inspirals is via the method of osculating geodesics (OG) which we employ along with a toy model for the gravitational self-force. Using this method requires resolving tens of thousands radial and polar orbital librations over the long duration of the signal which makes the inspiral trajectory expensive to compute. In this work we accelerate these calculations by employing Near-Identity (averaging) Transformations. However, this averaging technique breaks down at orbital resonances where the radial and polar frequencies are an integer ratio of each other. Thus, we switch to a partial averaging transformation in the vicinity of the resonance where the dynamics are characterised by the slow evolution of the so-called "resonant phase". Additionally, we develop an optimal switching criterion to minimise the computation time while maximising accuracy. We find the error in the waveform phase is improved from $\mathcal{O}(ε^{-1/2})$ in the fully averaged scheme to $\mathcal{O}(ε^{4/7})$ in the switching scheme. At the same time, this scheme improves the scaling of the computation time from being inversely proportional to $ε$ using OG, to a very weak scaling with $ε$. This results in a speed-up of at least two orders of magnitude for LISA EMRIs with room for further optimisation. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 48 pages, 11 figures

Worldtube puncture scheme for first- and second-order self-force calculations in the Fourier domain

Authors: Jeremy Miller, Benjamin Leather, Adam Pound, Niels Warburton

Abstract: 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 in the Fourier domain. In this paper we present that method, which involves dividing the… ▽ More 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 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. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Comments: 40 pages, 10 figures

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: 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 cosmologic… ▽ More 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. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 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

Accuracy Requirements: Assessing the Importance of First Post-Adiabatic Terms for Small-Mass-Ratio Binaries

Authors: Ollie Burke, Gabriel Andres Piovano, Niels Warburton, Philip Lynch, Lorenzo Speri, Chris Kavanagh, Barry Wardell, Adam Pound, Leanne Durkan, Jeremy Miller

Abstract: 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 $ε\gtrsim 10^{-5}$ for quasi circular… ▽ More 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 $ε\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 $ε\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. △ Less

Submitted 23 May, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 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

Terrestrial Very-Long-Baseline Atom Interferometry: Workshop Summary

Authors: Sven Abend, Baptiste Allard, Iván Alonso, John Antoniadis, Henrique Araujo, Gianluigi Arduini, Aidan Arnold, Tobias Aßmann, Nadja Augst, Leonardo Badurina, Antun Balaz, Hannah Banks, Michele Barone, Michele Barsanti, Angelo Bassi, Baptiste Battelier, Charles Baynham, Beaufils Quentin, Aleksandar Belic, Ankit Beniwal, Jose Bernabeu, Francesco Bertinelli, Andrea Bertoldi, Ikbal Ahamed Biswas, Diego Blas , et al. (228 additional authors not shown)

Abstract: This document presents a summary of the 2023 Terrestrial Very-Long-Baseline Atom Interferometry Workshop hosted by CERN. The workshop brought together experts from around the world to discuss the exciting developments in large-scale atom interferometer (AI) prototypes and their potential for detecting ultralight dark matter and gravitational waves. The primary objective of the workshop was to lay… ▽ More This document presents a summary of the 2023 Terrestrial Very-Long-Baseline Atom Interferometry Workshop hosted by CERN. The workshop brought together experts from around the world to discuss the exciting developments in large-scale atom interferometer (AI) prototypes and their potential for detecting ultralight dark matter and gravitational waves. The primary objective of the workshop was to lay the groundwork for an international TVLBAI proto-collaboration. This collaboration aims to unite researchers from different institutions to strategize and secure funding for terrestrial large-scale AI projects. The ultimate goal is to create a roadmap detailing the design and technology choices for one or more km-scale detectors, which will be operational in the mid-2030s. The key sections of this report present the physics case and technical challenges, together with a comprehensive overview of the discussions at the workshop together with the main conclusions. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: Summary of the Terrestrial Very-Long-Baseline Atom Interferometry Workshop held at CERN: https://indico.cern.ch/event/1208783/

Adiabatic inspirals under electromagnetic radiation reaction on Kerr spacetime

Authors: Ethan J German, Kevin Cunningham, Visakan Balakumar, Niels Warburton, Sam R Dolan

Abstract: A compact body in orbit about a black hole loses orbital energy and angular momentum through radiation-reaction processes, inspiralling towards the black hole until a final plunge. Here we consider a scenario with a charged compact body in which fluxes of electromagnetic radiation drive this inspiral. We calculate trajectories in the $(p,e)$ plane for inspirals in the equatorial plane of a rotatin… ▽ More A compact body in orbit about a black hole loses orbital energy and angular momentum through radiation-reaction processes, inspiralling towards the black hole until a final plunge. Here we consider a scenario with a charged compact body in which fluxes of electromagnetic radiation drive this inspiral. We calculate trajectories in the $(p,e)$ plane for inspirals in the equatorial plane of a rotating black hole within the adiabatic (orbital-averaged-dissipative) approximation. We make comparisons with a non-relativistic Keplerian approximation based on the Abraham-Lorentz force law, and with standard gravitational-wave driven scenarios. We find that EM-driven inspirals are less efficiently circularized (i.e. orbits remain more eccentric at the point of plunge) than their gravitational counterparts, and we quantify the effect of black hole spin. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 23 pages, 7 figures

Fast and Fourier: Extreme Mass Ratio Inspiral Waveforms in the Frequency Domain

Authors: Lorenzo Speri, Michael L. Katz, Alvin J. K. Chua, Scott A. Hughes, Niels Warburton, Jonathan E. Thompson, Christian E. A. Chapman-Bird, Jonathan R. Gair

Abstract: Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the firs… ▽ More Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the first ready-to-use Schwarzschild eccentric EMRI waveform implementation in the frequency domain for use with either graphics processing units (GPUs) or central processing units (CPUs). We present the overall waveform implementation and test the accuracy and performance of the frequency domain waveforms against the time domain implementation. On GPUs, the frequency domain waveform takes in median $0.044$ seconds to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses $\geq 2 \times 10^6 \,{\rm M_\odot}$ and initial eccentricities $e_0 > 0.2$. On CPUs, the median waveform evaluation time is $5$ seconds, and it is five times faster in the frequency domain than in the time domain. Using a sparser frequency array can further speed up the waveform generation, reaching up to $ 0.3$ seconds. This enables us to perform, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs. Future EMRI models which encompass wider source characteristics (particularly black hole spin and generic orbit geometries) will require significantly more harmonics. Frequency-domain models will be essential analysis tools for these astrophysically realistic and important signals. △ Less

Submitted 15 January, 2024; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 23 pages, 6 figures

Journal ref: Front. Appl. Math. Stat. 9 (2023)

Applying the effective-source approach to frequency-domain self-force calculations for eccentric orbits

Authors: Benjamin Leather, Niels Warburton

Abstract: Extreme mass-ratio inspirals (EMRIs) are expected to have considerable eccentricity when emitting gravitational waves (GWs) in the LISA band. Developing GW templates that remain phase accurate over these long inspirals requires the use of second-order self-force theory and practical second-order self-force calculations are now emerging for quasi-circular EMRIs. These calculations rely on effective… ▽ More Extreme mass-ratio inspirals (EMRIs) are expected to have considerable eccentricity when emitting gravitational waves (GWs) in the LISA band. Developing GW templates that remain phase accurate over these long inspirals requires the use of second-order self-force theory and practical second-order self-force calculations are now emerging for quasi-circular EMRIs. These calculations rely on effective-source regularization techniques in the frequency domain that presently are specialized to circular orbits. Here we make a first step towards more generic second-order calculations by extending the frequency domain effective-source approach to eccentric orbits. In order to overcome the slow convergence of the Fourier sum over radial modes, we develop a new extended effective-sources approach which builds upon the method of extended particular solutions. To demonstrate our new computational technique we apply it a toy scalar-field problem which is conceptually similar to the gravitational case. △ Less

Submitted 24 October, 2023; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 23 pages, 16 figures; updated to reflect published version

Journal ref: Phys. Rev. D 108, 084045 (2023)

Self-forced inspirals with spin-orbit precession

Authors: Philip Lynch, Maarten van de Meent, Niels Warburton

Abstract: We develop the first model for extreme mass-ratio inspirals (EMRIs) with misaligned angular momentum and primary spin, and zero eccentricity -- also known as quasi-spherical inspirals -- evolving under the influence of the first-order in mass ratio gravitational self-force. The forcing terms are provided by an efficient spectral interpolation of the first-order gravitational self-force in the outg… ▽ More We develop the first model for extreme mass-ratio inspirals (EMRIs) with misaligned angular momentum and primary spin, and zero eccentricity -- also known as quasi-spherical inspirals -- evolving under the influence of the first-order in mass ratio gravitational self-force. The forcing terms are provided by an efficient spectral interpolation of the first-order gravitational self-force in the outgoing radiation gauge. In order to speed up the calculation of the inspiral we apply a near-identity (averaging) transformation to eliminate all dependence of the orbital phases from the equations of motion while maintaining all secular effects of the first-order gravitational self-force at post-adiabatic order. The resulting solutions are defined with respect to `Mino time' so we perform a second averaging transformation so the inspiral is parametrized in terms of Boyer-Lindquist time, which is more convenient of LISA data analysis. We also perform a similar analysis using the two-timescale expansion and find that using either approach yields self-forced inspirals that can be evolved to sub-radian accuracy in less than a second. The dominant contribution to the inspiral phase comes from the adiabatic contributions and so we further refine our self-force model using information from gravitational wave flux calculations. The significant dephasing we observe between the lower and higher accuracy models highlights the importance of accurately capturing adiabatic contributions to the phase evolution. △ Less

Submitted 4 June, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 25 pages, 9 figures, published in Physical Review D

Journal ref: Phys. Rev. D 109, 084072 (2024)

Enhancing the SEOBNRv5 effective-one-body waveform model with second-order gravitational self-force fluxes

Authors: Maarten van de Meent, Alessandra Buonanno, Deyan P. Mihaylov, Serguei Ossokine, Lorenzo Pompili, Niels Warburton, Adam Pound, Barry Wardell, Leanne Durkan, Jeremy Miller

Abstract: 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 qu… ▽ More 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. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Comparing second-order gravitational self-force and effective one body waveforms from inspiralling, quasi-circular and nonspinning black hole binaries II: the large-mass-ratio case

Authors: Angelica Albertini, Alessandro Nagar, Adam Pound, Niels Warburton, Barry Wardell, Leanne Durkan, Jeremy Miller

Abstract: 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 d… ▽ More 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_ω\equiv ω^2/\dotω$, where $ω$ 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 △ Less

Submitted 2 November, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 14 pages, 8 figures, published on Phys. Rev. D

Journal ref: Phys.Rev.D 106 (2022) 8, 084062

Comparing second-order gravitational self-force, numerical relativity and effective one body waveforms from inspiralling, quasi-circular and nonspinning black hole binaries

Authors: Angelica Albertini, Alessandro Nagar, Adam Pound, Niels Warburton, Barry Wardell, Leanne Durkan, Jeremy Miller

Abstract: 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 impl… ▽ More 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_ω(ω)\equiv ω^2/\dotω$, where $ω$ 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. △ Less

Submitted 2 November, 2022; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: 31 pages, 19 figures, published on Phys. Rev. D

Journal ref: Phys.Rev.D 106 (2022) 8, 084061

Slow evolution of the metric perturbation due to a quasicircular inspiral into a Schwarzschild black hole

Authors: Leanne Durkan, Niels Warburton

Abstract: Extreme mass-ratio inspirals (EMRIs) are one of the most highly anticipated sources of gravitational radiation novel to detection by millihertz space-based detectors. To accurately estimate the parameters of EMRIs and perform precision tests of general relativity, their models should incorporate self-force theory through second-order in the small mass ratio. Due to their extreme mass ratio, EMRIs… ▽ More Extreme mass-ratio inspirals (EMRIs) are one of the most highly anticipated sources of gravitational radiation novel to detection by millihertz space-based detectors. To accurately estimate the parameters of EMRIs and perform precision tests of general relativity, their models should incorporate self-force theory through second-order in the small mass ratio. Due to their extreme mass ratio, EMRIs inspiral slowly when sufficiently far from merger. As such, the slow evolution of the first-order metric perturbation contributes to the source for the second-order metric perturbation, and must be accounted for in EMRI waveform models. In this paper we calculate the slow evolution of the first-order metric perturbation in the Lorenz gauge for quasicircular orbits on a Schwarzschild background in the frequency domain. Lorenz gauge solutions to the first-order metric perturbation and its slow evolution are obtained via a gauge transformation from Regge-Wheeler gauge solutions. The slow evolution of Regge-Wheeler and Zerilli master functions, in addition to a gauge field are determined using the method of partial annihilators. △ Less

Submitted 18 May, 2023; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 23 pages, (15 pages with appendices and references), 7 figures, 1 table

Hyperboloidal method for frequency-domain self-force calculations

Authors: Rodrigo Panosso Macedo, Benjamin Leather, Niels Warburton, Barry Wardell, Anıl Zenginoğlu

Abstract: Gravitational self-force theory is the leading approach for modeling gravitational wave emission from small mass-ratio compact binaries. This method perturbatively expands the metric of the binary in powers of the mass ratio. The source for the perturbations depends on the orbital configuration, calculational approach, and the order of the perturbative expansion. These sources fall into three broa… ▽ More Gravitational self-force theory is the leading approach for modeling gravitational wave emission from small mass-ratio compact binaries. This method perturbatively expands the metric of the binary in powers of the mass ratio. The source for the perturbations depends on the orbital configuration, calculational approach, and the order of the perturbative expansion. These sources fall into three broad classes: (i) distributional, (ii) worldtube, and (iii) unbounded support. The latter, in particular, is important for emerging second-order (in the mass ratio) calculations. Traditional frequency domain approaches employ the variation of parameters method and compute the perturbation on standard time slices with numerical boundary conditions supplied at finite radius from series expansions of the asymptotic behavior. This approach has been very successful, but the boundary conditions calculations are tedious, and the approach is not well suited to unbounded sources where homogeneous solutions must be computed at all radii. This work develops an alternative approach where hyperboloidal slices foliate the spacetime, and compactifying coordinates simplify the boundary treatment. We implement this approach with a multi-domain spectral solver with analytic mesh refinement and use the scalar-field self-force on circular orbits around a Schwarzschild black hole as an example problem. The method works efficiently for all three source classes encountered in self-force calculations and has distinct advantages over the traditional approach. For example, our code efficiently computes the perturbation for orbits with extremely large orbital radii ($r_{p}>10^5M$) or modes with very high spherical harmonic mode index ($\ell \ge 100$). Our results indicate that hyperboloidal methods can play an essential role in self-force calculations. △ Less

Submitted 1 June, 2022; v1 submitted 3 February, 2022; originally announced February 2022.

Comments: 25 pages, 19 figures; updated to reflect published version

Journal ref: Phys. Rev. D 105, 104033 (2022)

Gravitational waveforms for compact binaries from second-order self-force theory

Authors: Barry Wardell, Adam Pound, Niels Warburton, Jeremy Miller, Leanne Durkan, Alexandre Le Tiec

Abstract: 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 fr… ▽ More 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. △ Less

Submitted 29 March, 2023; v1 submitted 22 December, 2021; originally announced December 2021.

Eccentric self-forced inspirals into a rotating black hole

Authors: Philip Lynch, Maarten van de Meent, Niels Warburton

Abstract: We develop the first model for extreme mass-ratio inspirals (EMRIs) into a rotating massive black hole driven by the gravitational self-force. Our model is based on an action angle formulation of the method of osculating geodesics for eccentric, equatorial (i.e., spin-aligned) motion in Kerr spacetime. The forcing terms are provided by an efficient spectral interpolation of the first-order gravita… ▽ More We develop the first model for extreme mass-ratio inspirals (EMRIs) into a rotating massive black hole driven by the gravitational self-force. Our model is based on an action angle formulation of the method of osculating geodesics for eccentric, equatorial (i.e., spin-aligned) motion in Kerr spacetime. The forcing terms are provided by an efficient spectral interpolation of the first-order gravitational self-force in the outgoing radiation gauge. We apply a near-identity (averaging) transformation to eliminate all dependence of the orbital phases from the equations of motion, while maintaining all secular effects of the first-order gravitational self-force at post-adiabatic order. This implies that the model can be evolved without having to resolve all $\mathcal{O}(10^5)$ orbit cycles of an EMRI, yielding an inspiral model that can be evaluated in less than a second for any mass-ratio. In the case of a non-rotating central black hole, we compare inspirals evolved using self-force data computed in the Lorenz and radiation gauges. We find that the two gauges generally produce differing inspirals with a deviation of comparable magnitude to the conservative self-force correction. This emphasizes the need for including the (currently unknown) dissipative second order self-force to obtain gauge independent, post-adiabatic waveforms. △ Less

Submitted 3 August, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: 46 pages, 12 figures

Journal ref: Class. Quantum Grav. 39:145004 (2022)

Survey of gravitational wave memory in intermediate mass ratio binaries

Authors: Tousif Islam, Scott E. Field, Gaurav Khanna, Niels Warburton

Abstract: The non-linear gravitational wave (GW) memory effect is a distinct prediction in general relativity. While the effect has been well studied for comparable mass binaries, it has mostly been overlooked for intermediate mass ratio inspirals (IMRIs). We offer a comprehensive analysis of the phenomenology and detectability of memory effects, including contributions from subdominant harmonic modes, in h… ▽ More The non-linear gravitational wave (GW) memory effect is a distinct prediction in general relativity. While the effect has been well studied for comparable mass binaries, it has mostly been overlooked for intermediate mass ratio inspirals (IMRIs). We offer a comprehensive analysis of the phenomenology and detectability of memory effects, including contributions from subdominant harmonic modes, in heavy IMRIs consisting of a stellar mass black hole and an intermediate mass black hole. When formed through hierarchical mergers, for example when a GW190521-like remnant captures a stellar mass black hole, IMRI systems have a large total mass, large spin on the primary, and possibly residual eccentricity; features that potentially raise the prospect for memory detection. We compute both the displacement and spin non-linear GW memory from the $m \neq 0$ gravitational waveforms computed within a black hole perturbation theory framework that is partially calibrated to numerical relativity waveforms. We probe the dependence of memory effects on mass ratio, spin, and eccentricity and consider the detectability of a memory signal from IMRIs using current and future GW detectors. We find that (i) while eccentricity introduces additional features in both displacement and spin memory, it does not appreciatively change the prospects of detectability, (ii) including higher modes into the memory computation can increase singal-to-noise (SNR) values by about 7\% in some cases, (iii) the SNR from displacement memory dramatically increases as the spin approaches large, positive values, (iv) spin memory from heavy IMRIs would, however, be difficult to detect with future generation detectors even from highly spinning systems. Our results suggest that hierarchical binary black hole mergers may be a promising source for detecting memory and could favorably impact memory forecasts. △ Less

Submitted 24 July, 2023; v1 submitted 2 September, 2021; originally announced September 2021.

Comments: 18 pages, 17 figures, matches the version published in PRD

Journal ref: Phys. Rev. D 108, 024046 (2023)

Gravitational-wave energy flux for compact binaries through second order in the mass ratio

Authors: Niels Warburton, Adam Pound, Barry Wardell, Jeremy Miller, Leanne Durkan

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 9 October, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: 5 pages + supplemental material, 7 figures; Updated to reflect published version

Journal ref: Phys. Rev. Lett. 127, 151102 (2021)

The missing link in gravitational-wave astronomy: A summary of discoveries waiting in the decihertz range

Authors: Manuel Arca Sedda, Christopher P L Berry, Karan Jani, Pau Amaro-Seoane, Pierre Auclair, Jonathon Baird, Tessa Baker, Emanuele Berti, Katelyn Breivik, Chiara Caprini, Xian Chen, Daniela Doneva, Jose M Ezquiaga, K E Saavik Ford, Michael L Katz, Shimon Kolkowitz, Barry McKernan, Guido Mueller, Germano Nardini, Igor Pikovski, Surjeet Rajendran, Alberto Sesana, Lijing Shao, Nicola Tamanini, Niels Warburton , et al. (3 additional authors not shown)

Abstract: Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of m… ▽ More Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of massive black holes binaries. Between the $\sim 10$-$10^3~\mathrm{Hz}$ band of ground-based observatories and the $\sim10^{-4}$-$10^{-1}~\mathrm{Hz}$ band of LISA lies the uncharted decihertz gravitational-wave band. We propose a Decihertz Observatory to study this frequency range, and to complement observations made by other detectors. Decihertz observatories are well suited to observation of intermediate-mass ($\sim10^2$-$10^4 M_\odot$) black holes; they will be able to detect stellar-mass binaries days to years before they merge, providing early warning of nearby binary neutron star mergers and measurements of the eccentricity of binary black holes, and they will enable new tests of general relativity and the Standard Model of particle physics. Here we summarise how a Decihertz Observatory could provide unique insights into how black holes form and evolve across cosmic time, improve prospects for both multimessenger astronomy and multiband gravitational-wave astronomy, and enable new probes of gravity, particle physics and cosmology. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: 13 pages, 1 figure. Published in Experimental Astronomy. Summarising white paper arXiv:1908.11375

Journal ref: Experimental Astronomy; 51(3):1427-1440; 2021

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

Authors: Michael L. Katz, Alvin J. K. Chua, Lorenzo Speri, Niels Warburton, Scott A. Hughes

Abstract: We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI w… ▽ More We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI waveforms. Because the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, we also introduce an improved Augmented Analytic Kludge (AAK) model that describes generic Kerr inspirals. Both waveform models can be accelerated using graphics processing unit (GPU) hardware. With the GPU-accelerated waveforms in hand, a variety of studies are performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation. We find relativistic EMRI waveform templates can be generated with fewer harmonic modes ($\sim10-100$) without biasing signal extraction. However, we show for the first time that extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space. △ Less

Submitted 27 September, 2021; v1 submitted 9 April, 2021; originally announced April 2021.

Comments: 26 pages, 12 Figures, FastEMRIWaveforms Package: bhptoolkit.org/FastEMRIWaveforms/

Journal ref: Phys. Rev. D 104, 064047 (2021)

Adiabatic waveforms for extreme mass-ratio inspirals via multivoice decomposition in time and frequency

Authors: Scott A. Hughes, Niels Warburton, Gaurav Khanna, Alvin J. K. Chua, Michael L. Katz

Abstract: We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, a property which can be exploited to efficien… ▽ More We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, a property which can be exploited to efficiently produce waveform models that faithfully encode the properties of EMRI systems. We look at examples for a range of different orbital geometries: spherical orbits, equatorial eccentric orbits, and one example of generic (inclined and eccentric) orbits. To our knowledge, this is the first calculation of a generic EMRI waveform that uses strong-field radiation reaction. We examine waveforms in both the time and frequency domains. Although EMRIs evolve slowly enough that the stationary phase approximation (SPA) to the Fourier transform is valid, the SPA calculation must be done to higher order for some voices, since their instantaneous frequency can change from chirping forward ($\dot f > 0$) to chirping backward ($\dot f < 0$). The approach we develop can eventually be extended to more complete EMRI waveform models, for example to include effects neglected by the adiabatic approximation such as the conservative self force and spin-curvature coupling. △ Less

Submitted 28 March, 2024; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 31 pages, 17 figures. Update: corrected a sign error in a term in Eq. (B5). This factor was correct in the related code, it was purely a typo in the manuscript. An erratum to PRD has been submitted. Many thanks to Angelica Albertini for bringing this error to our attention. Update 2: corrected a missing factor of $m$ in Eq. (3.30)

Journal ref: Phys. Rev. D 103, 104014 (2021)

Rapid generation of fully relativistic extreme-mass-ratio-inspiral waveform templates for LISA data analysis

Authors: Alvin J. K. Chua, Michael L. Katz, Niels Warburton, Scott A. Hughes

Abstract: The future space mission LISA will observe a wealth of gravitational-wave sources at millihertz frequencies. Of these, the extreme-mass-ratio inspirals of compact objects into massive black holes are the only sources that combine the challenges of strong-field complexity with that of long-lived signals. Such signals are found and characterized by comparing them against a large number of accurate w… ▽ More The future space mission LISA will observe a wealth of gravitational-wave sources at millihertz frequencies. Of these, the extreme-mass-ratio inspirals of compact objects into massive black holes are the only sources that combine the challenges of strong-field complexity with that of long-lived signals. Such signals are found and characterized by comparing them against a large number of accurate waveform templates during data analysis, but the rapid generation of such templates is hindered by computing the $\sim10^3$-$10^5$ harmonic modes in a fully relativistic waveform. We use order-reduction and deep-learning techniques to derive a global fit for these modes, and implement it in a complete waveform framework with hardware acceleration. Our high-fidelity waveforms can be generated in under $1\,\mathrm{s}$, and achieve a mismatch of $\lesssim 5\times 10^{-4}$ against reference waveforms that take $\gtrsim 10^4$ times longer. This marks the first time that analysis-length waveforms with full harmonic content can be produced on timescales useful for direct implementation in LISA analysis algorithms. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 051102 (2021)

Dissipation in extreme-mass ratio binaries with a spinning secondary

Authors: Sarp Akcay, Sam R. Dolan, Chris Kavanagh, Jordan Moxon, Niels Warburton, Barry Wardell

Abstract: We present the gravitational-wave flux balance law in an extreme mass-ratio binary with a spinning secondary. This law relates the flux of energy (angular momentum) radiated to null infinity and through the event horizon to the local change in the secondary's orbital energy (angular momentum) for generic (non-resonant) bound orbits in Kerr spacetime. As an explicit example we compute these quantit… ▽ More We present the gravitational-wave flux balance law in an extreme mass-ratio binary with a spinning secondary. This law relates the flux of energy (angular momentum) radiated to null infinity and through the event horizon to the local change in the secondary's orbital energy (angular momentum) for generic (non-resonant) bound orbits in Kerr spacetime. As an explicit example we compute these quantities for a spin-aligned body moving on a circular orbit around a Schwarzschild black hole. We perform this calculation both analytically, via a high-order post-Newtonian expansion, and numerically in two different gauges. Using these results we demonstrate explicitly that our new balance law holds. △ Less

Submitted 25 February, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: Minor typos corrected

Journal ref: Phys. Rev. D 102, 064013 (2020)

The location of the last stable orbit in Kerr spacetime

Authors: Leo C. Stein, Niels Warburton

Abstract: 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 st… ▽ More 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. △ Less

Submitted 5 March, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: 11+4 pages, 6 figures, 1 ancillary Mathematica file. v2: Minor edits, matched version published in PRD. Code available at https://bhptoolkit.org/

Journal ref: Phys. Rev. D 101, 064007 (2020)

Unveiling the Gravitational Universe at μ-Hz Frequencies

Authors: Alberto Sesana, Natalia Korsakova, Manuel Arca Sedda, Vishal Baibhav, Enrico Barausse, Simon Barke, Emanuele Berti, Matteo Bonetti, Pedro R. Capelo, Chiara Caprini, Juan Garcia-Bellido, Zoltan Haiman, Karan Jani, Oliver Jennrich, Peter Johansson, Fazeel Mahmood Khan, Valeriya Korol, Astrid Lamberts, Alessandro Lupi, Alberto Mangiagli, Lucio Mayer, Germano Nardini, Fabio Pacucci, Antoine Petiteau, Alvise Raccanelli , et al. (9 additional authors not shown)

Abstract: We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $μ$-Hz frequency range. By the 2040s', the $μ$-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 astro… ▽ More We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $μ$-Hz frequency range. By the 2040s', the $μ$-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 $μ$-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. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: 28 pages, 8 figures; White Paper submitted to ESA's Voyage 2050 call for papers on behalf of the LISA Consortium 2050 Task Force

Probing the Nature of Black Holes: Deep in the mHz Gravitational-Wave Sky

Authors: Vishal Baibhav, Leor Barack, Emanuele Berti, Béatrice Bonga, Richard Brito, Vitor Cardoso, Geoffrey Compère, Saurya Das, Daniela Doneva, Juan Garcia-Bellido, Lavinia Heisenberg, Scott A. Hughes, Maximiliano Isi, Karan Jani, Chris Kavanagh, Georgios Lukes-Gerakopoulos, Guido Mueller, Paolo Pani, Antoine Petiteau, Surjeet Rajendran, Thomas P. Sotiriou, Nikolaos Stergioulas, Alasdair Taylor, Elias Vagenas, Maarten van de Meent , et al. (4 additional authors not shown)

Abstract: 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 tel… ▽ More 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. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: 24 pages, 4 figures. White paper submitted to ESA's Voyage 2050 on behalf of the LISA Consortium 2050 Task Force

Journal ref: Exp Astron (2021)

The Missing Link in Gravitational-Wave Astronomy: Discoveries waiting in the decihertz range

Authors: Manuel Arca Sedda, Christopher P. L. Berry, Karan Jani, Pau Amaro-Seoane, Pierre Auclair, Jonathon Baird, Tessa Baker, Emanuele Berti, Katelyn Breivik, Adam Burrows, Chiara Caprini, Xian Chen, Daniela Doneva, Jose M. Ezquiaga, K. E. Saavik Ford, Michael L. Katz, Shimon Kolkowitz, Barry McKernan, Guido Mueller, Germano Nardini, Igor Pikovski, Surjeet Rajendran, Alberto Sesana, Lijing Shao, Nicola Tamanini , et al. (5 additional authors not shown)

Abstract: 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. Betwe… ▽ More 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. △ Less

Submitted 27 July, 2020; v1 submitted 29 August, 2019; originally announced August 2019.

Comments: 52 pages, 5 figures, 4 tables. Submitted to Classical & Quantum Gravity. Based upon a white paper for ESA's Voyage 2050 on behalf of the LISA Consortium 2050 Task Force

Second-order self-force calculation of the gravitational binding energy in compact binaries

Authors: Adam Pound, Barry Wardell, Niels Warburton, Jeremy Miller

Abstract: 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-fo… ▽ More 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. △ Less

Submitted 19 November, 2019; v1 submitted 20 August, 2019; originally announced August 2019.

Comments: 6 pages + supplemental material, 3 figures, revised title and two added references

Journal ref: Phys. Rev. Lett. 124, 021101 (2020)

Factorization and resummation: A new paradigm to improve gravitational wave amplitudes. III: the spinning test-body terms

Authors: Alessandro Nagar, Francesco Messina, Chris Kavanagh, Georgios Lukes-Gerakopoulos, Niels Warburton, Sebastiano Bernuzzi, Enno Harms

Abstract: We present new calculations of the energy flux of a spinning test-body on circular orbits around a Schwarzschild black hole at linear order in the particle spin. We compute the multipolar fluxes up to $\ell=m=6$ using two independent numerical solvers of theTeukolsky equation, one in the time domain and the other in the frequency domain. After linearization in the spin of the particle, we obtain a… ▽ More We present new calculations of the energy flux of a spinning test-body on circular orbits around a Schwarzschild black hole at linear order in the particle spin. We compute the multipolar fluxes up to $\ell=m=6$ using two independent numerical solvers of theTeukolsky equation, one in the time domain and the other in the frequency domain. After linearization in the spin of the particle, we obtain an excellent agreement ($\sim 10^{-5}$) between the two numerical results.The calculation of the multipolar fluxes is also performed analytically (up to $\ell=7$) using the post-Newtonian (PN) expansion of the Teukolsky equation solution; each mode is obtained at 5.5PN order beyond the corresponding leading-order contribution. From the analytical fluxes we obtain the PN-expanded analytical waveform amplitudes. These quantities are then resummed using new procedures either based on the factorization of the orbital contribution (and resumming it independently from the spin-dependent factor) or on the factorization of the tail contribution solely for odd-parity multipoles. We compare these prescriptions and the resummation procedure proposed in Pan et al. [Phys. Rev. D 83 (2011) 064003] to the numerical data. We find that the new procedures significantly improve over the existing one that, notably, is inconsistent with the numerical data for $\ell+m=\text{odd}$ multipoles already at low orbital frequencies. Our study suggests that the approach to waveform resummation used in current effective-one-body-based waveform models should be modified to improve its robustness and accuracy all over the binary parameter space. △ Less

Submitted 30 July, 2019; v1 submitted 29 July, 2019; originally announced July 2019.

Comments: 33 pages, 7 figures. Submitted to Phys. Rev. D

Journal ref: Phys. Rev. D 100, 104056 (2019)

Gravitational waves from bodies orbiting the Galactic Center black hole and their detectability by LISA

Authors: Eric Gourgoulhon, Alexandre Le Tiec, Frederic H. Vincent, Niels Warburton

Abstract: 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 induc… ▽ More 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. △ Less

Submitted 9 June, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Comments: 27 pages, 22 figures, accepted for publication in Astronomy and Astrophysics; text slightly modified and figures added/changed with respect to v2; see related paper by Amaro-Seoane 2019, arXiv:1903.10871

Journal ref: A&A 627, A92 (2019)

High-order asymptotics for the Spin-Weighted Spheroidal Equation at large real frequency

Authors: Marc Casals, Adrian C. Ottewill, Niels Warburton

Abstract: The spin-weighted spheroidal eigenvalues and eigenfunctions arise in the separation by variables of spin-field perturbations of Kerr black holes. We derive a large, real-frequency asymptotic expansion of the spin-weighted spheroidal eigenvalues and eigenfunctions to high order. This expansion corrects and extends existing results in the literature and we validate it via a numerical calculation. The spin-weighted spheroidal eigenvalues and eigenfunctions arise in the separation by variables of spin-field perturbations of Kerr black holes. We derive a large, real-frequency asymptotic expansion of the spin-weighted spheroidal eigenvalues and eigenfunctions to high order. This expansion corrects and extends existing results in the literature and we validate it via a numerical calculation. △ Less

Submitted 6 February, 2019; v1 submitted 30 September, 2018; originally announced October 2018.

Comments: 11 pages, 2 figures

Journal ref: Proc. R. Soc. A 475:20180701 (2019)

Fast Self-forced Inspirals

Authors: Maarten van de Meent, Niels Warburton

Abstract: We present a new, fast method for computing the inspiral trajectory and gravitational waves from extreme mass-ratio inspirals that can incorporate all known (and future) self-force results. Using near-identity (averaging) transformations we formulate equations of motion that do not explicitly depend upon the orbital phases of the inspiral, making them fast to evaluate, and whose solutions track th… ▽ More We present a new, fast method for computing the inspiral trajectory and gravitational waves from extreme mass-ratio inspirals that can incorporate all known (and future) self-force results. Using near-identity (averaging) transformations we formulate equations of motion that do not explicitly depend upon the orbital phases of the inspiral, making them fast to evaluate, and whose solutions track the evolving constants of motion, orbital phases and waveform phase of a full self-force inspiral to $O(η)$, where $η$ is the (small) mass ratio. As a concrete example, we implement these equations for inspirals of non-spinning (Schwarzschild) binaries. Our code computes inspiral trajectories in milliseconds which is a speed up of 2-5 orders of magnitude (depending on the mass-ratio) over previous self-force inspiral models which take minutes to hours to evaluate. Computing two-year duration waveforms using our new model we find a mismatch better than $\sim 10^{-4}$ with respect to waveforms computed using the (slower) full self-force models. The speed of our new approach is comparable with kludge models but has the added benefit of easily incorporating self-force results which will, once known, allow the waveform phase to be tracked to sub-radian accuracy over an inspiral. △ Less

Submitted 26 July, 2018; v1 submitted 14 February, 2018; originally announced February 2018.

Comments: 33 pages, code available at http://bhptoolkit.org/

Journal ref: Class. Quant. Grav. 35:144003 (2018)

Accelerated motion and the self-force in Schwarzschild spacetime

Authors: Anna Heffernan, Adrian C. Ottewill, Niels Warburton, Barry Wardell, Peter Diener

Abstract: We provide expansions of the Detweiler-Whiting singular field for motion along arbitrary, planar accelerated trajectories in Schwarzschild spacetime. We transcribe these results into mode-sum regularization parameters, computing previously unknown terms that increase the convergence rate of the mode-sum. We test our results by computing the self-force along a variety of accelerated trajectories. F… ▽ More We provide expansions of the Detweiler-Whiting singular field for motion along arbitrary, planar accelerated trajectories in Schwarzschild spacetime. We transcribe these results into mode-sum regularization parameters, computing previously unknown terms that increase the convergence rate of the mode-sum. We test our results by computing the self-force along a variety of accelerated trajectories. For non-uniformly accelerated circular orbits we present results from a new 1+1D discontinuous Galerkin time-domain code which employs an effective-source. We also present results for uniformly accelerated circular orbits and accelerated bound eccentric orbits computed within a frequency-domain treatment. Our regularization results will be useful for computing self-consistent self-force inspirals where the particle's worldline is accelerated with respect to the background spacetime. △ Less

Submitted 17 July, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: 19 pages, 6 figures (accepted CQG special issue article version)

Journal ref: Class. Quantum Grav. 35, 194001 (2018)

Evolution of small-mass-ratio binaries with a spinning secondary

Authors: Niels Warburton, Thomas Osburn, Charles R. Evans

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2017; originally announced August 2017.

Comments: 17 pages, 6 figures

Journal ref: Phys. Rev. D 96, 084057 (2017)

Inspiral into Gargantua

Authors: Samuel E. Gralla, Scott A. Hughes, Niels Warburton

Abstract: 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 con… ▽ More 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. △ Less

Submitted 19 September, 2016; v1 submitted 3 March, 2016; originally announced March 2016.

Comments: 10 pages, 5 figures, updated to journal version

Journal ref: Class. Quant. Grav. 33 155002 (2016)

Highly eccentric inspirals into a black hole

Authors: Thomas Osburn, Niels Warburton, Charles R. Evans

Abstract: We model the inspiral of a compact stellar-mass object into a massive nonrotating black hole including all dissipative and conservative first-order-in-the-mass-ratio effects on the orbital motion. The techniques we develop allow inspirals with initial eccentricities as high as $e\sim0.8$ and initial separations as large as $p\sim 50$ to be evolved through many thousands of orbits up to the onset o… ▽ More We model the inspiral of a compact stellar-mass object into a massive nonrotating black hole including all dissipative and conservative first-order-in-the-mass-ratio effects on the orbital motion. The techniques we develop allow inspirals with initial eccentricities as high as $e\sim0.8$ and initial separations as large as $p\sim 50$ to be evolved through many thousands of orbits up to the onset of the plunge into the black hole. The inspiral is computed using an osculating elements scheme driven by a hybridized self-force model, which combines Lorenz-gauge self-force results with highly accurate flux data from a Regge-Wheeler-Zerilli code. The high accuracy of our hybrid self-force model allows the orbital phase of the inspirals to be tracked to within $\sim0.1$ radians or better. The difference between self-force models and inspirals computed in the radiative approximation is quantified. △ Less

Submitted 20 March, 2016; v1 submitted 4 November, 2015; originally announced November 2015.

Comments: Updated to reflect published version

Journal ref: Phys. Rev. D 93, 064024 (2016)

Particle on the Innermost Stable Circular Orbit of a Rapidly Spinning Black Hole

Authors: Samuel E. Gralla, Achilleas P. Porfyriadis, Niels Warburton

Abstract: 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 sho… ▽ More 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. △ Less

Submitted 23 September, 2015; v1 submitted 28 June, 2015; originally announced June 2015.

Comments: 15 pages, 4 figures, v2: reference added, minor changes, matches published version

Journal ref: Phys. Rev. D 92, 064029 (2015)

Applying the effective-source approach to frequency-domain self-force calculations: Lorenz-gauge gravitational perturbations

Authors: Barry Wardell, Niels Warburton

Abstract: 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 replac… ▽ More 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. △ Less

Submitted 11 June, 2021; v1 submitted 28 May, 2015; originally announced May 2015.

Comments: Several typos corrected, including some that were identified in the published version

Journal ref: Phys. Rev. D 92, 084019 (2015)

Octupolar invariants for compact binaries on quasi-circular orbits

Authors: Patrick Nolan, Chris Kavanagh, Sam R Dolan, Adrian C Ottewill, Niels Warburton, Barry Wardell

Abstract: We extend the gravitational self-force methodology to identify and compute new $O(μ)$ tidal invariants for a compact body of mass $μ$ on a quasi-circular orbit about a black hole of mass $M \gg μ$. In the octupolar sector we find seven new degrees of freedom, made up of 3+3 conservative/dissipative `electric' invariants and 3+1 `magnetic' invariants, satisfying 1+1 and 1+0 trace conditions. After… ▽ More We extend the gravitational self-force methodology to identify and compute new $O(μ)$ tidal invariants for a compact body of mass $μ$ on a quasi-circular orbit about a black hole of mass $M \gg μ$. In the octupolar sector we find seven new degrees of freedom, made up of 3+3 conservative/dissipative `electric' invariants and 3+1 `magnetic' invariants, satisfying 1+1 and 1+0 trace conditions. After formulating for equatorial circular orbits on Kerr spacetime, we calculate explicitly for Schwarzschild spacetime. We employ both Lorenz gauge and Regge-Wheeler gauge numerical codes, and the functional series method of Mano, Suzuki and Takasugi. We present (i) highly-accurate numerical data and (ii) high-order analytical post-Newtonian expansions. We demonstrate consistency between numerical and analytic results, and prior work. We explore the application of these invariants in effective one-body models, and binary black hole initial-data formulations, and conclude with a discussion of future work. △ Less

Submitted 17 May, 2015; originally announced May 2015.

Comments: 24 pages, 3 figures, 4 tables

Journal ref: Phys. Rev. D 92, 123008 (2015)

Comparison Between Self-Force and Post-Newtonian Dynamics: Beyond Circular Orbits

Authors: Sarp Akcay, Alexandre Le Tiec, Leor Barack, Norichika Sago, Niels Warburton

Abstract: The gravitational self-force (GSF) and post-Newtonian (PN) schemes are complementary approximation methods for modelling the dynamics of compact binary systems. Comparison of their results in an overlapping domain of validity provides a crucial test for both methods, and can be used to enhance their accuracy, e.g. via the determination of previously unknown PN parameters. Here, for the first time,… ▽ More The gravitational self-force (GSF) and post-Newtonian (PN) schemes are complementary approximation methods for modelling the dynamics of compact binary systems. Comparison of their results in an overlapping domain of validity provides a crucial test for both methods, and can be used to enhance their accuracy, e.g. via the determination of previously unknown PN parameters. Here, for the first time, we extend such comparisons to noncircular orbits---specifically, to a system of two nonspinning objects in a bound (eccentric) orbit. To enable the comparison we use a certain orbital-averaged quantity $\langle U \rangle$ that generalizes Detweiler's redshift invariant. The functional relationship $\langle U \rangle(Ω_r,Ω_φ)$, where $Ω_r$ and $Ω_φ$ are the frequencies of the radial and azimuthal motions, is an invariant characteristic of the conservative dynamics. We compute $\langle U \rangle(Ω_r,Ω_φ)$ numerically through linear order in the mass ratio $q$, using a GSF code which is based on a frequency-domain treatment of the linearized Einstein equations in the Lorenz gauge. We also derive $\langle U \rangle(Ω_r,Ω_φ)$ analytically through 3PN order, for an arbitrary $q$, using the known near-zone 3PN metric and the generalized quasi-Keplerian representation of the motion. We demonstrate that the $\mathcal{O}(q)$ piece of the analytical PN prediction is perfectly consistent with the numerical GSF results, and we use the latter to estimate yet unknown pieces of the 4PN expression at $\mathcal{O}(q)$. △ Less

Submitted 10 March, 2015; v1 submitted 4 March, 2015; originally announced March 2015.

Comments: 44 pages, 2 figures, 4 tables

Journal ref: Phys. Rev. D 91, 124014 (2015)

Self force on a scalar charge in Kerr spacetime: inclined circular orbits

Authors: Niels Warburton

Abstract: 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… ▽ More 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. △ Less

Submitted 2 February, 2015; v1 submitted 12 August, 2014; originally announced August 2014.

Comments: 17 pages, 5 figures. Updated to reflect published version

Journal ref: Phys. Rev. D 91, 024045 (2015)

Tidal invariants for compact binaries on quasi-circular orbits

Authors: Sam R. Dolan, Patrick Nolan, Adrian C. Ottewill, Niels Warburton, Barry Wardell

Abstract: We extend the gravitational self-force approach to encompass `self-interaction' tidal effects for a compact body of mass $μ$ on a quasi-circular orbit around a black hole of mass $M \gg μ$. Specifically, we define and calculate at $O(μ)$ (conservative) shifts in the eigenvalues of the electric- and magnetic-type tidal tensors, and a (dissipative) shift in a scalar product between their eigenbases.… ▽ More We extend the gravitational self-force approach to encompass `self-interaction' tidal effects for a compact body of mass $μ$ on a quasi-circular orbit around a black hole of mass $M \gg μ$. Specifically, we define and calculate at $O(μ)$ (conservative) shifts in the eigenvalues of the electric- and magnetic-type tidal tensors, and a (dissipative) shift in a scalar product between their eigenbases. This approach yields four gauge-invariant functions, from which one may construct other tidal quantities such as the curvature scalars and the speciality index. First, we analyze the general case of a geodesic in a regular perturbed vacuum spacetime admitting a helical Killing vector and a reflection symmetry. Next, we specialize to focus on circular orbits in the equatorial plane of Kerr spacetime at $O(μ)$. We present accurate numerical results for the Schwarzschild case for orbital radii up to the light-ring, calculated via independent implementations in Lorenz and Regge-Wheeler gauges. We show that our results are consistent with leading-order post-Newtonian expansions, and demonstrate the existence of additional structure in the strong-field regime. We anticipate that our strong-field results will inform (e.g.) effective one-body models for the gravitational two-body problem that are invaluable in the ongoing search for gravitational waves. △ Less

Submitted 25 January, 2015; v1 submitted 18 June, 2014; originally announced June 2014.

Comments: 29 pages, 5 figures, 3 tables. Corrected data in Table I (cf arXiv:1409.6933) to match published version

Journal ref: Phys. Rev. D 91, 023009 (2015)

Gravitational Self-Force Correction to the Innermost Stable Circular Equatorial Orbit of a Kerr Black Hole

Authors: Soichiro Isoyama, Leor Barack, Sam R. Dolan, Alexandre Le Tiec, Hiroyuki Nakano, Abhay G. Shah, Takahiro Tanaka, Niels Warburton

Abstract: For a self-gravitating particle of mass μin orbit around a Kerr black hole of mass M >> μ, we compute the O(μ/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-def… ▽ More For a self-gravitating particle of mass μin orbit around a Kerr black hole of mass M >> μ, we compute the O(μ/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. △ Less

Submitted 16 October, 2014; v1 submitted 24 April, 2014; originally announced April 2014.

Comments: 5 pages, 1 table, 1 figure, matches version published in PRL. Raw data of H_int/mu are available at http://link.aps.org/supplemental/10.1103/PhysRevLett.113.161101

Report number: KUNS-2490,YITP-14-33

Journal ref: Phys. Rev. Lett. 113, 161101 (2014)

Gravitational self-torque and spin precession in compact binaries

Authors: Sam R. Dolan, Niels Warburton, Abraham I. Harte, Alexandre Le Tiec, Barry Wardell, Leor Barack

Abstract: We calculate the effect of self-interaction on the "geodetic" spin precession of a compact body in a strong-field orbit around a black hole. Specifically, we consider the spin precession angle $ψ$ per radian of orbital revolution for a particle carrying mass $μ$ and spin $s \ll (G/c) μ^2$ in a circular orbit around a Schwarzschild black hole of mass $M \gg μ$. We compute $ψ$ through $O(μ/M)$ in pe… ▽ More We calculate the effect of self-interaction on the "geodetic" spin precession of a compact body in a strong-field orbit around a black hole. Specifically, we consider the spin precession angle $ψ$ per radian of orbital revolution for a particle carrying mass $μ$ and spin $s \ll (G/c) μ^2$ in a circular orbit around a Schwarzschild black hole of mass $M \gg μ$. We compute $ψ$ through $O(μ/M)$ in perturbation theory, i.e, including the correction $δψ$ (obtained numerically) due to the torque exerted by the conservative piece of the gravitational self-field. Comparison with a post-Newtonian (PN) expression for $δψ$, derived here through 3PN order, shows good agreement but also reveals strong-field features which are not captured by the latter approximation. Our results can inform semi-analytical models of the strong-field dynamics in astrophysical binaries, important for ongoing and future gravitational-wave searches. △ Less

Submitted 7 March, 2014; v1 submitted 3 December, 2013; originally announced December 2013.

Comments: 5 pages, 1 table, 1 figure. Minor changes to match published version

Journal ref: Phys. Rev. D 89, 064011 (2014)

Applying the effective-source approach to frequency-domain self-force calculations

Authors: Niels Warburton, Barry Wardell

Abstract: 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 effecti… ▽ More 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. △ Less

Submitted 26 October, 2014; v1 submitted 13 November, 2013; originally announced November 2013.

Comments: 13 pages, 1 figure. Correct minor typo and update references

Journal ref: Phys. Rev. D 89, 044046 (2014)

Frequency-domain algorithm for the Lorenz-gauge gravitational self-force

Authors: Sarp Akcay, Niels Warburton, Leor Barack

Abstract: State-of-the-art computations of the gravitational self-force (GSF) on massive particles in black hole spacetimes involve numerical evolution of the metric perturbation equations in the time-domain, which is computationally very costly. We present here a new strategy, based on a frequency-domain treatment of the perturbation equations, which offers considerable computational saving. The essential… ▽ More State-of-the-art computations of the gravitational self-force (GSF) on massive particles in black hole spacetimes involve numerical evolution of the metric perturbation equations in the time-domain, which is computationally very costly. We present here a new strategy, based on a frequency-domain treatment of the perturbation equations, which offers considerable computational saving. The essential ingredients of our method are (i) a Fourier-harmonic decomposition of the Lorenz-gauge metric perturbation equations and a numerical solution of the resulting coupled set of ordinary equations with suitable boundary conditions; (ii) a generalized version of the method of extended homogeneous solutions [Phys. Rev. D {\bf 78}, 084021 (2008)] used to circumvent the Gibbs phenomenon that would otherwise hamper the convergence of the Fourier mode-sum at the particle's location; and (iii) standard mode-sum regularization, which finally yields the physical GSF as a sum over regularized modal contributions. We present a working code that implements this strategy to calculate the Lorenz-gauge GSF along eccentric geodesic orbits around a Schwarzschild black hole. The code is far more efficient than existing time-domain methods; the gain in computation speed (at a given precision) is about an order of magnitude at an eccentricity of 0.2, and up to three orders of magnitude for circular or nearly circular orbits. This increased efficiency was crucial in enabling the recently reported calculation of the long-term orbital evolution of an extreme mass ratio inspiral [Phys. Rev. D {\bf 85}, 061501(R) (2012)]. Here we provide full technical details of our method to complement the above report. △ Less

Submitted 4 November, 2013; v1 submitted 23 August, 2013; originally announced August 2013.

Comments: 27 pages, 4 figures

Journal ref: Phys. Rev. D88:104009 (2013)

Isofrequency pairing of geodesic orbits in Kerr geometry

Authors: Niels Warburton, Leor Barack, Norichika Sago

Abstract: 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 i… ▽ More 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. △ Less

Submitted 5 April, 2013; v1 submitted 16 January, 2013; originally announced January 2013.

Comments: 11 pages, 8 figures. Updated to reflect published version

Journal ref: Phys. Rev. D 87, 084012 (2013)

Evolution of inspiral orbits around a Schwarzschild black hole

Authors: Niels Warburton, Sarp Akcay, Leor Barack, Jonathan R. Gair, Norichika Sago

Abstract: 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 o… ▽ More 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. △ Less

Submitted 10 December, 2013; v1 submitted 29 November, 2011; originally announced November 2011.

Comments: 4 pages, 3 figures, Fast GSF calculator code available at: http://www.nielswarburton.net/lib_Sch_GSF/

Journal ref: Phys. Rev. D 85, 061501(R) (2012)