Eccentricity estimation from initial data for Numerical Relativity Simulations

TL;DR

This paper introduces an instantaneous eccentricity measure for initial data in binary black hole simulations, validated against post-Newtonian approximations and numerical evolutions, aiding in predicting and characterizing eccentricity.

Contribution

It presents a new method to estimate eccentricity at the start of numerical relativity simulations, incorporating spin effects and validated against full numerical evolutions.

Findings

Excellent agreement between 3PN estimates and full orbit calculations.

The method effectively predicts eccentricity in spinning binary black hole simulations.

Initial eccentricity estimates correlate well with evolved eccentricity during simulations.

Abstract

We describe and study an instantaneous definition of eccentricity to be applied at the initial moment of full numerical simulations of binary black holes. The method consists of evaluating the eccentricity at the moment of maximum separation of the binary. We estimate it using up to third post-Newtonian (3PN) order, and compare these results with those of evolving (conservative) 3PN equations of motion for a full orbit and compute the eccentricity $e_r$ from the radial turning points, finding excellent agreement. We next include terms with spins up to 3.5PN, and then compare this method with the corresponding estimates of the eccentricity $e_r^{NR}$ during full numerical evolutions of spinning binary black holes, characterized invariantly by a fractional factor $0\leq f\leq1$ of the initial tangential momenta. It is found that our initial instantaneous definition is a very useful tool to predict and characterize even highly eccentric full numerical simulations.