An Efficient Numerical Method for Computing Gravitational Waves Induced by a Particle Moving on Eccentric Inclined Orbits around a Kerr Black Hole

Gravitational Waveforms for Compact Binaries from Second-Order Self-Force Theory

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 modeling extreme-mass-ratio inspirals for the LISA mission and intermediate-mass-ratio systems currently being observed by the LIGO-Virgo-KAGRA Collaboration.

Adiabatic Waveforms from Extreme-Mass-Ratio Inspirals: An Analytical Approach

Scientific analysis for the gravitational wave detector LISA will require theoretical waveforms from extreme-mass-ratio inspirals (EMRIs) that extensively cover all possible orbital and spin configurations around astrophysical Kerr black holes. However, on-the-fly calculations of these waveforms have not yet overcome the high dimensionality of the parameter space. To confront this challenge, we present a user-ready EMRI waveform model for generic (eccentric and inclined) orbits in Kerr spacetime, using an analytical self-force approach. Our model accurately covers all EMRIs with arbitrary inclination and black hole spin, up to modest eccentricity (≲0.3) and separation (≳2-10  M from the last stable orbit). In that regime, our waveforms are accurate at the leading "adiabatic" order, and they approximately capture transient self-force resonances that significantly impact the gravitational wave phase. The model fills an urgent need for extensive waveforms in ongoing data-analysis studies, and its individual components will continue to be useful in future science-adequate waveforms.

Rapid Generation of Fully Relativistic Extreme-Mass-Ratio-Inspiral Waveform Templates for LISA Data Analysis

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 templates is hindered by computing the ∼10^{3}-10^{5} harmonic modes in a fully relativistic waveform. We use order-reduction and deep-learning techniques to derive a global fit for the ≈4000 modes in the special case of an eccentric Schwarzschild orbit, and implement the fit in a complete waveform framework with hardware acceleration. Our high-fidelity waveforms can be generated in under 1 s, and achieve a mismatch of ≲5×10^{-4} against reference waveforms that take ≳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.

Intermediate Mass-Ratio Black Hole Binaries: Applicability of Small Mass-Ratio Perturbation Theory

The inspiral phasing of binary black holes at intermediate mass ratios (m_{2}/m_{1}∼10^{-3}) is important for gravitational wave observations, but not accessible to standard modeling techniques: The accuracy of the small mass-ratio (SMR) expansion is unknown at intermediate mass ratios, whereas numerical relativity simulations cannot reach this regime. This article assesses the accuracy of the SMR expansion by extracting the first three terms of the SMR expansion from numerical relativity data for nonspinning, quasicircular binaries. We recover the leading term predicted by SMR theory and obtain a robust prediction of the next-to-leading term. The influence of higher-order terms is bounded to be small, indicating that the SMR series truncated at next-to-leading order is quite accurate at intermediate mass ratios and even at nearly comparable mass binaries. We estimate the range of applicability for SMR and post-Newtonian series for nonspinning, quasicircular inspirals.

Exploring New Physics Frontiers Through Numerical Relativity

The demand to obtain answers to highly complex problems within strong-field gravity has been met with significant progress in the numerical solution of Einstein's equations - along with some spectacular results - in various setups. We review techniques for solving Einstein's equations in generic spacetimes, focusing on fully nonlinear evolutions but also on how to benchmark those results with perturbative approaches. The results address problems in high-energy physics, holography, mathematical physics, fundamental physics, astrophysics and cosmology.

Transient resonances in the inspirals of point particles into black holes

We show that transient resonances occur in the two-body problem in general relativity for spinning black holes in close proximity to one another when one black hole is much more massive than the other. These resonances occur when the ratio of polar and radial orbital frequencies, which is slowly evolving under the influence of gravitational radiation reaction, passes through a low order rational number. At such points, the adiabatic approximation to the orbital evolution breaks down, and there is a brief but order unity correction to the inspiral rate. The resonances cause a perturbation to orbital phase of order a few tens of cycles for mass ratios ∼10(-6), make orbits more sensitive to changes in initial data (though not quite chaotic), and are genuine nonperturbative effects that are not seen at any order in a standard post-Newtonian expansion. Our results apply to an important potential source of gravitational waves, the gradual inspiral of white dwarfs, neutron stars, or black holes into much more massive black holes. Resonances' effects will increase the computational challenge of accurately modeling these sources.

Analytic Black Hole Perturbation Approach to Gravitational Radiation

We review the analytic methods used to perform the post-Newtonian expansion of gravitational waves induced by a particle orbiting a massive, compact body, based on black hole perturbation theory. There exist two different methods of performing the post-Newtonian expansion. Both are based on the Teukolsky equation. In one method, the Teukolsky equation is transformed into a Regge-Wheeler type equation that reduces to the standard Klein Gordon equation in the flat-space limit, while in the other method (which was introduced by Mano, Suzuki, and Takasugi relatively recently, the Teukolsky equation is used directly in its original form. The former's advantage is that it is intuitively easy to understand how various curved space effects come into play. However, it becomes increasingly complicated when one goes to higher and higher post-Newtonian orders. In contrast, the latter's advantage is that a systematic calculation to higher post-Newtonian orders can be implemented relatively easily, but otherwise, it is so mathematical that it is hard to understand the interplay of higher order terms. In this paper, we review both methods so that their pros and cons may be seen clearly. We also review some results of calculations of gravitational radiation emitted by a particle orbiting a black hole.