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 (?,|m|)=(2,2) mode, specifically (?,|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 improvement in accuracy when 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 ??5 ranges from 90% to 99.9% for binaries with total mass 20–200 M? (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%. Starting from the complete inspiral-merger-ringdown model, we 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-precessional effects, and will be employed in upcoming observing runs of Advanced LIGO and Virgo.
As gravitational-wave detectors become more sensitive and broaden their frequency bandwidth, we will access a greater variety of signals emitted by compact binary systems, shedding light on their ...astrophysical origin and environment. A key physical effect that can distinguish among different formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such precessing signals, especially when masses and spins vary in the wide astrophysical range, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the entire coalescence stage (i.e., inspiral, merger and ringdown) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from black-hole perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to 118 SXS precessing NR waveforms, produced as part of this project, which span mass ratios 1-4 and (dimensionless) black-hole's spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for 94% (57%) of the cases, the maximum value, in the total mass range 20 − 200 M⊙, is below 3% (1%). Those numbers change to 83% (20%) when using the inspiral-merger-ringdown, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two Bayesian, parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models as a function of the mass and spin parameters to identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.
This Letter presents a publicly available catalog of 174 numerical binary black hole simulations following up to 35 orbits. The catalog includes 91 precessing binaries, mass ratios up to 8∶1, orbital ...eccentricities from a few percent to 10(-5), black hole spins up to 98% of the theoretical maximum, and radiated energies up to 11.1% of the initial mass. We establish remarkably good agreement with post-Newtonian precession of orbital and spin directions for two new precessing simulations, and we discuss other applications of this catalog. Formidable challenges remain: e.g., precession complicates the connection of numerical and approximate analytical waveforms, and vast regions of the parameter space remain unexplored.
Binary black hole simulations with black hole excision using spectral methods require a coordinate transformation into a co rotating coordinate system where the black holes are essentially at rest. ...This paper presents and discusses two coordinate transformations that are applicable to precessing binary systems, one based on Euler angles, the other on quaternions. Both approaches are found to work well for binaries with moderate precession, i.e., for cases where the orientation of the orbital plane changes by << 90degrees. For strong precession, performance of the Euler-angle parametrization deteriorates, eventually failing fora 90degrees change in orientation because of singularities in the parametrization ("gimbal lock"). In contrast, the quaternion representation is invariant under an overall rotation and handles any orientation of the orbital plane as well as the Euler-angle technique handles nonprecessing binaries.
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, covering the entire space in the quartic case, and an extensive portion of interest in the solitonic case. We find that the tidal deformability for boson stars with a quartic self-interaction is bounded below by Λmin≈280 and for those with a solitonic interaction by Λmin≈1.3. We summarize our results as ready-to-use fits for practical applications. 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 novel strategy to improve the distinguishability between black holes/neutrons stars and boson stars by combining tidal 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 (asymmetric) 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.
With the first direct detections of gravitational waves (GWs) from the coalescence of compact binaries observed by the advanced LIGO and VIRGO interferometers, the era of GW astronomy has begun. ...Whilst there is strong evidence that the observed GWs are connected to the merger of two black holes (BH) or two neutron stars (NS), future detections may present a less consistent picture. Indeed, the possibility that the observed GW signal was created by a merger of exotic compact objects (ECOs) such as boson stars (BS) or axion stars (AS) has not yet been fully excluded. For a detailed understanding of the late stages of the coalescence full 3D numerical relativity simulations are essential. In this paper, we extend the infrastructure of the numerical relativity code BAM, to permit the simultaneous simulation of baryonic matter with bosonic scalar fields, thus enabling the study of BS-BS, BS-NS, and BS-BH mergers. We present a large number of single star evolutions to test the newly implemented routines, and to quantify the numerical challenges of such simulations, which we find to partially differ from the default NS case. We also compare head-on BS-BS simulations with independent numerical relativity codes, namely the SpEC and the GRChombo codes, and find good general agreement. Finally, we present what are, to the best of our knowledge, the first full NR simulations of BS-NS mergers, a first step towards identifying the hallmarks of BS-NS interactions in the strong gravity regime, as well as possible GW and electromagnetic observables.
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 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 nonspinning systems with mass ratios up to q=10 and for the spin-orbit interaction for mass ratios up to q=8 obtaining accuracies of ≲0.1% and ≲6%, respectively.