Evolution of precessing binary black holes on eccentric orbits using orbit-averaged evolution equations

Abstract: The most general bound binary black hole (BBH) system has an eccentric orbit and precessing spins. The detection of such a system with significant eccentricity close to the merger would be a clear signature of dynamical formation. In order to study such systems, it is important to be able to evolve their spins and eccentricity from the larger separations at which the binary formed to the smaller separations at which it is detected, or vice versa. Knowledge of the precessional evolution of the binary's orbital angular momentum can also be used to twist up aligned-spin eccentric waveform models to create a spin-precessing eccentric waveform model. In this paper, we present a new publicly available code to evolve eccentric, precessing BBHs using orbit-averaged post-Newtonian (PN) equations from the literature. The spin-precession dynamics is 2PN accurate, i.e., with the leading spin-orbit and spin-spin corrections. The evolution of orbital parameters (orbital frequency, eccentricity, and periastron precession), which follow the quasi-Keplerian parametrization, is 3PN accurate in the point particle terms and includes the leading order spin-orbit and spin-spin effects. All the spin-spin terms include the quadrupole-monopole interaction. The eccentricity enhancement functions in the fluxes use the high-accuracy hyperasymptotic expansions from Loutrel and Yunes [Classical Quantum Gravity {\bf 34} 044003 (2017)]. We discuss various features of the code and study the evolution of the orbital and spin-precession parameters of eccentric, precessing BBHs. In particular, we study the dependence of the spin morphologies on eccentricity, where we find that the transition point from one spin morphology to another can depend nonmonotonically on eccentricity, and the fraction of binaries in a given morphology at a given point in the evolution of a population depends on the instantaneous eccentricity.