Search for: All records

The asymptotic nature of unbound binary-black-hole encounters provides a clean method for comparing different approaches for modeling the two-body problem in general relativity. In this work, we use numerical relativity simulations of black-hole scattering, generated using the Spectral Einstein Code, to explore the self-force and post-Minkowskian expansions of the scattering angle. First, we use a set of unequal-mass simulations to extract the self-force contributions to the scattering angle. Our main result is that using information up to second-order in the symmetric mass ratio (2SF) reproduces numerical relativity within the error bars across the full range of mass-ratios, including equal mass. Next, we compare our numerical relativity results to state-of-the-art post-Minkowskian predictions at larger impact parameters than previously explored. We find good agreement in the weak-field regime and discuss the relative importance of higher order terms. 

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 ( ${\chi }_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. 

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (FIL, SpEC), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. We also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors. 

Mounting evidence indicates that some of the gravitational wave signals observed by ground-based observatories might arise from eccentric compact object binaries, increasing the urgency for accurate waveform models for such systems. While for noneccentric binaries, surrogate models are efficient and accurate, the additional features due to eccentricity have posed a challenge. In this work, we present a method for decomposing eccentric numerical relativity waveforms, making them amenable to surrogate modeling techniques. We parametrize the inspiral in the radial-phase domain, factoring out eccentricity-induced dephasing and thus enhancing compressibility and accuracy. This is combined with a second surrogate for the merger ringdown in the time domain and a technique to take advantage of the approximate periodicity with radial oscillations during the inspiral. We apply this procedure to the (2, 2) mode for nonspinning black hole binaries and demonstrate that the resulting surrogate, , is able to faithfully reproduce the underlying numerical relativity waveforms, with maximum mismatches of $6\times {10}^{-4}$and median mismatches of $4\times {10}^{-5}$. This technique paves the way for high-accuracy parameter estimation with eccentric models, a key ingredient for astrophysical inference and tests of general relativity. 

Abstract When numerically solving Einstein’s equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ‘junk radiation.’ This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produced. Thus, an automated procedure for identifying junk radiation is required. To this end, we present a new algorithm based on the empirical mode decomposition (EMD) from the Hilbert–Huang transform. This approach allows us to isolate and measure the high-frequency oscillations present in the measured irreducible masses of the black holes. The decay of these oscillations allows us to estimate the time from which the junk radiation can be ignored. To make this procedure more precise, we propose three distinct threshold criteria that specify how small the contribution of junk radiation has to be before it can be considered negligible. We apply this algorithm to 3403 BBH simulations from the Simulating eXtreme Spacetime catalog to find appropriate values for the thresholds in the three criteria. We find that this approach yields reliable decay time estimates, i.e. when to consider the simulation physical, for $>$98.5% of the simulations studied. This demonstrates the efficacy of the EMD as a suitable tool to automatically isolate and characterize junk radiation in the simulation of binary black hole systems. 

Abstract The Simulating eXtreme Spacetimes Collaboration’s codeSpECcan now routinely simulate binary black hole mergers undergoing $\sim 25$orbits, with the longest simulations undergoing nearly $\sim 180$orbits. While this sounds impressive, the mismatch between the highest resolutions for this long simulation is $O\left({10}^{-1}\right)$. Meanwhile, the mismatch between resolutions for the more typical simulations tends to be $O\left({10}^{-4}\right)$, despite the resolutions being similar to the long simulations’. In this note, we explain why mismatch alone gives an incomplete picture of code—and waveform—quality, especially in the context of providing waveform templates for LISA and 3G detectors, which require templates with $O\left({10}^3\right)-O\left({10}^5\right)$orbits. We argue that to ready the GW community for the sensitivity of future detectors, numerical relativity groups must be aware of this caveat, and also run future simulations with at least three resolutions to properly assess waveform accuracy.