Search

We investigate the late-time tail behavior in gravitational waves from merging eccentric binary black holes (BBH) using black hole perturbation theory. For simplicity, we focus only on the dominant quadrupolar mode of the radiation. We demonstrate that such tails become more prominent as eccentricity increases. Exploring the phenomenology of the tails in both spinning and non-spinning eccentric binaries, with the spin magnitude varying from $\chi=-0.6$ to $\chi=+0.6$ and eccentricity as high as $e=0.98$, we find that these tails can be well approximated by a slowly decaying power law. We study the power law for varying systems and find that the power law exponent lies close to the theoretically expected value $-4$. Finally, using both plunge geodesic and radiation-reaction-driven orbits, we perform a series of numerical experiments to understand the origin of the tails in BBH simulations. Our results suggest that the late-time tails are strongly excited in eccentric BBH systems when the smaller black hole is in the neighborhood of the apocenter, as opposed to any structure in the strong field of the larger black hole. Our analysis framework is publicly available through the \texttt{gwtails} Python package., Comment: 16 pages, 18 figures

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.

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, ($\epsilon \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}(\epsilon^{-1/2})$ in the fully averaged scheme to $\mathcal{O}(\epsilon^{4/7})$ in the switching scheme. At the same time, this scheme improves the scaling of the computation time from being inversely proportional to $\epsilon$ using OG, to a very weak scaling with $\epsilon$. This results in a speed-up of at least two orders of magnitude for LISA EMRIs with room for further optimisation., Comment: 50 pages, 11 figures

In this work, we test an effective-one-body radiation-reaction force for eccentric planar orbits of a test mass in a Kerr background, which contains third-order post-Newtonian (PN) non-spinning and second-order PN spin contributions. We compare the analytical fluxes connected to two different resummations of this force, truncated at different PN orders in the eccentric sector, with the numerical fluxes computed through the use of frequency- and time-domain Teukolsky-equation codes. We find that the different PN truncations of the radiation-reaction force show the expected scaling in the weak gravitational-field regime, and we observe a fractional difference with the numerical fluxes that is $<5 \%$, for orbits characterized by eccentricity $0 \le e \le 0.7$, central black-hole spin $-0.99 M \le a \le 0.99 M$ and fixed orbital-averaged quantity $x=\langle M\Omega \rangle^{2/3} = 0.06$, corresponding to the mildly strong-field regime with semilatera recta $9 M

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

Astrophysical black holes do not exist in vacuum, and their motion is affected by the galactic environment. As a black hole moves it attracts stars and matter, creating a wake that, in turn, exerts an effective friction slowing down the black hole. This force is known as dynamical friction, and has significant consequences, ranging from the formation of supermassive black hole binaries to modifications in the phase of binary mergers. In this work we explore the motion of spinning black holes on a medium. We find that the classical drag along the velocity direction is modified and two novel forces appear: a rotational force, which in the context of fluid dynamics is dubbed the Magnus force, and a lift, orthogonal to the direction of motion. We develop a first-principles fully-relativistic treatment of these spin-induced aerodynamic forces in two types of environment: i) collisionless corpuscular matter and ii) a light scalar field, exploring the differences between both cases. In both cases we find that the total rotational force acts precisely in the opposite direction as compared to the classical set-up of a spinning ball moving through a fluid. Finally, we comment on the consequences of these new effects for astrophysics and gravitational wave observations., Comment: 15 pages, 8 figures, minor update to text and references

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

We 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., Comment: 25 pages, 9 figures, published in Physical Review D

The two-body problem is extensively studied in open systems and asymptotically flat spacetimes. However, there are many systems where radiation is trapped: they range from radiating charges in cavities to low-energy excitations of massive degrees of freedom, to anti-de Sitter spacetimes. Here, we study the problem of motion of a pointlike particle orbiting a massive compact object inside a cavity. We first show that - assuming circular motion - there are initial conditions for which the self-force vanishes and the binary is eternal. We then consider the evolution of the system under radiation reaction in a toy model which we argue captures the essentials of orbiting particles. We show that eternal circular binaries may exist. We also show that the presence of cavity modes leads to chaos in regimes of strong coupling or when the system is initialized close enough to a resonance. Our results have implications for physics in anti-de Sitter spacetimes and possibly for binaries evolving within dark matter haloes, if it consists on massive fundamental fields., Comment: 21 pages, 17 figures

We present SEOBNRv5HM, a more accurate and faster inspiral-merger-ringdown gravitational waveform model for quasi-circular, spinning, nonprecessing binary black holes within the effective-one-body (EOB) formalism. Compared to its predecessor, SEOBNRv4HM, the waveform model i) incorporates recent high-order post- Newtonian results in the inspiral, with improved resummations, ii) includes the gravitational modes (l, |m|) = (3, 2), (4, 3), in addition to the (2, 2), (3, 3), (2, 1), (4, 4), (5, 5) modes already implemented in SEOBNRv4HM, iii) is calibrated to larger mass-ratios and spins using a catalog of 442 numerical-relativity (NR) simulations and 13 additional waveforms from black-hole perturbation theory, iv) incorporates information from second-order gravitational self-force (2GSF) in the nonspinning modes and radiation-reaction force. Computing the unfaithfulness against NR simulations, we find that for the dominant (2, 2) mode the maximum unfaithfulness in the total mass range $10-300 M_{\odot}$ is below $10^{-3}$ for 90% of the cases (38% for SEOBNRv4HM). When including all modes up to l = 5 we find 98% (49%) of the cases with unfaithfulness below $10^{-2} (10^{-3})$, while these numbers reduce to 88% (5%) when using SEOBNRv4HM. Furthermore, the model shows improved agreement with NR in other dynamical quantities (e.g., the angular momentum flux and binding energy), providing a powerful check of its physical robustness. We implemented the waveform model in a high-performance Python package (pySEOBNR), which leads to evaluation times faster than SEOBNRv4HM by a factor 10 to 50, depending on the configuration, and provides the flexibility to easily include spin-precession and eccentric effects, thus making it the starting point for a new generation of EOBNR waveform models (SEOBNRv5) to be employed for upcoming observing runs of the LIGO-Virgo-KAGRA detectors.

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.

Eccentric compact binary mergers are significant scientific targets for current and future gravitational wave observatories. To detect and analyze eccentric signals, there is an increasing effort to develop waveform models, numerical relativity simulations, and parameter estimation frameworks for eccentric binaries. Unfortunately, current models and simulations use different internal parameterisations of eccentricity in the absence of a unique natural definition of eccentricity in general relativity, which can result in incompatible eccentricity measurements. In this paper, we adopt a standardized definition of eccentricity and mean anomaly based solely on waveform quantities, and make our implementation publicly available through an easy-to-use Python package, gw_eccentricity. This definition is free of gauge ambiguities, has the correct Newtonian limit, and can be applied as a postprocessing step when comparing eccentricity measurements from different models. This standardization puts all models and simulations on the same footing and enables direct comparisons between eccentricity estimates from gravitational wave observations and astrophysical predictions. We demonstrate the applicability of this definition and the robustness of our implementation for waveforms of different origins, including post-Newtonian theory, effective one body, extreme mass ratio inspirals, and numerical relativity simulations. We focus on binaries without spin-precession in this work, but possible generalizations to spin-precessing binaries are discussed., Comment: Python implementation available at https://pypi.org/project/gw-eccentricity/, matches published version

We present closed-form solutions for plunging geodesics in the extended Kerr spacetime using Boyer-Lindquist coordinates. Our solutions directly solve for the dynamics of generic timelike plunges, we also specialise to the case of test particles plunging from a precessing innermost stable circular orbit (ISSO). We find these solutions in the form of elementary and Jacobi elliptic functions parameterized by Mino time. In particular, we demonstrate that solutions for the ISSO case can be determined almost entirely in terms of elementary functions, depending only on the spin parameter of the black hole and the radius of the ISSO. This extends recent work on the case of equatorial plunges from the innermost stable circular orbit. Furthermore, we introduce a new equation that characterizes the radial inflow from the ISSO to the horizon, taking into account the inclination. For ease of application, our results have been implemented in the KerrGeodesics package in the Black Hole Perturbation Toolkit., Comment: 22 pages, 7 Figures, v3 updated minor typos in Eq.51 and Eq.53

The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a function of mass ratio for eccentric non-spinning binary black holes. We produce $52$ NR simulations with mass ratios between $1:10$ and $1:1$ and initial eccentricities up to $0.7$. From these we extract quantities like gravitational wave energy and angular momentum fluxes and periastron advance, and assess their accuracy. To facilitate comparison, we develop tools to map between NR and SMR inspiral evolutions of eccentric binary black holes. We derive post-Newtonian accurate relations between different definitions of eccentricity. Based on these analyses, we introduce a new definition of eccentricity based on the (2,2)-mode of the gravitational radiation, which reduces to the Newtonian definition of eccentricity in the Newtonian limit. From the comparison between NR simulations and SMR results, we quantify the unknown next-to-leading order SMR contributions to the gravitational energy and angular momentum fluxes, and periastron advance. We show that in the comparable mass regime these contributions are subdominant and higher order SMR contributions are negligible.

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., Comment: 46 pages, 12 figures

The next generation of ground-based gravitational-wave detectors will observe coalescences of black holes and neutron stars throughout the cosmos, thousands of them with exceptional fidelity. The Science Book is the result of a 3-year effort to study the science capabilities of networks of next generation detectors. Such networks would make it possible to address unsolved problems in numerous areas of physics and astronomy, from Cosmology to Beyond the Standard Model of particle physics, and how they could provide insights into workings of strongly gravitating systems, astrophysics of compact objects and the nature of dense matter. It is inevitable that observatories of such depth and finesse will make new discoveries inaccessible to other windows of observation. In addition to laying out the rich science potential of the next generation of detectors, this report provides specific science targets in five different areas in physics and astronomy and the sensitivity requirements to accomplish those science goals. This report is the second in a six part series of reports by the GWIC 3G Subcommittee: i) Expanding the Reach of Gravitational Wave Observatories to the Edge of the Universe, ii) The Next Generation Global Gravitational Wave Observatory: The Science Book (this report), iii) 3G R&D: R&D for the Next Generation of Ground-based Gravitational Wave Detectors, iv) Gravitational Wave Data Analysis: Computing Challenges in the 3G Era, v) Future Ground-based Gravitational-wave Observatories: Synergies with Other Scientific Communities, and vi) An Exploration of Possible Governance Models for the Future Global Gravitational-Wave Observatory Network., Comment: 69 pages, 18 figures

Extreme mass-ratio inspirals (EMRIs) detectable by the Laser Interferometer Space Antenna are unique probes of the nature of supermassive compact objects. We compute the gravitational-wave signal emitted by a stellar-mass compact object in a circular equatorial orbit around a Kerr-like horizonless supermassive object defined by an effective radius and a reflectivity coefficient. The Teukolsky equations are solved consistently with suitable (frequency-dependent) boundary conditions, and the modified energy and angular-momentum fluxes are used to evolve the orbital parameters adiabatically. The gravitational fluxes have resonances corresponding to the low-frequency quasinormal modes of the central object, which can contribute significantly to the gravitational-wave phase. Overall, the absence of a classical event horizon in the central object can affect the gravitational-wave signal dramatically, with deviations even larger than those previously estimated by a model-independent analysis of the tidal heating. We estimate that EMRIs could potentially place the most stringent constraint on the reflectivity of supermassive compact objects at the remarkable level of ${\cal O}(10^{-6})\%$ and would allow to constrain various models which are not ruled out by the ergoregion instability. In particular, an EMRI detection could allow one to rule out (or provide evidence for) signatures of quantum black-hole horizons with Boltzmann reflectivity. Our results provide motivation for performing rigorous parameter estimation to assess the detectability of these effects., Comment: 11+3 pages, 9 figures. v2: minor revisions to match the version in PRD

The dynamics of merging black hole binaries with small mass-ratios (SMR) can be expressed as a series in the symmetric mass-ratio, the post-adiabatic expansion. Using numerical relativity simulations, we compute the first three terms in this series for quasicircular non-spinning binaries. We recover the leading term predicted by SMR perturbation theory, and obtain a robust prediction of the 1-post-adiabatic term. The 2-post-adiabatic term is found to be small, indicating that the 1PA approximation may be quite accurate even for nearly comparable mass binaries. We also estimate the range of applicability for SMR and post-Newtonian series for non-spinning, quasi-circular inspirals., Comment: 6 pages, 4 figures, published version

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., Comment: 24 pages, 4 figures. White paper submitted to ESA's Voyage 2050 on behalf of the LISA Consortium 2050 Task Force

The self-force program aims at accurately modeling relativistic two-body systems with a small mass ratio (SMR). In the context of the effective-one-body (EOB) framework, current results from this program can be used to determine the effective metric components at linear order in the mass ratio, resumming post-Newtonian (PN) dynamics around the test-particle limit in the process. It was shown in [Akcay et al., Phys. Rev. D 86 (2012)] that, in the original (standard) EOB gauge, the SMR contribution to the metric component $g^\text{eff}_{tt}$ exhibits a coordinate singularity at the light-ring (LR) radius. In this paper, we adopt a different gauge for the EOB dynamics and obtain a Hamiltonian that is free of poles at the LR, with complete circular-orbit information at linear order in the mass ratio and non-circular-orbit and higher-order-in-mass-ratio terms up to 3PN order. We confirm the absence of the LR-divergence in such an EOB Hamiltonian via plunging trajectories through the LR radius. Moreover, we compare the binding energies and inspiral waveforms of EOB models with SMR, PN and mixed SMR-3PN information on a quasi-circular inspiral against numerical-relativity predictions. We find good agreement between NR simulations and EOB models with SMR-3PN information for both equal and unequal mass ratios. In particular, when compared to EOB inspiral waveforms with only 3PN information, EOB Hamiltonians with SMR-3PN information improves the modeling of binary systems with small mass ratios $q \lesssim 1/3$, with a dephasing accumulated in $\sim$30 gravitational-wave (GW) cycles being of the order of few hundredths of a radian up to 4 GW cycles before merger., Comment: 16 pages, 5 figures

If a small compact object orbits a black hole, it is known that it can excite the black hole's quasinormal modes (QNMs), leading to high-frequency oscillations ("wiggles") in the radiated field at $\mathcal{J}^+$, and in the radiation-reaction self-force acting on the object after its orbit passes through periapsis. Here we survey the phenomenology of these wiggles across a range of black hole spins and equatorial orbits. In both the scalar-field and gravitational cases we find that wiggles are a generic feature across a wide range of parameter space, and that they are observable in field perturbations at fixed spatial positions, in the self-force, and in radiated fields at $\mathcal{J}^+$. For a given charge or mass of the small body, the QNM excitations have the highest amplitudes for systems with a highly spinning central black hole, a prograde orbit with high eccentricity, and an orbital periapsis close to the light ring. The QNM amplitudes depend smoothly on the orbital parameters, with only very small amplitude changes when the orbit's (discrete) frequency spectrum is tuned to match QNM frequencies. The association of wiggles with QNM excitations suggest that they represent a situation where the \emph{nonlocal} nature of the self-force is particularly apparent, with the wiggles arising as result of QNM excitation by the compact object near periapsis, and then encountered later in the orbit. Astrophysically, the effects of wiggles at $\mathcal{J}^+$ might allow direct observation of Kerr QNMs in extreme-mass-ratio inspiral (EMRI) binary black hole systems, potentially enabling new tests of general relativity., Comment: 26 pages, 12 color pdf figures and 6 color png bitmap figures, full-resolution versions of png bitmap figures included as ancillary files in anc/ directory; revised in response to referee comments and to specifically note the occurence of wiggles for retrograde orbits

We provide analytical closed form solutions for the parallel transport along a bound geodesic in Kerr spacetime. This can be considered the lowest order approximation for the motion a spinning black hole in an extreme mass-ratio inspiral. As an illustration of the usefulness of our new found expressions we scope out the locations of spin-spin resonances in quasi-circular EMRIs. All solutions are given as functions of Mino time, which facilitates the decoupling of the equations of motion. To help physical interpretation, we also provide an analytical expression for the proper time along a geodesic as a function of Mino time.

Advanced methods for computing perturbative, quantum-gravitational scattering amplitudes show great promise for improving our knowledge of classical gravitational dynamics. This is especially true in the weak-field and arbitrary-speed (post-Minkowskian, PM) regime, where the conservative dynamics at 3PM order has been recently determined for the first time, via an amplitude calculation. Such PM results are most relevantly applicable to relativistic scattering (unbound orbits), while bound/inspiraling binary systems, the most frequent sources of gravitational waves for the LIGO and Virgo detectors, are most suitably modeled by the weak-field and slow-motion (post-Newtonian, PN) approximation. Nonetheless, it has been suggested that PM results can independently lead to improved modeling of bound binary dynamics, especially when taken as inputs for effective-one-body (EOB) models of inspiraling binaries. Here, we initiate a quantitative study of this possibility, by comparing PM, EOB and PN predictions for the binding energy of a two-body system on a quasi-circular inspiraling orbit against results of numerical relativity (NR) simulations. The binding energy is one of the two central ingredients (the other being the gravitational-wave energy flux) that enters the computation of gravitational waveforms employed by LIGO and Virgo detectors, and for (quasi-)circular orbits it provides an accurate diagnostic of the conservative sector of a model. We find that, whereas 3PM results do improve the agreement with NR with respect to 2PM (especially when used in the EOB framework), it is crucial to push PM calculations at higher orders if one wants to achieve better performances than current waveform models used for LIGO/Virgo data analysis., Comment: 13 pages, 7 figures. v2: published

We generalize to Kerr spacetime previous gravitational self-force results on gyroscope precession along circular orbits in the Schwarzschild spacetime. In particular we present high order post- Newtonian expansions for the gauge invariant precession function along circular geodesics valid for arbitrary Kerr spin parameter and show agreement between these results and those derived from the full post-Newtonian conservative dynamics. Finally we present strong field numerical data for a range of the Kerr spin parameter, showing agreement with the GSF-PN results, and the expected lightring divergent behaviour. These results provide useful testing benchmarks for self- force calculations in Kerr spacetime, and provide an avenue for translating self-force data into the spin-spin coupling in effective-one-body models., Comment: 11 pages, 3 figures

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

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(\eta)$, where $\eta$ 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., Comment: 33 pages, code available at http://bhptoolkit.org/

In this work we present the first calculation of the gravitational self-force on generic bound geodesics in Kerr spacetime to first order in the mass-ratio. That is, the local correction to equations of motion for a compact object orbiting a larger rotating black hole due to its own impact on the gravitational field. This includes both dissipative and conservative effects. Our method builds on and extends earlier methods for calculating the gravitational self-force on equatorial orbits. In particular we reconstruct the local metric perturbation in the outgoing radiation gauge from the Weyl scalar $\psi_4$, which in turn is obtained by solving the Teukolsky equation using semi-analytical frequency domain methods. The gravitational self-force is subsequently obtained using (spherical) $l$-mode regularization. We test our implementation by comparing the large $l$-behaviour against the analytically known regularization parameters. In addition we validate our results be comparing the long-term average changes to the energy, angular momentum, and Carter constant to changes to these constants of motion inferred from the gravitational wave flux to infinity and down the horizon.

We prove a key result regarding the mass and angular momentum content of linear vacuum perturbations of the Kerr metric obtained through the formalism developed by Chrzarnowski, Cohen, and Kegeles (CCK). More precisely we prove that the Abbott-Deser mass and angular momentum integrals of any such perturbation vanish, when that perturbation was obtained from a regular Fourier mode of the Hertz potential. As a corollary we obtain a generalization of previous results on the completion of the `no string' radiation gauge metric perturbation generated by a point particle. We find that for any bound orbit around a Kerr black hole, the mass and angular momentum perturbations completing the CCK metric are simply the energy and angular momentum of the particle "outside" the orbit and vanish "inside" the orbit., Comment: Version accepted be CQG

We employ a recently developed mode-sum regularization method to compute the renormalized stress-energy tensor of a quantum field in the Kerr background metric (describing a stationary spinning black hole). More specifically, we consider a minimally-coupled massless scalar field in the Unruh vacuum state, the quantum state corresponding to an evaporating black hole. The computation is done here for the case $a=0.7M$, using two different variants of the method: $t$-splitting and $\varphi$-splitting, yielding good agreement between the two (in the domain where both are applicable). We briefly discuss possible implications of the results for computing semiclassical corrections to certain quantities, and also for simulating dynamical evaporation of a spinning black hole., Comment: 5 pages, 4 figures

The linear in mass ratio correction to the periapsis advance of equatorial orbits around a spinning black hole is calculated for the first time and to very high precision, providing a key benchmark for different approaches modelling spinning binaries. The high precision of the calculation is leveraged to discriminate between two recent incompatible derivations of the 4PN equations of motion. Finally, the limit of the periapsis advance near the innermost stable orbit (ISCO) allows determination of the ISCO shift, validating previous calculations using the first law of binary mechanics. Calculation of the ISCO shift is further extended into the near extremal regime (with spins up to $1-a=10^{-20}$), revealing new unexpected phenomenology. In particular, we find that the shift of the ISCO does not have a well defined extremal limit, but instead continues to oscillate., Comment: 6 pages, 3 figures, Updated to match PRL accepted version

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

This paper presents the first calculation of the gravitational self-force on a small compact object on an eccentric equatorial orbit around a Kerr black hole to first order in the mass-ratio. That is the pointwise correction to the object's equations of motion (both conservative and dissipative) due to its own gravitational field treated as a linear perturbation to the background Kerr spacetime generated by the much larger spinning black hole. The calculation builds on recent advances on constructing the local metric and self-force from solutions of the Teukolsky equation, which led to the calculation of the Detweiler-Barack-Sago redshift invariant on eccentric equatorial orbits around a Kerr black hole in a previous paper. After deriving the necessary expression to obtain the self-force from the Weyl scalar $\psi_4$, we perform several consistency checks of the method and numerical implementation, including a check of the balance law relating orbital average of the self-force to average flux of energy and angular momentum out of the system. Particular attention is paid to the pointwise convergence properties of the sum over frequency modes in our method, identifying a systematic inherent loss of precision that any frequency domain calculation of the self-force on eccentric orbits must overcome., Comment: Various typos and correction to match version accepted for publication in PRD

The effective-one-body theory (EOB) describes the conservative dynamics of compact binary systems in terms of an effective Hamiltonian approach. The Hamiltonian for moderately eccentric motion of two non-spinning compact objects in the extreme mass-ratio limit is given in terms of three potentials: $a(v), \bar{d}(v), q(v)$. By generalizing the first law of mechanics for (non-spinning) black hole binaries to eccentric orbits, [\prd{\bf92}, 084021 (2015)] recently obtained new expressions for $\bar{d}(v)$ and $q(v)$ in terms of quantities that can be readily computed using the gravitational self-force approach. Using these expressions we present a new computation of the EOB potential $q(v)$ by combining results from two independent numerical self-force codes. We determine $q(v)$ for inverse binary separations in the range $1/1200 \le v \lesssim 1/6$. Our computation thus provides the first-ever strong-field results for $q(v)$. We also obtain $\bar{d}(v)$ in our entire domain to a fractional accuracy of $\gtrsim 10^{-8}$. We find to our results are compatible with the known post-Newtonian expansions for $\bar{d}(v)$ and $q(v)$ in the weak field, and agree with previous (less accurate) numerical results for $\bar{d}(v)$ in the strong field., Comment: 4 figures, numerical data at the end. Fixed the typos, added the journal reference

It is known that a near-extremal Kerr black hole can be spun up beyond its extremal limit by capturing a test particle. Here we show that overspinning is always averted once back-reaction from the particle's own gravity is properly taken into account. We focus on nonspinning, uncharged, massive particles thrown in along the equatorial plane, and work in the first-order self-force approximation (i.e., we include all relevant corrections to the particle's acceleration through linear order in the ratio, assumed small, between the particle's energy and the black hole's mass). Our calculation is a numerical implementation of a recent analysis by two of us [Phys.\ Rev.\ D {\bf 91}, 104024 (2015)], in which a necessary and sufficient "censorship" condition was formulated for the capture scenario, involving certain self-force quantities calculated on the one-parameter family of unstable circular geodesics in the extremal limit. The self-force information accounts both for radiative losses and for the finite-mass correction to the critical value of the impact parameter. Here we obtain the required self-force data, and present strong evidence to suggest that captured particles never drive the black hole beyond its extremal limit. We show, however, that, within our first-order self-force approximation, it is possible to reach the extremal limit with a suitable choice of initial orbital parameters. To rule out such a possibility would require (currently unavailable) information about higher-order self-force corrections., Comment: 13 pages, 3 figures

We present the first numerical calculation of the (local) metric perturbation produced by a small compact object moving on an eccentric equatorial geodesic around a Kerr black hole, accurate to first order in the mass ratio. The procedure starts by first solving the Teukolsky equation to obtain the Weyl scalar $\psi_4$ using semi-analytical methods. The metric perturbation is then reconstructed from $\psi_4$ in an (outgoing) radiation gauge, adding the appropriate non-radiative contributions arising from the shifts in mass and angular momentum of the spacetime. As a demonstration we calculate the generalized redshift $U$ as a function of the orbital frequencies $\Omega_r$ and $\Omega_\phi$ to linear order in the mass ratio, a gauge invariant measure of the conservative corrections to the orbit due to self-interactions. In Schwarzschild, the results surpass the existing result in the literature in accuracy, and we find new estimates for some of the unknown 4PN and 5PN terms in the post-Newtonian expansion of $U$. In Kerr, we provide completely novel values of $U$ for eccentric equatorial orbits. Calculation of the full self-force will appear in a forthcoming paper., Comment: Fix typo in eq 60

We calculate the kick generated by an eccentric black hole binary inspiral as it evolves through a resonant orbital configuration where the precession of the system temporarily halts. As a result, the effects of the asymmetric emission of gravitational waves build up coherently over a large number of orbits. Our results are calculate using black hole perturbation theory in the limit where the ratio of the masses of the orbiting objects $\epsilon=m/M$ is small. The resulting kick velocity scales as $\epsilon^{3/2}$, much faster than the $\epsilon^2$ scaling of the kick generated by the final merger. For the most extreme case of a very eccentric ($e\sim 1$) inspiral around a maximally spinning black hole, we find kicks close to $30,000\;\epsilon^{3/2}$~km/s, enough to dislodge a black hole from its host cluster or even galaxy. In reality, such extreme inspirals should be very rare. Nonetheless, the astrophysical impact of kicks in less extreme inspirals could be astrophysically significant., Comment: Updated to match published version

We investigate the possibility of sustained orbital resonances in extreme mass ratio inspirals. Using a near-identity averaging transformation, we reduce the equations of motion for a particle moving in Kerr spacetime with self-force corrections in the neighbourhood of a resonant geodesic to a one dimensional equation for a particle moving in an effective potential. From this effective equation we obtain the necessary and sufficient conditions that the self-force needs to satisfy to allow inspiralling orbits to be captured in sustained resonance. Along the way we also obtain the full non-linear expression for the jump in the adiabatic constants of motion incurred as an inspiral transiently evolves through a strong resonance to first-order in the mass ratio. Finally, we find that if the resonance is strong enough to allow capture in sustained resonance, only a small fraction (order of the square root of mass-ratio) of all inspirals will indeed be captured. This makes observation of sustained resonances in EMRIs -- if they exist -- very unlikely for space based observatories like eLisa., Comment: 14 pages, 6 figures, corrected a number of typesetting errors in Sec. IIID

We study the geometry generated by a massless cosmic string. We find that this is given by a Riemann flat spacetime with a conical singularity along the worldsheet of the string. The geometry of such a spacetime is completely fixed by the holonomy of a simple loop wrapping the conical singularity. In the case of a massless cosmic string, this holonomy is a null-rotation/parabolic Lorentz transformation with a parabolic angle given by the linear energy density of the cosmic string. This description explicitly shows that there is no gravitational shockwave accompanying the massless cosmic string as has been suggested in the past. To illustrate the non-singular nature of the surrounding geometry, we construct a metric for the massless cosmic string that is smooth everywhere outside the conical singularity., Comment: (13 pages, 4 figures)

We study a model for gravity in 3+1 dimensions, inspired in general relativity in 2+1 dimensions. In contrast regular general relativity in 3+1 dimensions, the model postulates that space in absence of matter is flat. The requirement that the Einstein equation still holds for the complete spacetime, implies that matter may only appear as conical defects of co-dimension 2, which may be interpreted as straight cosmic strings moving at a constant velocity. The study of collisions of these defects reveals that the dynamics of the model is incomplete. Certain highly energetic collisions of almost parallel defects suggest that no dynamic completion may exist that is fully compatible with the principles on which the model was based. We also study the phase space of the model in the continuum limit. We find that even though does not contain gravitational waves at the fundamental level, they do appear as an emergent feature in the continuum limit., Comment: PhD thesis, Utrecht University

We generalize our previous linear result [1] in obtaining gravitational waves from our piecewise flat model for gravity in 3+1 dimensions to exact piecewise flat configurations describing exact planar gravitational waves. We show explicitly how to construct a piecewise flat spacetime that describes an impulsive plane wavefront. From these wavefronts more general plane waves may be constructed.

We examine the continuum limit of the piecewise flat locally finite gravity model introduced by 't Hooft. In the linear weak field limit we find the energy--momentum tensor and metric perturbation of an arbitrary configuration of defects. The energy--momentum turns out to be restricted to satisfy certain conditions. The metric perturbation is mostly fixed by the energy--momentum except for its lightlike modes which reproduce linear gravitational waves, despite no such waves being present at the microscopic level., Comment: Accepted for publication in Classical and Quantum gravity

We consider the (3+1)-dimensional locally finite gravity model proposed by 't Hooft. In particular we revisit the problem of resolving collisions of string defects. We provide a new geometric description of the configurations of strings using piecewise flat manifolds, and use it to resolve a more general class of collisions. We argue that beyond certain bounds for the deficiency/surplus angles no resolutions may be found that satisfy the imposed causality conditions., Comment: 29 pages, 12 figures