Realistic consecutive galaxy mergers form eccentric PTA sources

TL;DR

This study uses high-resolution simulations of galaxy mergers to show that supermassive black hole binaries entering the PTA frequency band are highly eccentric, affecting gravitational wave background interpretations.

Contribution

It provides new insights into the dynamical evolution and high eccentricity of BHBs from successive galaxy mergers, emphasizing the importance of eccentricity in PTA gravitational wave signals.

Findings

Binaries enter PTA band with eccentricity > 0.85

Dynamical friction lasts a few hundred Myrs

Binary formation time is resolution independent

Abstract

Results from pulsar timing arrays (PTAs) show evidence of a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole binaries (BHBs). The observed spectrum shows a flattening at lower frequencies that can be explained by a population of eccentric BHBs. This study aims to determine the dynamical evolution and merger timescales of the most massive BHBs, which are potential sources of the GWB. We select successive galactic major mergers from the IllustrisTNG100-1 cosmological simulation and re-simulate them at high resolution with the N-body code Griffin, down to binary separations of the order of a parsec. Coalescence timescales are estimated using a semi-analytical model that incorporates gravitational wave emission and stellar hardening. Throughout our investigation, we consider the impact of prior mergers on the remnant galaxy in the form of core scouring and anisotropy, which can influence the subsequent formation and evolution of BHBs. We find that all the binaries in our sample enter the PTA band with an eccentricity e>0.85: such a large eccentricity can impact the shape of the PTA observed GWB spectrum, and it highlights the importance of including the eccentricity of binaries when interpreting the PTA signal. Furthermore, we find that: (i) starting from initial separations of a few tens of kpcs, the dynamical friction phase lasts for a few hundred Myrs; (ii) the binary formation time is not resolution dependent; (iii) the scatter on the eccentricity at binary formation decreases with increasing resolution; (iv) triple systems form whenever a third galaxy interacts with a binary which hasn't yet reached coalescence.