Hanging on the cliff: Extreme mass ratio inspiral formation with local two-body relaxation and post-Newtonian dynamics

TL;DR

This paper introduces a novel Monte Carlo method incorporating post-Newtonian dynamics and local two-body relaxation to better predict the formation of extreme mass ratio inspirals (EMRIs) around massive black holes, revealing the significance of cliffhangers especially for low-mass black holes.

Contribution

The study develops a new simulation approach that accounts for local relaxation and post-Newtonian effects, challenging previous assumptions and highlighting the importance of cliffhangers in EMRI formation.

Findings

Cliffhangers are more common than previously thought, especially for MBHs below 3×10^5 solar masses.

The post-Newtonian description significantly increases the predicted number of EMRIs.

Local relaxation treatment boosts cliffhanger formation for small MBHs.

Abstract

Extreme mass ratio inspirals (EMRIs) are anticipated to be primary gravitational wave sources for LISA (Laser Interferometer Space Antenna). They form in dense nuclear clusters when a compact object (CO) is captured by the central massive black holes (MBHs) due to frequent two-body interactions among orbiting objects. We present a novel Monte Carlo approach to evolve the post-Newtonian (PN) equations of motion of a CO orbiting an MBH accounting for two-body relaxation locally on the fly, without the assumption of orbit-averaging. We estimate the fraction $S(a_0)$ of EMRIs to total captures (including direct plunges, DPs) as a function of the initial semi-major axis $a_0$ for COs around MBHs of $M_\bullet\in[10^4\,{\rm M}_\odot,4\times10^6\,{\rm M}_\odot]$. Previous results indicate $S(a_0)\rightarrow 0$ at large $a_0$, with a sharp transition from EMRIs to DPs around a critical scale $a_{\rm c}$. This notion has been recently challenged for low-mass MBHs, with EMRIs forming at $a\gg a_{\rm c}$, the so-called "cliffhangers''. Our simulations confirm their existence, at larger numbers than previously expected. Cliffhangers start to appear for $M_\bullet\lesssim3\times 10^5\,{\rm M}_\odot$ and can account for up to 55% of the overall EMRIs formed. We find $S(a_0)\gg 0$ for $a\gg a_{\rm c}$, reaching values as high as 0.6 for $M_\bullet=10^4\,{\rm M}_\odot$, much larger than previously found. We find that the PN description of the system greatly enhances the number of EMRIs by shifting $a_{\rm c}$ to larger values at all MBH masses, and that the local treatment of relaxation significantly boosts the number of cliffhangers for small MBHs. Our work shows the limitations of standard assumptions for estimating EMRI formation rates, most importantly their dynamical models. Future estimates of rates and properties of EMRIs detectable by LISA should account for these improvements.