Numerical simulations of compact object binaries

Coalescing compact object binaries consisting of black holes and/or Neutron stars are a prime target for ground-based gravitational wave detectors. This article reviews the status of numerical simulations of these systems, with an emphasis on recent progress.Comment: Minor corrections, including typos and gramma

Measuring neutron star tidal deformability with Advanced LIGO: a Bayesian analysis of neutron star - black hole binary observations

The discovery of gravitational waves (GW) by Advanced LIGO has ushered us into an era of observational GW astrophysics. Compact binaries remain the primary target sources for LIGO, of which neutron star-black hole (NSBH) binaries form an important subset. GWs from NSBH sources carry signatures of (a) the tidal distortion of the neutron star by its companion black hole during inspiral, and (b) its potential tidal disruption near merger. In this paper, we present a Bayesian study of the measurability of neutron star tidal deformability $\Lambda_\mathrm{NS}\propto (R/M)^{5}$ using observation(s) of inspiral-merger GW signals from disruptive NSBH coalescences, taking into account the crucial effect of black hole spins. First, we find that if non-tidal templates are used to estimate source parameters for an NSBH signal, the bias introduced in the estimation of non-tidal physical parameters will only be significant for loud signals with signal-to-noise ratios $> 30$. For similarly loud signals, we also find that we can begin to put interesting constraints on $\Lambda_\mathrm{NS}$ (factor of 1-2) with individual observations. Next, we study how a population of realistic NSBH detections will improve our measurement of neutron star tidal deformability. For astrophysical populations of $disruptive$ NSBH mergers, we find 20-35 events to be sufficient to constrain $\Lambda_\mathrm{NS}$ within $\pm 25-50\%$, depending on the chosen equation of state. In this we also assume that LIGO will detect black holes with masses within the astrophysical $mass$-$gap$. If the mass-gap remains preserved in NSBHs detected by LIGO, we estimate that $25\%$ $additional$ detections will furnish comparable tidal measurement accuracy. In both cases, we find that the loudest 5-10 events to provide most of the tidal information, thereby facilitating targeted follow-ups of NSBHs in the upcoming LIGO-Virgo runs.Comment: 21 pages, 17 figure

Numerical Relativity Injection Infrastructure

This document describes the new Numerical Relativity (NR) injection infrastructure in the LIGO Algorithms Library (LAL), which henceforth allows for the usage of NR waveforms as a discrete waveform approximant in LAL. With this new interface, NR waveforms provided in the described format can directly be used as simulated GW signals ("injections") for data analyses, which include parameter estimation, searches, hardware injections etc. As opposed to the previous infrastructure, this new interface natively handles sub-dominant modes and waveforms from numerical simulations of precessing binary black holes, making them directly accessible to LIGO analyses. To correctly handle precessing simulations, the new NR injection infrastructure internally transforms the NR data into the coordinate frame convention used in LAL.Comment: 20 pages, 2 figures, technical repor

Compact Binary Waveform Center-of-Mass Corrections

We present a detailed study of the center-of-mass (c.m.) motion seen in simulations produced by the Simulating eXtreme Spacetimes (SXS) collaboration. We investigate potential physical sources for the large c.m. motion in binary black hole simulations and find that a significant fraction of the c.m. motion cannot be explained physically, thus concluding that it is largely a gauge effect. These large c.m. displacements cause mode mixing in the gravitational waveform, most easily recognized as amplitude oscillations caused by the dominant (2,$\pm$2) modes mixing into subdominant modes. This mixing does not diminish with increasing distance from the source; it is present even in asymptotic waveforms, regardless of the method of data extraction. We describe the current c.m.-correction method used by the SXS collaboration, which is based on counteracting the motion of the c.m. as measured by the trajectories of the apparent horizons in the simulations, and investigate potential methods to improve that correction to the waveform. We also present a complementary method for computing an optimal c.m. correction or evaluating any other c.m. transformation based solely on the asymptotic waveform data.Comment: 20 pages, 15 figure

IMEX evolution of scalar fields on curved backgrounds

Inspiral of binary black holes occurs over a time-scale of many orbits, far longer than the dynamical time-scale of the individual black holes. Explicit evolutions of a binary system therefore require excessively many time steps to capture interesting dynamics. We present a strategy to overcome the Courant-Friedrichs-Lewy condition in such evolutions, one relying on modern implicit-explicit ODE solvers and multidomain spectral methods for elliptic equations. Our analysis considers the model problem of a forced scalar field propagating on a generic curved background. Nevertheless, we encounter and address a number of issues pertinent to the binary black hole problem in full general relativity. Specializing to the Schwarzschild geometry in Kerr-Schild coordinates, we document the results of several numerical experiments testing our strategy.Comment: 28 pages, uses revtex4. Revised in response to referee's report. One numerical experiment added which incorporates perturbed initial data and adaptive time-steppin

Approximate initial data for binary black holes

We construct approximate analytical solutions to the constraint equations of general relativity for binary black holes of arbitrary mass ratio in quasicircular orbit. We adopt the puncture method to solve the constraint equations in the transverse-traceless decomposition and consider perturbations of Schwarzschild black holes caused by boosts and the presence of a binary companion. A superposition of these two perturbations then yields approximate, but fully analytic binary black hole initial data that are accurate to first order in the inverse of the binary separation and the square of the black holes' momenta.Comment: 13 pages, 4 figures, added comparison to numerical calculations, accepted to PR

Binary-black-hole initial data with nearly-extremal spins

There is a significant possibility that astrophysical black holes with nearly-extremal spins exist. Numerical simulations of such systems require suitable initial data. In this paper, we examine three methods of constructing binary-black-hole initial data, focusing on their ability to generate black holes with nearly-extremal spins: (i) Bowen-York initial data, including standard puncture data (based on conformal flatness and Bowen-York extrinsic curvature), (ii) standard quasi-equilibrium initial data (based on the extended-conformal-thin-sandwich equations, conformal flatness, and maximal slicing), and (iii) quasi-equilibrium data based on the superposition of Kerr-Schild metrics. We find that the two conformally-flat methods (i) and (ii) perform similarly, with spins up to about 0.99 obtainable at the initial time. However, in an evolution, we expect the spin to quickly relax to a significantly smaller value around 0.93 as the initial geometry relaxes. For quasi-equilibrium superposed Kerr-Schild (SKS) data [method (iii)], we construct initial data with \emph{initial} spins as large as 0.9997. We evolve SKS data sets with spins of 0.93 and 0.97 and find that the spin drops by only a few parts in 10^4 during the initial relaxation; therefore, we expect that SKS initial data will allow evolutions of binary black holes with relaxed spins above 0.99. [Abstract abbreviated; full abstract also mentions several secondary results.