Black Hole and Neutron Star Binary Mergers in Triple Systems: Merger Fraction and Spin-Orbit Misalignment

Black hole (BH) mergers driven by gravitational perturbations of external companions constitute an important class of formation channels for merging BH binaries detected by LIGO. We have studied the orbital and spin evolution of binary BHs in triple systems, where the tertiary companion excites large eccentricity in the inner binary through Lidov-Kozai oscillations, causing the binary to merge via gravitational radiation. Using the single-averaged and double-averaged secular dynamics of triples (where the equations of motion are averaged over the inner orbit and both orbits, respectively), we perform a large set of numerical integrations to determine the merger window (the range of companion inclinations that allows the inner binary to merge within similar to 10 Gyr) and the merger fraction as a function of various system parameters (e.g., the binary masses m(1), m(2) and initial semimajor axis a(0), the mass, semimajor axis, and eccentricity e(out) of the outer companion). For typical BH binaries (m(1,2) similar or equal to 20 M-circle dot-30 M-circle dot and a(0) greater than or similar to 10 au), the merger fraction increases rapidly with e(out) because of the octupole perturbation, ranging from similar to 1% at e(out) = 0 to 10%-20% at e(out) = 0.9. We derive analytical expressions and approximate scaling relations for the merger window and merger fraction for systems with negligible octupole effect, and apply them to neutron star binary mergers in triples. We also follow the spin evolution of the BHs during the companion-induced orbital decay, where de Sitter spin precession competes with Lidov-Kozai orbital precession/nutation. Starting from aligned spin axes (relative to the orbital angular momentum axis), a wide range of final spin-orbit misalignment angle theta(f)(sl) can be generated when the binary enters the LIGO sensitivity band. For systems where the octupole effect is small (such as those with m(1) similar or equal to m(2) or e(out) similar to 0), the distribution of theta(f)(sl) peaks around 90 degrees. As the octupole effect increases, a more isotropic distribution of final spin axis is produced. Overall, merging BH binaries produced by Lidov-Kozai oscillations in triples exhibit a unique distribution of the effective (mass-weighted) spin parameter chi(eff); this may be used to distinguish this formation channel from other dynamical channels.