Intermediate-mass-ratio black-hole binaries: numerical relativity meets perturbation theory

Detecting Subsolar-Mass Primordial Black Holes in Extreme Mass-Ratio Inspirals with LISA and Einstein Telescope

Primordial black holes possibly formed in the early Universe could provide a significant fraction of the dark matter and would be unique probes of inflation. A smoking gun for their discovery would be the detection of a subsolar mass compact object. We argue that extreme mass-ratio inspirals will be ideal to search for subsolar-mass black holes not only with LISA but also with third-generation ground-based detectors such as Cosmic Explorer and the Einstein Telescope. These sources can provide unparalleled measurements of the mass of the secondary object at a subpercent level for primordial black holes as light as O(0.01)  M_{⊙} up to luminosity distances around hundred megaparsec and few gigaparsec for LISA and Einstein Telescope, respectively, in a complementary frequency range. This would allow claiming, with very high statistical confidence, the detection of a subsolar-mass black hole, which would also provide a novel (and currently undetectable) family of sources for third-generation detectors.

Exploring the Small Mass Ratio Binary Black Hole Merger via Zeno's Dichotomy Approach

We perform a sequence of binary black hole simulations with increasingly small mass ratios, reaching to a 128:1 binary that displays 13 orbits before merger. Based on a detailed convergence study of the q=m_{1}/m_{2}=1/15 nonspinning case, we apply additional mesh refinement levels around the smaller hole horizon [30] to reach successively the q=1/32, q=1/64, and q=1/128 cases. Roughly a linear computational resources scaling with 1/q is observed on eight-nodes simulations. We compute the remnant properties of the merger: final mass, spin, and recoil velocity, finding precise consistency between horizon and radiation measures. We also compute the gravitational waveforms: their peak frequency, amplitude, and luminosity. We compare those values with predictions of the corresponding phenomenological formulas, reproducing the particle limit within 2%, and we then use the new results to improve their fitting coefficients.

Numerical relativity of compact binaries in the 21st century

We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star-neutron star mergers. Modeling such events requires numerical integration of the field equations of general relativity in three spatial dimensions, coupled, in the case of neutron-star containing binaries, with increasingly sophisticated treatment of fluids, electromagnetic fields, and neutrino radiation. However, it was not until 2005 that accurate long-term evolutions of binaries containing black holes were even possible (Pretorius 2005 Phys. Rev. Lett. 95 121101, Campanelli et al 2006 Phys. Rev. Lett. 96 111101, Baker et al 2006 Phys. Rev. Lett. 96 111102). Since then, there has been an explosion of new results and insights into the physics of strongly-gravitating system. Particular emphasis has been placed on understanding the gravitational wave and electromagnetic signatures from these extreme events. And with the recent dramatic discoveries of gravitational waves from merging black holes by the Laser Interferometric Gravitational Wave Observatory and Virgo, and the subsequent discovery of both electromagnetic and gravitational wave signals from a merging neutron star-neutron star binary, numerical relativity became an indispensable tool for the new field of multimessenger astronomy.

Extraction of gravitational waves in numerical relativity

A numerical-relativity calculation yields in general a solution of the Einstein equations including also a radiative part, which is in practice computed in a region of finite extent. Since gravitational radiation is properly defined only at null infinity and in an appropriate coordinate system, the accurate estimation of the emitted gravitational waves represents an old and non-trivial problem in numerical relativity. A number of methods have been developed over the years to "extract" the radiative part of the solution from a numerical simulation and these include: quadrupole formulas, gauge-invariant metric perturbations, Weyl scalars, and characteristic extraction. We review and discuss each method, in terms of both its theoretical background as well as its implementation. Finally, we provide a brief comparison of the various methods in terms of their inherent advantages and disadvantages.

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.

Gravitational self-force correction to the innermost stable circular equatorial orbit of a Kerr black hole

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 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 innermost stable circular equatorial orbit 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.

Orbital evolution of extreme-mass-ratio black-hole binaries with numerical relativity

We perform the first fully nonlinear numerical simulations of black-hole binaries with mass ratios 100∶1. Our technique is based on the moving puncture formalism with a new gauge condition and an optimal choice of the mesh refinement. The evolutions start with a small nonspinning black hole just outside the ISCO that orbits twice before plunging. We compute the gravitational radiation, as well as the final remnant parameters, and find close agreement with perturbative estimates. We briefly discuss the relevance of these simulations for Advanced LIGO, third-generation ground-based detectors, LISA observations, and self-force computations.