Towards long and accurate numerical relativity waveforms of binary black holes beyond general relativity

Abstract

Numerical relativity (NR) simulations of compact binaries in theories beyond general relativity (GR) will be pivotal for the continued development of future tests of gravity with gravitational waves (GWs). In this Letter, we show that the combination of spectral methods and the "fixing-the-equations" approach allows us to produce the longest waveforms in the literature for a genuine beyond-GR theory, thus bringing NR methods for alternative theories of gravity closer to the state-of-the-art in GR. For concreteness, we focus on the well-known shift-symmetric version of scalar Gauss-Bonnet gravity, a theory postulating the existence of an additional dynamical scalar and describing black holes (BHs) different from the Kerr solution. We extract the gravitational and scalar waveforms at future null infinity for equal-mass, nonspinning, eccentricity-reduced BH binaries, and quantify the phase errors to be $\lesssim$ 1 rad after 40+ GW cycles (20+ orbits). We also show that the GW phase corrections in this alternative theory are distinguishable from Einstein's theory and lead to an earlier coalescence time than in GR. Obtaining such waveforms is a stepping stone to perform precise comparisons with Post-Newtonian theory and to calibrate waveform models beyond GR.