Assessing the Energetics of Spinning Binary Black Hole Systems

In this work we study the dynamics of spinning binary black hole systems in the strong field regime. For this purpose we extract from numerical relativity simulations the binding energy, specific orbital angular momentum, and gauge-invariant orbital frequency. The goal of our work is threefold: First, we extract the individual spin contributions to the binding energy, in particular the spin-orbit, spin-spin, and cubic-in-spin terms. Second, we compare our results with predictions from waveform models and find that while post-Newtonian approximants are not capable of representing the dynamics during the last few orbits before merger, there is good agreement between our data and effective-one-body approximants as well as the numerical relativity surrogate models. Finally, we present phenomenological representations for the binding energy for non-spinning systems with mass ratios up to $q = 10$ and for the spin-orbit interaction for mass ratios up to $q = 8$ obtaining accuracies of $\lesssim 0.1\%$ and $\lesssim 6\%$, respectively

Distinguishing Boson Stars from Black Holes and Neutron Stars from Tidal Interactions in Inspiraling Binary Systems

Binary systems containing boson stars---self-gravitating configurations of a complex scalar field--- can potentially mimic black holes or neutron stars as gravitational-wave sources. We investigate the extent to which tidal effects in the gravitational-wave signal can be used to discriminate between these standard sources and boson stars. We consider spherically symmetric boson stars within two classes of scalar self-interactions: an effective-field-theoretically motivated quartic potential and a solitonic potential constructed to produce very compact stars. We compute the tidal deformability parameter characterizing the dominant tidal imprint in the gravitational-wave signals for a large span of the parameter space of each boson star model. We find that the tidal deformability for boson stars with a quartic self-interaction is bounded below by $\Lambda_{\rm min}\approx 280$ and for those with a solitonic interaction by $\Lambda_{\rm min}\approx 1.3$. Employing a Fisher matrix analysis, we estimate the precision with which Advanced LIGO and third-generation detectors can measure these tidal parameters using the inspiral portion of the signal. We discuss a new strategy to improve the distinguishability between black holes/neutrons stars and boson stars by combining deformability measurements of each compact object in a binary system, thereby eliminating the scaling ambiguities in each boson star model. Our analysis shows that current-generation detectors can potentially distinguish boson stars with quartic potentials from black holes, as well as from neutron-star binaries if they have either a large total mass or a large mass ratio. Discriminating solitonic boson stars from black holes using only tidal effects during the inspiral will be difficult with Advanced LIGO, but third-generation detectors should be able to distinguish between binary black holes and these binary boson stars.Comment: 18 pages, 8 figures. Submitted to Physical Review

Improvements to the construction of binary black hole initial data

Construction of binary black hole initial data is a prerequisite for numerical evolutions of binary black holes. This paper reports improvements to the binary black hole initial data solver in the Spectral Einstein Code, to allow robust construction of initial data for mass-ratio above 10:1, and for dimensionless black hole spins above 0.9, while improving efficiency for lower mass-ratios and spins. We implement a more flexible domain decomposition, adaptive mesh refinement and an updated method for choosing free parameters. We also introduce a new method to control and eliminate residual linear momentum in initial data for precessing systems, and demonstrate that it eliminates gravitational mode mixing during the evolution. Finally, the new code is applied to construct initial data for hyperbolic scattering and for binaries with very small separation.Comment: 28 pages, 13 figures, 1 tabl

Enriching the Symphony of Gravitational Waves from Binary Black Holes by Tuning Higher Harmonics

For the first time, we construct an inspiral-merger-ringdown waveform model within the effective-one-body formalism for spinning, nonprecessing binary black holes that includes gravitational modes beyond the dominant $(\ell,|m|) = (2,2)$ mode, specifically $(\ell,|m|)=(2,1),(3,3),(4,4),(5,5)$. Our multipolar waveform model incorporates recent (resummed) post-Newtonian results for the inspiral and information from 157 numerical-relativity simulations, and 13 waveforms from black-hole perturbation theory for the (plunge-)merger and ringdown. We quantify the improved accuracy including higher-order modes by computing the faithfulness of the waveform model against the numerical-relativity waveforms used to construct the model. We define the faithfulness as the match maximized over time, phase of arrival, gravitational-wave polarization and sky position of the waveform model, and averaged over binary orientation, gravitational-wave polarization and sky position of the numerical-relativity waveform. When the waveform model contains only the $(2,2)$ mode, we find that the averaged faithfulness to numerical-relativity waveforms containing all modes with $\ell \leq$ 5 ranges from $90\%$ to $99.9\%$ for binaries with total mass $20-200 M_\odot$ (using the Advanced LIGO's design noise curve). By contrast, when the $(2,1),(3,3),(4,4),(5,5)$ modes are also included in the model, the faithfulness improves to $99\%$ for all but four configurations in the numerical-relativity catalog, for which the faithfulness is greater than $98.5\%$. Using our results, we also develop also a (stand-alone) waveform model for the merger-ringdown signal, calibrated to numerical-relativity waveforms, which can be used to measure multiple quasi-normal modes. The multipolar waveform model can be extended to include spin-precession, and will be employed in upcoming observing runs of Advanced LIGO and Virgo.Comment: 28 page

Comparing Post-Newtonian and Numerical-Relativity Precession Dynamics

Binary black-hole systems are expected to be important sources of gravitational waves for upcoming gravitational-wave detectors. If the spins are not colinear with each other or with the orbital angular momentum, these systems exhibit complicated precession dynamics that are imprinted on the gravitational waveform. We develop a new procedure to match the precession dynamics computed by post-Newtonian (PN) theory to those of numerical binary black-hole simulations in full general relativity. For numerical relativity NR) simulations lasting approximately two precession cycles, we find that the PN and NR predictions for the directions of the orbital angular momentum and the spins agree to better than $\sim 1^{\circ}$ with NR during the inspiral, increasing to $5^{\circ}$ near merger. Nutation of the orbital plane on the orbital time-scale agrees well between NR and PN, whereas nutation of the spin direction shows qualitatively different behavior in PN and NR. We also examine how the PN equations for precession and orbital-phase evolution converge with PN order, and we quantify the impact of various choices for handling partially known PN terms

Theoretical groundwork supporting the precessing-spin two-body dynamics of the effective-one-body waveform models SEOBNRv5

Waveform models are essential for gravitational-wave (GW) detection and parameter estimation of coalescing compact-object binaries. More accurate models are required for the increasing sensitivity of current and future GW detectors. The effective-one-body (EOB) formalism combines the post-Newtonian (PN) and small mass-ratio approximations with numerical-relativity results, and produces highly accurate inspiral-merger-ringdown waveforms. In this paper, we derive the analytical precessing-spin two-body dynamics for the \texttt{SEOBNRv5} waveform model, which has been developed for the upcoming LIGO-Virgo-KAGRA observing run. We obtain an EOB Hamiltonian that reduces to the exact Kerr Hamiltonian in the test-mass limit. It includes the full 4PN precessing-spin information, and is valid for generic compact objects (i.e., for black holes or neutron stars). We also build an efficient and accurate EOB Hamiltonian that includes partial precessional effects, notably orbit-averaged in-plane spin effects for circular orbits, and derive 4PN-expanded precessing-spin equations of motion, consistent with such an EOB Hamiltonian. The results were used to build the computationally-efficient precessing-spin multipolar \texttt{SEOBNRv5PHM} waveform model.Comment: 35 page

Binary Neutron Stars with Arbitrary Spins in Numerical Relativity

We present a code to construct initial data for binary neutron star systems in which the stars are rotating. Our code, based on a formalism developed by Tichy, allows for arbitrary rotation axes of the neutron stars and is able to achieve rotation rates near rotational breakup. We compute the neutron star angular momentum through quasi-local angular momentum integrals. When constructing irrotational binary neutron stars, we find a very small residual dimensionless spin of $\sim 2\times 10^{-4}$. Evolutions of rotating neutron star binaries show that the magnitude of the stars' angular momentum is conserved, and that the spin- and orbit-precession of the stars is well described by post-Newtonian approximation. We demonstrate that orbital eccentricity of the binary neutron stars can be controlled to $\sim 0.1\%$. The neutron stars show quasi-normal mode oscillations at an amplitude which increases with the rotation rate of the stars.Comment: 20 pages, 22 figure

pySEOBNR: a software package for the next generation of effective-one-body multipolar waveform models

We present pySEOBNR, a Python package for gravitational-wave (GW) modeling developed within the effective-one-body (EOB) formalism. The package contains an extensive framework to generate state-of-the-art inspiral-merger-ringdown waveform models for compact-object binaries composed of black holes and neutron stars. We document and demonstrate how to use the built-in quasi-circular precessing-spin model SEOBNRv5PHM, whose aligned-spin limit (SEOBNRv5HM) has been calibrated to numerical-relativity simulations and the nonspinning sector to gravitational self-force data using pySEOBNR. Furthermore, pySEOBNR contains the infrastructure necessary to construct, calibrate, test, and profile new waveform models in the EOB approach. The efficiency and flexibility of pySEOBNR will be crucial to overcome the data-analysis challenges posed by upcoming and next-generation GW detectors on the ground and in space, which will afford the possibility to observe all compact-object binaries in our Universe.Comment: 21 pages, 4 figure