Search

The phase evolution of gravitational waves (GWs) can be modulated by the astrophysical environment surrounding the source, which provides a probe for the origin of individual binary black holes (BBHs) using GWs alone. We here study the evolving phase of the GW waveform derived from a large set of simulations of BBH mergers forming in dense stellar clusters through binary-single interactions. We uncover that a well-defined fraction of the assembled eccentric GW sources will have a notable GW phase shift induced by the remaining third object. The magnitude of the GW phase shift often exceeds conservative analytical estimates due to strong 3-body interactions, which occasionally results in GW sources with clearly shifted and perturbed GW waveforms. This opens up promising opportunities for current and future GW detectors, as observing such a phase shift can identify the formation environment of a BBH, as well as help to characterise the local properties of its surrounding environment.

Although planets have been found orbiting binary systems, whether they can survive binary interactions is debated. While the tightest-orbit binaries should host the most dynamically stable and long-lived circumbinary planetary systems, they are also the systems that are expected to experience mass transfer, common envelope evolution, or stellar mergers. In this study, we explore the effect of stable non-conservative mass transfer on the dynamical evolution of circumbinary planets. We present a new script that seamlessly integrates binary evolution data from the 1D binary stellar evolution code MESA into the N-body simulation code REBOUND. This integration framework enables a comprehensive examination of the dynamical evolution of circumbinary planets orbiting mass-transferring binaries, while simultaneously accounting for the detailed stellar structure evolution. In addition, we introduce a recalibration method to mitigate numerical errors from updates of binary properties during the system's dynamical evolution. We construct a reference binary model in which a $2.21\, M_\odot$ star loses its hydrogen-rich envelope through non-conservative mass transfer to the $1.76\, M_\odot$ companion star, creating a $0.38\,M_\odot$ subdwarf. We find the tightest stable orbital separation for circumbinary planets to be $\simeq 2.5$ times the binary separation after mass transfer. Accounting for tides by using the interior stellar structure, we find that tidal effects become apparent after the rapid mass transfer phase and start to fade away during the latter stage of the slow mass transfer phase. Our research provides a new framework for exploring circumbinary planet dynamics in interacting binary systems., Comment: Submitted to MNRAS, comments are welcome

Recent numerical results seem to suggest that in certain regimes of typical particle velocities the gravitational $N-$body problem (for $3\leq N\lesssim 10^3$) is intrinsically less chaotic when the post-Newtonian (PN) force terms are included, with respect to its classical counterpart that exhibits a slightly larger maximal Lyapunov exponent $\Lambda_{\rm max}$. In this work we explore the dynamics of wildly chaotic, regular and nearly regular configurations of the 3-body problem with and without the PN corrective terms aiming at shedding some light on the behaviour of the Lyapunov spectra under the effect of said corrections. Because the interaction of the tangent-space dynamics in gravitating systems, needed to evaluate the Lyapunov exponents, becomes rapidly computationally heavy due to the complexity of the higher order force derivatives involving multiple powers of $v/c$, we introduce a technique to compute a proxy of the Lyapunov spectrum based on the time-dependent diagonalization of the inertia tensor of a cluster of trajectories in phase-space. We find that, for a broad range of orbital configurations, the relativistic 3-body problem has a smaller $\Lambda_{\rm max}$ than its classical counterpart starting with the exact same initial condition. However, the rest of the Lyapunov spectrum can be either lower or larger in the classical case, suggesting that the relativistic precession effectively reduces chaos only along one (or few) directions in phase-space. As a general trend, the dynamical entropy of the relativistic simulations as function of the rescaled speed of light always has a regime in which falls below the classical value.} We observe that, the sole analysis of $\Lambda_{\rm max}$ could induce possibly misleading conclusions on the chaoticity of systems with small (and possibly large $N$., Comment: 9 pages, 8 figures. Submitted, comments welcome

This study explores the formation and implications of mini-active galactic nuclei (mAGN) disks around intermediate-mass black holes (IMBHs) embedded in gas-rich globular/nuclear clusters (GCs). We examine the parameter space for stable mAGN disks, considering the influence of IMBH mass, disk radius, and gas density on disk stability. The dynamics of stars and black holes within the mAGN disk are modeled, with a focus on gas-induced migration and gas dynamical friction. These dynamical processes can lead to several potentially observable phenomena, including the alignment of stellar orbits into the disk plane, the enhancement of gravitational wave mergers (particularly IMRIs and EMRIs), and the occurrence of mili/centi-tidal disruption events (mTDEs/cTDEs) with unique observational signatures. We find that gas hardening can significantly accelerate the inspiral of binaries within the disk, potentially leading to a frequency shift in the emitted gravitational waves. Additionally, we explore the possibility of forming accreting IMBH systems from captured binaries within the mAGN disk, potentially resulting in the formation of ultraluminous X-ray sources (ULXs). The observational implications of such accreting systems, including X-ray emission, optical signatures, and transient phenomena, are discussed. Furthermore, we investigate the possibility of large-scale jets emanating from gas-embedded IMBHs in GCs. While several caveats and uncertainties exist, our work highlights the potential for mAGN disks to provide unique insights into IMBH demographics, accretion physics, and the dynamics of GCs., Comment: Comments welcome!

We explore the stability and fate of gravitational triple systems comprising a central massive body and a tight binary of less massive pairs. Our present purpose is two fold; (1) to improve the Hill-type stability criterion for the binary in those systems, and (2) to examine the fate of the triple systems after the binary break-up, with particular attention to the effects of the eccentricities of the inner and outer orbits. We perform direct Newtonian N-body simulations over much longer integration times than previous studies, which is essential to determine the eventual fate of those systems statistically in a reliable fashion. We obtain an empirical fitting formula of the binary stability boundary that incorporates effects of the inner and outer eccentricities, the mutual inclination of the inner and outer orbits, and the mass ratios of the three bodies. We also find that those triple systems are stable for a much longer timescale after the binary break-up, and that their final fates (ejection of the outer body, merger to the central massive body, and collision of two less massive bodies) are very sensitive to the initial outer eccentricity., Comment: 22 pages, 14 figures, submitted to ApJ

We study the gravitational wave (GW) emission of sources perturbed by periodic dynamical forces which do not cause secular evolution in the orbital elements. We construct a corresponding post-Newtonian waveform model and provide estimates for the detectability of the resulting GW phase perturbations, for both space-based and future ground-based detectors. We validate our results by performing a set of Bayesian parameter recovery experiments with post-Newtonian waveforms. We find that, in stark contrast to the more commonly studied secular dephasing, periodic phase perturbations do not suffer from degeneracies with any of the tested vacuum binary parameters. We discuss the applications of our findings to a range of possible astrophysical scenarios, finding that such periodic perturbations may be detectable for massive black hole binaries embedded in circum-binary discs, extreme mass-ratio inspirals in accretion discs, as well as stellar-mass compact objects perturbed by tidal fields. We argue that modelling conservative sub-orbital dynamics opens up a promising new avenue to detect environmental effects in binary sources of GWs that should be included in state-of-the-art waveform templates., Comment: Accepted in PRD. Comments welcome

Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use the MIST stellar models and our open-source population synthesis code SEVN to evolve $5.6 \times 10^8$ binaries exploring 20 sets of parameters. We find that systems like Gaia BH3 form preferentially from binaries initially composed of a massive star ($40-60$ M$_\odot$) and a low mass companion ($<1$ M$_\odot$) in a wide ($P>10^3$ days) and eccentric orbit ($e>0.6$). Such progenitor binary stars do not undergo any Roche-lobe overflow episode during their entire evolution, so that the final orbital properties of the BH-star system are determined at the core collapse of the primary star. Low natal kicks ($\lesssim$ 10~km/s) significantly favour the formation of Gaia BH3-like systems, but high velocity kicks up to $\approx 220$ km/s are also allowed. We estimate the formation efficiency for Gaia BH3-like systems in old ($t>10$ Gyr) and metal-poor ($Z<0.01$) populations to be $\sim 4 \times 10^{-8}$ M$_\odot^{-1}$ (for our fiducial model), representing $\sim 3\%$ of the whole simulated BH-star population. We expect up to $\approx 4000$ BH-star systems in the Galactic halo formed through isolated evolution, of which $\approx 100$ are compatible with Gaia BH3-like. Given the density profile of the Galactic halo we do not expect more than one at the observed distance of Gaia BH3. Our models show that, even if it was born inside a stellar cluster, Gaia BH3 is compatible with a primordial binary star that escaped from its parent cluster without experiencing significant dynamical interactions., Comment: Accepted for publication in A&A. The produced datasets are publicly available in Zenodo. I Initiail conditions and run scripts: https://zenodo.org/records/11617742; IIA simulation outcomes: https://zenodo.org/records/11992265; IIB simulation outcomes, different metallicities: https://zenodo.org/records/12188369; IIC simulation outcomes, fiducial 1E8, https://zenodo.org/records/12189734

The three-body problem (3BP) poses a longstanding challenge in physics and celestial mechanics. Despite the impossibility of obtaining general analytical solutions, statistical theories have been developed based on the ergodic principle. This assumption is justified by chaos, which is expected to fully mix the accessible phase space of the 3BP. This study probes the presence of regular (i.e. non chaotic) trajectories within the 3BP and assesses their impact on statistical escape theories. Using numerical simulations, we establish criteria for identifying regular trajectories and analyse their impact on statistical outcomes. Our analysis reveals that regular trajectories occupy up to 32% of the phase space, and their outcomes defy the predictions of statistical escape theories. The coexistence of regular and chaotic regions at all scales is characterized by a multi-fractal behaviour. Integration errors manifest as numerical chaos, artificially enhancing the mixing of the phase space and affecting the reliability of individual simulations, yet preserving the statistical correctness of an ensemble of realizations. Our findings underscore the challenges in applying statistical escape theories to astrophysical problems, as they may bias results by excluding the outcome of regular trajectories. This is particularly important in the context of formation scenarios of gravitational wave mergers, where biased estimates of binary eccentricity can significantly impact estimates of coalescence efficiency and detectable eccentricity., Comment: 15 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Models of accretion discs surrounding active galactic nuclei (AGNs) find vast applications in high-energy astrophysics. The broad strategy is to parametrize some of the key disc properties such as gas density and temperature as a function of the radial coordinate from a given set of assumptions on the underlying physics. Two of the most popular approaches in this context were presented by Sirko & Goodman (2003) and Thompson et al. (2005). We present a critical reanalysis of these widely used models, detailing their assumptions and clarifying some steps in their derivation that were previously left unsaid. Our findings are implemented in the pAGN module for the Python programming language, which is the first public implementation of these accretion-disc models. We further apply pAGN to the evolution of stellar-mass black holes embedded in AGN discs, addressing the potential occurrence of migration traps., Comment: 18 pages, 8 figures, 2 tables, code available at https://github.com/DariaGangardt/pAGN

We explore three-body binary formation (3BBF), the formation of a bound system via gravitational scattering of three initially unbound bodies (3UB), using direct numerical integrations. For the first time, we consider systems with unequal masses, as well as finite-size and post-Newtonian effects. Our analytically derived encounter rates and numerical scattering results reproduce the 3BBF rate predicted by Goodman \& Hut (1993) for hard binaries in dense star clusters. We find that 3BBF occurs overwhelmingly through nonresonant encounters and that the two most massive bodies are never the most likely to bind. Instead, 3BBF favors pairing the two least massive bodies (for wide binaries) or the most plus least massive bodies (for hard binaries). 3BBF overwhelmingly favors wide binary formation with super-thermal eccentricities, perhaps helping to explain the eccentric wide binaries observed by Gaia. Hard binary formation is far rarer, but with a thermal eccentricity distribution. The semimajor axis distribution scales cumulatively as $a^3$ for hard and slightly wider binaries. Though mergers are rare between black holes when including relativistic effects, direct collisions occur frequently between main-sequence stars -- more often than hard 3BBF. Yet, these collisions do not significantly suppress hard 3BBF at the low velocity dispersions typical of open or globular clusters. Energy dissipation through gravitational radiation leads to a small probability of a bound, hierarchical triple system forming directly from 3UB., Comment: 28 pages, 16 figures, submitted to ApJ. Comments welcome

Dynamical encounters of stellar-mass 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 3-body dynamics. 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. We perform simulations of 4-body encounters (i.e. with the SMBH as fourth particle) using the highly-accurate, regularized code TSUNAMI, which includes post-Newtonian corrections up to order 3.5PN. 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 post-encounter 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, and one dominated by the Keplerian shear. 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 anti-aligned binary BHs., Comment: 12 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Theoretical predictions suggest that very massive stars have the potential to form through multiple collisions and eventually evolve into intermediate-mass black holes (IMBHs) within Population III star clusters embedded in mini dark matter haloes. In this study, we investigate the long-term evolution of Population III star clusters, including models with a primordial binary fraction of $f_{\rm b}=0$ and 1, using the $N$-body simulation code PETAR. We comprehensively examine the phenomenon of hierarchical triple black holes in the clusters, specifically focusing on their merging inner binary black holes (BBHs), with post-Newtonian correction, by using the TSUNAMI code. Our findings suggest a high likelihood of the inner BBHs containing IMBHs with masses on the order of $\mathcal{O}(100)M_{\odot}$, and as a result, their merger rate could be up to $0.1{\rm Gpc}^{-3}{\rm yr}^{-3}$. The orbital eccentricities of some merging inner BBHs oscillate over time periodically, known as the Kozai-Lidov oscillation, due to dynamical perturbations. Detectable merging inner BBHs for mHz GW detectors LISA/TianQin/Taiji concentrate within $z<3$. More distant sources would be detectable for CE/ET/LIGO/KAGRA/DECIGO, which are sensitive from $\mathcal{O}(0.1)$Hz to $\mathcal{O}(100)$Hz. Furthermore, compared with merging isolated BBHs, merging inner BBHs affected by dynamical perturbations from tertiary BHs tend to have higher eccentricities, with a significant fraction of sources with eccentricities closing to 1 at mHz bands. GW observations would help constrain formation channels of merging BBHs, whether through isolated evolution or dynamical interaction, by examining eccentricities., Comment: 19 pages, 22 figures, and 3 tables, accepted by MNRAS

Nearly a hundred of binary black holes (BBHs) have been discovered with gravitational-wave signals emitted at their merging events. Thus, it is quite natural to expect that significantly more abundant BBHs with wider separations remain undetected in the universe, or even in our Galaxy. We consider a possibility that star-BH binary candidates may indeed host an inner BBH, instead of a single BH. We present a detailed feasibility study of constraining the binarity of the currently available two targets, Gaia BH1 and Gaia BH2. Specifically, we examine three types of radial velocity (RV) modulations of a tertiary star in star-BBH triple systems; short-term RV modulations induced by the inner BBH, long-term RV modulations induced by the nodal precession, and long-term RV modulations induced by the von Zeipel-Kozai-Lidov oscillations. Direct three-body simulations combined with approximate analytic models reveal that Gaia BH1 system may exhibit observable signatures of the hidden inner BBH if it exists at all. The methodology that we examine here is quite generic, and is expected to be readily applicable to future star-BH binary candidates in a straightforward manner., Comment: 21 pages, 11 figures, 1 table, ApJ, in press

We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with $N$, based on the Multi-Particle Collision (MPC) scheme, previously applied in Plasma Physics. We simulate globular clusters with a realistic number of stellar particles (at least up to several times $10^6$) on a standard workstation. We simulate clusters hosting an intermediate mass black hole (IMBH), probing a broad range of BH-cluster and BH-average-star mass ratios, unrestricted by the computational constraints affecting direct N-body codes. We use either single mass models or models with a Salpeter mass function, with the IMBH initially sitting at the centre. The force exerted by and on the IMBH is evaluated with a direct scheme. We measure the evolution of the Lagrangian radii and core density and velocity dispersion over time. In addition, we study the evolution of the velocity anisotropy profiles. We find that models with an IMBH undergo core collapse at earlier times, the larger the IMBH mass the shallower, with an approximately constant central density at core collapse. The presence of an IMBH tends to lower the central velocity dispersion. These results hold independently of the mass function. For the models with Salpeter MF we observe that equipartition of kinetic energies is never achieved. Orbital anisotropy at large radii appears driven by energetic escapers on radial orbits. We measure the wander radius. Among the results we obtained, which mostly confirm or extend previously known trends that had been established over the range of parameters accessible to direct N-body simulations, we underline that the leptokurtic nature of the IMBH wander radius distribution might lead to IMBHs presenting as off-center more frequently than expected, with implications on observational IMBH detection., Comment: 11 pages, 10 figures. Matching the version of the ms accepted in A&A

The flux-based statistical theory of the non-hierarchical three-body system predicts that the chaotic outcome distribution reduces to the chaotic emissivity function times a known function, the asymptotic flux. Here, we measure the chaotic emissivity function (or equivalently, the absorptivity) through simulations. More precisely, we follow millions of scattering events only up to the point when it can be decided whether the scattering is regular or chaotic. In this way, we measure a tri-variate absorptivity function. Using it, we determine the flux-based prediction for the chaotic outcome distribution over both binary binding energy and angular momentum, and we find good agreement with the measured distribution. This constitutes a detailed confirmation of the flux-based theory, and demonstrates a considerable reduction in computation to determine the chaotic outcome distribution., Comment: 10 pages, 8 figures

Population-synthesis codes are an unique tool to explore the parameter space of massive binary star evolution and binary compact object (BCO) formation. Most population-synthesis codes are based on the same stellar evolution model, limiting our ability to explore the main uncertainties. Here, we present the new version of the code SEVN, which overcomes this issue by interpolating the main stellar properties from a set of pre-computed evolutionary tracks. We describe the new interpolation and adaptive time-step algorithms of SEVN, and the main upgrades on single and binary evolution. With SEVN, we evolved $1.2\times10^9$ binaries in the metallicity range $0.0001\leq Z \leq 0.03$, exploring a number of models for electron-capture, core-collapse and pair-instability supernovae, different assumptions for common envelope, stability of mass transfer, quasi-homogeneous evolution and stellar tides. We find that stellar evolution has a dramatic impact on the formation of single and binary compact objects. Just by slightly changing the overshooting parameter ($\lambda_{\rm ov}=0.4,0.5$) and the pair-instability model, the maximum mass of a black hole can vary from $\approx{60}$ to $\approx{100}\ \mathrm{M}_\odot$. Furthermore, the formation channels of BCOs and the merger efficiency we obtain with SEVN show significant differences with respect to the results of other population-synthesis codes, even when the same binary-evolution parameters are used. For example, the main traditional formation channel of BCOs is strongly suppressed in our models: at high metallicity ($Z\gtrsim{0.01}$) only $<20$% of the merging binary black holes and binary neutron stars form via this channel, while other authors found fractions $>70$%. The local BCO merger rate density of our fiducial models is consistent with the most recent estimates by the LIGO--Virgo--KAGRA collaboration., Comment: Accepted for publication in MNRAS. The SEVN code is available at https://gitlab.com/sevncodes/sevn.git. All the data underlying this article are available in Zenodo at the link https://doi.org/10.5281/zenodo.7794546. All the Jupyter notebooks used to produce the plots in the paper are available in the gitlab repository https://gitlab.com/iogiul/iorio22_plot.git

We present a novel, few-body computational framework designed to shed light on the likelihood of forming intermediate-mass (IM) and supermassive (SM) black holes (BHs) in nuclear star clusters (NSCs) through successive BH mergers, initiated with a single BH seed. Using observationally motivated NSC profiles, we find that the probability of a ${\sim}100 \, M_\odot$ BH to grow beyond ${\sim}1000 \, M_\odot$ through successive mergers ranges from ${\sim}0.1\%$ in low-density, low-mass clusters to nearly $90\%$ in high-mass, high-density clusters. However, in the most massive NSCs, the growth timescale can be very long ($\gtrsim 1\,$Gyr); vice versa, while growth is least likely in less massive NSCs, it is faster there, requiring as little as ${\sim}0.1\,$Gyr. The increased gravitational focusing in systems with lower velocity dispersions is the primary contributor to this behavior. We find that there is a simple "7-strikes-and-you're-in" rule governing the growth of BHs: our results suggest that if the seed survives 7 to 10 successive mergers without being ejected (primarily through gravitational wave recoil kicks), the growing BH will most likely remain in the cluster and will then undergo runaway, continuous growth all the way to the formation of an SMBH (under the simplifying assumption adopted here of a fixed background NSC). Furthermore, we find that rapid mergers enforce a dynamically-mediated "mass gap" between about ${50-300 \, M_\odot}$ in an NSC., Comment: 26 pages, 21 figures, and 3 tables. Submitted to MNRAS. Comments welcome

While there have been many studies examining the stability of hierarchical triple systems, the meaning of ``stability'' is somewhat vague and has been interpreted differently in previous literatures. The present paper focuses on ``Lagrange stability'', which roughly refers to the stability against the escape of a body from the system, or ``disruption'' of the triple system, in contrast to ``Lyapunov-like stability'' that is related to the chaotic nature of the system dynamics. We compute the evolution of triple systems using direct $N$-body simulations up to $10^7 P_\mathrm{out}$, which is significantly longer than previous studies (with $P_\mathrm{out}$ being the initial orbital period of the outer body). We obtain the resulting disruption timescale $T_\mathrm{d}$ as a function of the triple orbital parameters with particular attention to the dependence on the mutual inclination between the inner and outer orbits, $i_\mathrm{mut}$. By doing so, we have clarified explicitly the difference between Lagrange and Lyapunov stabilities in astronomical triples. Furthermore, we find that the von Zeipel-Kozai-Lidov oscillations significantly destabilize inclined triples (roughly with $60^\circ < i_\mathrm{mut} < 150^\circ$) relative to those with $i_\mathrm{mut}=0^\circ$. On the other hand, retrograde triples with $i_\mathrm{mut}>160^\circ$ become strongly stabilized with much longer disruption timescales. We show the sensitivity of the normalized disruption timescale $T_\mathrm{d}/P_\mathrm{out}$ to the orbital parameters of triple system. The resulting $T_\mathrm{d}/P_\mathrm{out}$ distribution is practically more useful in a broad range of astronomical applications than the stability criterion based on the Lyapunov divergence., Comment: 13 pages, 3 figures, ApJ, in press

We examine the stability of hierarchical triple systems using direct $N$-body simulations without adopting a secular perturbation approximation. We estimate their disruption timescales in addition to the mere stable/unstable criterion, with particular attention to the mutual inclination between the inner and outer orbits. First, we improve the fit to the dynamical stability criterion by \citet{Mardling1999,Mardling2001} widely adopted in the previous literature. Especially, we find that that the stability boundary is very sensitive to the mutual inclination; coplanar retrograde triples and orthogonal triples are much more stable and unstable, respectively, than coplanar prograde triples. Next, we estimate the disruption timescales of triples satisfying the stability condition up to $10^9$ times the inner orbital period. The timescales follow the scaling predicted by \citet{Mushkin2020}, especially at high $e_\mathrm{out}$ where their random walk model is most valid. We obtain an improved empirical fit to the disruption timescales, which indicates that the coplanar retrograde triples are significantly more stable than the previous prediction. We furthermore find that the dependence on the mutual inclination can be explained by the energy transfer model based on a parabolic encounter approximation. We also show that the disruption timescales of triples are highly sensitive to the tiny change of the initial parameters, reflecting the genuine chaotic nature of the dynamics of those systems., Comment: 33 pages, 19 figures, 1 table, ApJ, in press. Comments welcome

On 11 February 2016, the LIGO and Virgo scientific collaborations announced the first direct detection of gravitational waves, a signal caught by the LIGO interferometers on 14 September 2015, and produced by the coalescence of two stellar-mass black holes. The discovery represented the beginning of an entirely new way to investigate the Universe. The latest gravitational-wave catalog by LIGO, Virgo and KAGRA brings the total number of gravitational-wave events to 90, and the count is expected to significantly increase in the next years, when additional ground-based and space-born interferometers will be operational. From the theoretical point of view, we have only fuzzy ideas about where the detected events came from, and the answers to most of the five Ws and How for the astrophysics of compact binary coalescences are still unknown. In this work, we review our current knowledge and uncertainties on the astrophysical processes behind merging compact-object binaries. Furthermore, we discuss the astrophysical lessons learned through the latest gravitational-wave detections, paying specific attention to the theoretical challenges coming from exceptional events (e.g., GW190521 and GW190814)., Comment: Invited review paper for Special Issue "Present and Future of Gravitational Wave Astronomy'' - MDPI, Editor: Dr. Gabriele Vajente

Understanding the long-term evolution of hierarchical triple systems is challenging due to its inherent chaotic nature, and it requires computationally expensive simulations. Here we propose a convolutional neural network model to predict the stability of hierarchical triples by looking at their evolution during the first $5 \times 10^5$ inner binary orbits. We employ the regularized few-body code TSUNAMI to simulate $5\times 10^6$ hierarchical triples, from which we generate a large training and test dataset. We develop twelve different network configurations that use different combinations of the triples' orbital elements and compare their performances. Our best model uses 6 time-series, namely, the semimajor axes ratio, the inner and outer eccentricities, the mutual inclination and the arguments of pericenter. This model achieves an area under the curve of over $95\%$ and informs of the relevant parameters to study triple systems stability. All trained models are made publicly available, allowing to predict the stability of hierarchical triple systems $200$ times faster than pure $N$-body methods., Comment: 12 pages, 6 figures, accepted for publication in ApJ

In recent years, an increasing amount of attention is being paid to the gravitational few-body problem and its applications to astrophysical scenarios. Among the main reasons for this renewed interest there is large number of newly discovered exoplanets and the detection of gravitational waves. Here, we present two numerical codes to model three- and few-body systems, called TSUNAMI and OKINAMI. The TSUNAMI code is a direct few-body code with algorithmic regularization, tidal forces and post-Newtonian corrections. OKINAMI is a secular, double-averaged code for stable hierarchical triples. We describe the main methods implemented in our codes, and review our recent results and applications to gravitational-wave astronomy, planetary science and statistical escape theories., Comment: 6 pages, 2 figures

We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory. Our formalism can model the eccentricity during the common envelope inspiral, and it gives results consistent with smoothed particles hydrodynamical simulations. We apply our formalism to common envelope events from binary population synthesis models and find that the final eccentricity distribution resembles the observed distribution of post-common-envelope binaries. Our model can be used for time-resolved common-envelope evolution in population synthesis calculations or as part of binary interactions in direct N-body simulations of star clusters., Comment: 14 Pages, 8 Figures, accepted for publication in PRD

Triple stars and compact objects are ubiquitously observed in nature. Their long-term evolution is complex; in particular, the von-Zeipel-Lidov-Kozai (ZLK) mechanism can potentially lead to highly eccentric encounters of the inner binary. Such encounters can lead to a plethora of interacting binary phenomena, as well as stellar and compact-object mergers. Here we find explicit analytical formulae for the maximal eccentricity, $e_{\rm max}$, of the inner binary undergoing ZLK oscillations, where both the test particle limit (parametrised by the inner-to-outer angular momentum ratio $\eta$) and the double-averaging approximation (parametrised by the period ratio, $\epsilon_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $\eta$ or $\epsilon_{\rm SA} \to 0$) and verify the validity of our model using numerical simulations. We test our results with two accurate numerical N-body codes, $\texttt{Rebound}$ for Newtonian dynamics and $\texttt{Tsunami}$ for general-relativistic (GR) dynamics, and find excellent correspondence. We discuss the implications of our results for stellar triples and both stellar and supermassive triple black hole mergers.

In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission (using post-Newtonian (PN) corrections) and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH-star mergers by up to 75 per cent but it does not affect the number of BH-BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH-BH mergers during binary-single encounters. However, BH-star mergers are largely unaffected by the inclusion of these terms., Comment: 14 pages, 6 figures, 12 tables, accepted for publication in MNRAS

Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass ratio ($10^{-3}-1$), the ratio of the binary semimajor axis to the hardening radius ($10^{-4}-1$), the binary eccentricity ($0.0-0.9$) and the stellar mass ($0.3-3M_{\odot}$). This is a significant extension of the parameter space compared to previous work. We show that the encounter probabilities of all three events are well-described by the encounter cross section, which is proportional to the pericenter distance. The probability of full disruptions by supermassive black hole binaries can be enhanced by up to a factor of $40-50$ or suppressed by up to a factor of $10$, relative to that by single black holes, depending on the binary parameters. Relativistic effects are not important for the primary black hole mass $\leq 10^{7}M_{\odot}$, but can provide an additional enhancement of the full disruption probability by less than a factor of $2-3$ for higher primary black hole masses. We provide a fitting formula for the full disruption probability by the hard supermassive black hole binaries that works for a wide range of parameters. We also find that partial disruption events can occur multiple times before full disruptions or direct captures, and their probabilities can be greater than that of full disruption events by a factor of three. Because partial disruption events can induce stellar spins and mass loss and change the orbits, it can significantly affect the overall full disruption event rate and the shape of the light curves., Comment: 16 pages, 14 figures and 2 tables, accepted for publication in MNRAS

A binary star orbited by an outer companion constitutes a hierarchical triple system. The outer body may excite the eccentricity of the inner binary through the von~Zeipel-Lidov-Kozai (ZLK) mechanism, triggering the gravitational wave (GW) coalescence of the inner binary when its members are compact objects. Here, we study a sample of hierarchical triples with an inner black hole (BH) -- BH binary, BH -- neutron star (NS) binary, and BH -- white dwarf (WD) binary, formed via dynamical interactions in low-mass young star clusters. Our sample of triples was obtained self-consistently from direct N-body simulations of star clusters which included up-to-date stellar evolution. We find that the inner binaries in our triples cannot merge via GW radiation alone, and the ZLK mechanism is essential to trigger their coalescence. Contrary to binaries assembled dynamically in young star clusters, binary BHs merging in triples have preferentially low mass ratios (q ~ 0.3) and higher primary masses (m_p > 40 MSun). We derive a local merger rate density of 0.60, 0.11 and 0.5 yr^-1 Gpc^-3 for BH-BH, BH-NS and BH-WD binaries, respectively. Additionally, we find that merging binaries have high eccentricities across the GW spectrum, including the LIGO-Virgo-KAGRA (LVK), LISA, and DECIGO frequencies. About 7% of BH-BH and 60% of BH-NS binaries will have detectable eccentricities in the LVK band. Our results indicate that the eccentricity and the mass spectrum of merging binaries are the strongest features for the identification of GW mergers from triples., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitational wave (GW) observations. Notably, Pop III and EMP ($< 10^{-2}$ $Z_\odot$) binary stars yield most of the pair instability (PI) mass gap events with $m_1 = 65$--$130$ $M_\odot$. Pop III binary stars contribute more to the PI mass gap events with increasing redshift, and all the PI mass gap events have the Pop III origin at redshifts $\gtrsim 8$. Our result can be assessed by future GW observations in the following two points. First, there are no binary BHs with $m_1=100$--$130$ $M_\odot$ in our result, and thus the $m_1$ distribution should suddenly drop in the range of $m_1=100$--$130$ $M_\odot$. Second, the PI mass gap event rate should increase toward higher redshift up to $\sim 11$, since those events mainly originate from the Pop III binary stars. We find that the following three assumptions are needed to reproduce the current GW observations: a top-heavy stellar initial mass function and the presence of close binary stars for Pop III and EMP binary stars, and inefficient convective overshoot in the main-sequence phase of stellar evolution. Without any of the above, the number of PI mass gap events becomes too low to reproduce current GW observations., Comment: 26 pages, 13 figures, 2 tables, ApJ accepted. BSEEMP, a binary population synthesis code used in this paper, is published on github: https://github.com/atrtnkw/bseemp

We present numerical simulations of the gravitational three-body problem, in which three particles lie at rest close to the vertices of an equilateral triangle. In the unperturbed problem, the three particles fall towards the center of mass of the system to form a three-body collision, or singularity, where the particles overlap in space and time. By perturbing the initial positions of the particles, we are able to study chaos in the vicinity of the singularity. Here we cover both the singular region close to the unperturbed configuration and the binary-single scattering regime where one side of the triangle is very short compared to the other two. We make phase space plots to study the regular and ergodic subsets of our simulations and compare them with the outcomes expected from the statistical escape theory of the three-body problem. We further provide fits to the ergodic subset to characterize the properties of the left-over binaries. We identify the discrepancy between the statistical theory and the simulations in the regular subset of interactions, which only exhibits weak chaos. As we decrease the scale of the perturbations in the initial positions, the phase space becomes entirely dominated by regular interactions, according to our metric for chaos. Finally, we show the effect of general relativity corrections by simulating the same scenario with the inclusion of post-Newtonian corrections to the equations of motion., Comment: 23 pages, 9 figures, 1 table, accepted in SciPost Physics Core

We performed direct N-body simulations of open clusters with four different metallicities. To investigate the effective spins of merging binary black holes (BBHs) originated from these open clusters, we calculated the spin evolution of Wolf-Rayet (WR) stars with close companion stars (BBH progenitors), taking into account stellar wind mass loss and tidal spin-up of the WR stars. We found that BBH progenitors with smaller semi-major axes evolve to merging BBHs with greater effective spins because of strong tidal forces. In the local Universe, about 16% of merging BBHs get effective spins larger than 0.1 even if BHs and their progenitors do not get spin angular momenta other than tidal forces exerted by their companion stars. If we assume that WR stars have flat and isotropic distribution of dimensionless spins just after common envelope phases, the effective spin distribution of merging BBHs is similar to that inferred from gravitational wave observations with LIGO and Virgo., Comment: 8 pages, 10 figures, submitted to MNRAS

Recent studies indicate that the progenitors of merging black hole (BH) binaries from young star clusters can undergo a common envelope phase just like isolated binaries. If the stars emerge from the common envelope as naked cores, tidal interactions can efficiently synchronize their spins before they collapse into BHs. Contrary to the isolated case, these binary BHs can also undergo dynamical interactions with other BHs in the cluster before merging. The interactions can tilt the binary orbital plane, leading to spin-orbit misalignment. We estimate the spin properties of merging binary BHs undergoing this scenario by combining up-to-date binary population synthesis and accurate few-body simulations. We show that post-common envelope binary BHs are likely to undergo only a single encounter, due to the high binary recoil velocity and short coalescence times. Adopting conservative limits on the binary-single encounter rates, we obtain a local BH merger rate density of ~6.6 yr^-1 Gpc^-3. Assuming low (<0.2) natal BH spins, this scenario reproduces the trends in the distributions of effective spin Xeff and precession parameters Xp inferred from GWTC-2, including the peaks at (Xeff, Xp) ~ (0.1, 0.2) and the tail at negative Xeff values., Comment: 10 pages, 6 figures, accepted for publication in MNRAS

We present an extensive comparison between the statistical properties of non-hierarchical three-body systems and the corresponding three-body theoretical predictions. We perform and analyze 1 million realizations for each different initial condition considering equal and unequal mass three-body systems to provide high accuracy statistics. We measure 4 quantities characterizing the statistical distribution of ergodic disintegrations: escape probability of each body, the characteristic exponent for escapes by a narrow margin, predicted absorptivity as a function of binary energy and binary angular momentum, and, finally, the lifetime distribution. The escape probabilities are shown to be in agreement down to the 1% level with the emissivity-blind, flux-based theoretical prediction. This represents a leap in accuracy compared to previous three-body statistical theories. The characteristic exponent at the threshold for marginally unbound escapes is an emissivity-independent flux-based prediction, and the measured values are found to agree well with the prediction. We interpret both tests as strong evidence for the flux-based three-body statistical formalism. The predicted absorptivity and lifetime distributions are measured to enable future tests of statistical theories., Comment: 17 pages, 11 figures

All the giant planets in the solar system host a large number of natural satellites. Moons in extrasolar systems are difficult to detect, but a Neptune-sized exomoon candidate has been recently found around a Jupiter-sized planet in the Kepler-1625bsystem. Due to their relative ease of detection, hot Jupiters (HJs), which reside in close orbits around their host stars with a period of a few days, may be very good candidates to search for exomoons. It is still unknown whether the HJ population can host (or may have hosted) exomoons. One suggested formation channel for HJs is high-eccentricity migration induced by a stellar binary companion combined with tidal dissipation. Here, we investigate under which circumstances an exomoon can prevent or allow high-eccentricity migration of a HJ, and in the latter case, if the exomoon can survive the migration process. We use both semianalytic arguments, as well as direct N-body simulations including tidal interactions. Our results show that massive exomoons are efficient at preventing high-eccentricity migration. If an exomoon does instead allow for planetary migration, it is unlikely that the HJ formed can host exomoons since the moon will either spiral onto the planet or escape from it during the migration process. A few escaped exomoons can become stable planets after the Jupiter has migrated, or by tidally migrating themselves. The majority of the exomoons end up being ejected from the system or colliding with the primary star and the host planet. Such collisions might nonetheless leave observable features, such as a debris disc around the primary star or exorings around the close-in giant., Comment: 12 Pages, 6 Figures, accepted for publication in MNRAS

We present the merger rate density of Population (Pop.) III binary black holes (BHs) by means of a widely-used binary population synthesis code BSE with extensions to very massive and extreme metal-poor stars. We consider not only low-mass BHs (lBHs: $5-50 M_\odot$) but also high-mass BHs (hBHs: $130-200 M_\odot$), where lBHs and hBHs are below and above the pair-instability mass gap ($50-130 M_\odot$), respectively. Pop. III BH-BHs can be categorized into three subpopulations: BH-BHs without hBHs (hBH0s: $m_{\rm tot} \lesssim 100 M_\odot$), with one hBH (hBH1s: $m_{\rm tot} \sim 130-260 M_\odot$), and with two hBHs (hBH2s: $m_{\rm tot} \sim 270-400 M_\odot$), where $m_{\rm tot}$ is the total mass of a BH-BH. Their merger rate densities at the current universe are $\sim 0.1$ yr$^{-1}$ Gpc$^{-3}$ for hBH0s, and $\sim 0.01$ yr$^{-1}$ Gpc$^{-3}$ for the sum of hBH1s and hBH2s, provided that the mass density of Pop. III stars is $\sim 10^{13} M_\odot$ Gpc$^{-3}$. These rates are modestly insensitive to initial conditions and single star models. The hBH1 and hBH2 mergers can dominate BH-BHs with hBHs discovered in near future. They have low effective spins $\lesssim 0.2$ in the current universe. The number ratio of the hBH2s to the hBH1s is high, $\gtrsim 0.1$. We also find BHs in the mass gap (up to $\sim 85 M_\odot$) merge. These merger rates can be reduced to nearly zero if Pop. III binaries are always wide ($\gtrsim 100 R_\odot$), and if Pop. III stars always enter into chemically homogeneous evolution. The presence of close Pop. III binaries ($\sim 10 R_\odot$) are crucial for avoiding the worst scenario., Comment: 38 pages, 33 figures, accepted for publication in ApJ

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $\rho_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively., Comment: 14 pages, 10 figures, 5 tables, accepted for publication in MNRAS

We study chaos and L\'evy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the ergodic interactions and L\'evy flight interactions. This enables us to study the three-body lifetime distribution in more detail by decomposing it into the individual distributions from the different kinds of interactions. We observe that ergodic interactions follow an exponential decay distribution similar to that of radioactive decay. Meanwhile, L\'evy flight interactions follow a power-law distribution. L\'evy flights in fact dominate the tail of the general three-body lifetime distribution, providing conclusive evidence for the speculated connection between power-law tails and L\'evy flight interactions. We propose a new physically-motivated model for the lifetime distribution of three-body systems and discuss how it can be used to extract information about the underlying ergodic and L\'evy flight interactions. We discuss mass ejection probabilities in three-body systems in the ergodic limit and compare it to previous ergodic formalisms. We introduce a novel mechanism for a three-body relaxation process and discuss its relevance in general three-body systems., Comment: 20 pages, 14 figures

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 \sigma$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo)., Comment: accepted for publication in AJ

High-density cusps of compact remnants are expected to form around supermassive black holes (SMBHs) in galactic nuclei via dynamical friction and two-body relaxation. Due to the high density, binaries in orbit around the SMBH can frequently undergo close encounters with compact remnants from the cusp. This can affect the gravitational wave merger rate of compact binaries in galactic nuclei. We investigated this process by means of high accuracy few-body simulations, performed with a novel Monte Carlo approach. We find that, around a SgrA*-like SMBH, three-body encounters increase the number of mergers by a factor of 3. This occurs because close encounters can reorient binaries with respect to their orbital plane around the SMBH, increasing the number of Kozai-Lidov induced mergers. We obtain a binary black hole merger rate of 1.6x10^-6/yr per Milky Way-like nucleus., Comment: Proceedings of the IAU symposium 351 (MODEST-19)

The origin of the molecular clouds orbiting SgrA* is one of the most debated questions about our Galactic Centre. Here, we present the hydrodynamic simulation of a molecular cloud infalling towards SgrA*, performed with the adaptive-mesh-refinement code RAMSES. Through such simulation, we propose that the circumnuclear disc and the +20 km/s cloud originated from the same tidal disruption episode, occurred less than 1 Myr ago. We also show that recent star formation is to be expected in the +20 km/s cloud, as also suggested by recent observations., Comment: 10 pages, 10 figures. To be published in MNRAS

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs., Comment: 13 pages, 11 figures, Accepted for publication in MNRAS

We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects the breakup of triples by tidally limiting the system, so that the triples break earlier with lower breakup velocity, leaving behind slightly larger binaries (relative to the isolated case). The breakup direction becomes anisotropic and tends to follow the shape of the Hill region of the triple, favouring breakups in the radial direction. Furthermore, the tidal field can torque the system, leading to angular momentum exchanges between the triple and its orbit about the SMBH. This process changes the properties of the final binary, depending on the initial angular momentum of the triple. Finally, the tidal field also affects binary-single encounters: binaries tend to become both harder and more eccentric with respect to encounters that occur in isolation. Consequently, single-binary scattering in a deep Keplerian potential produces binaries with shorter gravitational wave merger timescales., Comment: 15 pages, 11 Figures, 4 tables, expanded and revised version, accepted for publication in ApJ

Recent spectroscopic analysis has set an upper limit to the age of the S-stars, the ~30 B-type stars in highly eccentric orbits around the supermassive black hole (SMBH) in the Galactic center. The inferred age (<15 Myr) is in tension with the binary break-up scenario proposed to explain their origin. However, the new estimate is compatible with the age of the disk of O-type stars that lies at a farther distance from the SMBH. Here we investigate a new formation scenario, assuming that both S-stars and the O-type stars were born in the same disk around SgrA*. We simulate encounters between binaries of the stellar disk and stellar black holes from a dark cusp around SgrA*. We find that B-type binaries can be easily broken up by the encounters and their binary components are kicked into highly eccentric orbits around the SMBH. In contrast, O-type binaries are less frequently disrupted and their members remain in low eccentricity orbits. This mechanism can reproduce 12 S-stars just by assuming that the binaries initially lie within the stellar disk as observed nowadays. To reproduce all the S-stars, the original disk must have been extended down to 0.006 pc. However in this case many B- and O-type stars remain in low eccentricity orbits below 0.03 pc, in contrast with the observations. Therefore, some other mechanism is necessary to disrupt the disk below 0.03 pc. This scenario can also explain the high eccentricity of the G-objects, if they have a stellar origin., Comment: 9 pages, 8 figures, accepted in ApJ

Studying the formation and evolution of black hole binaries (BHBs) is essential for the interpretation of current and forthcoming gravitational wave (GW) detections. We investigate the statistics of BHBs that form from isolated binaries, by means of a new version of the SEVN population-synthesis code. SEVN integrates stellar evolution by interpolation over a grid of stellar evolution tracks. We upgraded SEVN to include binary stellar evolution processes and we used it to evolve a sample of $1.5\times{}10^8$ binary systems, with metallicity in the range $\left[10^{-4};4\times 10^{-2}\right]$. From our simulations, we find that the mass distribution of black holes (BHs) in double compact-object binaries is remarkably similar to the one obtained considering only single stellar evolution. The maximum BH mass we obtain is $\sim 30$, $45$ and $55\, \mathrm{M}_\odot$ at metallicity $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. A few massive single BHs may also form ($\lesssim 0.1\%$ of the total number of BHs), with mass up to $\sim 65$, $90$ and $145\, \mathrm{M}_\odot$ at $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and $10^{-4}$, respectively. These BHs fall in the mass gap predicted from pair-instability supernovae. We also show that the most massive BHBs are unlikely to merge within a Hubble time. In our simulations, merging BHs like GW151226 and GW170608, form at all metallicities, the high-mass systems (like GW150914, GW170814 and GW170104) originate from metal poor ($Z\lesssim{}6\times 10^{-3}$) progenitors, whereas GW170729-like systems are hard to form, even at $Z = 10^{-4}$. The BHB merger rate in the local Universe obtained from our simulations is $\sim 90 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$, consistent with the rate inferred from LIGO-Virgo data., Comment: 26 pages, 14 figures, MNRAS accepted

Planet host stars with well-constrained ages provide a rare window to the time domain of planet formation and evolution. The NASA K2 mission has enabled the discovery of the vast majority of known planets transiting stars in clusters, providing a valuable sample of planets with known ages and radii. We present the discovery of two planets transiting K2-264, an M2 dwarf in the intermediate age (600-800 Myr) Praesepe open cluster (also known as the Beehive Cluster, M44, or NGC 2632), which was observed by K2 during Campaign 16. The planets have orbital periods of 5.8 and 19.7 days, and radii of $2.2 \pm 0.2 $ and $2.7 \pm 0.2$ $R_\oplus$, respectively, and their equilibrium temperatures are $496 \pm 10$ and $331 \pm 7$ $K$, making this a system of two warm sub-Neptunes. When placed in the context of known planets orbiting field stars of similar mass to K2-264, these planets do not appear to have significantly inflated radii, as has previously been noted for some cluster planets. As the second known system of multiple planets transiting a star in a cluster, K2-264 should be valuable for testing theories of photoevaporation in systems of multiple planets. Follow-up observations with current near-infrared (NIR) spectrographs could yield planet mass measurements, which would provide information about the mean densities and compositions of small planets soon after photoevaporation is expected to have finished. Follow-up NIR transit observations using Spitzer or large ground-based telescopes could yield improved radius estimates, further enhancing the characterization of these interesting planets., Comment: 12 pages, 11 figures, 2 tables, published in MNRAS; see also secular dynamics at https://youtu.be/MSLodufHRes

We investigate the formation of circumnuclear gas structures from the tidal disruption of molecular clouds in galactic nuclei, by means of smoothed particle hydrodynamics simulations. We model galactic nuclei as composed of a supermassive black hole (SMBH) and a nuclear star cluster (NSC) and consider different mass ratios between the two components. We find that the relative masses of the SMBH and the NSC have a deep impact on the morphology of the circumnuclear gas. Extended disks form only inside the sphere of influence of the SMBH. In contrast, compact rings naturally form outside the SMBH's sphere of influence, where the gravity is dominated by the NSC. This result is in agreement with the properties of the Milky Way's circumnuclear ring, which orbits outside the SMBH sphere of influence. Our results indicate that compact circumnuclear rings can naturally form outside the SMBH sphere of influence., Comment: Accepted for publication in ApJ. 12 pages, 6 figures, 3 tables. Comments welcome

We present a drag force model for evolving chaotic few-body interactions with the inclusion of orbital energy losses, such as tidal dissipation and gravitational wave (GW) emission. The main effect from such losses is the formation of two-body captures, that for compact objects result in GW mergers, and for stars lead to either compact binaries, mergers or disruptions. Studying the inclusion of energy loss terms in few-body interactions is therefore likely to be important for modeling and understanding the variety of transients that soon will be observed by current and upcoming surveys. However, including especially tides in few-body codes has been shown to be technically difficult and computationally heavy, which has lead to very few systematic tidal studies. In this paper we derive a drag force term that can be used to model the effects from tidal, as well as other, energy losses in few-body interactions, if the two-body orbit averaged energy loss is known a priori. This drag force model is very fast to evolve, and gives results in agreement with other approaches, including the impulsive and affine tide approximations., Comment: accepted version. MNRAS

Recent observations suggest ongoing planet formation in the innermost parsec of the Galactic center (GC). The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. Photoevaporation by the ultraviolet field of young stars, combined with ongoing tidal disruption, could enhance the near-infrared luminosity of such starless planets, making their detection possible even with current facilities. In this paper, we investigate the chance of planet tidal captures by means of high-accuracy N-body simulations exploiting Mikkola's algorithmic regularization. We consider both planets lying in the clockwise (CW) disk and planets initially bound to the S-stars. We show that tidally captured planets remain on orbits close to those of their parent star. Moreover, the semi-major axis of the planet orbit can be predicted by simple analytic assumptions in the case of prograde orbits. We find that starless planets that were initially bound to CW disk stars have mild eccentricities and tend to remain in the CW disk. However, we speculate that angular momentum diffusion and scattering with other young stars in the CW disk might bring starless planets on low-angular momentum orbits. In contrast, planets initially bound to S-stars are captured by the SMBH on highly eccentric orbits, matching the orbital properties of the G1 and G2 clouds. Our predictions apply not only to planets but also to low-mass stars initially bound to the S-stars and tidally captured by the SMBH., Comment: 12 pages, 12 figures, 1 table, accepted for publication in The Astrophysical Journal