Binary mergers in strong gravity background of Kerr black hole

Binary-black-hole (BBH) mergers can take place close to a supermassive black hole (SMBH) while being in a bound orbit around the SMBH. In this paper, we study such bound triple systems and show that including the strong gravity effects of describing the SMBH with a Kerr metric can significantly modify the dynamics, as compared to a Newtonian point particle description of the SMBH. We extract the dynamics of the system, using a quadrupole approximation to the tidal forces due to the SMBH. We exhibit how the gyroscope precession is built into this dynamics, and find the secular Hamiltonian by both averaging over the inner and outer orbits, the latter being the orbit of the BBH around the SMBH. We study the long-time-scale dynamics, including the periastron precession and GW radiation-reaction of the binary system, finding that the strong gravity effects of the SMBH can enhance the von Zeipel-Lidov-Kozai mechanism, resulting in more cycles, higher maximum eccentricity, and thereby a shorter merger time, particularly when the binary is close to, or at, the innermost stable orbit of the SMBH. We end with an analysis of the peak frequency of the GW emission from the binary system, highlighting possible observable signatures in the ET and LISA frequency bands.

So far, nearly all of the literature has studied BBHs in vicinity of a SMBH using the approximation that all three BHs are particles.This is valid provided the BHs only affect each other through weak gravity and with relative velocities much smaller than the speed of ★ camilloni@itp.uni-frankfurt.de† harmark@nbi.ku.dk ‡ gianluca.grignani@unipg.it§ orselli@pg.infn.it¶ daniele.pica@nbi.ku.dk light, such that the effects due to General Relativity (GR) can be added through Post-Newtonian (PN) corrections.
However, some BBH mergers can conceivably take place very close to the SMBH, perhaps even near its innermost stable circular orbit (ISCO) Peng & Chen (2021).In this case the particle/PN approximation of the supermassive black hole breaks down since it is only valid if the distance between the BBH and the SMBH is much larger than the Schwarzschild radius of the SMBH.The consequence of this breakdown is that one can only study the system using a metric for the supermassive BH.Indeed, the BBH is exposed to strong gravity effects since its orbit around the SMBH can have a relative velocity comparable to the speed of light and have effects of the gravitational potential that can be considered non-perturbative in comparison with the PN expansion.
In this paper we study a BBH exposed to such strong gravity effects with the SMBH described via the Kerr metric. 1 The only necessary assumption is that the size of the binary system is very small compared to the curvature scale of the SMBH.This alone means that we can model the interaction with the SMBH through tidal forces.
1 As we comment further on below, Maeda et al. (2023a) studied a Newtonian binary system in circular orbit around a large Schwarzschild BH along with results on ZLK oscillations.Our work differs in several ways, but mainly that we consider a large Kerr BH and we include post-Newtonian effects in the binary system such as periastron precession and radiation-reaction effects which have a major effect at long time-scales.Note also footnote 8 regarding how the gyroscope precession enters in the secular Hamiltonian.
The main goal of this paper is to study the long-time-scale dynamics of the BBH-SMBH system, and to understand how this is modified by strong gravity effects arising from a metric description of the SMBH, as opposed to a Newtonian point particle description.
To accomplish this goal, we use the first part of the paper (Secs. 2 to 6) to set up and derive the dynamical equations of the BBH-SMBH system, ending with the secular Hamiltonian that describes the long-time-scale dynamics.In the second part (Secs.7 and 8) we use the secular Hamiltonian to study the long-time-scale dynamics, including the interplay between the ZLK mechanism, the periastron precession, and the GW radiation-reaction forces in the binary system.As we shall see, the strong gravity effects of the SMBH can significantly change the dynamics of the binary system when it is close to the SMBH.
We begin the first part of the paper, in Secs. 2 and 3, with deriving a Hamiltonian for a BBH system in the vicinity of an SMBH described by the Kerr metric.The center of mass motion of the BBH is in our approximation moving on a geodesic of the Kerr metric.For simplicity, we consider a circular equatorial geodesic.The effect of the SMBH on the BBH is described via quadrupole tidal forces, using the work of Marck Marck (1983); Camilloni et al. (2023).At this point, we take the Newtonian approximation of the internal dynamics of the binary system, while the tidal forces are the full GR result from the Kerr metric, which includes strong gravity effects.
In Sec. 4 we discuss the gyroscope precession effect on the SMBH-BBH system.This warrants changing the frame of reference, from a local inertial frame of reference, as found by Marck in Marck (1983), to a non-inertial frame of reference that we call distant-star frame as found in Bini et al. (2016).The distant-star frame provides a global point of view on the orientation of the SMBH-BBH system since the frame orientation can be seen as aligning with the orientation of an asymptotic observer.In this frame, the gyroscope precession angle is directly built into the local dynamics.
In Sec. 5 we exhibit the appropriate Euler angles as well as action-angle variables in the distant-star frame, thereby setting up our study of the dynamics.
In Sec. 6 we reach the goal of the first part of the paper, with the derivation of the long-time-scale dynamics of the BBH-SMBH system including the strong gravity effects arising from a metric description of the SMBH In detail, we find the secular Hamiltonian in Sec.6 by averaging both over the inner orbit of the BBH, as well as the outer orbit, the latter corresponding to the orbit of the BBH around the SMBH.An important consistency check on the secular Hamiltonian is that we obtain the same equivalent dynamics from both Marck's frame of reference as well as the distant-star frame of reference.This is non-trivial as the outer orbit average is taken over different period angles.As we explain, this is all tied together by the gyroscope precession.
We begin the second part of the paper with Sec. 7 in which we take the first steps towards understanding the long-time-scale dynamics.This type of dynamics is interesting due to the ZLK mechanism Lidov & Ziglin (1976); Kozai (1962); Naoz (2016) which is an astrophysically relevant effect observed in triple systems consisting of an exchange between the angular momentum of the inner and outer binary which has an impact on the eccentricity and inclination angle of the inner binary. 2 We show that a change in the ZLK frequency Ω ZLK is the sole difference between the purely Newtonian triple system and our triple BBH-SMBH system in which we consider the full metric of the SMBH.Therefore, we make a careful study of Ω ZLK , including both the frequency Ω (GR) ZLK that one obtains using the local proper time  of the BBH, but also the frequency Ω (∞) ZLK measured by a faraway observer, the difference being a redshift factor.Both of these frequencies are relevant to understand the dynamics of the BBH.We compare these frequencies to Ω (N) ZLK being the frequency obtained by describing all three black holes as Newtonian point particles.In particular, we find the important result that Ω ZLK when the binary system moves on the ISCO, no matter what the SMBH spin is.This is a key to understand many of the further results of this paper.
In Sec. 7 we also consider the large radius limit of the outer orbit, thus making contact with previous results in the literature for which one studied PN corrections to the long-time-scale dynamics of the triple system Liu et al. (2019).We, furthermore, predict new higher-order terms.
The paper culminates with Sec. 8, in which we add the PN corrections of periastron precession and radiation-reaction effect of GW emission to the long-time-scale evolution equations of the binary system, following a similar procedure as Randall & Xianyu (2018), but now including the effect of using the Kerr metric to describe the SMBH.The added PN corrections are highly important to the dynamics of the binary system.The periastron precession is important as it can counteract the ZLK mechanism, narrowing the window of possible initial inclination angles for which it can occur.Furthermore, the radiation-reaction effect of GW emission from the binary system models how the binary system loses angular momentum and energy, ending with the binary merger.
The main result of Sec. 8 is that the strong gravity effect, arising from treating the SMBH via a metric, significantly changes the dynamics of the binary system: • The strong gravity effects are enhanced by a closer proximity of the BBH to the SMBH.We study in particular when the BBH is moving along the ISCO of the SMBH, both with and without spin, corresponding to the largest possible strong gravity effects of the SMBH in our model.
• The strong gravity effects can enhance the ZLK mechanism, resulting in more ZLK cycles, a larger peak eccentricity, and thereby a shorter merger time, since more energy is emitted by GWs when the eccentricity is high.
• The spin of the SMBH can significantly impact the dynamics of the binary system.We find that it plays a significant role in whether the SMBH has zero spin, or whether it counter-rotates or co-rotates with respect to the outer orbit of the binary system around the SMBH.
• Another important aspect of the binary dynamics is the redshift factor, which is used to translate the local time scales of the binary system to the asymptotic time scales as measured by a far away observer.The redshift factor is always greater than 1, increases with counter-rotating spin and decreases with co-rotating spin.
• The various effects have an interesting interplay, with the ZLK mechanism enhanced by strong gravity, accelerating the merger of the binary system, while the redshift factor in a sense slows it down, since the local clock near the binary goes slower than that of an asymptotic observer.
To look for potentially observable signatures for the third-generation interferometer Einstein Telescope (ET), as well as for the Laser Interferometer Space Antenna (LISA), we examine the peak frequency  GW of the spectrum of GWs emitted from the binary system.We find that it is conceivable that one can observe the ZLK cycles of the pre-merger phase of the binary system.Compared to a Newtonian point particle description of the SMBH, the strong gravity effects can result in a larger number of ZLK cycles as well as shorter merger times.Moreover, the maximum peak frequency is lowered by the redshift factor.
We end the paper with our conclusions and give an outlook on future directions in Sec. 9.
Note added: While the preparation of this paper was in its final stages, the preprint of Maeda et al. (2023b) appeared with some overlap to this work.

BINARY SYSTEM IN A TIDAL FORCE BACKGROUND
In this section, we review the basic dynamics of a binary system of black holes in a general tidal force background up to the quadrupole approximation, where the tidal forces originate from the curvature of the space-time that the binary system moves in.For simplicity, we take the Newtonian approximation of the black holes for which we can regard them as particles.Instead, the quadrupole tidal force background is kept completely general.We assume the binary system to be freely falling, hence to leading order its center of mass moves on a geodesic of the background space-time.See Maeda et al. (2023a) for a similar analysis.
In general, we can describe the metric near a geodesic of a background space-time using the Thorne-Hartle version of the Fermi-normal coordinates Thorne & Hartle (1985) where ,  = 1, 2, 3, the square of the geodesic distance to the geodesic is  2 =     and R is the curvature length scale.To use the tidal force description of the background we need that the curvature length scale R of the background is significantly larger than the size  of the binary system, hence  ≪ R Poisson & Vlasov (2010).This is known as the "small-tide approximation".
That we keep terms up to order  2 /R 2 in (2.1) encaptures the quadrupole approximation of the tidal forces, for which the electric and magnetic tidal moments E   and B   are related to the Riemann curvature tensor as with , , ,  = 1, 2, 3. Note that these quantities depend in general on the proper time  of the geodesic.We assume that the background space-time has   = 0 which is indeed true for the Kerr space-time.
Below in Section 3 we will restrict ourselves to tidal moments E   and B   for geodesics around a Kerr black hole, as computed in Marck (1983); Camilloni et al. (2023), but for the moment we keep our considerations general.
We now model the two BHs in the binary system as two particles moving in the background metric (2.1) with the center of mass of the binary located at the geodesic.The Lagrangian for particle 1 is Here   (1) is the velocity of particle 1.The metric   is evaluated at the position ì  1 of particle 1 and it consists of the background metric   with tidal forces (2.1) plus the gravitational potentials generated by particle 2. For the latter, we use the PN expansion of Sections 5.1.4and 5.1.5 in the book Maggiore (2007).Explicitly, we have up to 1PN where the standard PN gravitational potentials  (2) ,  (2) and (2)  arise solely from the gravitational interaction with particle 2, see Maggiore (2007) for details.Here we are neglecting possible terms that correspond to mixed couplings between the tidal forces and the PN expansion.Note that the first appearance would be terms of the type  (2) times the tidal moments.
One can now perform a PN expansion of the Lagrangian (2.3) for particle 1, with the restriction that the tidal moments provide a small perturbation that only should be included up to the first order.Schematically, this gives Here the first term  (1) | E =B=0 is the Lagrangian that arises solely from the gravitational interaction with particle 2, as well as the kinetic energy of particle 1.
The second term in (2.5), coupling to E   , is the leading coupling to the quadrupole tidal moments of the background space-time.
The third term in (2.5), coupling to B   , is subleading as it is suppressed by   (1) /.Note that the coupling between the gravitational potentials and the tidal moments would be of order  2 (1) / 2 or higher, thus not affecting this term.
Restricting for a moment ourselves to the Newtonian limit (with respect to particles 1 and 2) of all of the terms, and adding the analogous Lagrangian for particle 2, we get the Lagrangian for the Newtonian limit of the binary system Introducing now the center of mass quantities (2.7) the Lagrangian (2.6) becomes Since the center of mass motion of the binary is decoupled from the relative motion of the binary system we can consistently set This means we put the center of mass of the binary system at the geodesic, which is what one physically would imagine as well (for large distance scales the binary system should be seen as one particle of mass  placed at the geodesic).Using moreover the general formula (Poisson & Vlasov (2010)) where E  is the quadrupole tidal potential, we get the binary system Lagrangian (2.11) One can now easily Legendre transform this to the Hamiltonian (2.12) with   =   .This Hamiltonian describes the general situation of a binary system of particles moving on a geodesic of a background space-time subject to quadrupole tidal forces from the curvature of the background space-time.This is valid in a local inertial system as set by the Fermi-normal coordinates (2.1).Below in Section 4 we shall explore the same dynamics in a non-inertial system as well.

BINARY SYSTEM IN THE BACKGROUND OF LARGE KERR BLACK HOLE
We consider now a system of three BHs of masses  1 ,  2 and  3 .The two BHs with masses  1 and  2 are in a bound motion and constitute the BBH system.We will often refer to it as the inner binary, and their motion as the inner orbit.Their masses are taken to be much smaller than the mass of the third BH, namely  1 ,  2 ≪  3 , thus identifying  3 as the SMBH.
We will assume two independent separations of scales, in order to gain analytical control of the dynamics.The first one is that the two BHs in the inner binary are in the particle/PN regime for which their Schwarzschild radii are much smaller than their separation : (3.1) In particular, their relative velocity is much smaller than the speed of light.
Secondly, we assume the small-tide approximation, i.e. that the size of the binary system, , is very small compared to the curvature length scale R of the SMBH.This enables us to treat the influence of the SMBH on the BBH via a tidal force approximation, even when the BBH is close to the SMBH.Under these assumptions one has that the motion of the BBH, to leading order, is a geodesic of the Kerr metric describing the SMBH.This we will refer to as the outer orbit.On the geodesic the small-tide approximation  ≪ R can be written explicitly as with r the radial coordinate of the Kerr metric that we introduce below.
We will restrict ourselves in this paper to the leading quadrupole effect, arising from the Riemann curvature tensor of the Kerr metric evaluated on the geodesic.Below we shall include PN effects in the binary dynamics, in the form of the periastron precession and the GW radiation-reaction, thus we assume that the octupole tidal forces are smaller than these PN effects.

Kerr black hole background
In our setup we consider the SMBH to be a Kerr black hole, whose line element, represented in Boyer-Lindquist coordinates x = (t, r, θ, φ), reads (3.3) In the metric above,  3 is the black hole mass,  3 =  3 /(  3 ) is the specific angular momentum and we defined In the following it is convenient to introduce the dimensionless spin parameter The location of the event horizon is given by the major root of the equation Δ = 0 and the cosmic censorship conjecture requires 0 ≤  3 ≤ ( 3 / 2 ), ..0 ≤  ≤ 1.
A generic geodesic x () in the Kerr space-time is parametrized by three constants of motion, respectively representing the energy Ê, the angular momentum L and the Carter constant  per unit of rest energy Carter (1968).
Upon restricting to the circular ( r constant) and equatorial ( θ = /2) case, the tangent vector   , representing the four-velocity of a particle moving along the geodesic can be written as where the redshift factor   and the coordinate angular velocity Ω are defined as and where we use  = ±1 to respectively distinguish orbits that are co-rotating and counter-rotating relative to the angular momentum of the Kerr black hole. 3For later convenience, we also report the orbital angular velocity defined with respect to the proper time, given by Ω φ =   Ω, namely (3.8) In this case it is possible to obtain explicit expressions for the constants of motion in terms of the orbital radius r, the black hole mass  3 , and the specific spin value  3 .More specifically, one has the energy, orbital angular momentum, and Carter constant per unit of rest energy (3.9) For later use, we introduce here the ISCO, which identifies the innermost stable circular orbit.In the equatorial plane of a Kerr black hole there exist two ISCOs, one of which is co-rotating ( = +1) and one which is counter-rotating ( = −1) with the black hole.The value of the radial coordinate at which the ISCOs are located in the equatorial plane of the Kerr space-time is given by Bardeen et al. (1972a) r where (3.11) One can see immediately that for  = 0 eq.(3.10) reduces to rISCO = 6 3 / 2 , i.e. the usual ISCO position for a Schwarzschild black hole.

Binary system on a Kerr geodesic
As explained above, we can regard the black holes of the BBH system as particles in a Newtonian limit since we assume the condition (3.1) holds.The BBH system moves along a circular equatorial orbit in the background of a Kerr black hole describing the SMBH.As seen in Section 2, for a sufficiently small binary system (3.2) we can approximate the influence of the SMBH through quadrupole tidal forces acting on the BBH, while the center of mass of the BBH moves approximately on a geodesic of the Kerr metric.This means we can employ the results of Section 2.
For the Kerr metric (3.3), we can use the Marck's tetrad Marck (1983) to describe the local Fermi-normal coordinates (2.1) of Section 2. Marck's tetrad is given by the four vectors    which provide an orthonormal basis for the vector space at each point of the Kerr-geodesic since         =   and         =   , where in particular   0 =   is the four-velocity.Note that  = 0, 1, 2, 3 are the flat tetrad indices.One can equivalently represent Marck's tetrad as the orthonormal one-forms    =        .One can now employ the standard map of the Fermi-Normal coordinates between vectors on the Kerr geodesic and events in the neighborhood of the geodesic.For Marck's tetrad, we have in addition to the time coordinate  0 = , which is the proper time on the geodesic, also the spatial coordinates   parametrizing an orthogonal vector 3 =1      at the geodesic.In this way the coordinates   describe a neighborhood of the geodesic.We will call this coordinate system the Marck's frame of reference.Since Marck's tetrad is parallel-transported along the geodesic, it provides an inertial frame, meaning it is characterized by a vanishing acceleration and vanishing angular velocity of rotation of spatial basis vectors Misner et al. (1973) 4 For circular geodesics in the equatorial plane ( θ = /2), Marck's tetrad can be conveniently written in the one-forms basis as where r is the constant radius of the circular equatorial geodesic and the angle Ψ is introduced to ensure that the tetrad is parallel transported along the geodesic Marck (1983) Misner et al. (1973).
2, given in Eqs.(2.11) and (2.12), respectively, we record that the electric tidal moments in the equatorial plane are Marck (1983); Camilloni et al. ( 2023) (3.16) Using these with Eq. ( 2.10), one gets an explicit expression for the scalar quadrupole electric tidal moment induced by the Kerr black hole, as measured by an observer using Fermi-normal coordinates, namely (3.17) With this, we can describe the dynamics of the binary system in the approximations (3.1) and (3.2) via the Lagrangian (2.11) and Hamiltonian (2.12) of the BBH system.

GYROSCOPE PRECESSION AND THE DISTANT-STAR FRAME
In this section, we consider the gyroscope precession of Marck's parallel transported frame of reference by introducing a non-inertial frame of reference that we dub the distant-star frame.5 The BBH has an orbital angular momentum that undergoes a gyroscope precession in its motion along an equatorial circular geodesic in Kerr space-time around the SMBH.This is induced by the curvature of the background, like in the case of the Earth-Moon binary system in orbit around the Sun Rindler & Perlick (1990).In the Schwarzschild space-time this gyroscope precession is known as the Fokker-de Sitter precession de Sitter (1917), whereas in the equatorial plane of a Kerr black hole, it takes the name of Schiff's precession Schiff (1960).
The origin of this precession is the difference between the local and global points of view for our binary system moving on a geodesic.Marck's frame represents the local view point, where we have an approximate inertial system close to the center of mass of the binary that moves on the geodesic.However, there is also a global point of view, in which the global properties of the Kerr space-time are taken into account.This is clear in our case of an equatorial circular geodesic.For this motion, the only spatial coordinate in the BL coordinates of Kerr that changes is the angle φ as where Ω φ is given in (3.8).A period of motion is obviously when φ changes with 2.However, Marck's frame is not the same after one period, since the Ψ angle has changed with ΔΨ = 2(Ω Ψ /Ω φ − 1), which gives the gyroscope precession.
Since the precession is not observable in Marck's frame of reference by itself, it is useful to define a non-inertial frame of reference in which the gyroscope precession is manifest.Such a frame, here called the distant-star frame of reference, is constructed in Bini et al. (2016), simply by rotation of Marck's frame with angle φ − Ψ such that the distant-star frame is periodic under rotations with respect to the φ angle.Specifically, seeing it as a tetrad    , it is defined by the following rotation of Marck's tetrad where we introduced the gyroscope angular velocity Ω g as It is easy to check that Ω  = 0 at r = 4/9 ( 3 / 2 )  2 , and that this location never corresponds to a stable orbital radius, since r★ < r  ISCO for 0 ≤  ≤ 1.For completeness, we report that at the ISCO one has √ 6 in the Schwarzschild case  = 0, whereas in the extreme Kerr case,  = 1, the gyroscope precession diverges as Ω  ≈ 2/ √ 3( r − r+ ISCO ) −1 for co-rotating orbits, and ) for the counter-rotating ones.As we shall see below, the distant-star frame provides a local coordinate system close to the circular equatorial geodesic in which one can directly observe the precession as a fictitious force in the Lagrangian description.
In general, a local observation of a precession angle is not possible, since one cannot compare angles between two events in space-time in a path-independent manner.However, the construction of the distant-star tetrad is based on the global structure of the Kerr space-time, being stationary and axisymmetric, which gives a natural definition of angular and radial directions in the equatorial plane through Carter's tetrad Bardeen et al. (1972b).Thus, in this sense, one can meaningfully claim the distant-star frame is fixed with respect to the asymptotic definition of the rotating angle, justifying its name as an angle with respect to distant stars.I.e. the distant-star frame provides a Cartesian frame that keeps a fixed orientation with respect to distant stars Bini et al. (2016).Hence, the non-inertial distant-star frame of reference provides a global point of view, contrary to the local inertial Marck's frame of reference.
One has that the distant-star tetrad    is characterized by a vanishing acceleration, being anchored to the geodesic, but a non-vanishing angular velocity of rotation relative to the Marck's tetrad Misner et al. (1973 where we used the vector tetrad basis The spatial coordinates   associated with the distant-star tetrad are given by6 In the following, we shall use the Cartesian vector notation where we defined the unit vectors The consequence of going to this non-inertial frame for the binary system is the introduction of fictitious forces.Indeed, the local Lagrangian now becomes where we introduced the Cartesian vectors  =  / and  g = Ω g ẑ in agreement with Eq. (4.4).Now Eq. (2.10) reads (4.9) To find the Hamiltonian we define the canonical momentum as and the canonical angular momentum as where we adopt the subscript "in" for later convenience to distinguish the angular momentum of the inner BBH system and the angular momentum associated with the outer orbit.The Hamiltonian, thus, reads with E  given in Eq.( 4.9).The extra term in (4.12) with respect to (2.12) is responsible for the gyroscope precession of  in .

EULER ANGLES AND ACTION-ANGLE VARIABLES
In Sections 5.1 and 5.2, we review certain standard definitions of angular coordinates and momenta that are highly useful to describe the inner binary, and to derive the secular Hamiltonian in Section 6.
All of these definitions are with respect to the distant-star frame of reference introduced in Section 4. Subsequently, in Section 5.3 we introduce for later use the action-angle variables in Marck's frame of reference as well, and find the canonical transformation between the distant-star and Marck's frame Hamiltonians.

Euler angles
While the center of mass for the inner Newtonian BBH system moves along the circular equatorial geodesic in the Kerr space-time, the orientation of its inner orbital plane can vary with respect to the outer orbital plane of the Kerr SMBH.
The vector (4.6) describing the relative position of a body in a Newtonian elliptic orbit can be represented as where where  and  are respectively the semi-major axis and the eccentricity of the orbit, whereas  is the angular coordinate that keeps track of the body motion along the orbit, namely the true anomaly.
The two directions û and v have a precise geometrical meaning, with û ( = 0) identifying the periapsis direction, and v ( = /2) lying along the direction of the ascending nodes.The space spanned by these two vectors specifies the inner orbital plane.Including the direction of the angular momentum for the inner binary Lin = û × v, one obtains a triad of orthonormal vectors.
Since the inner binary is assumed to be in a Newtonian regime, it is possible to unambiguously introduce the eccentric anomaly  and the mean anomaly , defined by means of (5. 3) The first relation defines  in terms of the true anomaly , whereas the second is the Kepler equation.Kepler's equation is a transcendental equation and no closed-form solution is known that allows the expression of the eccentric anomaly  in terms of the mean anomaly .
The mean anomaly  is particularly useful as it represents the angle that a fictitious body moving in a circular orbit would span if it had the same orbital frequency as the actual body moving along the elliptic orbit.In other words its motion is uniform in time,  = √︁  / 3 .
One can obtain a generic orientation of the orbit using the Euler angles defined through the following rotation matrices where the angle  is called the longitude of ascending nodes,  is the orbital inclination, and  is the argument of the periapsis.We refer the reader to Fig. 1 for an illustration of the orbital parameters.
A Newtonian orbit with arbitrary orientation can be obtained by performing a rotation around a reference plane.In our setup, where For later use, we define in addition the eccentricity vector as which corresponds to the dimensionless version of the Laplace-Runge-Lenz vector Goldstein (1980).(5.8) Since the reference plane in our setup corresponds to the Kerr equatorial plane, θ = /2, and the outer orbit is a fully relativistic equatorial circular geodesic, the magnitude of the outer orbit angular momentum is L, as given in (3.9), and its direction in the distant-star frame is (5.9) Notice that, by definition, the inclination angle quantifies the projection of the inner orbit angular momentum on the outer orbit one, Lin • Lout = cos .

Action-angle variables
For any libration of periodic motion, we can introduce action-angle variables to describe the momenta, which is advantageous since they are constants of motion.This is useful for our description of the BBH-SMBH system, as it means that the non-tidal part of the Hamiltonian (2.12) is expressed purely via constants of motion.
We shall use the action-angle variables known as Delauney variables, with the position given by the three angles (, , ) , (5.10) each periodic with period 2, as well as the corresponding action-angle variables (5.11) with  and  defined in Eq. (2.7) as the total mass and the reduced mass of the BBH system respectively.
It is worth noting that the Delauney action variables are related to the magnitude and the orientation of the angular momentum of the inner binary with respect to the reference plane.In particular (5.12) The total Hamiltonian of the inner binary (4.12) in the distant-star frame of reference can therefore be simply expressed as (5.13) and the quadrupole tidal part is expressed in terms of Euler angles according to (5.14) We notice that the non-tidal part of (5.13) indeed is given by an action-angle variable.

Action-angle variables for Marck's frame
For use below, we also introduce the action-angle variables for the Hamiltonian in Marck's frame of reference given by Eqs.(2.12) and (3.17).
It follows from Section 4 that the longitude of ascending nodes for Marck's frame of reference is (5.15) in terms of the corresponding angle  of the distant-star frame.We can also write this as  ′ − Ψ =  − φ.Instead, the other angles  and  remain the same.
The Hamiltonian along with the action-angle variables conjugate to the angles (, ,  ′ ) are most efficiently found using a canonical transformation of the second type Goldstein (1980).This reveals that the momenta (  ,   ,   ) given in the distant-star frame (5.11) are the same for Marck's frame, for it to be a canonical transformation.
In detail, we have the generating function (5.16) with (, ) and (, ) being the distant-star and the Marck phase space variables, respectively, and with the identifications (5.17) This gives   =  2 /  and   =  2 /  as needed.The transformed Hamiltonian therefore becomes (5.18) From this we get the Hamiltonian in Marck's frame as (5.19) with quadrupole tidal part (5.20) One can easily check that if one introduces Delauney/action-angle variables directly for Marck's frame Hamiltonian given by Eqs.(2.12) and (3.17) one would get the same result as above.

SECULAR HAMILTONIAN
In this section, we obtain the secular dynamics of the BBH-SMBH system.The secular dynamics describes the system at timescales that are larger than both the inner and outer orbit periods, i.e. both the time scale of internal motion within the BBH as well as the time scale of the geodesic motion of the BBH around the Kerr black hole.
The secular Hamiltonian is obtained by taking the average both over the inner orbit motion as well as the outer orbit motion.As explained above, these two motions can be separated to leading order in our regime (3.2), in that the outer orbit motion corresponds to the center of mass of the binary system moving on a circular geodesic in the equatorial plane of the supermassive Kerr black hole.
The inner binary is described by a Newtonian elliptic motion, which is perturbed by tidal forces. 7To take the average, we need an angle that grows uniformly with time in the elliptic motion.This is provided by the mean anomaly  defined in Eq. (5.3).However, since the tidal part of the Hamiltonian (5.14) is a function of the true anomaly  instead, we translate since this follows from Eq. (5.3).
For the outer binary, the averaging procedure is more subtle, reflecting the fact that general relativistic effects play a role.In terms of the action-angle variables introduced in Section 5, we have found the Hamiltonian in both the non-inertial distant-star frame with Eqs.(5.13)-(5.14)as well as in the inertial Marck's frame with Eqs.(5.19)-(5.20).
We begin by considering the outer orbit average in the distant-star frame since this bears the closest resemblance to the Newtonian case (see for example Randall & Xianyu (2018)) and therefore is more intuitive.Indeed, considering the Hamiltonian (5.13)-(5.14),we notice that it is periodic in the angle φ with period 2.This periodicity is precisely associated with one outer orbit cycle of motion.Moreover, the angle grows linearly with proper time, as one can infer from Eq. (4.1), making it the relativistic analog of the outer orbit angle, that one for instance uses in Randall & Xianyu (2018).Therefore, the outer orbit average in this frame is simply performed 7 In Section 8 we shall include the 1PN effect of the periastron precession as well as the leading GW radiation-reaction effect.
as (2) −1 ∫ 2  0  φ.The secular Hamiltonian is thus computed as the following doubleaverage of the Hamiltonian (5.13) in the distant-star frame with where the subscript ZLK refers to the ZLK effect which will be extensively discussed in the remaining Sec. 7 and in Sec.8.8 The expressions above make manifest the fact that the GR effects associated with the Kerr perturber enter the dynamics through the overall prefactor Ω (GR) ZLK in the averaged tidal Hamiltonian (6.4).It is immediate to see from (6.4) that in the weak field regime r → ∞ one recovers the Newtonian secular Hamiltonian of Ref. Randall & Xianyu (2018).These GR effects are completely accounted for by the term in Ω However, there are other GR effects as well, some of which are not immediately apparent from the above secular Hamiltonian of Eqs.(6.2) and (6.4).The most obvious one is the gyroscope precession of the binary system introduced in Section 4, here arising from the term Ω g   in the total Hamiltonian (6.2).Another less explicit effect is the time dilation of the proper time used above, relative to the asymptotic time t.This we shall include later in Section 8.Both of these effects are related to how an asymptotic observer will view the binary system, i.e. the global point of view, rather than the local point of view.Furthermore, in Section 8 we shall add further relativistic effects to the binary dynamics, to describe gravitational backreaction due to gravitational waves.
So far we have considered the averaging procedure only in the distantstar frame.It is important to check that one can obtain the same secular average in Marck's frame of reference.We notice immediately that the Hamiltonian (5.19)-(5.20) is periodic in Ψ with period 2.Also, the angle grows linearly with time Ψ = Ω Ψ .Therefore, we conclude that, in Marck's frame, one should compute the outer orbit average as with It is now straightforward to see that double-average over the tidal part of the Hamiltonian in the two different frames agree This means the only difference between the secular Hamiltonians in the two frames is the constant term Ω g   , accounting for the fictitious forces.
However, as explained in Section 4, the gyroscope precession essentially measures the difference between φ and Ψ when they have gone through one cycle in the outer orbit motion.So how can the outer orbit average over the tidal contribution give the same result in the two different frames, as we are averaging over two different angles?
The answer lies in the formula φ −  = Ψ −  ′ .In the distant-star frame one should keep fixed  in taking the average, as this angle is fixed during the motion.But, for Marck's frame, it is instead the angle  ′ that one should keep fixed.Thus, the reason that the outer orbit averages give the same result in the two frames is that the difference between  and  ′ precisely accounts for the difference φ − Ψ, which is the gyroscope precession.
Notice that the canonical transformation detailed in Sec.5.3 can also be directly used to relate the secular Hamiltonian in the distant-star frame with the secular Hamiltonian in Marck's frame.Since neither depend on the respective longitude of ascending nodes angle ( and  ′ ), this transformation simply relates the secular Hamiltonians as ⟨⟩ = ⟨H ⟩ − Ω g   .

ZLK MECHANISM IN A STRONG GR BACKGROUND
In this section, we apply the result for the secular Hamiltonian in the distant-star reference frame, as derived in the previous section, to study the long timescale dynamics of the BBH system moving on an equatorial circular geodesic of the external Kerr SMBH.
An important result of Sec.7.1 is that we can quantify to what extent the ZLK frequency departs from its Newtonian value when one takes into account strong gravity effects associated to the general relativistic description that we adopt for the outer orbit.
Moreover, we derive in Section 7.2 the equations of motion for the inner orbital parameters, from which one can study the evolution of the ZLK mechanism.This is used in Section 7.3 to compare the weak-gravity limit of our results to the PN corrections found in the literature.
In Section 8 we build on the results of this section by refining the equations of motion for the inner orbital parameters found in Section 7.3 to include the periastron precession and GW emission.This is used to study the ZLK mechanism and its influence on the binary merger time.

ZLK frequency in the vicinity of SMBH
From the distant-star frame secular Hamiltonian, given by Eqs.(6.2) and (6.4), it is immediate to observe that two main effects govern the secular dynamics of the BBH system: the ZLK mechanism, which manifests due to the tidal interaction with the external SMBH, and the gyroscope precession, which is present in the distant-star frame of reference whenever a general relativistic description for the SMBH is adopted.
We begin this section by discussing the ZLK mechanism which in the distant-star secular Hamiltonian is modeled by the quadrupole tidal contribution ⟨H  ⟩ in Eq. (6.4).A well-known result in the literature concerning the ZLK mechanism Lidov & Ziglin (1976); Kozai (1962); Naoz (2016) is that the set of parameters characterizing the outer orbit only enters in ⟨H  ⟩ through the frequency of the eccentricity/inclination oscillations.We observe from (6.4) that this remains true in our case, as all the information concerning the outer orbit enter through the frequency Ω (GR)

ZLK
. Thus, all new tidal force effects that arise from an exact metric description of the SMBH as a Kerr black hole enters through this frequency.
Therefore, the main aim of the following is to show that BBH systems close enough to an external SMBH to probe the strong gravity regime can manifest substantial deviations in the frequency of the ZLK oscillations, compared to the frequency one gets from employing a Newtonian point particle approximation.
We begin by comparing Ω (GR) ZLK to the Newtonian frequency Ω (N) ZLK , i.e. the frequency that would have resulted from Newtonian quadrupole tidal forces induced by a particle of mass  3 .We find where we defined for convenience the dimensionless radius  for the equatorial circular orbit as Remarkably, at the ISCO r = r  ISCO the ratio (7.1) takes the universal value 9 ZLK .
(7.3)This result will be highly important in Section 8 where we consider the evolution of the BBH-SMBH system in detail.One can check that (7.3) gives the maximal value of the ratio (7.1) that the binary system can attain.Instead, for large  the ratio goes to one.Both of these statements are illustrated in Fig. 2. The result (7. 3) shows that one has an order one difference between the weak-field Newtonian result Ω (N) ZLK and our novel strong field result Ω (GR) ZLK when close to the SMBH.However, it is important to note here that the frequency Ω (GR) ZLK is measured with respect to the proper time of the BBH orbit.Thus, this is not the frequency that an asymptotic observer would measure.To find the corresponding asymptotic ZLK frequency we need to incorporate the redshift factor as follows where the redshift factor can be written as which one can check is always greater than 1, and it is a decreasing function of   for fixed .Using (7.5) we get the following ratio between the asymptotically measured ZLK frequency Ω ZLK is plotted in the counter-rotating case  = −1 as a function of the dimensionless radius  for several values of the spin.The figure shows that the maximum value for the ratio is always at the ISCO.
now takes into account all the GR effects, and the corresponding Newtonian frequency Ω The inclusion of the redshift factor gives a more refined difference in the ratio of the frequencies.One finds that in the counter-rotating case  = −1 the maximal value of the ratio (7.6) is at the ISCO, as illustrated in Fig. 3.For the co-rotating case  = 1, the same is true for the range 0 ≤  ≤ 0.69.However, as illustrated in Fig. 4, this behavior starts changing in the range 0.69 ≤  ≤ 0.7, so that for  ≥ 0.7 the maximal value of the ratio (7.6) is no longer reached at the ISCO.We see from Fig. 5 that for  ≥ 0.75 it is instead the minimal value that one reaches at the ISCO.For all three Figs.3, 4, and 5 we note that the ratio (7.6) approaches 1 for  going to infinity, as one would expect.Finally, we have also plotted the value of the ratio (7.6) at the ISCO in Fig. 6.
The above results for ZLK show that the GR effects that arise from being in close vicinity to the SMBH are highly significant.As mentioned in the Introduction, this is particularly relevant in the case of bound systems of BBHs situated in the GC.Indeed, we see the importance of including strong-gravity effects as they significantly alter the frequency, and therefore the timescale, of the ZLK is plotted in the co-rotating case  = 1 as function of the dimensionless radius  for the range 0.61 ≤  ≤ 0.79.
ZLK is plotted in the co-rotating case  = 1 as a function of the dimensionless radius  for several values of the spin.ZLK oscillations.One can also see that the spin of the SMBH, as modeled by a Kerr black hole, can significantly alter the dynamics.

Evolution equations for ZLK mechanism
Using the secular Hamiltonian (6.4) we can now derive the evolution equations for the orbital variables describing the inner BBH system.It is possible to derive an evolution equation for the orbital inclination  by exploiting the fact that   =   cos , Eq. (5.11), which yields This general equation will be useful when including the loss of angular momentum associated with the emission of GWs, which we postpone to Sec. 8.In this section instead, we only focus on the effect of tidal deformations resulting in the ZLK mechanism.
For the secular Hamiltonian (6.4) the Euler angle  is a cyclic variable so that its conjugate momentum   is a constant of motion,   / = −  ⟨H ⟩ = 0, and the last term in (7.7) does not contribute if we ignore the emission of GWs.
Similarly, the equation of motion for the eccentricity of the inner binary follows from the definition of the Delauney variable   in Eq. (5.11).In the following we use the fact that no variation of the semi-major axis  exists in the absence of GW emission, so that Therefore, the evolution equations for the orbital elements can be derived from the equation of motion of   .This follows from Hamilton's equations and only involves the tidal Hamiltonian From the relations (7.7), (7.8), and (7.9) one can derive the ZLK contributions to the evolution equations for the orbital inclination  and for the eccentricity From these equations of motion, it is immediate to notice that the stationary points for minimum and maximum inclination and eccentricity correspond to  = 0, /2.As we mentioned earlier, as long as the gravitational backreaction is neglected, ) cos  is conserved, meaning that the orbital inclination has a maximum  max when the eccentricity is minimal  min , and viceversa.In particular, by computing the second derivative in (7.10) and (7.11), it is easy to show that the pair ( max ,  min ) occur for  = 0, whereas one has ( min ,  max ) for  = /2.
The evolution equation for the longitude of ascending nodes can be derived from Eq. (6.4), upon making in ⟨H ⟩ the substitution cos  =   /  .One has where it is evident the contribution of the gyroscope precession in the distant-star frame.Using the canonical transformation detailed in Sec.5.3 it is immediate to derive an analogous equation for the shifted angle  ′ =  − Ω   in Marck's frame, which indeed lacks the gyroscope precession contribution.We stress that in passing from the distant-star to Marck's frame of reference only  changes in  ′ , and therefore all other equations of motion written before remain unaltered.
Finally, the equation of motion for the argument of the periapsis  can be found by trading the eccentricity  for the angular momentum   and using again that cos  =   /  .We get (7.13)

General-relativistic effects in the ZLK mechanism
The ZLK mechanism is an effect for which a binary system under the influence of the tidal forces of an outer third body can exhibit a periodic exchange of eccentricity and orbital inclination, with a timescale much larger than its orbital period Kozai (1962); Lidov & Ziglin (1976); Naoz (2016).
To understand under which circumstances this mechanism can operate, one can start by recalling that, in the Newtonian approximation for the inner binary, the secularly-averaged Hamiltonian ⟨H ⟩ and the angular momentum projection   are conserved quantities.Their values can therefore be estimated by fixing initial conditions for the eccentricity  0 and the inclination  0 .The corresponding values for ⟨H ⟩ and   will be labeled as ⟨H ⟩ 0 and (  ) 0 .By restricting to the case in which the inner orbit is initially circular,  0 = 0, from the conservation of energy and angular momentum, explicitly Notice that, being  min independent of the initial inclination  0 , it does not only constitute the minimum inclination reached in the ZLK oscillation but also the critical angle for the onset of the ZLK effect Naoz (2016).For the ZLK mechanism to work, one has that the initial inclination  0 should obey |cos  0 | < √︁ 3/5.For an initial inclination  0 ≈ /2 of the inner binary system in circular orbit ( 0 = 0), the system exhibits high eccentricity  max ≈ 1.We thus find that including full general relativistic effects for the outer orbit does not alter the way the ZLK mechanism is triggered, it only affects the frequencies associated with the ZLK oscillations, as discussed in Sec.7.1.In Sec. 8 we will include PN effects in the inner binary dynamics that will modify the condition (7.14)Will (2014).
A more exhaustive way to see how the ZLK mechanism manifests itself together with the general-relativistic gyroscope precession consists of considering the interchange between the eccentricity vector  of the inner binary and its angular momentum  in relative to the orbital plane of the outer binary.More specifically, from Eqs. (5.7),(5.8)and (5.12), one has where we used   = Ω  Lout and the results obtained in Eqs.(7.8) and (7.9).
Eqs. (7.18) and (7.19) are written in terms of the proper time  related to the inner binary system.From the point of view of an asymptotic observer, they become where we defined Notice that the redshifted gyroscope precession frequency above is always finite: at the ISCO, for instance, one has Ω 3) in the non-spinning case  = 0, whereas at extremality,  = 1, one finds Ω (∞)  = 1/2(/ + ISCO ) and Ω  = (/ r − ISCO ) (4/ √ 3 − 9)/26 respectively for the co-rotating and counter-rotating orbits.Eq.s (7.20) and (7.21) include the gyroscope precession and extend the results previously known in the literature for the ZLK effect Liu et al. (2015); Liu et al. (2019); Kuntz et al. (2023) to the case in which the external body, in this case a supermassive Kerr black hole, is described using the Kerr metric (3.3), thus being in a strong GR regime.
In the case of a circular orbit for the inner binary,  = 0, which constitutes a solution for the equation for the eccentricity, Eq. (7.20) becomes which describes the precession of the angular momentum of the inner binary  in around the direction of the angular momentum of the outer binary Lout .In the weak-field limit r → ∞, one has where in the first line, the first term represents the 1PN contribution due to the gyroscope precession, the second term is a relativistic effect related to the spin of the SMBH and the third term represents the precession generated by the standard ZLK mechanism.This is consistent with the results presented in Ref. Liu et al. (2019), in the hierarchical regime for circular outer orbits and where the contribution proportional to the spin of the SMBH comes from the Lense-Thirring precession.The second line is instead a new result and represents higher-order contributions which we predict using the result in Eq. (7.23).
As pointed out in the past literature (see for instance Petrovich & Antonini (2017), Liu & Lai (2017) and Liu et al. (2019) ) the interplay of the ZLK mechanism with additional precessing effects can lead to significant alteration in the BBH dynamics and exhibit chaotic features.The dynamical behaviour can be identified by means of an adiabatic parameter, that we define as where the ratio is evaluated at the ISCO, since r  ISCO marks the scale at which the strong-gravity effects are more relevant in our setup.Fig. 7 shows plots for the ratio  in a specific configuration for the BBH-SMBH system as a function of the inner BBH system eccentricity  and the SMBH spin parameter .It would be interesting to study the interplay between the gyroscope precession and the ZLK mechanism further.

BINARY MERGER TIME CLOSE TO A SUPERMASSIVE BH
In this section, we refine the dynamics of the inner binary system by adding the periastron precession and GW emission, so that we can study the interplay of these effects together with the ZLK mechanism.This enables us to study how treating the SMBH in strong gravity can alter the dynamics of the BBH system compared to if one included only the Newtonian gravity effect of the SMBH Blaes et al. (2002); Randall & Xianyu (2018).
We recall that to treat the presence of the SMBH as a perturbation of the BBH system the condition (3.2) must be satisfied.After introducing the parametrization for an elliptic orbit as in (5.2), the tidal condition can be rewritten in terms of the semi-major axis  of the binary system as This ensures that we can safely neglect the GW backreaction of the outer orbit and consider therefore only the GWs emitted by the inner BBH Randall & Xianyu (2018).

Post-Newtonian dynamics of the binary system
GWs emitted by the inner BBH system reduce its energy and angular momentum and consequently its semi-major axis  and eccentricity .Peters's equations Peters & Mathews (1963); Peters (1964) keep into account this gravitational backreaction by providing the orbitaveraged evolution of  and  for an isolated binary.The average variations read These equations were obtained to model the gravitational backreaction for an isolated binary, but one can include the influence of a third external body, that interacts with the binary through the ZLK mechanism, by adding these PN contributions to the system of equations describing the time evolution of the orbital variables (, , , ), as given in Eqs.(7.10), (7.11), (7.13).In particular, the first of the two PN contributions above is responsible for triggering the inspiral phase for the binary system, by decreasing the relative distance  between the two masses.As opposed to the ZLK mechanism, which can lead to an increase of the orbital eccentricity , the second equation shows how the emission of GWs tends to circularise the orbit.
We recall that   represents the magnitude of the angular momentum for a Newtonian binary system.Hence, by combining Eqs.(8.2) with the definition of   , according to Eq. (5.11), it is immediate to estimate the loss of angular momentum associated to the GW emission, namely Landau & Lifschits (1975) It is important to stress that the gravitational backreaction only reduces the magnitude of the angular momentum, without affecting its direction.In other words, if we ignore the presence of the SMBH, the binary system remains in the same plane when including the GW emission.Thus, we have ⟨/⟩ GW = 0 which from Eq. (7.7) gives It is immediate to notice, using Eq. ( 8.3), that the loss of angular momentum (or, analogously, the decreasing of the orbital distance and eccentricity, according to Eq. (8.2)) due to GW emission becomes extremely efficient when the condition  ≈ 1 is met.
As we discuss this in more detail in the following subsection, the enhancement in the eccentricity due to the ZLK mechanism can therefore boost the merger of eccentric binary systems.While such ZLK-boosted mergers were already noticed in previous works that analyzed the combined ZLK dynamics with PN effects (see for instance Refs.Blaes et al. (2002); Randall & Xianyu (2018)), in the next subsection we show that strong gravity effects associated to the external SMBH can lead to further significant changes in the frequency of the ZLK-oscillations and in the merger time.
Finally, following Randall & Xianyu (2018), we include also the effect of the periastron precession of elliptic orbits as an additional contribution to the evolution equations, since this 1PN effect plays a significant role in the ZLK mechanism.Indeed, as we shall review later, it is known to limit the range in which the ZLK mechanism is valid.Notice that the periastron precession preserves the angular momentum vector   in / = 0, meaning for instance that the orbital plane remains unchanged and only manifests itself as an apsidal advance  êin / = (/) v (see Eqs. (7.16) and (7.17)).The periastron precession can therefore be written in terms of a first derivative for the argument of periapsis , according to As discussed below, these PN effects alter the dynamics of the ZLK mechanism for BBH systems influenced by the presence of an external mass.
We can now write down a closed set of equations by combining the evolution equations Eqs.(7.10), (7.11), (7.13) of Sec.7.2, which considers the strong gravity effects of an external spinning SMBH, with the effects of GW emission and periastron precession introduced above.This gives the evolution equations where the equations have to be supplemented with the definition , as in Eq. ( 5.11), and with the generalrelativistic definition of the ZLK frequency Ω (GR) ZLK , that we derived in Eq. (7.1).These evolution equations are written with respect to the proper time , being the time associated with the inner binary system.However, one can use   t = 1      to translate these equations into evolution equations with respect to the asymptotic time t to obtain the dynamical description as seen by an asymptotic observer.
The system of equations (8.6) allows us to study how the orbital parameters (, , , ) of the inner binary, evolve in time.To solve them, one has to specify initial conditions for the inner BBH system, ( 0 ,  0 ,  0 ,  0 ), a set of parameters that characterize the outer circular equatorial orbit around the Kerr SMBH ( r,  3 ,  = ±1), and the three masses ( 1 ,  2 ,  3 ).
In the next subsection, we solve numerically the system of equations (8.6), for some specific case of interest.We will work in astronomical units, i.e. (AU,  ⊙ and years), where  = 4 2 and  = 63072 AU/years.
As previously explained, we need to impose the condition (8.1) for the tidal approximation to be valid.This allows us to neglect the GW backreaction of the outer orbit over the timescale that governs the inner binary dynamics.Furthermore, we are working in the near-Newtonian regime (3.1) in which we can treat the black holes in the inner binary approximately as particles.
The ZLK oscillations take place on a timescale much longer than the orbital period of the binary system around the SMBH.Thus, it becomes important to ensure that the triple system is stable, allowing the ZLK mechanism to enhance the eccentricity of the BBH system.When the binary system is too close to the SMBH, the interaction between the three bodies can result in a tidal breakup.To ensure that the presence of an SMBH does not lead to a tidal break up of the binary system, the following condition must be satisfied Miller et al. (2005) where we recall that  =  1 +  2 is the total mass of the BBH system and  3 the mass of the SMBH perturbing the binary.In our analysis, we carefully chose the radial distance between the binary system and the SMBH so as to always guarantee the stability of the triple system, more specifically in all our examples, we will consider r ≫ rtide .As an illustrative example, in Fig. 8 we plot the regions of the parameter space for the SMBH which are consistent with the tidal breakup condition (8.7), in the case of a stellar-mass BBH system with masses  1 =  2 = 10  ⊙ and separation  = 0.1 AU.To prepare for the analysis we perform in the remaining of this section, we now briefly introduce respectively the GW and PN timescales at play in our investigation and the peak frequency of GW emitted by the BBH system.These will play a key role in the next subsections where we show how the GR effects in the strong field regime affect the merger time of the binary system and the emission of GW waves.

Time Scales
For an isolated BBH system the merger time due to the emission of GWs is Peters (1964); Maggiore ( 2007) where . (8.9) The function  ( 0 ) ∈ [0.979, 1.81] for  0 ∈ [0, 1], and for an order of magnitude estimate it can be ignored.For all practical purposes, a good approximation for an isolated BBH system merger time is thus given by The prefactor (1 −  2 0 ) 7/2 shows that for very high values of the eccentricity, the merger time can be drastically reduced.Tidal interactions generated by an external SMBH, through the ZLK mechanism, can produce high eccentricities, thus catalyzing the coalescence for highly-inclined BBHs and speeding up the merger to timescales much shorter than those characteristic of isolated binaries with the same masses and relative distance Antonini & Perets (2012a).
From Eq.(8.5) we can compute also the time scale associated with the periastron precession as In this regard, it is interesting to notice that, when the periastron timescale becomes comparable with the ZLK timescale  ZLK = 2/Ω ZLK , the periastron precession destroys the ZLK resonance and the binary system begins to evolve as if it was isolated.To see the competing effects of the periastron precession and the ZLK mechanism, it is convenient to neglect the GW backreaction and analyze under which conditions ⟨/⟩ = 0 and ⟨/⟩ = 0 in Eq. (8.6), where the second equation is automatically satisfied by setting  = /2.One finds which replaces Eq.(7.14) when the periastron precession is taken into account.Due to the presence of Ω (GR) ZLK , this generalizes the Newtonian formula found in Blaes et al. (2002).The window of values for the critical inclination that triggers the ZLK resonance is reduced by the presence of the periastron effect.This makes it complicated to find a condition on r that would result in the condition cos 2  0 > 0. To illustrate this, consider the simpler case of a nonspinning SMBH for which we find the following condition r3 In the limit in which this bound is saturated, the critical inclination value is close to  0 ≈ 90 • , so that the minimum and maximum values for the inclination almost coincide.It is clear from (8.12) that, consistently, the maximum value of the eccentricity is lowered compared to the case in which the periastron precession is neglected.

Peak Frequency
For eccentric binaries, the GW spectrum is spread across an infinite number of harmonics Peters & Mathews (1963), with frequencies that are integer multiples  of the fundamental Keplerian frequency 1/(2) √︁  / 3 , and peaked approximately at Wen ( 2003) (1 + ) 1.1954 .(8.14)In the evolution of an isolated binary the gravitational backreaction, encoded in Eqs.(8.2), would contribute to circularize the orbit and move the peak frequency towards the usual  = 2 harmonic long before the merger takes place Maggiore (2007).Under the influence of an external body, however, the ZLK mechanism can provide large eccentricity oscillations in highly-inclined binaries and can enhance the peak frequency to values high enough to enter in the sensitivity band of future space-based GW detectors Antonini & Perets (2012b).

Binary merger time in the weak field limit
In Sec. 4 we studied how GR effects induced by the presence of a spinning SMBH can lead to a significant enhancement of the frequency for the ZLK resonance in the case of a Newtonian binary system moving on a circular geodesic.By superimposing the ZLK mechanism with the PN dynamics of the inner BBH system, according to the discussion made in the previous subsection, we are now in a position to study how the BBH merger time is influenced by the presence of an external SMBH.
We will start the study of the evolution equations (8.6) by considering a case in which we are in the weak gravity regime, regarding the influence of the SMBH on the binary system.This regime has already been studied previously, e.g. in Randall & Xianyu (2018), but here we use it to provide a baseline for Section 8.3, where we find novel results for the strong field regime, also described by Eqs.(8.6).
We investigate the evolution equations (8.6) by solving them numerically.For the weak-field example of this section we choose fixed masses  1 =  2 = 10  ⊙ and  3 = 4 × 10 6  ⊙ for the three black holes. 10We use these values to connect to previous literature (e.g.Ref. Randall & Xianyu (2018)) and since  3 is roughly the mass of Sagittarius A*.Moreover, we choose the initial conditions for the inner binary to be  0 = 0.1,  0 = 0.1 AU ,  0 = 0 • and  0 = 89.9• , so that we only need to specify the outer orbit parameters.11 For a non-spinning SMBH with the above-mentioned choice for  3 we have the ISCO radius rISCO = 6  3 / 2 ≈ 0.24 AU.Naturally, the strong gravity effects from the presence of the SMBH are most significant at the ISCO.We choose here the outer orbit radius to be r = 120 AU ≈ 500 rISCO which means we are within the weak-field regime (concerning the gravitational influence of the SMBH on the binary system).In this regime, the redshift factor is   ≈ 1 according to Eq. (3.7).At the same time, r = 120 AU ensures we are well below the tidal breakup limit (8.7).
From Eq. (8.10) one can observe that, with our choice for the orbital parameters, it would take approximately ∼ 10 10 years for an isolated binary to coalesce.In Fig. 9 we show our numerical solution, which presents the characteristic quasi-periodic behavior of the BBH system influenced by an external SMBH.From Fig. 9 it is immediate to see that the time for the merger to occur is drastically reduced to order ∼ 10 2 years by the presence of the eccentricity oscillations associated with the weak-gravity ZLK mechanism.
It is possible to distinguish three different phases that characterize the merger process in the example we consider.At earlier stages ( PN ≫  ZLK ) the ZLK mechanism dominates the dynamics and, due to the inner binary's high initial inclination, leads to large-amplitude oscillations for the eccentricity and inclination.More specifically one can see from Fig. 9 that, in agreement with the ZLK mechanism, the maximum eccentricity  (minima of the black curve) corresponds to a minimum of the inclination  (maxima of the cyan curve).
The temporary increase of the eccentricity up to  ≈ 0.999 leads to an efficient loss of angular momentum via GW emission in a small amount of time, which translates into the typical step-like monotonic decreasing of the semi-major axis  (see Fig. 10), and in a progressive reduction of the ZLK oscillations amplitude.
In the intermediate stage, the ZLK timescale becomes comparable with the periastron precession timescale ( ZLK ∼  PN ), and the eccentricity/inclination oscillations are hampered by the PN contributions in the BBH dynamics.
In the final phase ( PN ≪  ZLK ), the ZLK resonance is suppressed.In .Semi-major axis as a function of the asymptotic time t in the non-spinning case for r = 500 rISCO ≈ 120 AU.The step-like behavior is a consequence of the loss of angular momentum, which is maximized when the eccentricity oscillations are close to their maximum values (the minima of the black curve in Fig. 9).other words, the BBH evolves as if the system were isolated and characterized by a high initial eccentricity inherited by the ZLK mechanism.The Peter's contributions in Eqs.(8.6) then become dominant, causing a fast orbital circularisation and a prompt decreasing in the semi-major axis of the BBH system via GW emission.Even though a realistic description of the final stages of the inspiral phase would require much more advanced analytic or numerical frameworks, here we limit to convene that the BBH merger takes place when both  and  vanish.

Binary merger time in the strong field limit
We are now ready to study the evolution equations (8.6) for configurations in which one can see the strong field effects in full, regarding the influence of the SMBH on the binary system.Indeed, we shall consider configurations for which the radius r of the outer orbit is comparable with the ISCO radius r  ISCO of the SMBH.Since we are in a strong-gravity regime, it means that we will also be able to account for the dependence on the spin of the SMBH and in the following we will describe both the case where the BBH system moves on a co-rotating ( = +1) circular equatorial orbit or a counter-rotating ( = −1) circular equatorial orbit.
For the configurations of this subsection, we choose the mass of the SMBH to be  3 = 2 × 10 9  ⊙ which happens to be roughly the mass of M87*.This choice allows us to explore the situation where a stellar BBH system with masses  1 =  2 = 10  ⊙ and semi-major axis  = 0.1 AU is near or at the ISCO of the SMBH, while at the same time being sufficiently far away from the tidal breakup condition (8.7), see Fig. 8.
We consider three cases for the spin of the SMBH: the non-spinning case  = 0, and the highly spinning cases  = 0.95 with co-rotation ( = 1) and counter-rotation ( = −1).12Finally, we choose the outer orbit radius to be r = 180 AU, respectively corresponding to r ≈ 4.68 r+ ISCO for co-rotating and r ≈ 1.02 r − ISCO for counter-rotating orbits.Notice that this value is also consistent in the case of a non-spinning SMBH of the same mass, since it would correspond to r ≈ 1.51 rISCO .We emphasize that the chosen values of the parameters ensure that the BBH system is not subjected to a tidal breakup and can experience strong gravity effects due to the tidal interaction with the external SMBH.
We have solved the evolution equations (8.6) for these choices of configurations.The result is depicted in Fig. 11 where we show the evolution of 1 −  2 as a function of the asymptotic time t for four different initial inclination angles  0 = 88.5 • , 89 • , 89.3 • , 89.5 • .Moreover, we choose the remaining initial conditions for the inner binary to be  0 = 0.1,  0 = 0.1 AU and  0 = 0 • .As one can see, in Fig. 11 we compare the BBH eccentricity evolution for the case where the inner binary is moving on a co-rotating and a counter-rotating circular equatorial Kerr geodesic (respectively purple and orange curves), with the case in which the external body is a non-spinning SMBH (black curve) described in the Newtonian point particle approximation (grey curve).This is to highlight the effect of the SMBH spin.
We point out that when the BBH is at the ISCO of the SMBH, or close to it, the GR effects become quite significant and the Newtonian point particle approximation is no longer valid to describe the dynamics of the system.Here we use the point particle The picture shows the evolution of the eccentricity for a BBH system perturbed by an external SMBH of mass  3 = 2 × 10 9  ⊙ .The four panels correspond to different values of the initial inclination  0 .The gray curve represents the case in which the external black hole is treated in the Newtonian point particle approximation, the black curve includes general-relativistic effects associated with a non-spinning black hole whereas the purple and orange curves describe a Kerr black hole with spin parameter  = 0.95 in the co-rotating (purple line) and counter-rotating (orange line) case.Here we choose r = 180 AU.In terms of the ISCO radius it means r ∼ 1.5 rISCO for  = 0 while for  = 0.95 we have r ∼ 4.7 r+ ISCO in the co-rotating case and r ∼ 1.02 r − ISCO in the counter-rotating case.The four panels show that small variations in the initial inclination  0 correspond to huge variations in the merger time, both for the GR and the Newtonian case.The amplitude of the ZLK oscillations shows that the maximum eccentricity grows as  0 → 90 • .In particular, note that when the the ZLK mechanism is most efficient, i.e. when  0 ∼ 90 • , strong-gravity effects contribute to accelerate the BBH merger compared to the Newtonian case (bottom right panel).Instead, for smaller values of  0 , the ZLK mechanism is less efficient and the merger time is thus slowed down by the redshift factor (top panels and bottom left panel).description for the SMBH to illustrate the comparison between our novel results valid in the full GR regime and the results obtained in the standard point particle approximation.We can conclude that describing the SMBH with the Kerr metric, instead of as a point particle, leads to effects that significantly impact the dynamics of the BBH system.
The four different panels with four different initial inclination angles in Fig 11 aim to illustrate the strong dependence of the merger dynamics on the initial inclination angle  0 : it is possible to observe that the more the BBH system is inclined, the more the ZLK mechanism is effective in reaching high values for the eccentricity and, thus, in reducing the merger time.This feature of the ZLK mechanism remains true regardless of whether the SMBH is described in a Newtonian or in a general-relativistic manner.Indeed, in all cases, one can observe that the maximum eccentricity grows as  0 increases approaching 90 • and, consequently, a smaller number of ZLK oscillations is needed in order to boost the merger.
In particular, by making use of the terminology adopted in Ref. Randall & Xianyu (2018), the two top panels and the bottom left panel in Fig. 11 describe a slow-merger dynamics, i.e. when the BBH system undergoes more than one ZLK cycle before merging, whereas the bottom right panel depicts a fast-merger dynamics, in which the BBH system coalesces in just one cycle of the ZLK mechanism, and the merger time t★ is directly given by t★ ≈    ZLK Randall & Xianyu (2018).
Another feature associated with strong gravity effects that emerges from Fig. 11 is related to the impact of the redshift factor on the merger time.We recall that in the evolution equations of the system, Eqs.(8.6), the details of the description adopted for the SMBH enter through the ZLK frequency Ω (GR) ZLK which now, from the point of view of an observer in the asymptotic region of the Kerr space-time, translates into the quantity Ω (∞) ZLK defined in Eq. (7.4), where the role of the redshift factor is evident.When the SMBH is described using the Kerr metric (or the Schwarzschild metric for non-spinning black holes), it is the redshift factor which is responsible for slowing down the evolution of the system and causing an increase in the merger time compared to the case where one uses the standard Newtonian point particle description for the SMBH.The effect of the gravitational redshift is particularly evident in the two top panels of Fig. 11, where the BBH evolution is characterized by a large number of ZLK cycles.Here one sees that for example the orange curve, with a redshift factor   ≈ 1.24, corresponds to a slower merger time compared to the black curve with redshift factor   ≈ 1.22 and so on.
Finally, Fig. 11 highlights another consequence of being in a stronggravity regime.In Sec.7.1 we discussed how a general-relativistic description of the external SMBH can lead to a significant increase in the local ZLK frequency when compared to the Newtonian description.As opposed to the gravitational redshift, thus, the stronggravity effects in the ZLK frequency contribute to catalyzing the BBH merger at earlier times than in the Newtonian case.In particular, from Fig. 2 it is immediate to see that for a fixed value of the radial coordinate r and of the spin parameter , the deviation of Ω ZLK is greater for counter-rotating orbits.The interplay between the GR enhancement of the local ZLK frequency, which accelerates the merger, and the gravitational redshift, which instead tends to slow down the merger, is particularly evident in the bottom right panel of Fig. 11.For the co-rotating case We depict here the eccentricity of a BBH system at the ISCO of an SMBH with  = 0.The solid black and red lines represent the system in terms of the local time  and the asymptotic time t, respectively.The masses are chosen as  1 =  2 = 10  ⊙ and  3 = 2 × 10 9  ⊙ .In contrast, the dashed black line is the result obtained by using the Newtonian point particle approximation for the SMBH with the same outer orbit radius 119 AU.The initial conditions for the BBH orbital parameters are  0 = 0.1,  0 = 0.1 AU,  0 = 0 • and  0 = 89.9• .
(purple), the gravitational redshift   ≈ 1.2 almost entirely compensates the GR frequency enhancement ZLK ≈ 1.This makes the purple curve almost indistinguishable from the Newtonian one (the gray curve).In the non-spinning (black) and counter-rotating (orange) case, instead, one respectively has Ω ZLK ≈ 1.57 respectively so that the strong-gravity effects in the ZLK frequency dominate over the gravitational redshift, thus accelerating the merger.

GW peak frequency at the ISCO
As already remarked, the strong-gravity effects should be maximal when the BBH is closest to the SMBH, which, for equatorial orbits, is when the outer orbit is at the ISCO.In this section we analyze this scenario in detail.The scenario is realistic because interactions with the surrounding material inside the accretion disc can lead compact objects to migrate in the neighborhood of the last stable orbit and remain confined there for long enough timescales, that allow binaries to form and coalesce.For further details about these migration-trap mechanisms we refer the reader to Ref.We begin with the case of a high initial inclination,  0 = 89.9• , which maximizes the ZLK mechanism.We have depicted this in Fig. 12, showing the evolution of the BBH eccentricity in terms of the proper time  (solid black line) and the asymptotic time t (red line) compared with the case in which the external SMBH is treated in a Newtonian point particle approximation (dashed curve).The masses are chosen as in Sec.8.3, with  1 =  2 = 10 ⊙ for the BBH and  3 = 2 × 10 9  ⊙ for the SMBH, ensuring that it is possible to approach the ISCO while avoiding tidal breakup.For simplicity we have chosen to focus on the non-spinning case  = 0, for which r ≈ 119 AU, and   = √ 2.
We see from Fig. 12 that the combination of a high initial inclination angle and the strong-gravity effects at the ISCO speeds up the merger, which occurs in just a single ZLK oscillation.Accordingly, Fig.
12 also shows that the maximum eccentricity is enhanced when including GR effects, reaching extremely high values 1 −  2 max < 10 −4 , and further contributing to reduce the merger time (t ≈ 0.025 yrs in terms of the asymptotic time) compared to the examples presented before in Fig. 9. Finally, it is important to observe that, in terms of the proper time  (black curve in Fig. 9), the merger time is halved, as compared to the merger time when using the Newtonian point particle approximation for the SMBH (dashed curve in Fig. 9), in perfect agreement with the analytical prediction we derived in Sec.7.1 that at the ISCO of an SMBH one has Ω (GR) ZLK /Ω (N) ZLK = 2.Other than the merger time, it is of obvious importance to understand the frequency spectrum of the GWs emitted from the BBH, to see further observational signatures of placing the BBH at the ISCO of the SMBH.A measure of this is the peak frequency  GW of Eq. (8.14) that marks the dominant frequency in the GW spectrum.
As pointed out in previous works Antonini & Perets (2012a); Randall & Xianyu (2018), an exhaustive characterization of ZLK-boosted triple systems via GW signal observations would require combining future space-based interferometers such as LISA, to probe the earlier phase of the merger when the two black holes are widely separated and the ZLK mechanism dominates, together with current ground-based facilities (such as LIGO-Virgo-KAGRA), which can detect the merger of an already-circularized BBH system at late stages.
As is evidenced by Figs.11 and 12, highly inclined ( 0 > 89.4 • ) binaries in the strong-gravity regime are characterized by fast-merger dynamics.The GW signal in these cases would result in a single pulse observed by both ET and LISA, promptly followed by a chirp signal detected by LIGO-Virgo-KAGRA.
Lowering the initial inclination angle  0 should instead provide more ZLK oscillations, as seen in Fig. 11, since the ZLK mechanism is less hampered by the periastron precession.This could in turn provide a longer signal to be observed by ET and LISA.Indeed, one would expect the detection of a BBH merger by LIGO-Virgo-KAGRA to be preceded by a series of repeated pulses detected by ET and LISA at earlier times, that mark the presence of the ZLK cycles.
To show this, we plot in Figs. 13 and 14 the GW peak frequency  GW of the GW emission from the BBH as a function of time given a slightly lower initial inclination angle  0 = 89.4• .We consider two cases, one for zero spin  = 0 of the SMBH (Fig. 13) and one with non-zero co-rotating spin  = 0.3 (Fig. 14).In both cases the BBH moves along the ISCO of the SMBH.In the figures we compare the results of our novel GR description for the SMBH to what one would have obtained by treating the SMBH as a Newtonian point particle.In detail, the black curves in Figs. 13 and 14 mark the peak frequency  GW measured with proper time  on the left while the red curve marks the asymptotically measured peak frequency  GW with asymptotic time t of the Kerr space-time on the right.These curves are then compared to the peak frequency  GW that one obtains from SMBH as Newtonian point particle, marked by the black dashed curves.
While  = 0 in Fig. 13, we have turned on the spin of the SMBH in Fig. 14 with  = 0.3 corresponding to a co-rotating spin  = 1.Turning on the spin parameter for the SMBH allows binary systems to access regions of space-time which would be prohibited in the non-spinning case.More specifically, the co-rotating ISCO for Kerr lies closer to the black hole than in the non-spinning case.Left Panel: we compare the peak frequency when treating the SMBH by a full GR description using the Schwarzschild metric (solid black) to what one obtains from a Newtonian point particle description (dashed black).One notices that there is twice the number of peaks in the GR description compared to the Newtonian approximation, consistent with the fact that Ω (GR) ZLK at the ISCO.The gravitational redshift is not considered in this plot.Right Panel: the same comparison as in the left panel, but in terms of the asymptotic time t (red curve).This includes the gravitational redshift   .Notice that the inclusion of the redshift factor makes the merger time quite different in the two cases (4.8 months in the Newtonian case compared to 6.9 months in the GR case) and reduces the maximum value reached by the frequency.The blue curve reproduces the ratio between the PN precession and the ZLK timescales.Finally, the pink and green bands in the background respectively show the range of frequency detectable for the upcoming interferometers ET (1 Hz <  GW < 10 kHz) and LISA (0.001 Hz <  GW < 0.1 Hz) .13 between the Newtonian and the GR case but now for a spinning SMBH with  3 = 2 × 10 9  ⊙ and  = 0.3 (with  = 1).In this case, the co-rotating ISCO is located at r ≈ 99 AU.As in Fig. 13 the left and the right panels differ for the inclusion of the redshift factor, which in this case amounts to   ≈ 1.41 .
Comparing the spinning case to the non-spinning case, we observe that the binary system undergoes a larger number of ZLK cycles before the merger occurs.
For the left panels of Figs. 13 and 14 we have depicted the evolution of the peak frequency in terms of the proper time, which is the local time of the BBH.These panels highlight the fact that in the early phases of the merger when the period of ZLK oscillations is much shorter than the one of the periastron precession ( PN ≫  ZLK , as marked with the blue curve), the future interferometers ET and LISA would detect twice the number of pulses entering their frequency band (depicted with pink and green stripes respectively) compared to a Newtonian point particle description of the SMBH.This is consistent with our analytic prediction that at the ISCO of any SMBH, we have Ω ZLK .While the LISA sensitivity band (10 −3 Hz <  GW < 10 −1 ) detects GWs emitted only in the earlier phases of the inspiral, it is interesting to notice that the ET frequency band (1 Hz <  GW < 10 kHz) allows to observe GWs up to the final stages of the BBH merger.
In the right panels of Figs. 13 and 14 we depict instead the evolution of the peak frequency in terms of the asymptotic time t, i.e. the time appropriate for the GW detector.The redshift factor included in the asymptotic time is not only responsible for reducing the maximum peak frequency reached by the GW emitted by the binary, but it also shifts the positions of the pulses and significantly affects the merger time.
Finally, Figs. 13 and 14 show how in the last part of the merger ( PN ≪  ZLK ) the GW peak frequency enters the detectable band of LIGO/Virgo/Kagra (  GW > 10 Hz).

CONCLUSIONS
In this paper, we investigated the dynamics of a BBH in point particle approximation subject to the strong gravitational field of an SMBH as described by the Kerr metric.We analyze this in the hierarchical regime where the mass of the SMBH is much bigger than the masses of the black holes in the binary system and where the size of the binary system is much smaller than the curvature scale of the SMBH.In this way, the BBH system experiences the presence of the SMBH through tidal forces which we include up to the quadrupole order.

Figure 1 .
Figure 1.Illustration of orbital parameters.The green plane represents the reference plane while the intersection between the red line and the orbit constitutes the periapsis and the intersection between the blue line and the orbital plane provides the ascending node.
Figure 2. Diagram of the ratio Ω (GR) ZLK /Ω (N) ZLK versus the distance  from the SMBH.Different colors label various values of the black hole spin  where the solid and dashed lines represent co-rotating ( = +1) and counter-rotating ( = −1) orbits, respectively.Each curve terminates at the ISCO.

Figure 6 .
Figure6.Here the binary system is always placed at the ISCO.The blue line describes the counter-rotating case while the red line describes the co-rotating case.

Figure 7 .
Figure 7.In this figure we set  1 =  2 = 10  ⊙ ,  3 = 2 × 10 9  ⊙ and  = 0.1 AU. (Left panel) Absolute value of the ratio between the gyroscope precession frequency and the ZLK frequency as a function of the inner binary eccentricity .The different curves are obtained for different values of the dimensionless spin parameter  of the SMBH.(Right panel) Absolute value of the ratio between the gyroscope precession frequency and the ZLK frequency as a function of the dimensionless spin parameter  of the SMBH, for different values of the eccentricity of the inner binary .The solid and dashed coloured curves are respectively for co-rotating ( = +1) and counter-rotating ( = −1) outer geodesics.The value  = 1, depicted with the dot-dashed line, marks the trans-adiabatic regime in which Ω (∞)  ≈ Ω (∞) ZLK .The rightmost curve in the left panel, corresponding to  = 0.99999, is reported to show that  → ∞ for  → 1, in agreement with the behaviour presented in the right panel and the discussion made below Eq. (4.3).

Figure 8 .
Figure8.The colored regions represent regions of the parameter space for the SMBH which are consistent with the tidal breakup condition r/ > (3 3 / ) 1/3 .The green, orange and blue colors respectively identify r = r  ISCO , r = 10 r  ISCO and r = 100 r  ISCO , whereas the left and right panels distinguish co-rotating ( = +1) and counter-rotating orbits ( = −1).The plots are obtained by fixing  1 =  2 = 10  ⊙ and  = 0.1 AU, as representative values for a stellar-mass BBH system.

Figure 9 .
Figure9.We focus on the spinless case for the SMBH, i.e.  = 0.Here r = 120 AU, which corresponds to approximately r ≈ 500 rISCO for a non-spinning black hole with mass  3 = 4 × 10 6  ⊙ .We combine the effect of the ZLK mechanism with the periastron precession and GW emission.The picture highlights that the maximum values of the BBH eccentricity  corresponds to the minimum values for its orbital inclination  and viceversa.
Figure12.We depict here the eccentricity of a BBH system at the ISCO of an SMBH with  = 0.The solid black and red lines represent the system in terms of the local time  and the asymptotic time t, respectively.The masses are chosen as  1 =  2 = 10  ⊙ and  3 = 2 × 10 9  ⊙ .In contrast, the dashed black line is the result obtained by using the Newtonian point particle approximation for the SMBH with the same outer orbit radius 119 AU.The initial conditions for the BBH orbital parameters are  0 = 0.1,  0 = 0.1 AU,  0 = 0 • and  0 = 89.9• .

Figure 13 .
Figure13.GW peak frequency  GW emitted by the BBH system as a function of the proper time  (left panel) and asymptotic time t (right panel).The BBH is placed on the ISCO of a non-spinning SMBH with mass  3 = 2 × 10 9  ⊙ , corresponding to r ≈ 120 AU.The two black holes in the binary system have masses of  1 =  2 = 10  ⊙ with an initial separation of  0 = 0.1 AU, an initial inclination angle  0 = 89.4• , initial eccentricity  0 = 0.1 and  0 = 0 • .Left Panel: we compare the peak frequency when treating the SMBH by a full GR description using the Schwarzschild metric (solid black) to what one obtains from a Newtonian point particle description (dashed black).One notices that there is twice the number of peaks in the GR description compared to the Newtonian approximation, consistent with the fact that Ω

𝑇Figure 14 .
Figure14.Same peak frequency comparison as in Fig.13between the Newtonian and the GR case but now for a spinning SMBH with  3 = 2 × 10 9  ⊙ and  = 0.3 (with  = 1).In this case, the co-rotating ISCO is located at r ≈ 99 AU.As in Fig.13the left and the right panels differ for the inclusion of the redshift factor, which in this case amounts to   ≈ 1.41 .
, with the quantity Ω  =     −     + /  = −Ω it is natural to choose this to be the equatorial plane of the Kerr black hole.An arbitrary orientation of the inner binary is therefore derived by performing the rotation  = ( cos  x+  sin  ŷ), and by fixing  = 0, /2 one gets û =(cos  cos  − cos  sin  sin ) x + (cos  sin  + cos  sin  cos ) ŷ (5.5)