Impact of eccentricity and mean anomaly in numerical relativity mergers

Abstract: Accurate modelling of black hole binaries is critical to achieve the science goals of gravitational-wave detectors. Modelling such configurations relies strongly on calibration to numerical-relativity (NR) simulations. Binaries on quasi-circular orbits have been widely explored in NR, however, coverage of the broader 9-dimensional parameter space, including orbital eccentricity, remains sparse. This article develops a new procedure to control orbital eccentricity of binary black hole simulations that enables choosing initial data parameters with precise control over eccentricity and mean anomaly of the subsequent evolution, as well as the coalescence time. We then calculate several sequences of NR simulations that nearly uniformly cover the 2-dimensional eccentricity--mean anomaly space for equal mass, non-spinning binary black holes. We demonstrate that, for fixed eccentricity, many quantities related to the merger dynamics of binary black holes show an oscillatory dependence on mean anomaly. The amplitude of these oscillations scales nearly linearly with the eccentricity of the system. We find that for the eccentricities explored in this work, deviations in various quantities such as the merger amplitude and peak luminosity can approach $\sim5\%$ of their quasi-circular value. We use our findings to explain eccentric phenomena reported in other studies. We also show that methods for estimating the remnant mass employed in the effective-one-body approach exhibit similar deviations, roughly matching the amplitude of the oscillations we find in NR simulations. This work is an important step towards a complete description of eccentric binary black hole mergers, and demonstrates the importance of considering the entire 2-dimensional parameter subspace related to eccentricity.