Title: Accounting for numerical-relativity calibration uncertainty in gravitational-wave modeling and inference
The increasing sensitivity of current and upcoming gravitational-wave (GW) detectors poses stringent requirements on the accuracy of the GW models used for data analysis. If these requirements are not met, systematic errors could dominate over statistical uncertainties, hindering our ability to extract astrophysical and cosmological information, and conduct precise tests of general relativity. In this work, we present a novel method to mitigate the impact of waveform-systematic errors, by incorporating and marginalizing over waveform-uncertainty estimates, which are modeled as probability distributions for the numerical-relativity calibration parameters of effective-one-body waveform models. By analyzing simulated GW signals of loud “golden” binary-black-hole systems, we show that our method significantly reduces biases in the recovered parameters, highlighting its potential to improve the robustness of GW parameter estimation with upcoming observing runs and next-generation ground-based facilities, such as the Einstein Telescope and Cosmic Explorer.  more » « less
Award ID(s):
2512902 2407453
PAR ID:
10683314
Author(s) / Creator(s):
; ;
Publisher / Repository:
APS
Date Published:
Journal Name:
Physical Review D
Volume:
112
Issue:
4
ISSN:
2470-0010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. While the majority of gravitational wave (GW) events observed by the LIGO and Virgo detectors are consistent with mergers of binary black holes (BBHs) on quasicircular orbits, some events are also consistent with nonzero orbital eccentricity, indicating that the binaries could have formed via dynamical interactions. Moreover, there may be GW events which show support for spin precession, eccentricity, or both. In this work, we study the interplay of spins and eccentricity on the parameter estimation of GW signals from BBH mergers. We inject eccentric signals with no spins, aligned spins, and precessing spins using hybrids, -, and new numerical relativity simulations, respectively, and study the biases in the posteriors of source parameters when these signals are recovered with a quasicircular precessing-spin waveform model, as opposed to an aligned-spin eccentric waveform model. We find significant biases in the source parameters, such as chirp mass and spin precession ( χp ), when signals from highly eccentric BBHs are recovered with a quasicircular waveform model. Moreover, we find that for signals with both eccentricity and spin-precession effects, Bayes factor calculations confirm that an eccentric, aligned-spin model is preferred over a quasicircular precessing-spin model. Our study highlights the complex nature of GW signals from eccentric, precessing-spin binaries and the need for readily usable inspiral-merger-ringdown eccentric, spin-precessing waveform models for unbiased parameter estimation. 
    more » « less
  2. ABSTRACT Tidal interactions in coalescing binary neutron stars modify the dynamics of the inspiral and hence imprint a signature on their gravitational wave (GW) signals in the form of an extra phase shift. We need accurate models for the tidal phase shift in order to constrain the supranuclear equation of state from observations. In previous studies, GW waveform models were typically constructed by treating the tide as a linear response to a perturbing tidal field. In this work, we incorporate non-linear corrections due to hydrodynamic three- and four-mode interactions and show how they can improve the accuracy and explanatory power of waveform models. We set up and numerically solve the coupled differential equations for the orbit and the modes and analytically derive solutions of the system’s equilibrium configuration. Our analytical solutions agree well with the numerical ones up to the merger and involve only algebraic relations, allowing for fast phase shift and waveform evaluations for different equations of state over a large parameter space. We find that, at Newtonian order, non-linear fluid effects can enhance the tidal phase shift by $$\gtrsim 1\, {\rm radian}$$ at a GW frequency of 1000 Hz, corresponding to a $$10{{\%}}-20{{\%}}$$ correction to the linear theory. The scale of the additional phase shift near the merger is consistent with the difference between numerical relativity and theoretical predictions that account only for the linear tide. Non-linear fluid effects are thus important when interpreting the results of numerical relativity and in the construction of waveform models for current and future GW detectors. 
    more » « less
  3. The measurement of orbital eccentricity in gravitational-wave (GW) signals will provide unique insights into the astrophysical origin of binary systems, while ignoring eccentricity in waveform models could introduce significant biases in parameter estimation and tests of general relativity. Upcoming LIGO-Virgo-KAGRA observing runs are expected to detect a subpopulation of eccentric signals, making it vital to develop accurate waveform models for eccentric orbits. Here, employing recent analytical results through the third post-Newtonian order, we develop v5: a new time-domain, effective-one-body, multipolar waveform model for eccentric binary black holes with spins aligned (or antialigned) with the orbital angular momentum. Besides the dominant (2, 2) mode, the model includes the (2, 1), (3, 3), (3, 2), (4, 4), and (4, 3) modes. We validate the model’s accuracy by computing its unfaithfulness against 99 (28 public and 71 private) eccentric numerical-relativity (NR) simulations, produced by the Simulating eXtreme Spacetimes Collaboration. Importantly, for NR waveforms with initial GW eccentricities below 0.5, the maximum (2, 2)-mode unfaithfulness across the total mass range 20–200M⊙ is consistently below or close to 1%, with a median value of ∼0.02% , reflecting an accuracy improvement of approximately an order of magnitude compared to the previous-generation v4 and the state-of-the-art esumalí eccentric model. In the quasi-circular-orbit limit, v5 is in excellent agreement with the highly accurate v5 model. The accuracy, robustness, and speed of v5 make it suitable for data analysis and astrophysical studies. We demonstrate this by performing a set of recovery studies of synthetic NR-signal injections, and parameter-estimation analyses of the events GW150914 and GW190521, which we find to have no eccentricity signatures. 
    more » « less
  4. Gravitational Wave (GW) signals from binary neutron star (BNS) mergers provide critical insights into the properties of matter under extreme conditions. Due to the scarcity of observational data, numerical relativity (NR) simulations are indispensable for exploring these phenomena, without replacing the need for observational confirmation. However, simulating BNS mergers is a formidable challenge, and ensuring the consistency, reliability or convergence, especially in the post-merger, remains a work in progress. In this paper we assess the performance of current BNS merger simulations by analyzing open-source GW waveforms from five leading NR codes –SACRA,BAM,THC,Whisky, andSpEC. We focus on the accuracy of these simulations and on the effect of the equation of state on waveform predictions. We first check if different codes give similar results for similar initial data, then apply two methods to calculate convergence and quantify discretization errors. Lastly, we perform a thorough investigation into the effect of tidal interactions on key frequencies in the GW spectrum. We introduce a novel quasi-universal relation for the transient post-merger time, enhancing our understanding of remnant dynamics in this region. This detailed analysis clarifies agreements and discrepancies between these leading NR codes, and highlights necessary improvements for the advanced accuracy requirements of future GW detectors. 
    more » « less
  5. null (Ed.)
    Abstract Among the potential milliHz gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range 105M⊙ ≲ M ≲ 107M⊙ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN), where torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio q ≡ M2/M1, the disc viscosity (α), and gas temperature (or equivalently Mach number; $$\mathcal {M}$$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters. The gas imprint is detectable by LISA for binaries embedded in AGN discs with surface densities above $$\Sigma _0\ge 10^{4-6} \rm \, g cm^{-2}$$, depending on q, α and $$\mathcal {M}$$. Deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction. 
    more » « less