-
Parameter estimates of GW150914 were obtained using Bayesian inference, based
on three semi-analytic waveform models for binary black hole coalescences.
These waveform models differ from each other in their treatment of black hole
spins, and all three models make some simplifying assumptions, notably to
neglect sub-dominant waveform harmonic modes and orbital eccentricity.
Furthermore, while the models are calibrated to agree with waveforms obtained
by full numerical solutions of Einstein's equations, any such calibration is
accurate only to some non-zero tolerance and is limited by the accuracy of the
underlying phenomenology, availability, quality, and parameter-space coverage
of numerical simulations. This paper complements the original analyses of
GW150914 with an investigation of the effects of possible systematic errors in
the waveform models on estimates of its source parameters. To test for
systematic errors we repeat the original Bayesian analyses on mock signals from
numerical simulations of a series of binary configurations with parameters
similar to those found for GW150914. Overall, we find no evidence for a
systematic bias relative to the statistical error of the original parameter
recovery of GW150914 due to modeling approximations or modeling inaccuracies.
However, parameter biases are found to occur for some configurations disfavored
by the data of GW150914: for binaries inclined edge-on to the detector over a
small range of choices of polarization angles, and also for eccentricities
greater than $\sim$0.05. For signals with higher signal-to-noise ratio than
GW150914, or in other regions of the binary parameter space (lower masses,
larger mass ratios, or higher spins), we expect that systematic errors in
current waveform models may impact gravitational-wave measurements, making more
accurate models desirable for future observations.
-
This paper presents updated estimates of source parameters for GW150914, a
binary black-hole coalescence event detected by the Laser Interferometer
Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference
presented parameter estimation [2] of the source using a 13-dimensional,
phenomenological precessing-spin model (precessing IMRPhenom) and a
11-dimensional nonprecessing effective-one-body (EOB) model calibrated to
numerical-relativity simulations, which forces spin alignment (nonprecessing
EOBNR). Here we present new results that include a 15-dimensional
precessing-spin waveform model (precessing EOBNR) developed within the EOB
formalism. We find good agreement with the parameters estimated previously [2],
and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and
$30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric
credible intervals). We also present slightly tighter constraints on the
dimensionless spin magnitudes of the two black holes, with a primary spin
estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability.
Reference [2] estimated the systematic parameter-extraction errors due to
waveform-model uncertainty by combining the posterior probability densities of
precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two
precessing-spin models are in closer agreement, suggesting that these
systematic errors are smaller than previously quoted.
-
We compare GW150914 directly to simulations of coalescing binary black holes
in full general relativity, accounting for all the spin-weighted quadrupolar
modes, and separately accounting for all the quadrupolar and octopolar modes.
Consistent with the posterior distributions reported in LVC_PE[1] (at 90%
confidence), we find the data are compatible with a wide range of nonprecessing
and precessing simulations. Followup simulations performed using
previously-estimated binary parameters most resemble the data. Comparisons
including only the quadrupolar modes constrain the total redshifted mass Mz \in
[64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin
\chi_eff \in [-0.3, 0.2], where \chi_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M.
Including both quadrupolar and octopolar modes, we find the mass ratio is even
more tightly constrained. Simulations with extreme mass ratios and effective
spins are highly inconsistent with the data, at any mass. Several nonprecessing
and precessing simulations with similar mass ratio and \chi_{eff} are
consistent with the data. Though correlated, the components' spins (both in
magnitude and directions) are not significantly constrained by the data. For
nonprecessing binaries, interpolating between simulations, we reconstruct a
posterior distribution consistent with previous results. The final black hole's
redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the
final black hole's dimensionless spin parameter is consistent with af = 0.62 -
0.73. As our approach invokes no intermediate approximations to general
relativity and can strongly reject binaries whose radiation is inconsistent
with the data, our analysis provides a valuable complement to LVC_PE[1].
-
The Numerical INJection Analysis (NINJA) project is a collaborative effort
between members of the numerical relativity and gravitational-wave astrophysics
communities. The purpose of NINJA is to study the ability to detect
gravitational waves emitted from merging binary black holes and recover their
parameters with next-generation gravitational-wave observatories. We report
here on the results of the second NINJA project, NINJA-2, which employs 60
complete binary black hole hybrid waveforms consisting of a numerical portion
modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian
portion modelling the early inspiral. In a "blind injection challenge" similar
to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid
waveforms to two months of data recolored to predictions of Advanced LIGO and
Advanced Virgo sensitivity curves during their first observing runs. The
resulting data was analyzed by gravitational-wave detection algorithms and 6 of
the waveforms were recovered with false alarm rates smaller than 1 in a
thousand years. Parameter estimation algorithms were run on each of these
waveforms to explore the ability to constrain the masses, component angular
momenta and sky position of these waveforms. We also perform a large-scale
monte-carlo study to assess the ability to recover each of the 60 hybrid
waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our
results predict that early Advanced LIGO and Advanced Virgo will have a
volume-weighted average sensitive distance of 300Mpc (1Gpc) for
$10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole
coalescences. We demonstrate that neglecting the component angular momenta in
the waveform models used in matched-filtering will result in a reduction in
sensitivity for systems with large component angular momenta. [Abstract
abridged for ArXiv, full version in PDF]
-
This document proposes data formats to exchange numerical relativity results,
in particular gravitational waveforms. The primary goal is to further the
interaction between gravitational-wave source modeling groups and the
gravitational-wave data-analysis community. We present a simple and extendable
format which is applicable to various kinds of gravitational wave sources
including binaries of compact objects and systems undergoing gravitational
collapse, but is nevertheless sufficiently general to be useful for other
purposes.
-
The spin of the final black hole in the coalescence of nonspinning black
holes is determined by the ``residual'' orbital angular momentum of the binary.
This residual momentum consists of the orbital angular momentum that the binary
is not able to shed in the process of merging. We study the angular momentum
radiated, the spin of the final black hole and the gravitational bursts in a
series of orbits ranging from almost direct infall to numerous orbits before
infall that exhibit multiple bursts of radiation in the merger process. We show
that the final black hole gets a maximum spin parameter $a/M_h \le 0.78$, and
this maximum occurs for initial orbital angular momentum $L \approx M^2_h$.
-
We discuss results that have been obtained from the implementation of the
initial round of testbeds for numerical relativity which was proposed in the
first paper of the Apples with Apples Alliance. We present benchmark results
for various codes which provide templates for analyzing the testbeds and to
draw conclusions about various features of the codes. This allows us to sharpen
the initial test specifications, design a new test and add theoretical insight.
-
We present results from fully nonlinear simulations of unequal mass binary
black holes plunging from close separations well inside the innermost stable
circular orbit with mass ratios q = M_1/M_2 = {1,0.85,0.78,0.55,0.32}, or
equivalently, with reduced mass parameters $\eta=M_1M_2/(M_1+M_2)^2 = {0.25,
0.248, 0.246, 0.229, 0.183}$. For each case, the initial binary orbital
parameters are chosen from the Cook-Baumgarte equal-mass ISCO configuration. We
show waveforms of the dominant l=2,3 modes and compute estimates of energy and
angular momentum radiated. For the plunges from the close separations
considered, we measure kick velocities from gravitational radiation recoil in
the range 25-82 km/s. Due to the initial close separations our kick velocity
estimates should be understood as a lower bound. The close configurations
considered are also likely to contain significant eccentricities influencing
the recoil velocity.
Parameter estimates of GW150914 were obtained using Bayesian inference, based on three semi-analytic waveform models for binary black hole coalescences. These waveform models differ from each other in their treatment of black hole spins, and all three models make some simplifying assumptions, notably to neglect sub-dominant waveform harmonic modes and orbital eccentricity. Furthermore, while the models are calibrated to agree with waveforms obtained by full numerical solutions of Einstein's equations, any such calibration is accurate only to some non-zero tolerance and is limited by the accuracy of the underlying phenomenology, availability, quality, and parameter-space coverage of numerical simulations. This paper complements the original analyses of GW150914 with an investigation of the effects of possible systematic errors in the waveform models on estimates of its source parameters. To test for systematic errors we repeat the original Bayesian analyses on mock signals from numerical simulations of a series of binary configurations with parameters similar to those found for GW150914. Overall, we find no evidence for a systematic bias relative to the statistical error of the original parameter recovery of GW150914 due to modeling approximations or modeling inaccuracies. However, parameter biases are found to occur for some configurations disfavored by the data of GW150914: for binaries inclined edge-on to the detector over a small range of choices of polarization angles, and also for eccentricities greater than $\sim$0.05. For signals with higher signal-to-noise ratio than GW150914, or in other regions of the binary parameter space (lower masses, larger mass ratios, or higher spins), we expect that systematic errors in current waveform models may impact gravitational-wave measurements, making more accurate models desirable for future observations.
This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [2], and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and $30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability. Reference [2] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted.
We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin \chi_eff \in [-0.3, 0.2], where \chi_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and \chi_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1].
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]
This document proposes data formats to exchange numerical relativity results, in particular gravitational waveforms. The primary goal is to further the interaction between gravitational-wave source modeling groups and the gravitational-wave data-analysis community. We present a simple and extendable format which is applicable to various kinds of gravitational wave sources including binaries of compact objects and systems undergoing gravitational collapse, but is nevertheless sufficiently general to be useful for other purposes.
The spin of the final black hole in the coalescence of nonspinning black holes is determined by the ``residual'' orbital angular momentum of the binary. This residual momentum consists of the orbital angular momentum that the binary is not able to shed in the process of merging. We study the angular momentum radiated, the spin of the final black hole and the gravitational bursts in a series of orbits ranging from almost direct infall to numerous orbits before infall that exhibit multiple bursts of radiation in the merger process. We show that the final black hole gets a maximum spin parameter $a/M_h \le 0.78$, and this maximum occurs for initial orbital angular momentum $L \approx M^2_h$.
We discuss results that have been obtained from the implementation of the initial round of testbeds for numerical relativity which was proposed in the first paper of the Apples with Apples Alliance. We present benchmark results for various codes which provide templates for analyzing the testbeds and to draw conclusions about various features of the codes. This allows us to sharpen the initial test specifications, design a new test and add theoretical insight.