Issue 
A&A
Volume 683, March 2024



Article Number  A135  
Number of page(s)  12  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/202347920  
Published online  15 March 2024 
Threebody encounters in black hole discs around a supermassive black hole
The disc velocity dispersion and the Keplerian tidal field determine the eccentricity and spinorbit alignment of gravitational wave mergers
^{1}
Niels Bohr International Academy, Niels Bohr Institute, Blegdamsvej 17, 2100 Copenhagen, Denmark
email: aatrani@gmail.com
^{2}
Research Center for the Early Universe, School of Science, The University of Tokyo, Tokyo 1130033, Japan
^{3}
Okinawa Institute of Science and Technology, 19191 Tancha, Onnason, Okinawa 9040495, Japan
^{4}
Dipartimento di Fisica ‘G. Occhialini’, Università degli Studi di MilanoBicocca, Piazza della Scienza 3, 20126 Milano, Italy
^{5}
INFN, Sezione di MilanoBicocca, Piazza della Scienza 3, 20126 Milano, Italy
Received:
8
September
2023
Accepted:
18
December
2023
Context. Dynamical encounters of stellarmass black holes (BHs) in a disc of compact objects around a supermassive BH (SMBH) can accelerate the formation and coalescence of BH binaries. It has been proposed that binary–single encounters among BHs in such discs can lead to an excess of highly eccentric BH mergers. However, previous studies have neglected how the disc velocity dispersion and the SMBH’s tidal field affect the threebody dynamics.
Aims. We investigate the outcomes of binary–single encounters considering different values of the disc velocity dispersion, and examine the role of the SMBH’s tidal field. We then demonstrate how their inclusion affects the properties of merging BH binaries.
Methods. We performed simulations of fourbody encounters (i.e. with the SMBH as the fourth particle) using the highly accurate, regularised code TSUNAMI, which includes postNewtonian corrections up to order 3.5PN. To isolate the effect of the SMBH’s tidal field, we compared these simulations with those of threebody encounters in isolation.
Results. The disc velocity dispersion controls how orbits in the disc are aligned and circular, and determines the relative velocity of the binary–single pair before the encounter. As the velocity dispersion decreases, the eccentricity of postencounter binaries transitions from thermal to superthermal, and binaries experience enhanced hardening. The transition between these two regimes happens at disc eccentricities and inclinations of order e ∼ i ∼ 10^{−4}. These distinct regimes correspond to a disc dominated by random motions (e ∼ i ≳ 10^{−4}) and one dominated by the Keplerian shear (e ∼ i ≲ 10^{−4}). The effect of the SMBH’s tidal field depends on the velocity dispersion of the disc. When the velocity dispersion is low, the resulting binaries are less eccentric compared to isolated encounters. Conversely, binaries become less eccentric compared to isolated encounters at high velocity dispersion. This also affects the number of BH mergers.
Conclusions. The inclusion of the SMBH’s tidal field and the disc velocity dispersion can significantly affect the number of GW mergers, and especially the number of highly eccentric inspirals. These can be up to ∼2 times higher at low velocity dispersion, and ∼12 times lower at high velocity dispersions. The spin–orbit alignment is influenced by the tidal field exclusively at high velocity dispersions, effectively inhibiting the formation of antialigned binary BHs. Nonetheless, encounters in randommotiondominated discs around a SMBH are still more effective in producing GW mergers compared to those occurring in spherically symmetric nuclear star clusters without an SMBH.
Key words: black hole physics / gravitation / gravitational waves / methods: numerical / celestial mechanics / stars: kinematics and dynamics
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Discs of stellarmass black holes (BHs) around supermassive BHs (SMBHs) may be common in the Universe, as suggested by both observations (albeit indirectly) and theoretical analysis. On the observational side, WolfRayet and O/B stars have been observed orbiting SgrA* in one or more discs (Genzel et al. 2003; Levin & Beloborodov 2003; Paumard et al. 2006; Lu et al. 2009; Bartko et al. 2009; Yelda et al. 2014; von Fellenberg et al. 2022), which, at the end of their stellar lifecycle, will leave a disc of compact objects. Additionally, a disc of earlytype stars has been observed around the SMBH in M31, which suggests that ubiquitous mechanisms may form nuclear discs (Lauer et al. 1993, 2005; Tremaine 1995; Peiris & Tremaine 2003; Bender et al. 2005; Brown & Magorrian 2013).
Numerical studies have proposed a variety of ways to form nuclear discs around SMBHs, such as fragmentation of tidally disrupted molecular clouds (Mapelli et al. 2012; Trani et al. 2016, 2018; MastrobuonoBattisti et al. 2019), gas funnelling during galaxy mergers (Hopkins & Quataert 2010a,b), anisotropic mass segregation (Szölgyén et al. 2019; Kamlah et al. 2022; Máthé et al. 2023), globular cluster infall and disruption (ArcaSedda et al. 2015, 2018), star formation in active galactic nuclei (AGN) discs (Cantiello et al. 2021; Jermyn et al. 2022), or gasfriction capture in AGN discs (Yang et al. 2019; Fabj et al. 2020; MacLeod & Lin 2020; Nasim et al. 2023; Generozov & Perets 2023). In the last two cases, the nuclear disc is embedded into the AGN’s gaseous disc, which is in contrast to dry nuclear discs formed by the other processes.
Discs of compact objects around an SMBH are excellent candidate environments for producing gravitational wave (GW) sources for both groundbased and spaceborne GW detectors (Abbott et al. 2019, 2021a,b; The LIGO Scientific Collaboration & The Virgo Collaboration 2021; The LIGO Scientific Collaboration 2023a,b). In particular, the high stellar density in galactic nuclei makes them rich in dynamical interactions, which can give rise to unique signatures in the GW waveform that can be used to infer their origin (see e.g. Mapelli 2021; Spera et al. 2022). For instance, the gravitational interaction of BH binaries with the SMBH may give rise to oscillations in the binary eccentricity, called von ZeipelKozaiLidov (ZKL) oscillations, which can accelerate and trigger GW coalescence of the binary (von Zeipel 1910; Lidov 1962; Kozai 1962; Blaes et al. 2002; Miller & Lauburg 2009; Antonini & Perets 2012; Liu et al. 2015; Naoz 2016; Shevchenko 2017; Hoang et al. 2018; Martinez et al. 2022; Trani et al. 2022; Camilloni et al. 2023). In addition to ZKL oscillations, gravitational interactions might occur in the disc itself in the form of fewbody encounters amongst single objects and binaries. Such interactions can lead to single–single GW captures or accelerate the formation and merger of BH binaries (Stone et al. 2017; Bartos et al. 2017; Leigh et al. 2018; Secunda et al. 2019, 2021; Li et al. 2022, 2023; Boekholt et al. 2023; DeLaurentiis et al. 2023; Rowan et al. 2023; Rom et al. 2023; Qian et al. 2024). Samsing et al. (2022) argued that binary–single encounters in AGN discs may give rise to highly eccentric mergers because of the low relative velocity between bodies in the disc. However, their study did not consider the effect of the tidal field of the SMBH on the threebody dynamics in the disc. The role of this tidal field was first studied by Trani et al. (2019b) in the context of binary encounters with BHs from an isotropic cusp, but a detailed study on how it affects binary–single encounters in a disc is still missing.
In this paper, we focus on the outcome of threebody encounters between single and binary BHs in a dry nuclear disc around an SMBH by means of highly accurate fewbody simulations. Particularly, we describe (i) the effect of different disc velocity dispersions on the binary outcomes (Sect. 3.1) and (ii) the effect of the SMBH tidal field compared to threebody encounters in its absence (Sect. 3.2). We put our findings into the context of GW sources in Sect. 3.3. Our numerical methods and simulation setup are described in Sect. 2.
2. Numerical setup
2.1. Initial conditions
We run fourbody simulations of an encounter between a BH binary and a single BH around an SMBH, assuming that both objects are part of the same nuclear disc. The initial setup of the system is outlined in Fig. 1. The binary orbits around the SMBH with a semimajor axis and eccentricity a_{out} and e_{out} (hereafter we refer to this as the outer orbit). Likewise, the binary is characterised by an ‘inner’ semimajor axis and inner eccentricity a_{inn} and e_{inn}, respectively. The orbit of the single BH around the SMBH has a semimajor axis a_{sin} and eccentricity e_{sin} and is inclined by i_{sin} with respect to the orbital plane of the binary.
Fig. 1. Scheme of the initial setup of our simulations. The centre of mass of the binary and the single are set to encounter at a location R_{enc} with impact parameter b ∈ (0, 2a_{inn}). The eccentricities and mutual inclinations of the two orbits around the SMBH are drawn from a Rayleigh distribution (ℛ) with scale parameter σ and σ/2, respectively. 
The eccentricities e_{out} and e_{sin} are drawn from the same distribution, which is consistent with the binary and the single being part of the same disc, while i_{sin} relates to the opening angle of the disc. More generally, eccentricity and inclination express the deviation from an aligned circular orbit, and parametrise the random motion of Keplerian orbits in a disc in the same way as the velocity dispersion describes the random motion of stars in an isotropic stellar cluster. Specifically, the random velocity is related to e and i as (Lissauer & Stewart 1993)
where v_{K} is the Keplerian circular velocity. The random motions of particles in a disc will grow because of twobody gravitational interactions, keeping the ratio ⟨e⟩^{1/2}/⟨i⟩^{1/2} ≈ 2 approximatively constant, with eccentricities and inclinations growing with time as e, i ∝ t^{1/4}. In this regime, both the inclination and eccentricity distributions are accurately modelled by a Rayleigh distribution, a finding corroborated by both simulations and theoretical analyses (Ida & Makino 1992; Stewart & Ida 2000; Kokubo & Ida 2002). Specifically, a Rayleigh distribution in Keplerian orbital elements results from assuming that the random velocity field is described by a triaxial Gaussian distribution. This is only valid in the regime of small inclinations and eccentricities, where the orbits are well approximated by linear epicycles.
For this reason, we parametrise the disc velocity dispersion with a single dimensionless parameter σ ≡ σ_{e} = 2σ_{i}, where σ_{e} and σ_{i} are the scale parameters of the Rayleigh distributions (ℛ) of eccentricity and inclination. We choose σ = 10^{−1}, 10^{−2}, 10^{−3}, 10^{−4}, 10^{−5}, 10^{−6}, 10^{−7}, 0, corresponding to simulation sets σn, where n goes from −1 to −7, and then to 0. The sets correspond to eccentricity and inclinations of order (e, i)∼(1.25, 36° ) × σ, which is going from e, i ∼ 0.125, 3.6° for σ1, to e, i ∼ 1.25 × 10^{−7}, 3.6° ×10^{−7} for σ7. Values sampled from a Rayleigh distribution can in principle go above the physical upper limits of eccentricity and inclination (1 and π, respectively). We make sure that we discard and resample such values whenever they occur, which is extremely unlikely even for the set with σ = 0.1 (the probability density at e = 1 for σ_{e} = 0.1 is 10^{−20}).
These eight sets represent discs with increasing velocity dispersion, and can be interpreted as different evolutionary stages of the same disc, whose orbits become more eccentric and inclined due to mutual gravitational interactions. Our σ parameter can be thought of as an analogue to the W_{0} parameter of MichieKing globular clusters models (Michie 1963; King 1966), which controls their central density, and for which higher values of W_{0} correspond to more evolved clusters due to twobody relaxation (Spitzer 1987).
We assume the binary inner orbit to be initially circular (e_{inn} = 0). As the outer orbit of the binary defines a plane and therefore the direction of the orbital angular momentum of the binary relative to the SMBH, here we assume the inner orbit to be prograde and coplanar^{1}. The binary semimajor axis is drawn from a loguniform distribution (log 𝒰) between 0.1 au and half of its Hill radius R_{Hill}, which is given by
where m_{out} = m_{1} + m_{2} is the total mass of the binary and M_{SMBH} is the mass of the SMBH. We sample the mass of the stellarmass BHs from a loguniform distribution between 10 and 50 M_{⊙} and set the mass of the SMBH to 4.31 × 10^{6} M_{⊙}, which is the inferred mass of SgrA* (Gillessen et al. 2009, 2017). The true anomaly of the inner BH binary is uniformly sampled between 0 and 2π.
The binary and the single orbits cross at a distance R_{enc} from the SMBH. We set this distance by drawing it from a f(R_{enc})∝R_{enc} distribution between 0.03 pc and 0.1 pc; this latter range is inspired by the stellar disc observed around SgrA*. Given R_{enc}, e_{out}, and e_{sin}, we sample a_{sin} and a_{out} from the same R_{enc} distribution, adding the constraint that the orbits need to cross at R_{enc}. The orbit of the binary is fixed so that it passes through R_{enc}, and we add an impact parameter with respect to this location in order to determine where the orbit of the single will pass. The impact parameter is set to be orthogonal to the tangential velocity of the single at the encounter location, and we sample its magnitude from a f(b)∝b distribution between 0 and 2a_{inn}.
Once R_{enc}, e_{sin}, e_{out}, a_{sin}, a_{out}, and i_{sin} are fixed, the orbital parameters of the single and the centre of mass of the binary are uniquely determined. We then start the simulations with the binary and the single some time ΔT_{enc} before the encounter along their orbits. Here, we chose ΔT_{enc} = min(P_{out}, P_{sin})/16, where P_{out} and P_{sin} are the orbital periods of the binary and the single around the SMBH. We increase ΔT_{enc} if the binary and the single are already within the sum of their Hill radii at the beginning of the simulations, or if they are closer than 100 times the binary semimajor axis. In set σ0, the orbits are perfectly circular and aligned, meaning that the impact parameter is also the difference in semimajor axis between the two orbits; for more details, see Trani et al. (2019a,b). We run a total of 5 × 10^{4} realisations per set. The initial conditions are summarised in Table 1.
Initial conditions for our sets of simulations.
2.2. The TSUNAMI code
We run the simulations with TSUNAMI, an implementation of Mikkola’s algorithmic regularisation (Mikkola & Tanikawa 1999a,b; Trani & Spera 2023). TSUNAMI is ideally suited for integrating the dynamical evolution of fewbody systems in which strong gravitational encounters are frequent and the mass ratio between the interacting bodies is large. It implements the logarithmic Hamiltonian and the timetransformed leapfrog, along with the nonregularised leapfrog as described in Mikkola (2020). In the present paper, we employ the logarithmic Hamiltonian scheme. TSUNAMI uses a BulirschStoer algorithm with stepsize control to improve the accuracy of the integration, which would otherwise be of second order (Stoer & Bulirsch 1980). Finally, TSUNAMI solves the equations of motion in a system of relative coordinates based on a chain of interparticle vectors in order to reduce roundoff errors when calculating small distances between bodies (Mikkola & Aarseth 1993). TSUNAMI includes velocitydependent forces, such as postNewtonian (PN) corrections and tidal forces (Mikkola & Merritt 2006; Mikkola & Tanikawa 2013). In this work, we enable PN terms of order 1PN, 2PN, 2.5PN, and 3.5PN (Blanchet 2014).
2.3. Stopping criteria
The initial time of the simulation is defined when the single comes within 2a_{inn} of the centre of mass of the binary. Since exchanges and ionisations can occur during the simulation, we identify the particles forming the bound binary at any given moment in time. Specifically, we check for the most bound stellarmass pair in the system and label it as the binary, while the remaining body is labelled as the single. We flag an interaction as terminated if the single is unbound and has positive radial velocity with respect to the binary. We then terminate the simulation once the binary–single distance exceeds 20 binary semimajor axes. We also stop a simulation when there is no bound binary, and all stellarmass bodies are unbound with respect to each other and their common centre of mass. Finally, we define a merger when two bodies become closer than the sum of their collision radii, which we set to ten times their Schwarzschild radius.
2.4. Classification of interactions
A threebody interaction can be thought of as a succession of two microstates: (1) ‘democratic interactions’ where all three bodies are close to each other and freely exchange angular momentum and energy, and (2) ‘excursions’, where the binary and the single recoil on a bound orbit, temporarily forming an unstable hierarchical triple (e.g. Monaghan 1976a,b; Valtonen & Karttunen 2005; Stone & Leigh 2019; Kol 2021; Ginat & Perets 2021a). Eventually, in Newtonian gravity, an unstable threebody system will inevitably decay into an unbound binary–single pair; that is, the single escapes the system.
Throughout the paper, we employ the following language to describe the qualitative nature of a threebody interaction. If the escape of the single happens after only one democratic interaction without any excursions, we label the interaction as ‘prompt’. On the other hand, if there is at least one excursion, we label the interaction as ‘resonant’ (using the terminology established in Hut & Bahcall 1983). A flyby, where a single passes close to the binary without breaking it and escapes immediately, is necessarily a prompt interaction. On the other hand, an exchange, where the single takes the place of a binary member, can result from both prompt and resonant interactions (see e.g. Manwadkar et al. 2020, 2021). Analogously, if the final binary is the same as the initial one, it might have resulted from a prompt interaction –a flyby– or from a resonant interaction –possibly after several exchanges with the single. For this reason, we label those binaries formed by the same bodies that made up the initial binary as ‘original’, regardless of whether they were formed via a resonant or a prompt interaction. Figure 2 displays the trajectories of the stellarmass BHs during a resonant exchange.
Fig. 2. Trajectories of the three stellarmass BHs during a threebody encounter around the SMBH. We highlight the different scale of the z axis with respect to the x–y coordinates, used in order to better visualise the trajectories. The interaction is resonant, because there are multiple episodes of democratic interactions and excursions. 
3. Results and discussion
3.1. Varying disc velocity dispersions
As explained in Sect. 2.1, the disc velocity dispersion σ controls the relative velocity (v_{rel}) between the binary and the single at the encounter location, which can be thought of as the velocity at infinity in a binary–single encounter in isolation, before the encounter. To better understand the impact of the disc velocity dispersion, we show in Fig. 3 the distribution of v_{rel} corresponding to σ = 10^{−1}, 10^{−2}, 10^{−3}, and 10^{−4}. For lower disc velocity dispersion, the difference in velocity only arises from the Keplerian shear, and is of the order of 10^{−2} km s^{−1}. For comparison, we also show two MaxwellBoltzmann distributions with a dispersion of σ_{ℳℬ}, which accurately describe the velocity dispersion in an isotropic cluster.
Fig. 3. Distributions of relative velocity between the centre of mass of the binary and that of the single at the encounter location R_{enc}, for different disc velocity dispersions σ. The blue, orange, and red histograms correspond to σ = 10^{−1}, 10^{−2}, 10^{−3}, and 10^{−4}, respectively. The grey dotted line is a MaxwellBoltzmann distribution, which is characteristic of isotropic clusters, and matches the velocity dispersion in σ2. The grey dashed line indicates the typical velocity distribution for a NSC that does not host a central SMBH. 
The velocity dispersion in a Keplerian disc is qualitatively different from that of an isotropic star cluster. A MaxwellBoltzmann distribution with σ_{ℳℬ} = 0.6 km s^{−1} is too peaked and sharp to efficiently approximate the v_{rel} distribution of the σ = 10^{−2} set. Most importantly, σ_{ℳℬ} = 0.6 km s^{−1} is typically associated with lowmass open clusters and not with a nuclear star cluster (NSC). A typical NSC, like that of the Milky Way, would have a velocity dispersion of σ_{ℳℬ} = 50 km s^{−1} if it did not have an SMBH at its centre (e.g. Atallah et al. 2023), which is comparable only with the most eccentric and misaligned disc we consider (σ = 10^{−1}, e, i ∼ 0.125, 3.6°).
The low velocity dispersion in a Keplerian disc has major implications for the outcome of threebody encounters. Binaries are more likely to be hard in a lowvelocitydispersion environment, which makes them more likely to harden, that is, to shrink in separation (Heggie 1975). The chance that a binary is hard is also increased by the dynamical stability constraint of Eq. (2), which enforces a strict upper limit on the size of binaries, beyond which they will be disrupted by the tidal field of the SMBH. Given our choice of initial conditions, we find that only 73.22% of the binaries in set σ1 have a binding energy greater than the kinetic energy of the incoming single. In every other set, all binaries are hard. This explains why set σ1 displays qualitatively different outcomes with respect to sets σ2 and σ3, as we show below.
For decreasing velocity dispersion, the centre of mass of the binary is close to that of the single for a longer period of time before the encounter, which makes them interact gravitationally with each other even before the encounter happens. This breaks our assumption that the two bodies keep on their original Keplerian before the encounter, and causes them to miss each other upon a number of encounters; that is, the single never comes within 2a_{inn} of the binary. We exclude these flybys from the following statistical analysis because the encounter does not actually happen. The frequency of such cases is negligible for σ1 and σ2, but it increases to ∼34% for σ3 and to ∼38% for σ ≤ 10^{−4}.
Table 2 summarises the simulation outcomes. The only set that presents breakups of the binary is set σ1, which is consistent with it being the one with 27% of the binaries being soft. Both the number of mergers during the simulations and the number of binaries with a GW coalescence time (t_{GW}) of less than a Hubble time (t_{Hubble} ≡ 14 Gyr) increase with decreasing velocity dispersion. Similarly, the number of prompt interactions is higher in σ1, as expected from its high velocity dispersion. All the mergers happen between the stellarmass BHs, and do not involve the SMBH. In principle, the escaping single might be kicked onto a highly eccentric orbit around the SMBH, triggering an extrememassratio inspiral (AmaroSeoane et al. 2007). We calculated the postencounter orbital properties of the single and estimated its coalescence time with the SMBH, but we find no inspirals that would occur within t_{Hubble}.
Interestingly, the fraction of resonant interactions changes in a nontrivial way as a function of the velocity dispersion. The highest percentage (∼59%) of resonant encounters occurs at σ = 10^{−2} and decreases for lower velocity dispersions to ∼54% at σ = 0 (we note that Table 2 shows this classification only for encounters with a postencounter binary).
The top two panels of Fig. 4 show the final binary properties for those simulations that do not end in a breakup or merger. The semimajor axes tend to be larger for σ1, whose binaries in the upper tail of the distribution are close to being disrupted by the tidal field of the SMBH. All the sets with σ ≤ 10^{−4} do not display significant differences in the semimajor axes distributions.
Fig. 4. Distributions of binary semimajor axes (top panel) and eccentricities (middle panel) of postencounter binaries, and encounter duration (bottom panel) for different disc velocity dispersions σ. The blue, orange, and green histograms correspond to σ = 10^{−1}, 10^{−2}, and 10^{−3}, respectively. The encounter duration is shown in units of the initial inner binary period, which is different for each realisation. The eccentricity distributions become more skewed towards e ∼ 1 as the disc velocity dispersion decreases, so that for σ ≥ 10^{−3} they are roughly thermal, while for σ ≤ 10^{−4} they are superthermal. The distribution of encounter duration shows two peaks, corresponding to prompt and resonant interactions, with the former peak being more pronounced in set σ1. 
The eccentricity distribution of the binaries after the encounter is also strongly dependent on the disc velocity dispersion. For σ ≤ 10^{−4}, the eccentricity distribution is remarkably superthermal, with an excess of eccentric binaries at e ∼ 1. The superthermal eccentricity of lowangular momentum triples has often been described by the distribution (e.g. Valtonen & Karttunen 2005), but in our case it is rather a mix of thermal eccentricities at e ≲ 0.9 and a superthermal component at e ≳ 0.9. For σ ≥ 10^{−3}, the eccentricity distribution is very close to a thermal distribution (f(e)∝e), with just a slight excess of highly eccentric binaries at e ∼ 1, which is more prominent in σ3.
This result can be explained in terms of the geometry of the system. As inclination and eccentricity decrease, the binary–single interaction retains less angular momentum and becomes increasingly confined to a 2D plane. This restriction to the phasespace dimensionality of the interacting triple leads to higher eccentricities in the final binary. Although this never reaches the extremes observed in systems with zero angular momentum or those fully confined to two dimensions (as in the equilateral triangle experiments of Parischewsky et al. 2023), the trend is evident, and one can observe the eccentricities shifting to a superthermal distribution. The eccentricity has a strong impact on the GW coalescence timescale (Peters 1964), which makes the number of binaries with t_{GW} < t_{Hubble} decrease with increasing velocity dispersion.
The bottom plot in Fig. 4 shows the distribution of encounter durations in units of the initial binary period. All the distributions are clearly bimodal, with the left and right peaks corresponding to prompt and resonant interactions, respectively.
Overall, these results agree with what we would expect from basic considerations. The disc velocity dispersion affects the eccentricity and inclination of orbits, which in turn determine the relative velocity between the binary and the single. An aligned and coplanar disc will have the lowest relative velocity during encounters, which will lead to highly eccentric binaries and more GW mergers.
However, upon closer examination of the merger fraction concerning velocity dispersion, a notable trend emerges. The merger fraction rises more rapidly as the velocity dispersion decreases than anticipated by analytic expectations. Ginat & Perets (2023) establish semianalytically that the merger probability – as a function of the total energy of the triple – is well approximated by a power law with an exponent of ≈3.5, as illustrated in their Fig. 1. Repeating a parallel analysis using our simulated data, we observe that the merger fraction is only broadly consistent with a powerlaw exponent of ≈4.0. Moreover, it exhibits a more pronounced increase at low energies and a more gradual decrease at higher energies compared to a powerlaw trend. The observed discrepancy is ascribed to the tidal field of the SMBH, which introduces qualitative changes in the behaviour of dynamical interactions. This aspect is further investigated in the following section.
3.2. The role of the tidal field of the SMBH
3.2.1. Expectations from analytic theories and previous results
Our simulations fully include the effect of the tidal field of the SMBH, which is modelled as a point particle. Based on the results of Trani et al. (2019b), we can expect the tidal field to have two main effects on the threebody dynamics. First, the encounter will be shorter in duration than encounters in isolation because the system is tidally limited. Consequently, a threebody interaction may end prematurely with the first excursion that extends beyond the tidal radius of the system. Because the interaction would then go through fewer democratic interactions, we can expect the final binaries to be less eccentric and softer than if the same encounter had happened in isolation. Second, the tidal field will exert a torque on the threebody system, which will make the system gain or lose angular momentum. Lowangularmomentum triples are the ones that result in a superthermal distribution of final binaries (e.g. Mikkola 1986; Valtonen et al. 2005; Manwadkar et al. 2024; Leigh et al. 2022; Parischewsky et al. 2023; Ginat & Perets 2023), and therefore have the highest chance of producing a GW coalescence. Trani et al. (2019b) found that lowangularmomentum triples in a Keplerian tidal field will tend to gain angular momentum, so that the final binary eccentricity distribution will be thermal rather than superthermal.
It is easy to estimate the truncating effect of the tidal field using simple calculations. As shown in Trani et al. (2019b), the tidal field of the binary limits how extended the outer orbit of the temporarily bound triple can be during an excursion. This limit is represented by the Hill radius, R_{Hill}, of the triple system. From energy conservation, this implies that the upper limit of the semimajor axis that the inner binary can have increases for decreasing R_{Hill}:
where E_{0} is the total energy, m = m_{1} + m_{2} is the inner binary mass, and m_{3} is the tertiary body (see also Eq. (6) of Trani et al. 2019b). This result was confirmed by numerical experiments and later extended by Ginat & Perets (2021b) who derived the resulting distribution of the binary semimajor axes from statistical escape theories using an arbitrary distance cutoff instead of R_{Hill} (see their Fig. 1).
Based on these results, we expect a shift in the distribution of semimajor axes towards larger values as the effect of the tidal field is stronger. In our simulations, all the sets have comparable R_{Hill} values, which range from ∼50 au to ∼200 au, with a median value of ∼134 au. We can then estimate the strength of the tidal field by comparing Eq. (3) –using the values of E_{0} and R_{Hill} from the simulations– with the same expression for R_{Hill} → +∞. Even though the distribution of a^{max} shifts to smaller values compared to that of lim_{RHill → +∞}a^{max}, both distributions have similar median values of a^{max} ≃ 0.6 au.
This suggests that the effect of the tidal field will be relatively small, given our set of initial conditions. However, it is worth noting that the estimate above only considers the truncating effect of the tidal field, neglecting the impact due to exchanges of energy and angular momentum with the triple system. Instead of investigating the impact of varying tidal field strengths on the outcomes of threebody interactions, which was previously addressed in Trani et al. (2019b), this section focuses on a different comparison. In the following subsection, we introduce a new set of simulations specifically designed to isolate the influence of the tidal field. We compare threebody interactions occurring in the vicinity of the SMBH with those taking place in an isolated environment.
3.2.2. Comparison with isolated encounters
To quantify the role of the tidal field, we performed simulations without the SMBH using the relative velocity (Fig. 3) derived from the simulations with the SMBH as the velocity at infinity in a hyperbolic binary–single encounter in isolation. For the other initial conditions, we used exactly the same properties (masses, inclination of the single, etc.) as those in the sets with the SMBH. We indicate these sets with the suffix i to specify that the threebody encounter occurred in isolation; for example σni, where n goes from −7 to 0. In this way, we still use the characteristic velocity dispersion associated to a Keplerian disc, but isolate the effect of the SMBH’s tidal field. In these simulations, the inclination between the binary and the hyperbolic orbit is the same as i_{sin} in the SMBH simulations. Our setup is therefore similar to the that considered in Samsing et al. (2022), although we introduce a velocity dispersion characteristic similar to that found in a disc, rather than assuming that the bodies are initially at rest.
Sets σn and σni are meant to isolate the effect of the SMBH’s tidal field in the same environmental conditions, that is, they both have a velocity dispersion associated with a Keplerian disc. However, we can also compare them to the environmental condition proper of a spherically symmetric, isotropic NSC without a central SMBH. Therefore, we performed an additional set of simulations, named σℳℬ, with the same initial conditions as sets σ1i, but drawing the velocity at infinity from a MaxwellBoltzmann distribution with σ_{ℳℬ} = 50 km s^{−1}. However, unlike in set σ1i, in σℳℬ, the inclination between the binary and the hyperbolic orbit is drawn uniformly in cos i, because in an NSC we expect no preferential direction.
Table 3 shows the frequency of final outcomes in the isolated sets, while Fig. 5 compares the outcome fractions between the two sets. The number of resonant interactions becomes smaller in the absence of tidal field, and the number of prompt interactions increases accordingly. This outcome is both puzzling and unexpected because one would anticipate the tidal field to prematurely terminate the encounters, reducing the likelihood of the single body returning to interact with the binary. Consequently, one would expect a decrease in the frequency of resonant interactions. Instead, depending on the velocity dispersion, the sets without the SMBH exhibit between 8% and 17% less resonant interactions.
Fig. 5. Fraction of outcomes for each set of simulations. The velocity dispersion decreases from left to right, and goes from σ = 10^{−1} to σ = 0. The solid lines with circle markers are the simulation in isolation (σni). The dashed lines with cross markers indicate the simulations with the SMBH (σn). The horizontal dotted lines indicate the results for σℳℬ. The top panel shows the fraction of mergers (blue lines) and the fraction of mergers with high eccentricity at f_{GW} = 10 Hz (orange lines). The bottom panel shows the fraction of prompt (green lines) and resonant (red lines) interactions, and the fraction of encounters ending with the original binary (purple lines) or with an exchanged member (brown lines). 
We find that this unexpected result is strictly tied to the definition of a resonant interaction, which we define as interactions with one or more excursions. If we use a stricter definition, allowing only interactions with two or more excursions, the results are reversed. In this case, the number of resonant interactions in isolation is from 8% to 60% larger compared to simulations with the SMBH, which agrees with our naive expectations.
The source of this discrepancy lies in the fact that interactions around the SMBH have a significantly higher likelihood of experiencing one excursion compared to none, while the opposite is true for encounters in isolation. Thorough inspection and validation of the excursiondetection algorithm in our simulations revealed no issues.
We propose the following hypothesis to explain this behaviour. When the triple system disintegrates, the binary–single components are ejected on unbound orbits whether they are in isolation or around the SMBH. However, in the presence of the SMBH, the binary and the single continue on Keplerian orbits, which remain in close proximity unless the escape velocity of the binary–single is comparable to the orbital Keplerian velocity. Given the close proximity of the two orbits, they can more easily gravitationally focus and interact again, potentially leading to subsequent interactions. This contrasts with encounters in isolation, where the binary–single components remain unbound as they escape to infinity, precluding further gravitational interactions.
In practice, this mechanism requires that the escape velocity between the binary–single at breakup be sufficiently high that the two escape orbits become adequately separated. Otherwise, gravitational focusing may lead to their interaction. This can explain why interactions around the SMBH show a greater probability of having a single excursion compared to none, while the opposite is true for those in isolation.
Another way to see this mechanism is within the framework of the circular restricted threebody problem, where a test particle moves under the influence of two massive bodies, where one is significantly more massive than the other. If the test particle is within the Hill region of the less massive body, having negative binding energy with respect to it is not enough to guarantee escape; the particle may still remain trapped in the Hill region, ‘bouncing’ against the zerovelocity surface before finding the Lagrange points, where the Hill region opens.
Even though this mechanism exists, we find that the simulations in isolation have a longer duration with respect to those around the SMBH, confirming our initial hypothesis (and the findings of Trani et al. 2019b) that the tidal field can interrupt excursions and end the encounters prematurely. The average encounter duration in isolation is between approximately two and ten times longer than the average encounter duration of encounters around the SMBH, depending on the velocity dispersion.
More interesting is the evolution of the merger fraction for different velocity dispersions. In isolation, the merger fraction saturates at ∼8% for σ ≤ 10^{−5}, while for the sets with the SMBH it rapidly increases in the same range. Consequently, for σ ≥ 10^{−4}, the sets in isolation have a consistently larger fraction of mergers, that is, up to ∼12 times larger at σ = 10^{−2}. This reverses for σ ≤ 10^{−5}, for which the sets with the SMBH have a higher fraction of mergers, up to ∼2 times more at σ = 0.
The dependency on the disc velocity dispersion is further exemplified in Fig. 6, which compares the distributions of semimajor axis and eccentricity of the final binaries. The postencounter binaries have smaller semimajor axes in the sets with the SMBH for σ ≤ 10^{−3}. The distribution of eccentricity also has a strong dependency on σ, which mimics the trend in the merger fraction. For σ ≥ 10^{−3}, the eccentricity distribution for encounters with the SMBH is thermal, while it is superthermal for the sets in isolation. At σ = 10^{−4}, the trend begins to reverse, with the distribution of σ4i becoming more thermal and that of σ4 becoming more superthermal. For σ ≤ 10^{−5}, the eccentricity distribution is superthermal for both sets, although it is consistently more peaked at e ∼ 1 in the sets with the SMBH.
Fig. 6. Cumulative distributions of binary semimajor axes (top panels) and eccentricities (bottom panels) of postencounter binaries. Solid thin lines indicate the sets with the SMBH (σn), while the thick dashed line indicates the sets in isolation (σni). The disc velocity dispersion for each set of curves decreases from left to right and from top to bottom. For σ ≥ 10^{−4}, the number of highly eccentric binaries is higher in the absence of tidal field. For σ ≤ 10^{−5}, the situation is the opposite, and encounters around the SMBH result in slightly higher eccentricities. 
This trend is consistent with both the results of Trani et al. (2019b), who find that the SMBH’s tidal field causes the superthermal eccentricity distribution of binaries from lowangular momentum triples to become thermal, and those of Samsing et al. (2022), who instead find a superthermal eccentricity even without the presence of the SMBH. Here we show that the disc velocity dispersion is a critical parameter in determining the outcome of threebody encounters in discs, as it determines the shape of the final eccentricity distribution.
The comparison between set σℳℬ and set σ1i is particularly intriguing. Despite both sets exhibiting similar velocity dispersions, as illustrated in Fig. 3, the binaries in set σℳℬ have an isotropic orientation, in contrast to the ordered orientation with ⟨i⟩∼⟨i exp( − i^{2}/2σ)⟩∼3.6° observed in set σ1i. Consequently, the triple systems in set σℳℬ possess higher angular momentum overall, leading to a final eccentricity distribution of binaries that is thermal. This in turn brings it closer to set σ1, which includes the SMBH. As a result, the merger fraction in set σℳℬ is significantly lower than in set σ1i. This highlights the importance of the binary orientation with respect to the incoming single.
In conclusion, including the presence of the SMBH is crucial to correctly estimate the outcome properties of threebody encounters around an SMBH. At high velocity dispersions (σ ≳ 10^{−4}), not including the SMBH results in overestimation of resonant encounters and mergers, while at low velocity dispersions it results in underestimation of the same quantities.
3.3. Implications for compact object mergers: Eccentricity and spin–orbit misalignment
Does the tidal field of the SMBH enhance or impair the formation of highly eccentric GW inspirals? As expected, the answer critically depends on the disc velocity dispersion. We checked for eccentric inspirals that enter the LIGOVirgoKagra band (f_{GW} = 10 Hz) for the mergers that take place during the simulations by calculating the GW peak frequency (f_{GW}) using the numerical fit of Hamers (2021). The results are shown in the tables and in the top panel of Fig. 5.
The number of eccentric mergers closely traces the number of mergers that occur during the simulation, with about 40% of the mergers retaining a significant eccentricity at f_{GW} = 10 Hz. For σ ≳ 10^{−4}, the fraction of eccentric mergers is up to ∼12 times greater in the isolated sets than in the sets with the SMBH, and at σ ≤ 10^{−5} the fraction can be up to ∼2 times larger in the sets with the SMBH. This is consistent with the results of the previous section, where we observe a transition in the eccentricity distributions at σ = 10^{−4}.
In light of this, the simulations of Samsing et al. (2022), which were performed with zero velocity dispersion and neglect the SMBH, have likely underestimated the number of eccentric mergers by a factor of approximately two. However, the most significant comparison is with the set σℳℬ, which represents the initial conditions in NSC without an SMBH. The fraction of mergers, eccentric or not, is significantly smaller, from ∼0.5 times smaller at σ = 10^{−1} to 1.4 × 10^{−2} smaller at σ = 0. Therefore, we can safely conclude that the outcome of threebody encounters in a BH disc around an SMBH strongly favours the production of highly eccentric mergers compared to encounters occurring in an NSC without an SMBH.
In addition to eccentricity, another property that can be used to discriminate GW formation channels is the spin–orbit misalignment. Depending on the formation pathways of binary BHs, their spin might be aligned with the total angular momentum of the orbit. This can be the case for binaries formed through gas capture in AGN discs, which are spunup and aligned by gas accretion (Tagawa et al. 2020). Another possible pathway to spin–orbit alignment is via the tidal spinup of postcommonenvelope binaries (e.g. Kushnir et al. 2016; Eldridge & Stanway 2016; Kruckow et al. 2018; Spera et al. 2019; Bavera et al. 2020; Belczynski et al. 2020; Piran & Piran 2020; Tanikawa et al. 2021).
While it is generally assumed that binaries in dynamically active environments have isotropic spin–orbit misalignments, Trani et al. (2021) showed that spin–orbit aligned BH binaries undergoing encounters retain memory of their initial configuration, resulting in nonisotropic tilt angle upon merger. Samsing et al. (2022) found similar results for BHs undergoing encounters in an AGN disc.
The bottom panel of Fig. 7 shows the distribution of the spin–orbit tilt angles Δi of the final binaries, assuming that the initial binary had spins aligned with the orbit. In the case of threebody encounters without the SMBH’s tidal field, the tilt angle distribution is bimodal and favours both aligned and antialigned spins, which is consistent with the results of Samsing et al. (2022). However, once we introduce the tidal field, we begin to see deviations at high disc velocity dispersion.
Fig. 7. Distribution of spin–orbit tilt angle (Δi) caused by the tilting of the binary plane during the encounter. An isotropic distribution would appear as flat in cos θ. Solid thin lines indicate the sets with the SMBH (σn), while thick dashed lines indicate the sets in isolation (σni). The disc velocity dispersion for each distribution decreases from left to right and from top to bottom. For σ ≥ 10^{−4}, the bimodality of the tilt angle distribution is strongly damped by the tidal field of the SMBH, leaving a majority of spin–orbit aligned binaries. 
Binary encounters in isolation are strongly bimodal, with two sharp peaks at Δi ∼ 0° and Δi ∼ 180°. This is especially true at low velocity dispersions, which is easily explained by noting that σ controls the inclination between binaries and singles. As σ tends to zero, the system is increasingly constrained into a plane, meaning that the only two possible inclinations are aligned or antialigned.
On the other hand, in the presence of the SMBH, the tidal field can amplify the small differences in inclination between the binary and the single. Consequently, for σ ≥ 10^{−4}, the sets with the SMBH have flatter distributions in cos(Δi), and the bimodality is entirely absent for σ ≥ 10^{−3}. As with the eccentricity distribution, the binaries from the isotropic set σℳℬ have a Δi distribution that is closer to that of set σ1 than to that of set σ1i.
Therefore, our results show that we should only expect a completely bimodal distribution of spin–orbit misalignment in the most aligned and circular discs. Even a small disc velocity dispersion will prevent the formation of antialigned Δi ∼ 180° binaries, which is because of the presence of the SMBH’s tidal field.
3.4. Caveats
Here, we consider nuclear discs around a MilkyWaylike SMBH with different velocity dispersion. The velocity dispersion we chose is consistent with that arising from twobody interactions in a disc, and can be expressed as a single, dimensionless parameter σ (Stewart & Ida 2000). However, the degree of velocity dispersion in a disc of BHs around an SMBH will depend on its formation mechanism. For example, a disc formed by the infall and disruption of a molecular cloud might be eccentric and lopsided following its formation (Mapelli et al. 2012; Trani et al. 2016, 2018). This configuration would likely lead to an excess of encounters at pericenter, where most of the orbits would cross, and at apocenter, where the bodies would spend most of their time. We plan to investigate this specific configuration in a future work.
For discs of black holes formed in an AGN disc, σ is comparable to the disctoheight ratio (h/R) of the accretion disc. The typical h/R for an accretion disc around a 10^{6} M_{⊙} SMBH (Sirko & Goodman 2003, see also D. Gangardt, in prep.) is 0.1–10^{−2}, which corresponds to our σ1 and σ2 sets.
On the other hand, if the disc has formed through capture in an AGN gaseous disc, many BHs could have aligned and circular orbits – depending on the thickness and lifetime of the disc – because of gas drag forces (Nasim et al. 2023; Generozov & Perets 2023). Our set σ0 best represents the environmental conditions in a perfectly aligned and circular disc. Our results point out that encounters in discs with σ ≤ 10^{−5} produce outcomes that are approximately equivalent to those of a circular disc, because their velocity dispersion is dominated by the Keplerian shear only. A transition to a randommotiondominated disc begins to emerge only at σ ∼ 10^{−4}, which corresponds to eccentricities and inclinations of order e, i ∼ 0.000125, 0.0036°.
These small values of e and i can only be achieved in an AGN disc if the alignment and subsequent circularisation processes are exceptionally efficient. However, such efficiency is deemed unlikely due to the limited lifespan of AGNs, which allows only a fraction of orbiting bodies to align within the disc. Moreover, those bodies that do become aligned often retain substantial eccentricities (Generozov & Perets 2023). Additionally, the turbulence within the disc, which is responsible for driving the viscosity and mass accretion, introduces stochastic forces on embedded bodies, hindering complete circularisation. Depending on the density of the BH population, mutual gravitational interactions between BHs may prevail over gas drag, leading to an increase in velocity dispersion. Based on these considerations, we argue that BH discs in AGNs should have large enough eccentricities and inclinations to make them randommotion dominated.
Additionally, if the disc was formed in situ, it will be composed of both stars and compact objects. Encounters between main sequence stars and BHs might lead to interesting phenomena, such as tidal disruption events, but we reserve the study of these for future work.
Throughout this work, we focused on comparing the outcome properties of threebody encounters, without considering the encounter rates between binaries and singles. While we find that threebody encounters in a disc result in more eccentric mergers than encounters in isolation, the frequency of these encounters remains to be determined. Migration traps in AGN discs have been invoked as a catalyst for such encounters (e.g. Bellovary et al. 2016), but the existence of such traps still needs to be fully ascertained (e.g. Grishin et al. 2023). In addition, the feasibility of aligning and embedding compact objects from around the SMBH within the lifetime of the AGN still needs to be explored using detailed AGN disc models (D. Gangardt, in prep.).
Another crucial assumption in this work is that binaries lie in the same plane as their orbit around the SMBH. Recent hydrodynamical studies suggest that misaligned binaries embedded in AGNs may gradually align with the disc over time (Dittmann et al. 2023), but the authors caution that this assumption may not hold for all types of discs around the SMBH.
Finally, the extent to which the dynamics of BHs in a dry disc apply to BHs embedded in an AGN disc remains unclear. We will explore this aspect in a followup work.
4. Conclusion
We investigated the outcomes of binary–single encounters in a disc of BHs around an SMBH by means of highly accurate simulations including up to 3.5 postNewtonian corrections. We examined the effect of the disc velocity dispersion (Sect. 3.1), the impact of the SMBH’s tidal field (Sect. 3.2), and how these factors affect the properties of merging BH binaries (Sect. 3.3). Our results can be summarised as follows:

The disc velocity dispersion σ, which parametrises the eccentricity and inclination of its constituent orbits, has a significant impact on the outcome of threebody encounters. For decreasing velocity dispersion, the eccentricity of postencounter binaries goes from subthermal to superthermal, and binaries become harder (Fig. 4). This transition from thermal to superthermal eccentricity happens sharply at around σ ≈ 10^{−4}. Encounters in discs with σ ≤ 10^{−5} yield outcomes that are almost the same as in perfectly circular and aligned discs, because for σ ≤ 10^{−5} the relative velocity of the binary–single pair is dominated by the Keplerian shear rather than by the velocity dispersion (Fig. 3). The main reason for this is that a higher velocity dispersion will enlarge the hard–soft boundary of the binaries embedded in the disc. Consequently, assessing the proper velocity dispersion of nuclear discs is key to correctly estimating the merger properties of their binary population.

The impact of the SMBH’s tidal field also depends on the disc velocity dispersion. By comparing with encounters in isolation, we find that neglecting the SMBH overestimates the number of mergers (by up to a factor of ∼12) for σ ≥ 10^{−4} and underestimates it (by up to a factor of ∼2) for σ < 10^{−4} (Fig. 5). This occurs because, if we neglect the SMBH, the final binaries are harder and more eccentric relative to encounters in the presence of an SMBH, but only at high velocity dispersion (Fig. 6).
Highly eccentric mergers are frequent during encounters in lowvelocitydispersion environments, when the binary and the single lie in the same plane, as proposed by Samsing et al. (2022). However, we find that this simple picture is complicated when considering the tidal field of the SMBH on encounters in discs with a velocity dispersion.
For high velocity dispersions σ > 10^{−4}, our results are in agreement with the outcomes of the disintegration of triples in a Keplerian potential (Trani et al. 2019b), with these latter authors finding that the SMBH’s tidal field hampers the production of a highly eccentric binary. A low disc velocity dispersion (i.e. dominated by the Keplerian shear only) reverses the role of the SMBH, and helps in the production of a superthermal distribution of final binaries.
The interplay between disc velocity dispersion and SMBH tidal field is also visible in the distributions of spin–orbit misalignment (Fig. 7). For dispersiondominated discs, (σ > 10^{−4}), the SMBH’s tidal field tends to isotropise the orientation of the final binaries, hampering the formation of binaries with antialigned spins.
Regardless of the disc velocity dispersion, we show that a nuclear disc can still produce a larger number of mergers than NCSs that do not host an SMBH, as long as the binaries lie in the plane of the disc. Even the disc with the highest velocity dispersion considered here (σ = 10^{−1}, e, i ∼ 0.125, 3.6°, equivalent to σ_{ℳℬ} ∼ 50 km s^{−1}) has approximately twice the number of mergers compared to an NSC without a SMBH, where the encounters show no preference in regard to the orientation of the binaries with respect to the incoming singles.
Our work points out the importance of characterising the velocity dispersions of nuclear BH discs – whether they are formed via AGN interaction or via other mechanisms – in order to constrain the properties of GW sources in galactic nuclei. In our next work, we will perform a demographic study of merging BH binaries in a nuclear disc, taking into account binary–single encounter rates and the interplay with ZKL oscillations, even after the encounters have taken place.
Acknowledgments
We would like to thank the anonymous referee for multiple useful comments and suggestions that helped to improve this work. A.A.T. wishes to thank Evgeni Grishin and Rosemary Mardling for helfpul comments and criticism. A.A.T. acknowledges support from JSPS KAKENHI Grant Number 21K13914 and from the European Union’s Horizon 2020 and Horizon Europe research and innovation programs under the Marie SkłodowskaCurie grant agreements no. 847523 and 101103134. All the simulations were performed on the awamori computer cluster at The University of Tokyo.
References
 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2019, Phys. Rev. X, 9, 031040 [Google Scholar]
 Abbott, R., Abbott, T. D., Abraham, S., et al. 2021a, Phys. Rev. X, 11, 021053 [Google Scholar]
 Abbott, R., Abbott, T. D., Abraham, S., et al. 2021b, ApJ, 913, L7 [NASA ADS] [CrossRef] [Google Scholar]
 AmaroSeoane, P., Gair, J. R., Freitag, M., et al. 2007, CQG, 24, R113 [Google Scholar]
 Antonini, F., & Perets, H. B. 2012, ApJ, 757, 27 [NASA ADS] [CrossRef] [Google Scholar]
 ArcaSedda, M., CapuzzoDolcetta, R., Antonini, F., & Seth, A. 2015, ApJ, 806, 220 [NASA ADS] [CrossRef] [Google Scholar]
 ArcaSedda, M., Kocsis, B., & Brandt, T. D. 2018, MNRAS, 479, 900 [NASA ADS] [CrossRef] [Google Scholar]
 Atallah, D., Trani, A. A., Kremer, K., et al. 2023, MNRAS, 523, 4227 [NASA ADS] [CrossRef] [Google Scholar]
 Bartko, H., Martins, F., Fritz, T. K., et al. 2009, ApJ, 697, 1741 [Google Scholar]
 Bartos, I., Kocsis, B., Haiman, Z., & Márka, S. 2017, ApJ, 835, 165 [Google Scholar]
 Bavera, S. S., Fragos, T., Qin, Y., et al. 2020, A&A, 635, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belczynski, K., Klencki, J., Fields, C. E., et al. 2020, A&A, 636, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bellovary, J. M., Mac Low, M.M., McKernan, B., & Ford, K. E. S. 2016, ApJ, 819, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Bender, R., Kormendy, J., Bower, G., et al. 2005, ApJ, 631, 280 [NASA ADS] [CrossRef] [Google Scholar]
 Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 [NASA ADS] [CrossRef] [Google Scholar]
 Blanchet, L. 2014, Liv. Rev. Relativ., 17, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Boekholt, T. C. N., Rowan, C., & Kocsis, B. 2023, MNRAS, 518, 5653 [Google Scholar]
 Brown, C. K., & Magorrian, J. 2013, MNRAS, 431, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Camilloni, F., Grignani, G., Harmark, T., Orselli, M., & Pica, D., 2023, arXiv eprints [arXiv:2310.06894] [Google Scholar]
 Cantiello, M., Jermyn, A. S., & Lin, D. N. C. 2021, ApJ, 910, 94 [NASA ADS] [CrossRef] [Google Scholar]
 DeLaurentiis, S., EpsteinMartin, M., & Haiman, Z. 2023, MNRAS, 523, 1126 [NASA ADS] [CrossRef] [Google Scholar]
 Dittmann, A. J., Dempsey, A. M., & Li, H. 2023, arXiv eprints [arXiv:2310.03832] [Google Scholar]
 Eldridge, J. J., & Stanway, E. R. 2016, MNRAS, 462, 3302 [Google Scholar]
 Fabj, G., Nasim, S. S., Caban, F., et al. 2020, MNRAS, 499, 2608 [NASA ADS] [CrossRef] [Google Scholar]
 Generozov, A., & Perets, H. B. 2023, MNRAS, 522, 1763 [NASA ADS] [CrossRef] [Google Scholar]
 Genzel, R., Schödel, R., Ott, T., et al. 2003, ApJ, 594, 812 [Google Scholar]
 Gillessen, S., Eisenhauer, F., Trippe, S., et al. 2009, ApJ, 692, 1075 [Google Scholar]
 Gillessen, S., Plewa, P. M., Eisenhauer, F., et al. 2017, ApJ, 837, 30 [Google Scholar]
 Ginat, Y. B., & Perets, H. B. 2021a, Phys. Rev. X, 11, 031020 [NASA ADS] [Google Scholar]
 Ginat, Y. B., & Perets, H. B. 2021b, MNRAS, 508, 190 [NASA ADS] [CrossRef] [Google Scholar]
 Ginat, Y. B., & Perets, H. B. 2023, MNRAS, 519, L15 [Google Scholar]
 Grishin, E., Gilbaum, S., & Stone, N. C. 2023, MNRAS, submitted [arXiv:2307.07546] [Google Scholar]
 Hamers, A. S. 2021, Res. Notes Am. Astron. Soc., 5, 275 [Google Scholar]
 Heggie, D. C. 1975, MNRAS, 173, 729 [NASA ADS] [CrossRef] [Google Scholar]
 Hoang, B.M., Naoz, S., Kocsis, B., Rasio, F. A., & Dosopoulou, F. 2018, ApJ, 856, 140 [NASA ADS] [CrossRef] [Google Scholar]
 Hopkins, P. F., & Quataert, E. 2010a, MNRAS, 407, 1529 [Google Scholar]
 Hopkins, P. F., & Quataert, E. 2010b, MNRAS, 405, L41 [NASA ADS] [Google Scholar]
 Hut, P., & Bahcall, J. N. 1983, ApJ, 268, 319 [NASA ADS] [CrossRef] [Google Scholar]
 Ida, S., & Makino, J. 1992, Icarus, 96, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Jermyn, A. S., Dittmann, A. J., McKernan, B., Ford, K. E. S., & Cantiello, M. 2022, ApJ, 929, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Kamlah, A. W. H., Spurzem, R., Berczik, P., et al. 2022, MNRAS, 516, 3266 [CrossRef] [Google Scholar]
 King, I. R. 1966, AJ, 71, 276 [NASA ADS] [CrossRef] [Google Scholar]
 Kokubo, E., & Ida, S. 2002, ApJ, 581, 666 [Google Scholar]
 Kol, B. 2021, Celest. Mech. Dyn. Astron., 133, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Kozai, Y. 1962, AJ, 67, 591 [Google Scholar]
 Kruckow, M. U., Tauris, T. M., Langer, N., Kramer, M., & Izzard, R. G. 2018, MNRAS, 481, 1908 [CrossRef] [Google Scholar]
 Kushnir, D., Zaldarriaga, M., Kollmeier, J. A., & Waldman, R. 2016, MNRAS, 462, 844 [NASA ADS] [CrossRef] [Google Scholar]
 Lauer, T. R., Faber, S. M., Groth, E. J., et al. 1993, AJ, 106, 1436 [Google Scholar]
 Lauer, T. R., Faber, S. M., Gebhardt, K., et al. 2005, AJ, 129, 2138 [NASA ADS] [CrossRef] [Google Scholar]
 Leigh, N. W. C., Geller, A. M., McKernan, B., et al. 2018, MNRAS, 474, 5672 [Google Scholar]
 Leigh, N. W. C., Stone, N. C., Webb, J. J., & Lyra, W. 2022, MNRAS, 517, 3838 [CrossRef] [Google Scholar]
 Levin, Y., & Beloborodov, A. M. 2003, ApJ, 590, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Li, J., Lai, D., & Rodet, L. 2022, ApJ, 934, 154 [NASA ADS] [CrossRef] [Google Scholar]
 Li, J., Dempsey, A. M., Li, H., Lai, D., & Li, S. 2023, ApJ, 944, L42 [NASA ADS] [CrossRef] [Google Scholar]
 Lidov, M. L. 1962, Planet. Space Sci., 9, 719 [Google Scholar]
 Lissauer, J. J., & Stewart, G. R. 1993, in Protostars and Planets III, eds. E. H. Levy, & J. I. Lunine, 1061 [Google Scholar]
 Liu, B., Muñoz, D. J., & Lai, D. 2015, MNRAS, 447, 747 [NASA ADS] [CrossRef] [Google Scholar]
 Lu, J. R., Ghez, A. M., Hornstein, S. D., et al. 2009, ApJ, 690, 1463 [NASA ADS] [CrossRef] [Google Scholar]
 MacLeod, M., & Lin, D. N. C. 2020, ApJ, 889, 94 [NASA ADS] [CrossRef] [Google Scholar]
 Manwadkar, V., Trani, A. A., & Leigh, N. W. C. 2020, MNRAS, 497, 3694 [CrossRef] [Google Scholar]
 Manwadkar, V., Kol, B., Trani, A. A., & Leigh, N. W. C. 2021, MNRAS, 506, 692 [NASA ADS] [Google Scholar]
 Manwadkar, V., Trani, A. A., & Kol, B. 2024, Celest. Mech. Dyn. Astron., 136, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Mapelli, M. 2021, in Handbook of Gravitational Wave Astronomy, eds. C. Bambi, S. Katsanevas, & K. D. Kokkotas (Living Reference Work. Springer), 16 [Google Scholar]
 Mapelli, M., Hayfield, T., Mayer, L., & Wadsley, J. 2012, ApJ, 749, 168 [NASA ADS] [CrossRef] [Google Scholar]
 Martinez, M. A. S., Rodriguez, C. L., & Fragione, G. 2022, ApJ, 937, 78 [NASA ADS] [CrossRef] [Google Scholar]
 MastrobuonoBattisti, A., Perets, H. B., Gualandris, A., Neumayer, N., & Sippel, A. C. 2019, MNRAS, 490, 5820 [Google Scholar]
 Máthé, G., Szölgyén, Á., & Kocsis, B. 2023, MNRAS, 520, 2204 [CrossRef] [Google Scholar]
 Michie, R. W. 1963, MNRAS, 125, 127 [Google Scholar]
 Mikkola, S. 1986, MNRAS, 223, 757 [NASA ADS] [CrossRef] [Google Scholar]
 Mikkola, S. 2020, Gravitational FewBody Dynamics: A Numerical Approach (Cambridge University Press) [CrossRef] [Google Scholar]
 Mikkola, S., & Aarseth, S. J. 1993, Celest. Mech. Dyn. Astron., 57, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Mikkola, S., & Merritt, D. 2006, MNRAS, 372, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Mikkola, S., & Tanikawa, K. 1999a, MNRAS, 310, 745 [NASA ADS] [CrossRef] [Google Scholar]
 Mikkola, S., & Tanikawa, K. 1999b, Celest. Mech. Dyn. Astron., 74, 287 [NASA ADS] [CrossRef] [Google Scholar]
 Mikkola, S., & Tanikawa, K. 2013, MNRAS, 430, 2822 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, M. C., & Lauburg, V. M. 2009, ApJ, 692, 917 [Google Scholar]
 Monaghan, J. J. 1976a, MNRAS, 176, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Monaghan, J. J. 1976b, MNRAS, 177, 583 [NASA ADS] [CrossRef] [Google Scholar]
 Naoz, S. 2016, ARA&A, 54, 441 [Google Scholar]
 Nasim, S. S., Fabj, G., Caban, F., et al. 2023, MNRAS, 522, 5393 [NASA ADS] [CrossRef] [Google Scholar]
 Parischewsky, H. D., Trani, A. A., & Leigh, N. W. C. 2023, SciPost Phys. Core, 6, 016 [NASA ADS] [CrossRef] [Google Scholar]
 Paumard, T., Genzel, R., Martins, F., et al. 2006, ApJ, 643, 1011 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H. V., & Tremaine, S. 2003, ApJ, 599, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Peters, P. C. 1964, Phys. Rev., 136, B1224 [NASA ADS] [CrossRef] [Google Scholar]
 Piran, Z., & Piran, T. 2020, ApJ, 892, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Qian, K., Li, J., & Lai, D. 2024, ApJ, 962, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Rom, B., Sari, R., & Lai, D. 2023, ApJ, submitted [arXiv:2310.03801] [Google Scholar]
 Rowan, C., Boekholt, T., Kocsis, B., & Haiman, Z. 2023, MNRAS, 524, 2770 [NASA ADS] [CrossRef] [Google Scholar]
 Samsing, J., Bartos, I., D’Orazio, D. J., et al. 2022, Nature, 603, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Secunda, A., Bellovary, J., Mac Low, M.M., et al. 2019, ApJ, 878, 85 [NASA ADS] [CrossRef] [Google Scholar]
 Secunda, A., Hernandez, B., Goodman, J., et al. 2021, ApJ, 908, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Shevchenko, I. I. 2017, The LidovKozai Effect  Applications in Exoplanet Research and Dynamical Astronomy (Springer International Publishing), 441 [CrossRef] [Google Scholar]
 Sirko, E., & Goodman, J. 2003, MNRAS, 341, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Spera, M., Mapelli, M., Giacobbo, N., et al. 2019, MNRAS, 485, 889 [Google Scholar]
 Spera, M., Trani, A. A., & Mencagli, M. 2022, Galaxies, 10, 76 [NASA ADS] [CrossRef] [Google Scholar]
 Spitzer, L. 1987, Dynamical Evolution of Globular Clusters (Princeton, N.J. : Princeton University Press) [Google Scholar]
 Stewart, G. R., & Ida, S. 2000, Icarus, 143, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Stoer, J., & Bulirsch, R. 1980, Introduction to Numerical Analysis (Springer Verlag: New york), 430 [Google Scholar]
 Stone, N. C., Metzger, B. D., & Haiman, Z. 2017, MNRAS, 464, 946 [Google Scholar]
 Stone, N. C., & Leigh, N. W. C. 2019, Nature, 576, 406 [NASA ADS] [CrossRef] [Google Scholar]
 Szölgyén, Á., Meiron, Y., & Kocsis, B. 2019, ApJ, 887, 123 [CrossRef] [Google Scholar]
 Tagawa, H., Haiman, Z., & Kocsis, B. 2020, ApJ, 898, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Tanikawa, A., Susa, H., Yoshida, T., Trani, A. A., & Kinugawa, T. 2021, ApJ, 910, 30 [NASA ADS] [CrossRef] [Google Scholar]
 The LIGO Scientific Collaboration & The Virgo Collaboration (Abbott, R., et al.) 2021, arXiv eprints [arXiv:2108.01045] [Google Scholar]
 The LIGO Scientific Collaboration, The Virgo Collaboration,& The KAGRA Collaboration 2023a, Phys. Rev. X, 13, 041039 [Google Scholar]
 The LIGO Scientific Collaboration, The Virgo Collaboration,& The KAGRA Collaboration 2023b, Phys. Rev. X, 13, 011048 [NASA ADS] [Google Scholar]
 Trani, A. A., & Spera, M. 2023, IAU Symp., 362, 404 [NASA ADS] [Google Scholar]
 Trani, A. A., Mapelli, M., Bressan, A., et al. 2016, ApJ, 818, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Trani, A. A., Mapelli, M., & Ballone, A. 2018, ApJ, 864, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Trani, A. A., Fujii, M. S., & Spera, M. 2019a, ApJ, 875, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Trani, A. A., Spera, M., Leigh, N. W. C., & Fujii, M. S. 2019b, ApJ, 885, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Trani, A. A., Tanikawa, A., Fujii, M. S., Leigh, N. W. C., & Kumamoto, J. 2021, MNRAS, 504, 910 [NASA ADS] [CrossRef] [Google Scholar]
 Trani, A. A., Rastello, S., Carlo, U. N. D., et al. 2022, MNRAS, 511, 1362 [NASA ADS] [Google Scholar]
 Tremaine, S. 1995, AJ, 110, 628 [NASA ADS] [CrossRef] [Google Scholar]
 Valtonen, M., & Karttunen, H. 2005, The ThreeBody Problem (Cambridge University Press) [Google Scholar]
 Valtonen, M., Mylläri, A., Orlov, V., & Rubinov, A. 2005, MNRAS, 364, 91 [NASA ADS] [CrossRef] [Google Scholar]
 von Fellenberg, S. D., Gillessen, S., Stadler, J., et al. 2022, ApJ, 932, L6 [NASA ADS] [CrossRef] [Google Scholar]
 von Zeipel, H. 1910, Astron. Nachr., 183, 345 [Google Scholar]
 Yang, Y., Bartos, I., Haiman, Z., et al. 2019, ApJ, 876, 122 [NASA ADS] [CrossRef] [Google Scholar]
 Yelda, S., Ghez, A. M., Lu, J. R., et al. 2014, ApJ, 783, 131 [Google Scholar]
All Tables
All Figures
Fig. 1. Scheme of the initial setup of our simulations. The centre of mass of the binary and the single are set to encounter at a location R_{enc} with impact parameter b ∈ (0, 2a_{inn}). The eccentricities and mutual inclinations of the two orbits around the SMBH are drawn from a Rayleigh distribution (ℛ) with scale parameter σ and σ/2, respectively. 

In the text 
Fig. 2. Trajectories of the three stellarmass BHs during a threebody encounter around the SMBH. We highlight the different scale of the z axis with respect to the x–y coordinates, used in order to better visualise the trajectories. The interaction is resonant, because there are multiple episodes of democratic interactions and excursions. 

In the text 
Fig. 3. Distributions of relative velocity between the centre of mass of the binary and that of the single at the encounter location R_{enc}, for different disc velocity dispersions σ. The blue, orange, and red histograms correspond to σ = 10^{−1}, 10^{−2}, 10^{−3}, and 10^{−4}, respectively. The grey dotted line is a MaxwellBoltzmann distribution, which is characteristic of isotropic clusters, and matches the velocity dispersion in σ2. The grey dashed line indicates the typical velocity distribution for a NSC that does not host a central SMBH. 

In the text 
Fig. 4. Distributions of binary semimajor axes (top panel) and eccentricities (middle panel) of postencounter binaries, and encounter duration (bottom panel) for different disc velocity dispersions σ. The blue, orange, and green histograms correspond to σ = 10^{−1}, 10^{−2}, and 10^{−3}, respectively. The encounter duration is shown in units of the initial inner binary period, which is different for each realisation. The eccentricity distributions become more skewed towards e ∼ 1 as the disc velocity dispersion decreases, so that for σ ≥ 10^{−3} they are roughly thermal, while for σ ≤ 10^{−4} they are superthermal. The distribution of encounter duration shows two peaks, corresponding to prompt and resonant interactions, with the former peak being more pronounced in set σ1. 

In the text 
Fig. 5. Fraction of outcomes for each set of simulations. The velocity dispersion decreases from left to right, and goes from σ = 10^{−1} to σ = 0. The solid lines with circle markers are the simulation in isolation (σni). The dashed lines with cross markers indicate the simulations with the SMBH (σn). The horizontal dotted lines indicate the results for σℳℬ. The top panel shows the fraction of mergers (blue lines) and the fraction of mergers with high eccentricity at f_{GW} = 10 Hz (orange lines). The bottom panel shows the fraction of prompt (green lines) and resonant (red lines) interactions, and the fraction of encounters ending with the original binary (purple lines) or with an exchanged member (brown lines). 

In the text 
Fig. 6. Cumulative distributions of binary semimajor axes (top panels) and eccentricities (bottom panels) of postencounter binaries. Solid thin lines indicate the sets with the SMBH (σn), while the thick dashed line indicates the sets in isolation (σni). The disc velocity dispersion for each set of curves decreases from left to right and from top to bottom. For σ ≥ 10^{−4}, the number of highly eccentric binaries is higher in the absence of tidal field. For σ ≤ 10^{−5}, the situation is the opposite, and encounters around the SMBH result in slightly higher eccentricities. 

In the text 
Fig. 7. Distribution of spin–orbit tilt angle (Δi) caused by the tilting of the binary plane during the encounter. An isotropic distribution would appear as flat in cos θ. Solid thin lines indicate the sets with the SMBH (σn), while thick dashed lines indicate the sets in isolation (σni). The disc velocity dispersion for each distribution decreases from left to right and from top to bottom. For σ ≥ 10^{−4}, the bimodality of the tilt angle distribution is strongly damped by the tidal field of the SMBH, leaving a majority of spin–orbit aligned binaries. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.