The inspiral and merger of eccentric binaries leads to gravitational waveforms distinct from those generated by circularly merging binaries. Dynamical environments can assemble binaries with high eccentricity and peak frequencies within the LIGO band. In this paper, we study binary-single stellar scatterings occurring in dense stellar systems as a source of eccentrically-inspiraling binaries. Many interactions between compact binaries and single objects are characterized by chaotic resonances in which the binary-single system undergoes many exchanges before reaching a final state. During these chaotic resonances, a pair of objects has a non-negligible probability of experiencing a very close passage. Significant orbital energy and angular momentum are carried away from the system by gravitational wave (GW) radiation in these close passages and in some cases this implies an inspiral time shorter than the orbital period of the bound third body. We derive the cross section for such dynamical inspiral outcomes through analytical arguments and through numerical scattering experiments including GW losses. We show that the cross section for dynamical inspirals grows with increasing target binary semi-major axis, a, and that for equal-mass binaries it scales as a 2/7 . Thus, we expect wide target binaries to predominantly contribute to the production of these relativistic outcomes. We estimate that eccentric inspirals account for approximately one percent of dynamically assembled non-eccentric merging binaries. While these events are rare, we show that binary-single scatterings are a more effective formation channel than single-single captures for the production of eccentrically-inspiraling binaries, even given modest binary fractions.

We derive the probability for a newly formed binary black hole (BBH) to undergo an eccentric gravitational wave (GW) merger during binary-single interactions inside a stellar cluster. By integrating over the hardening interactions such a BBH must undergo before ejection, we find that the observable rate of BBH mergers with eccentricity > 0.1 at 10 Hz relative to the rate of circular mergers can be as high as ∼ 5% for a typical globular cluster (GC). This further suggests that BBH mergers forming through GW captures in binary-single interactions, eccentric or not, are likely to constitute ∼ 10% of the total BBH merger rate from GCs. Such GW capture mergers can only be probed with an N -body code that includes General Relativistic corrections, which explains why recent Newtonian cluster studies not have been able to resolve this population. Finally, we show that the relative rate of eccentric BBH mergers depends on the compactness of their host cluster, suggesting that an observed eccentricity distribution can be used to probe the origin of BBH mergers.

Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc −3 yr −1 to the binary merger rate in the local universe, increasing to 10 Gpc −3 yr −1 at z ∼ 3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.

We present the first systematic study of strong binary-single and binary-binary black hole (BH) interactions with the inclusion of general relativity. By including general relativistic effects in the equations of motion during strong encounters, the dissipation of orbital energy from the emission of gravitational waves (GWs) can lead to captures and subsequent inspirals with appreciable eccentricities when entering the sensitive frequency ranges of the LIGO and Virgo GW detectors. It has been shown that binary-single interactions significantly contribute to the rate of eccentric mergers, but no studies have looked exclusively into the contribution from binary-binary interactions. To this end, we perform binary-binary and binary-single scattering experiments with general relativistic dynamics up through the 2.5 post-Newtonian order included, both in a controlled setting to gauge the importance of non-dissipative post-Newtonian terms and derive scaling relations for the cross section of GW captures, as well as experiments tuned to the strong interactions from state-of-the art globular cluster (GC) models to assess the relative importance of the binary-binary channel in facilitating GW captures and the resultant eccentricity distributions of inspiral from channel. Although binary-binary interactions are 10-100 times less frequent in GCs than binary-single interactions, their longer lifetime and more complex dynamics leads to a higher probability for GW captures to occur during the encounter. We find that binarybinary interactions contribute 25-45% of the eccentric mergers that occur during strong BH encounters in GCs, regardless of the properties of the cluster environment. The inclusion of higher multiplicity encounters in dense star clusters therefore have major implications on the predicted rates of highly eccentric binaries potentially detectable by the LIGO/Virgo network. Becse gravitational waveforms of eccentric inspirals are distinct from those generated by merging binaries that have circularized, measurements of eccentricity in such systems would highly constrain their formation scenario.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters N GW needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for N GW to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically N GW 10.We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation-and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.