Drat version September 5, 217 Preprint typeset using L A TEX style emulateapj v. 12/16/11 SPIN-ORBIT MISALIGNMENT OF MERGING BLACK-HOLE BINARIES WITH TERTIARY COMPANIONS Bin Liu 1,2 and Dong Lai 2 1 Shanghai Astronomical Observatory, Chinese Academy o Sciences, 8 Nandan Road, Shanghai 2, China and 2 Cornell Center or Astrophysics and Planetary Science, Cornell University, Ithaca, NY 1485, USA Drat version September 5, 217 arxiv:176.29v2 [astro-ph.he] Sep 217 ABSTRACT We study the eect o external companion on the orbital and spin evolution o merging blackhole (BH) binaries. An suiciently close by and inclined companion can excite Lidov-Kozai (LK) eccentricity oscillations in the binary, thereby shortening its merger time. During such LK-enhanced orbital decay, the spin axis o the BH generally exhibits chaotic evolution, leading to a wide range ( -18 ) o inal spin-orbit misalignment angle rom an initially aligned coniguration. For systems that do not experience eccentricity excitation, only modest ( 2 ) spin-orbit misalignment can be produced, and we derive an analytic expression or the inal misalignment using the principle o adiabatic invariance. The spin-orbit misalignment directly impacts the gravitational waveorm, and can be used to constrain the ormation scenarios o BH binaries and dynamical inluences o external companions. 1. INTRODUCTION The recent breakthrough in the detection o gravitational waves (GWs) rom merging black hole (BH) binaries by advanced LIGO (Abbott et al. 216a,b, 217) has generated renewed interest in understanding the ormation mechanisms o compact BH binaries, rom the evolution o massive stellar binaries (Lipunov et al. 1997; Belczynski et al. 21, 216; Mandel & de Mink 216; Lipunov et al. 217) and triples (Silsbee & Tremaine 217; Antonini et al. 217) in the galactic ields, to dynamical interactions in galactic nuclei (Antonini & Perets 212; Petrovich & Antonini 217) and in the dense core o globular clusters (Miller & Hamilton 22; Rodriguez et al. 215; Chatterjee et al. 217). Because o the uncertainties associated with various ormation channels (e.g., common envelope evolution in the standard binary channel), it is diicult to distinguish the dierent ormation mechanisms based on BH mass measurement alone. The detection o eccentric systems would obviously indicates some dynamical processes at work (e.g., Antonini & Perets (212); Silsbee & Tremaine (217)). However, because o the eicient eccentricity damping by gravitational radiation, the vast majority o compact binaries will likely be circular when entering the LIGO sensitivity band regardless o the ormation channels. It has been suggested that the BH spin and the spin-orbit misalignment angle may be an important discriminant. In particular, the spin-orbit misalignment directly impacts the projected spin parameter o the merging binaries, χ e = m 1a 1 + m 2 a 2 m 1 + m 2 ˆL, (1) (where m 1,2 are the BH masses, a 1,2 = cs 1,2 /(Gm 2 1,2) are the dimensionless BH spins, and ˆL is the unit orbital angular momentum vector), which can be measured rom the phase evolution o GWs (Abbott et al. 216b, 217). In this paper, we study the merger and spin-orbit misalignment o BH binaries in the presence o tertiary companion. Such triple BH systems could be a direct product o massive triple stars in the galactic ield (Silsbee & Tremaine 217; Antonini et al. 217), or could be produced dynamically in a dense cluster (Miller & Hamilton 22; Rodriguez et al. 215; Antonini & Rasio 216). For binaries ormed near the center o a galaxy, the third body could be a supermassive BH (Antonini & Perets 212; Petrovich & Antonini 217). It is well known that a tertiary body on an inclined orbit can accelerate the decay o an inner binary by inducing Lidov-Kozai (LK) eccentrcity/inclination oscillations (Lidov 1962; Kozai 1962). This has been studied beore in the contexts o supermassive BH binary merger (Blaes et al. 22) and stellar mass BH binaries (e.g., Miller & Hamilton (22); Thompson (211); Antonini et al. (214); Silsbee & Tremaine (217)). We ocus on the latter in this paper. We show that as the BH binary undergoes LK-enhanced decay rom a wide orbit and eventually enters the LIGO band, the spin axis o individual BHs can experience chaotic evolution, so that a signiicant spin-orbit misalignment can be produced prior to merger even or binaries ormed with zero initial misalignment. We derive relevant analytic relations and quantiy how the inal spin-orbit misalignment angle depends on various parameters o the system (binary and external companion). Note that in this paper we ocus on triple systems with relatively small binary separations ( AU or the inner binaries), so that the inner binary can merger within 1 1 years either by itsel or through modest (e.99) LK eccentricity excitation. Such compact triple systems likely have gone through a complex (and highly uncertain) sequence o common envelop or mass transer evolution we do not study such evolution in this paper and thereore do not address issues related to the occurrence rate o compact triples. Our goal is to use such triple systems to illustrate the complex spin dynamics o the individual BHs. We expect that similar spin dynamics may take place in other types o triple systems, e.g., those with much larger initial separations, but experience extreme eccentricity excitation due to non-secular orcing rom tertiary companions (Silsbee & Tremaine 217;
2 Antonini et al. 217). 2. LIDOV-KOZAI CYCLES IN BH TRIPLES WITH GRAVITATIONAL RADIATION 2.1. Setup and Orbital Evolution Consider a hierarchical triple system, consisting o an inner BH binary with masses m 1, m 2 and a relatively distant companion o mass m. The reduced mass or the inner binary is µ in m 1 m 2 /m 12, with m 12 m 1 + m 2. Similarly, the outer binary has µ out (m 12 m )/m 12 with m 12 m 12 + m. The orbital semimajor axes and eccentricities are denoted by a in,out and e in,out, respectively. The orbital angular momenta o the inner and outer binaries are L in = L in ˆLin = µ in Gm 12 a in (1 e 2 in ) ˆL in, (2) L out = L out ˆLout = µ out Gm 12 a out (1 e 2 out) ˆL out,() where ˆL in,out are unit vectors. The relative inclination between ˆL in and ˆL out is denoted by I. For convenience, we will requently omit the subscript in. The merger time due to GW radiation o an isolated binary with initial a and e = is given by T m, = 5c 5 a 4 256G m 2 12 µ 1 1 ( 6M m 12 ) 2 ( )( 15M µ a 2AU ) 4 yrs. (4) A suiciently inclined external companion can raise the binary eccentricity through Lidov-Kozai oscillations, thereby reducing the merger time or making an otherwise non-merging binary merge within 1 1 years. To study the evolution o merging BH binary under the inluence o a companion, we use the secular equations to the octupole level in terms o the angular momentum L and eccentricity e vectors: dl dt = dl dt + dl LK dt, (5) GW de dt = de dt + de LK dt + de GR dt. (6) GW Here the Lidov-Kozai (LK) terms are given explicitly in (Liu et al. 215a) (we also evolve L out and e out ), and the associated timescale o LK oscillation is t LK = 1 n m 12 m ( aout,e a ), (7) where n = Gm 12 /a is the mean motion o the inner binary and a out,e a out 1 e 2 out. General Relativity (1-PN correction) induces pericenter precession de dt = Ω GR ˆL e, ΩGR = Gnm 12 GR c 2 a(1 e 2 ). (8) We include GW emission (2.5-PN eect) that causes orbital decay and circularization (Shapiro & Teukolsky 198), but not the extreme eccentricity excitation due to non-secular eects (Antonini et al. 214; Silsbee & Tremaine 217). aau 1.1 1 1 1..6.4 9 8 7 6 e IDeg ΘSLDeg 15 1 5 I 81,Θ SL 1 I 81,Θ SL I 8,Θ SL 1 1 1 5 1 7 tyrs Fig. 1. Sample orbital and spin evolution o a BH binary system with a tertiary companion. The top three panels show the semimajor axis, eccentricity and inclination (relative to ˆL out) o the inner BH binary, and the bottom panel shows the spin-orbit misalignment (the angle between S 1 and L). The parameters are m 1 = m 2 = m = M, a out = AU, e out =, and the initial a =.1AU, e = 1, I = 8 and θsl =. For this example, the octupole eect is absent. In the bottom panel, the results or slightly dierent values o I and θsl (as indicated) are plotted, showing a strong dependence o the inal θ SL on the initial conditions. The top three panels o Figure 1 show an example o the orbital evolution o a BH binary with an inclined companion (initial I = 8 ). We see that the inner binary undergoes cyclic excursions to maximum eccentricity e max, with accompanying oscillations in the inclination I. As the binary decays, the range o eccentricity oscillations shrinks. Eventually the oscillations reeze and the binary experiences pure orbital decay/circularization governed by GW dissipation. 2.2. Eccentricity Excitation and Merger Time In the quadrupole approximation, the maximum eccentricity e max attained in the LK oscillations (starting rom an initial I and e ) can be calculated analytically using energy and angular momentum conservation, according to the equation (Anderson et al. 217a) 8 j 2 min 1 j 2 min +η 2 j 4 min [ ( 5 cos I + η ) 2 ( + 4η cos I + 9 2 4 η2) jmin 2 ] + ε GR ( 1 j 1 min) =, (9)
elim I- I+ emax.6.4 m = 2 M, aout = H2 L1 AU m = M, aout = AU m = 2 M, aout = 2 H2 L m = M, aout = 2 AU 2 4 1 AU I @DegD 6 8 1 12 14 16 18 Fig. 2. The maximum eccentricity o the inner BH binary versus the initial inclination I o the tertiary companion, calculated using Equation (9). The inner binary has m1 = m2 = M, a =.1AU, and initial e '. The companion has a circular orbit and its mass and semimajor axis are as labeled. The emax (I ) curve depends mainly on m /aout. The horizontal (elim ) and vertical (I± ) lines are given by Equations (11) and (12), respectively. where jmin given by p 1 e2max, η (L/Lout )e= and εgr is Gm212 aout,e εgr = tlk ΩGR = c2 a4 m e= 2 1 4 m12 m aout,e a '.96.(1) 6M M AU.1AU Note that in the limit o η and εgr, Equation (9) yields the well-known relation emax = p 1 (5/) cos2 I. The maximum possible emax or all values o I, called elim, is given by 2 η2 4 2 1 (jlim 1) + jlim 1 + εgr 1 jlim =, 8 4 5 (11) 2 and is reached at cos I = (η/1)(4jlim 5). Eccentricity excitation (emax ) occurs within a window o inclinations (cos I ) 6 cos I 6 (cos I )+, where (Anderson et al. 217a) r 1 8 (cos I )± = η ± η 2 + 6 εgr. (12) 1 This window vanishes when 9 εgr + η 2 (no eccentricity excitation). (1) 4 8 Figure 2 shows some examples o the emax (I ) curves. For η. 1, these curves depend mainly on m /aout,e (or given inner binary parameters). As εgr increases (with decreasing m /aout,e ), the LK window shrinks and elim decreases. Eccentricity excitation is suppressed (or all I s) when Equation (1) is satisied. For systems with m1 6= m2 and eout 6=, so that m1 m2 a eout εoct (14) m12 aout 1 e2out TmTm,.6.4 elim Α 2 emax.1 Α 1 Α 2. 5 1.5 1 1 4 5.9.95 1. Tm Tm, 1. 1..4.6 emax 1. Fig.. The binary merger time Tm (in units o Tm, ; Equation 4) as a unction o the maximum eccentricity induced by a tertiary companion. The BH masses are m1 = m2 = m = M, and the inner and outer orbits are initially circular. We consider two sets o initial semimajor axes: a =.1AU, aout = AU (orange) and a = AU, aout = 5AU (green). Here Tm is computed numerically by integrating the secular evolution equations with dierent I s, and emax is calculated using Equation (9). The decreasing trend o Tm /Tm, as a unction o emax can be approximated by (1 e2max )α (dashed lines). The inset shows a zoomed-in portion o large emax s. is non-negligible, the octupole eect may become important (e.g. Ford et al. (2); Naoz (216)). This tends to widen the inclination window or large eccentricity excitation. However, the analytic expression or elim given by Equation (11) remains valid even or εoct 6= (Liu et al. 215a; Anderson & Lai 217b). Thereore Equation (1) still provides a good criterion or negligible eccentricity excitation (note that when εoct 6=, the inner binary cannot be exactly circular; e.g., Anderson & Lai (217b); Liu et al. (215b)). Eccentricity excitation leads to a shorter binary merger time Tm compared to the circular merger time Tm, (Figure ). We compute Tm by integrating the secular evolution equations o the BH triples with a range o I. Each I run has a corresponding emax, which can be calculated using Equation (9). We see that in general Tm /Tm, can be approximated by (1 e2max )α, with α ' 1.5, 2 and 2.5 or emax = (,.6), (.6, ) and (,.95), respectively.. SPIN-ORBIT DYNAMICS IN MERGING BH BINARY WITH AN EXTERNAL COMPANION.1. Spin-Orbit Coupling We now study how the BH spin evolves during the ˆ 1 (where S ˆ 1 is binary merger, considering only S1 = S1 S ˆ the unit vector). The de Sitter precession o S1 around ˆ is (Barker & O Connell 1975) L ˆ1 ds ˆ S ˆ 1, = ΩdS L dt ΩdS = Gn(m2 + µ/). 2c2 a(1 e2 ) (15) Note that there is a back-reaction torque rom S1 on L; this can be saely neglected since L S1. Beore presenting our numerical results, it is useul to note the dierent regimes or the evolution o the spinorbit misalignment angle θsl (the angle between S1 and ˆ precesses (and L). In general, the inner binary axis L nutates when e 6= ) around the total angular momentum J = L + Lout (recall that S1, S2 L, and J is constant in the absence o GW dissipation). The related
4 precession rate Ω L is o order t 1 LK at e (Equation 17), but increases with e. Depending on the ratio o Ω ds and Ω L, we expect three possible spin behaviors: (i) For Ω L Ω ds ( nonadiabatic ), the spin axis Ŝ1 cannot keep up with the rapidly changing ˆL, and thus eectively precesses around Ĵ, keeping θ SJ cos 1 (Ŝ1 Ĵ) constant. (ii) For Ω ds Ω L ( adiabatic ), Ŝ1 closely ollows ˆL, maintaining an approximately constant θ SL. (iii) For Ω ds Ω L ( trans-adiabatic ), the spin evolution can be chaotic due to overlapping resonances. Since both Ω ds and Ω L depend on e during the LK cycles, the precise transitions between these regimes can be uzzy (Storch et al. 214; Storch & Lai 215; Anderson et al. 216; Storch et al. 217). For circular orbits (e = ), the precession o ˆL is governed by the equation dˆl dt = Ω L ˆLout ˆL = Ω LĴ ˆL, (16) LK,e= where Ĵ is the unit vector along J = L + L out, and Ω L = 4t LK (ˆL ˆLout ), Ω L = Ω L J L out. (17) We can deine an adiabaticity parameter ( ) [ ]( ) ΩdS (m2 + µ/) m12 A.7 Ω L e,i= 5M 6M ( ) 1 ( ) ( ) 4 m aout,e a. M AU.1AU (18) As the binary orbit decays, the system may transition rom non-adiabatic (A 1) at large a s to adiabatic (A 1) at small a s, where the inal spin-orbit misalignment angle θsl is rozen. Note that A is directly related to ε GR by A = 2 m 2 + µ/. (19) ε GR m 12 Thus, when the initial value o ε GR (at a = a ) satisies ε GR, 9/4 (a necessary condition or LK eccentricity excitation; see Equation 1), we also have A (m 2 + µ)/(2m 12 ) 1. This implies that any system that experiences enhanced orbital decay due to LK oscillations must go through the trans-adiabatic regime and thereore possibly chaotic spin evolution. The bottom panel o Figure 1 shows a representative example o the evolution o the misalignment angle as the BH binary undergoes LK-enhanced orbital decay. We see that the BH spin axis can exhibit complex evolution even though the orbital evolution is regular. In particular, θ SL evolves in a chaotic way, with the inal value θsl depending sensitively on the initial conditions (the precise initial θ SL and I ). Also note that retrograde spin (θsl > 9 ) can be produced even though the binary always remains prograde with respect to the outer companion (I < 9 ). These behaviors are qualitatively similar to the chaotic evolution o stellar spin driven by Newtonian spin-orbit coupling with a giant planet undergoing high-eccentricity migration (Storch et al. 214; Storch & Lai 215; Anderson et al. 216; Storch et al. 217). We carry out a series o numerical integrations, evolving the orbit o the merging BH binary with a tertiary companion, along with spin-orbit coupling, to determine θsl or various triple parameters. In our population synthesis study, we consider initial conditions such that Ŝ 1 is parallel to ˆL, the binary inclinations are isotropically distributed (uniorm distribution in cos I ), and the orientations o e and e out (or systems with e out ) are random. All initial systems satisy the criterion o dynamical stability or triples (Mardling & Aarseth 21). Figure 4 shows our results or systems with equal masses, e out = (so that the octupole eect vanishes), and several values o a out. We see that when the eccentricity o the inner binary is excited (I lies in the LK window), a wide range o θsl is generated, including appreciable raction o retrograde (θsl > 9 ) systems (see the a out = AU case, or which e lim = 4). The memory o chaotic spin evolution is evident, as slightly dierent initial inclinations lead to vastly dierent θsl. The regular behavior o θsl around I = 9 (again or the a out = AU case) is intriguing, but may be understood using the theory developed in Re. (Storch et al. 217). For systems with no eccentricity excitation, θsl varies regularly as a unction o I this can be calculated analytically (see below). Figure 5 shows our results or systems with m 1 m 2 and e out, or which octupole terms may aect the orbital evolution. We see that eccentricity excitation and the corresponding reduction in T m occur outside the analytic (quadupole) LK window (see the e out = case, or which e lim =.66). As in the equal-mass case (Figure 4), a wide range o θsl values are produced whenever eccentricity excitation occurs. A larger raction (2%) o systems attain retrograde spin (θsl > 9 ). Again, or systems with negligible eccentricity excitation, θsl behaves regularly as a unction o I and agrees with the analytic result (the uzziness o the numerical result in this regime is likely due to the very small eccentricity o the inner binary; see Liu et al. (215b))..2. Analytical Calculation o θsl or Circular Binaries I the inner binary experiences no eccentricity excitation and remains circular throughout the orbital decay, the inal spin-orbit misalignment can be calculated analytically using the principle o adiabatic invariance. Equation (16) shows that ˆL rotates around the Ĵ axis at the rate ( Ω L ). In this rotating rame, the spin evolution equation (15) transorms to ( ) d Ŝ 1 = Ω e dt Ŝ1, Ω e Ω ds ˆL + Ω L Ĵ. (2) rot Note that in the absence o GW dissipation, ˆL and ˆL out are constants (in the rotating rame), and thus Ŝ1 precesses with a constant Ω e. The relative inclination between Ω e and ˆL is given by tan θ e,l = Ω L sin I Ω ds + ( η + cos I ) Ω L. (21)
5 Deg I Deg 18 12 9 6 15 1 m 1 M, m 2 M m M a.1 AU e 1 a out AU, A.7 a out 5AU, A 1.7 a out 1 AU, A 1.82 Deg I Deg 18 12 9 6 15 1 m 1 M m 2 1 M m M a 5 AU a out AU e out, A.1 e out.4, A 1.1 e out 1, A 1.42 5 5 Tm 1 8 yr 5 1 5.5 1 I I.1 1..5.5 1. cosi Tm 1 8 yr 5 1.5 I I.1 1..5.5 1. cosi N Nmax. 5.15.1 5 ΘSL 9 1.7 6 9 12 15 18 Deg N Nmax. 5.15.1 5 ΘSL 9 2.2 6 9 12 15 18 Deg Fig. 4. Final spin-orbit misalignment angle (top panel) and merger time (middle panel) as a unction o the initial inclination or equal mass triple systems with dierent outer semimajor axes (as labeled). The inner binary has ixed initial a =.1 AU and e = 1, and e out = or the outer binary. In the top panel, the dots are the result o numerical integration or the a out = AU system (a total o 4 runs on a uniorm cos I grid), and the dashed curves are the analytical results or circular orbits, as given by Equations (2) and (24). The initial value o the adiabaticity parameter A (Equation 18 with a = a ) is also given. The vertical lines (I ± ) shown in the middle panel correspond to the LK window o eccentricity excitation (Equation 12). The bottom panel shows the distribution o the inal spin-orbit misalignment angle or the system with a out = AU. Now i we include GW dissipation, ˆL ˆLout = cos I is exactly conserved, and Ω e becomes a slowly changing vector. When the rate o change o Ω e is much smaller than Ω e, the angle between Ω e and Ŝ1 is adiabatic invariant, i.e. θ e,s1 constant (adiabatic invariant). (22) This adiabatic invariance requires dω e /dt / Ω e Tm, 1 Ω e, or Ω e T m, 1, which is easily satisied. Suppose Ŝ1 and ˆL are aligned initially, we have θe,s 1 = θe,l (the superscript denotes initial value). Equation (22) then implies θ e,s1 θe,l at all times. Ater the binary has decayed, η, Ω ds Ω L, and thereore Ω e Ω ds ˆL, which implies θ SL θ e,s1. Thus we ind θsl θe,l. (2) Fig. 5. Same as Figure4, but or unequal masses and noncircular outer orbits (as labeled). The dots in the top panel are numerical results or the e out = system (a total o 4 runs on a uniorm cos I grid), and the bottom panel shows the distribution o θsl or such a system. That is, the inal spin-orbit misalignment angle is equal to the initial inclination angle between ˆL and Ω e, obtained by evaluating Equation (21) at a = a : tan θ e,l = sin I (A / cos I ) + η + cos I. (24) This analytic expression agrees with the numerical results shown in Figures 4-5 in the appropriate regime. Note that or systems that experience no eccentricity excitation or all I s, A (m 2 + µ)/(2m 12 ) 1 (see Equations 1 and 19), and thus θ SL θ e,l 2, i.e., only modest spin-orbit misalignment can be generated. For systems with A 1 (e.g., very distant companion), we have θ SL 1. The above analytical result can be easily generalized to the situation o non-zero initial spin-orbit misalignment. It shows that θ SL θ SL or A 1. 4. SUMMARY AND DISCUSSION We have studied the eect o external companion on the orbital and spin evolution o merging BH binaries due to gravitational radiation. A suiciently close by and inclined companion can excite Lidov-Kozai eccentricity oscillation in the binary, shortening its merger time compared to circular orbits [see Fig.]. We ind
6 that during the LK-enhanced orbital decay, the spin axis o the BH generally experiences complex, chaotic evolution, with the inal spin-orbit misalignment angle θsl depending sensitively on the initial conditions. A wide range o θsl (including θ SL > 9 ) can be produced rom an initially aligned (θsl = ) coniguration (see Figs.4-5). For systems that do not experience eccentricity excitation (because o relatively low orbital inclinations o the companion or/and suppression by GR-induced precession), modest ( 2 ) spin-orbit misalignment can be produced we have derived an analytic expression or θsl or such systems (Eqs.2-24). Note that while our numerical results reer to stellar-mass companions, our analysis is not restricted to any speciic binary ormation scenarios, and can be easily adapted to other types o systems (e.g. when the tertiary is a supermassive BH) by applying appropriate scaling relations. The key dimensionless parameter that determines the spin-orbit evolution is A (see Eq.18). The BH binaries detected by aligo so ar (Abbott et al. 216b, 217) have relatively small χ e (6 +.14.14 or GW15914, 1 +.1 or GW151226, and.12+1. or GW1714). These small values could be due to the slow rotation o the BHs (Zaldarriaga et al. 217) or spinorbit misalignments. The latter possibility would imply a dynamical ormation channel o the BH binaries (such as exchange interaction in globular clusters (Rodriguez et al. 215; Chatterjee et al. 217)) or, as our calculations indicate, dynamical inluences o external companions. BL thanks Natalia Storch, Feng Yuan and Xing- Hao Liao or discussions. This work is supported in part by grants rom the National Postdoctoral Program and NSFC (No.BX216179, No.216M6167 and No.11768). DL is supported by NASA grants NNX14AG94G and NNX14AP1G, and a Simons Fellowship in theoretical physics. This work made use o the Resource in the Core Facility or Advanced Research Computing at SHAO. Abbott, B. P., et al. (LIGO Scientiic and Virgo Collaboration) 216a, PhRvL, 116, 6112 Abbott, B. P., et al. (LIGO Scientiic and Virgo Collaboration) 216b, PhRvX, 6, 4115 Abbott, B. P., et al. (LIGO Scientiic and Virgo Collaboration) 217, PhRvL, 118, 22111 Anderson, K. R., Storch, N. I., & Lai, D. 216, MNRAS, 456, 671 Anderson, K. R., Lai, D., & Storch, N. I. 217, MNRAS, 467, 66 Anderson, K. R., & Lai, D. 217, arxiv:176.84 Antonini, F., & Perets, H. B. 212, ApJ, 757, 27 Antonini, F., Murray, N., & Mikkola, S. 214, ApJ, 781, 45 Antonini, F., & Rasio, F. A. 216, ApJ, 81, 187 Antonini, F., Toonen, S., & Hamers, A. S. 217, ApJ, 841, 77 Barker, B. M., & O Connell, R. F. 1975, PhRvD, 12, 29 Belczynski, K., Benacquista, M., & Bulik, T. 21, ApJ, 725, 816 Belczynski, K., Holz, D. E., Bulik, T., & O Shaughnessy, R. 216, Nature, 54, 512 Blaes, O., Lee, M. H., & Socrates, A. 22, ApJ, 578, 775 Chatterjee, S., Rodriguez, C. L., Kalogera, V., & Rasio, F. A. 217, ApJL, 86, L26 Ford, E. B., Kozinsky, B., & Rasio, F. A. 2, ApJ, 55, 85 Kozai, Y. 1962, AJ, 67, 591 Lidov, M. L. 1962, Planetary and Space Science, 9, 719 REFERENCES Lipunov, V. M., Postnov, K. A., & Prokhorov, M. E 1997, AstL, 2, 492 Lipunov, V. M., et al. 217, MNRAS, 465, 656 Liu, B., Muñoz, D. J., & Lai, D. 215a, MNRAS, 447, 747 Liu, B., Lai, D., & Yuan, Y.-F. 215b, PhRvD, 92, 12448 Mandel, I., & de Mink, S. E. 216, MNRAS, 458, 264 Mardling, R. A., & Aarseth, S. J. 21, MNRAS, 21, 98 Miller, M. C., & Hamilton, D. P. 22, ApJ, 576, 894 Naoz, S. 216, ARA&A, 54, 441 Petrovich, C., & Antonini, F. 217, arxiv:175.5848 Rodriguez, C. L., et al. 215, PhRvL, 115, 5111 Shapiro, S. L., & Teukolsky, S. A. 198, Black Holes, White Dwars, and Neutron Stars: The Physics o Compact Objects, (Wiley, New York, 198) Silsbee, K., & Tremaine, S. 217, ApJ, 86, 9 Storch, N. I., Anderson, K. R., & Lai, D. 214, Science 45, 117 Storch, N. I., & Lai, D. 215, MNRAS, 448, 1821 Storch, N. I., Lai, D., & Anderson, K. R. 217, MNRAS, 465, 927 Thompson, T. A. 211, ApJ, 741, 82 Zaldarriaga, M., Kushnir, D., & Kollmeier, J. A. 217, arxiv:172.885