Showing 1–50 of 53 results for author: Trani, A A

Swift-BAT GUANO follow-up of gravitational-wave triggers in the third LIGO-Virgo-KAGRA observing run

Authors: Gayathri Raman, Samuele Ronchini, James Delaunay, Aaron Tohuvavohu, Jamie A. Kennea, Tyler Parsotan, Elena Ambrosi, Maria Grazia Bernardini, Sergio Campana, Giancarlo Cusumano, Antonino D'Ai, Paolo D'Avanzo, Valerio D'Elia, Massimiliano De Pasquale, Simone Dichiara, Phil Evans, Dieter Hartmann, Paul Kuin, Andrea Melandri, Paul O'Brien, Julian P. Osborne, Kim Page, David M. Palmer, Boris Sbarufatti, Gianpiero Tagliaferri , et al. (1797 additional authors not shown)

Abstract: We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wav… ▽ More We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 50 pages, 10 figures, 4 tables

Detecting environmental effects in gravitational waves from binaries perturbed by periodic forces

Authors: Lorenz Zwick, Christopher Tiede, Alessandro A. Trani, Andrea Derdzinski, Zoltan Haiman, Daniel J. D'Orazio, Johan Samsing

Abstract: 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 se… ▽ More 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. △ Less

Submitted 15 May, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: Submitted to PRD. Comments welcome

The boring history of Gaia BH3 from isolated binary evolution

Authors: Giuliano Iorio, Stefano Torniamenti, Michela Mapelli, Marco Dall'Amico, Alessandro A. Trani, Sara Rastello, Cecilia Sgalletta, Stefano Rinaldi, Guglielmo Costa, Bera A. Dhal-Lahtinen, Gaston J. Escobar, Erika Korb, M. Paola Vaccaro, Elena Lacchin, Benedetta Mestichelli, Ugo Niccolò di Carlo, Mario Spera, Manuel Arca Sedda

Abstract: 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 our open-source population synthesis code SEVN to evolve… ▽ More 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 our open-source population synthesis code SEVN to evolve $3.2 \times 10^8$ binaries exploring 16 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 progenitors do not undergo any Roche-lobe overflow episode during their entire evolution, so that the final orbital properties of the BH-star system are essentially determined at the core collapse of the primary star. Low natal kicks ($\approx$ 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 $4 \times 10^{-8}$ M$_\odot^{-1}$, representing $\approx 3\%$ of the whole simulated BH-star population. We expect up to $\sim 3000$ BH-star systems in the Galactic halo formed through isolated evolution, of which $\sim 100$ are compatible with Gaia BH3 alike. Gaia BH3-like systems represent a common product of isolated binary evolution at low metallicity ($Z<0.01$), but given the steep density profile of the Galactic halo we do not expect more than one in the halo at the observed distance of Gaia BH3. Considering the estimated formation efficiency for both the isolated and dynamical formation channel, we conclude that they are almost equally likely to explain the origin of Gaia BH3. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Submitted to A&A, comments welcome

Isles of regularity in a sea of chaos amid the gravitational three-body problem

Authors: Alessandro Alberto Trani, Nathan W. C. Leigh, Tjarda C. N. Boekholt, Simon Portegies Zwart

Abstract: 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 c… ▽ More 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. △ Less

Submitted 3 July, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 15 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

pAGN: the one-stop solution for AGN disc modeling

Authors: Daria Gangardt, Alessandro Alberto Trani, Clément Bonnerot, Davide Gerosa

Abstract: 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 &… ▽ More 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. △ Less

Submitted 14 May, 2024; v1 submitted 29 February, 2024; originally announced March 2024.

Comments: 18 pages, 8 figures, 2 tables, code available at https://github.com/DariaGangardt/pAGN

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 4, (2024), pp.3689-3705

On Binary Formation from Three Initially Unbound Bodies

Authors: Dany Atallah, Newlin C. Weatherford, Alessandro A. Trani, Frederic Rasio

Abstract: 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 pred… ▽ More 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. △ Less

Submitted 22 May, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: 28 pages, 16 figures, submitted to ApJ. Comments welcome

Three-body encounters in black hole discs around a supermassive black hole: The disc velocity dispersion and the Keplerian tidal field determine the eccentricity and spin-orbit alignment of gravitational wave mergers

Authors: Alessandro Alberto Trani, Stefano Quaini, Monica Colpi

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 12 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Active learning meets fractal decision boundaries: a cautionary tale from the Sitnikov three-body problem

Authors: Nicolas Payot, Mario Pasquato, Alessandro Alberto Trani, Yashar Hezaveh, Laurence Perreault-Levasseur

Abstract: Chaotic systems such as the gravitational N-body problem are ubiquitous in astronomy. Machine learning (ML) is increasingly deployed to predict the evolution of such systems, e.g. with the goal of speeding up simulations. Strategies such as active Learning (AL) are a natural choice to optimize ML training. Here we showcase an AL failure when predicting the stability of the Sitnikov three-body prob… ▽ More Chaotic systems such as the gravitational N-body problem are ubiquitous in astronomy. Machine learning (ML) is increasingly deployed to predict the evolution of such systems, e.g. with the goal of speeding up simulations. Strategies such as active Learning (AL) are a natural choice to optimize ML training. Here we showcase an AL failure when predicting the stability of the Sitnikov three-body problem, the simplest case of N-body problem displaying chaotic behavior. We link this failure to the fractal nature of our classification problem's decision boundary. This is a potential pitfall in optimizing large sets of N-body simulations via AL in the context of star cluster physics, galactic dynamics, or cosmology. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 7+3 pages, 4 figures, Machine Learning and the Physical Sciences Workshop, NeurIPS 2023

The search for the lost attractor

Authors: Mario Pasquato, Syphax Haddad, Pierfrancesco Di Cintio, Alexandre Adam, Pablo Lemos, Noé Dia, Mircea Petrache, Ugo Niccolò Di Carlo, Alessandro Alberto Trani, Laurence Perreault-Levasseur, Yashar Hezaveh

Abstract: N-body systems characterized by inverse square attractive forces may display a self similar collapse known as the gravo-thermal catastrophe. In star clusters, collapse is halted by binary stars, and a large fraction of Milky Way clusters may have already reached this phase. It has been speculated -- with guidance from simulations -- that macroscopic variables such as central density and velocity d… ▽ More N-body systems characterized by inverse square attractive forces may display a self similar collapse known as the gravo-thermal catastrophe. In star clusters, collapse is halted by binary stars, and a large fraction of Milky Way clusters may have already reached this phase. It has been speculated -- with guidance from simulations -- that macroscopic variables such as central density and velocity dispersion are governed post-collapse by an effective, low-dimensional system of ODEs. It is still hard to distinguish chaotic, low dimensional motion, from high dimensional stochastic noise. Here we apply three machine learning tools to state-of-the-art dynamical simulations to constrain the post collapse dynamics: topological data analysis (TDA) on a lag embedding of the relevant time series, Sparse Identification of Nonlinear Dynamics (SINDY), and Tests of Accuracy with Random Points (TARP). △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted by ML4PS workshop at NeurIPS 2023

Merging Hierarchical Triple Black Hole Systems with Intermediate-mass Black Holes in Population III Star Clusters

Authors: Shuai Liu, Long Wang, Yi-Ming Hu, Ataru Tanikawa, Alessandro A. Trani

Abstract: 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 that are 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… ▽ More 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 that are 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 \textsc{petar}. We comprehensively examine the phenomenon of hierarchical triple black holes in the clusters, specifically focusing on the merging of their inner binary black holes (BBHs), with post-Newtonian correction, by using the \textsc{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}$. In the model with $f_{\rm b}=0$, the evolution of these merging inner BBHs is predominantly driven by their gravitational wave radiation at an early time, but their evolutionary dynamics are dominated by the interaction between them and tertiary BHs in the case with $f_{\rm b}=1$. The orbital eccentricities of some merging inner BBHs oscillate over time periodically, known as the Kozai-Lidov oscillation, due to dynamical perturbations. Furthermore, the merging inner BBHs tend to have highly eccentric orbits at low frequency range, some of them with low redshift would be detected by LISA/TianQin. In the higher frequency range, the merging inner BBHs across a wider redshift range would be detected by DECIGO/ET/CE/LIGO/KAGRA. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 15 pages, 16 figures, and 3 tables

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Constraining the binarity of black hole candidates: a proof-of-concept study of Gaia BH1 and Gaia BH2

Authors: Toshinori Hayashi, Yasushi Suto, Alessandro A. Trani

Abstract: 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… ▽ More 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. △ Less

Submitted 29 August, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: 21 pages, 11 figures, 1 table, ApJ, in press

Measurement of three-body chaotic absorptivity predicts chaotic outcome distribution

Authors: Viraj Manwadkar, Alessandro Alberto Trani, Barak Kol

Abstract: 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 i… ▽ More 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. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 10 pages, 8 figures

Dynamics of intermediate mass black holes in globular clusters. Wander radius and anisotropy profiles

Authors: Pierfrancesco Di Cintio, Mario Pasquato, Luca Barbieri, Alessandro A. Trani, Ugo N. Di Carlo

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 11 pages, 10 figures. Matching the version of the ms accepted in A&A

Journal ref: A&A 673, A8 (2023)

Compact object mergers: exploring uncertainties from stellar and binary evolution with SEVN

Authors: Giuliano Iorio, Michela Mapelli, Guglielmo Costa, Mario Spera, Gastón J. Escobar, Cecilia Sgalletta, Alessandro A. Trani, Erika Korb, Filippo Santoliquido, Marco Dall'Amico, Nicola Gaspari, Alessandro Bressan

Abstract: 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… ▽ More 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 ($λ_{\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. △ Less

Submitted 29 June, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 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

Growing Black Holes through Successive Mergers in Galactic Nuclei: I. Methods and First Results

Authors: Dany Atallah, Alessandro A. Trani, Kyle Kremer, Newlin C. Weatherford, Giacomo Fragione, Mario Spera, Frederic A. Rasio

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: 26 pages, 21 figures, and 3 tables. Submitted to MNRAS. Comments welcome

Lagrange vs. Lyapunov stability of hierarchical triple systems: dependence on the mutual inclination between inner and outer orbits

Authors: Toshinori Hayashi, Alessandro A. Trani, Yasushi Suto

Abstract: 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 ``Lyapun… ▽ More 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. △ Less

Submitted 16 December, 2022; v1 submitted 18 September, 2022; originally announced September 2022.

Comments: 13 pages, 3 figures, ApJ, in press

Dynamical disruption timescales and chaotic behavior of hierarchical triple systems

Authors: Toshinori Hayashi, Alessandro A. Trani, Yasushi Suto

Abstract: 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{Mar… ▽ More 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. △ Less

Submitted 5 September, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 33 pages, 19 figures, 1 table, ApJ, in press. Comments welcome

Compact Binary Coalescences: Astrophysical Processes and Lessons Learned

Authors: Mario Spera, Alessandro Alberto Trani, Mattia Mencagli

Abstract: 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, Vir… ▽ More 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). △ Less

Submitted 30 June, 2022; originally announced June 2022.

Comments: Invited review paper for Special Issue "Present and Future of Gravitational Wave Astronomy'' - MDPI, Editor: Dr. Gabriele Vajente

Journal ref: Galaxies 2022, 10(4), 76

Predicting the Stability of Hierarchical Triple Systems with Convolutional Neural Networks

Authors: Florian Lalande, Alessandro Alberto Trani

Abstract: 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… ▽ More 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. △ Less

Submitted 1 September, 2022; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 12 pages, 6 figures, accepted for publication in ApJ

Modeling gravitational few-body problems with TSUNAMI and OKINAMI

Authors: Alessandro Alberto Trani, Mario Spera

Abstract: 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. T… ▽ More 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. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 6 pages, 2 figures

Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes

Authors: Alessandro Alberto Trani, Steven Rieder, Ataru Tanikawa, Giuliano Iorio, Riccardo Martini, Georgii Karelin, Hila Glanz, Simon Portegies Zwart

Abstract: 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.… ▽ More 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. △ Less

Submitted 23 August, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 14 Pages, 8 Figures, accepted for publication in PRD

Extreme eccentricities of triple systems: Analytic results

Authors: Abhi Mangipudi, Evgeni Grishin, Alessandro A. Trani, Ilya Mandel

Abstract: 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… ▽ More 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 $η$) and the double-averaging approximation (parametrised by the period ratio, $ε_{\rm SA}$) are relaxed, for circular outer orbits. We recover known results in both limiting cases (either $η$ or $ε_{\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. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Influence of tidal dissipation on outcomes of binary-single encounters between stars and black holes in stellar clusters

Authors: Lucas Hellström, Abbas Askar, Alessandro A. Trani, Mirek Giersz, Ross P. Church, Johan Samsing

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 14 pages, 6 figures, 12 tables, accepted for publication in MNRAS

Tidal Disruption Events by Compact Supermassive Black Hole Binaries

Authors: Taeho Ryu, Alessandro A. Trani, Nathan W. C. Leigh

Abstract: 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 rat… ▽ More 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. △ Less

Submitted 2 August, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 16 pages, 14 figures and 2 tables, accepted for publication in MNRAS

Multiparticle collision simulations of dense stellar systems and plasmas

Authors: P. Di Cintio, M. Pasquato, L. Barbieri, H. Bufferand, L. Casetti, G. Ciraolo, U. N. Di Carlo, P. Ghendrih, J. P. Gunn, S. Gupta, H. Kim, S. Lepri, R. Livi, A. Simon-Petit, A. A. Trani, S. -J. Yoon

Abstract: We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range… ▽ More We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range fields, ensuring that simulation time scales with $N\log N$ in the number of particles, as opposed to $N^2$ for direct $N$-body. The collisional relaxation effects are modelled by computing particle interactions based on a collision operator approach that ensures rigorous conservation of energy and momenta and depends only on particles velocities and cell-based integrated quantities. △ Less

Submitted 11 February, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 7 pages, 2 figures. To appear in the proceedings of the 362 IAU symposium "Predictive Power of Computational Astrophysics as a Discovery Tool". Chamonix, France 8-12 Nov. 2021

Compact Object Mergers in Hierarchical Triples from Low-Mass Young Star Clusters

Authors: Alessandro Alberto Trani, Sara Rastello, Ugo N. Di Carlo, Filippo Santoliquido, Ataru Tanikawa, Michela Mapelli

Abstract: 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… ▽ More 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. △ Less

Submitted 17 January, 2022; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Merger rate density of binary black holes through isolated Population I, II, III and extremely metal-poor binary star evolution

Authors: Ataru Tanikawa, Takashi Yoshida, Tomoya Kinugawa, Alessandro A. Trani, Takashi Hosokawa, Hajime Susa, Kazuyuki Omukai

Abstract: 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 gravitat… ▽ More 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. △ Less

Submitted 2 March, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: 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

Journal ref: The Astrophysical Journal, 926, 83, 2022 February

Chaos in the vicinity of a singularity in the Three-Body Problem: The equilateral triangle experiment in the zero angular momentum limit

Authors: Hugo D. Parischewsky, Alessandro A. Trani, Nathan W. C. Leigh

Abstract: 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 particl… ▽ More 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. △ Less

Submitted 7 November, 2023; v1 submitted 13 August, 2021; originally announced August 2021.

Comments: 23 pages, 9 figures, 1 table, accepted in SciPost Physics Core

37 New Validated Planets in Overlapping K2 Campaigns

Authors: J. P. de Leon, J. Livingston, M. Endl, W. D. Cochran, T. Hirano, R. A. Garcia, S. Mathur, K. W. F. Lam, J. Korth, A. A. Trani, F. Dai, E. Diez Alonso, A. Castro-Gonzalez, M. Fridlund, A. Fukui, D. Gandolfi, P. Kabath, M. Kuzuhara, R. Luque, A. B. Savel, H. Gill, C. Dressing, S. Giacalone, N. Narita, E. Palle , et al. (2 additional authors not shown)

Abstract: We analysed 68 candidate planetary systems first identified during Campaigns 5 and 6 (C5 and C6) of the NASA \textit{K2} mission. We set out to validate these systems by using a suite of follow-up observations, including adaptive optics, speckle imaging, and reconnaissance spectroscopy. The overlap between C5 with C16 and C18, and C6 with C17, yields lightcurves with long baselines that allow us t… ▽ More We analysed 68 candidate planetary systems first identified during Campaigns 5 and 6 (C5 and C6) of the NASA \textit{K2} mission. We set out to validate these systems by using a suite of follow-up observations, including adaptive optics, speckle imaging, and reconnaissance spectroscopy. The overlap between C5 with C16 and C18, and C6 with C17, yields lightcurves with long baselines that allow us to measure the transit ephemeris very precisely, revisit single transit candidates identified in earlier campaigns, and search for additional transiting planets with longer periods not detectable in previous works. Using \texttt{vespa}, we compute false positive probabilities of less than 1\% for 37 candidates orbiting 29 unique host stars and hence statistically validate them as planets. These planets have a typical size of $2.2R_{\oplus}$ and orbital periods between 1.99 and 52.71 days. We highlight interesting systems including a sub-Neptune with the longest period detected by \textit{K2}, sub-Saturns around F stars, several multi-planetary systems in a variety of architectures. These results show that a wealth of planetary systems still remains in the \textit{K2} data, some of which can be validated using minimal follow-up observations and taking advantage of analyses presented in previous catalogs. △ Less

Submitted 6 September, 2021; v1 submitted 12 August, 2021; originally announced August 2021.

Comments: 25 pages, 12 figures, accepted for publication in MNRAS

Spin distribution of binary black holes formed in open clusters

Authors: Jun Kumamoto, Michiko S. Fujii, Alessandro A. Trani, Ataru Tanikawa

Abstract: 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… ▽ More 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. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 8 pages, 10 figures, submitted to MNRAS

Spin misalignment of black hole binaries from young star clusters: implications for the origin of gravitational waves events

Authors: Alessandro Alberto Trani, Ataru Tanikawa, Michiko S. Fujii, Nathan W. C. Leigh, Jun Kumamoto

Abstract: 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 int… ▽ More 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. △ Less

Submitted 2 April, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

Testing the Flux-based statistical prediction of the Three-Body Problem

Authors: Viraj Manwadkar, Barak Kol, Alessandro A. Trani, Nathan W. C. Leigh

Abstract: 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 dist… ▽ More 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. △ Less

Submitted 10 January, 2021; originally announced January 2021.

Comments: 17 pages, 11 figures

The ominous fate of exomoons around hot Jupiters in the high-eccentricity migration scenario

Authors: Alessandro A. Trani, Adrian S. Hamers, Aaron Geller, Mario Spera

Abstract: 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,… ▽ More 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. △ Less

Submitted 25 October, 2020; v1 submitted 31 August, 2020; originally announced August 2020.

Comments: 12 Pages, 6 Figures, accepted for publication in MNRAS

Merger rate density of Population III binary black holes below, above, and in the pair-instability mass gap

Authors: Ataru Tanikawa, Hajime Susa, Takashi Yoshida, Alessandro A. Trani, Tomoya Kinugawa

Abstract: 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 (… ▽ More 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. △ Less

Submitted 28 March, 2021; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 38 pages, 33 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 910:30 (31pp), 2021 March

Binary black holes in young star clusters: the impact of metallicity

Authors: Ugo N. Di Carlo, Michela Mapelli, Nicola Giacobbo, Mario Spera, Yann Bouffanais, Sara Rastello, Filippo Santoliquido, Mario Pasquato, Alessandro Ballone, Alessandro A. Trani, Stefano Torniamenti, Francesco Haardt

Abstract: 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 (d… ▽ More 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 $ρ_{\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. △ Less

Submitted 26 July, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: 14 pages, 10 figures, 5 tables, accepted for publication in MNRAS

Chaos and Lévy Flights in the Three-Body Problem

Authors: Viraj Manwadkar, Alessandro A. Trani, Nathan W. C. Leigh

Abstract: We study chaos and Lévy 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… ▽ More We study chaos and Lévy 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évy 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évy flight interactions follow a power-law distribution. Lévy 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évy 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évy 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. △ Less

Submitted 19 June, 2020; v1 submitted 11 April, 2020; originally announced April 2020.

Comments: 20 pages, 14 figures

LHS 1815b: The First Thick-Disk Planet Detected By TESS

Authors: Tianjun Gan, Avi Shporer, John H. Livingston, Karen A. Collins, Shude Mao, Alessandro A. Trani, Davide Gandolfi, Teruyuki Hirano, Rafael Luque, Keivan G. Stassun, Carl Ziegler, Steve B. Howell, Coel Hellier, Jonathan M. Irwin, Jennifer G. Winters, David R. Anderson, César Briceño, Nicholas Law, Andrew W. Mann, Xavier Bonfils, Nicola Astudillo-Defru, Eric L. N. Jensen, Guillem Anglada-Escudé, George R. Ricker, Roland Vanderspek , et al. (13 additional authors not shown)

Abstract: 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… ▽ More 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 σ$ 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). △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: accepted for publication in AJ

Do three-body encounters in galactic nuclei affect compact binary merger rates?

Authors: Alessandro Alberto Trani

Abstract: 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 inv… ▽ More 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. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Comments: Proceedings of the IAU symposium 351 (MODEST-19)

A common origin for the circumnuclear disc and the nearby molecular clouds in the Galactic Centre

Authors: Alessandro Ballone, Michela Mapelli, Alessandro Alberto Trani

Abstract: 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,… ▽ More 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. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: 10 pages, 10 figures. To be published in MNRAS

The evolution of kicked stellar-mass black holes in star cluster environments II. Rotating star clusters

Authors: Jeremy J. Webb, Nathan W. C. Leigh, Roberto Serrano, Jillian Bellovary, K. E. Saavik Ford, Barry McKernan, Mario Spera, Alessandro A. Trani

Abstract: 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 fra… ▽ More 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. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 13 pages, 11 figures, Accepted for publication in MNRAS

The Keplerian three-body encounter II. Comparisons with isolated encounters and impact on gravitational wave merger timescales

Authors: Alessandro A. Trani, Mario Spera, Nathan W. C. Leigh, Michiko S. Fujii

Abstract: 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 t… ▽ More 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. △ Less

Submitted 24 September, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: 15 pages, 11 Figures, 4 tables, expanded and revised version, accepted for publication in ApJ

The Keplerian three-body encounter I. Insights on the origin of the S-stars and the G-objects in the Galactic center

Authors: Alessandro A. Trani, Michiko S. Fujii, Mario Spera

Abstract: 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 f… ▽ More 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. △ Less

Submitted 1 April, 2019; v1 submitted 19 September, 2018; originally announced September 2018.

Comments: 9 pages, 8 figures, accepted in ApJ

Merging black hole binaries with the SEVN code

Authors: Mario Spera, Michela Mapelli, Nicola Giacobbo, Alessandro Alberto Trani, Alessandro Bressan, Guglielmo Costa

Abstract: 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 up… ▽ More 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. △ Less

Submitted 10 April, 2019; v1 submitted 12 September, 2018; originally announced September 2018.

Comments: 26 pages, 14 figures, MNRAS accepted

K2-264: A transiting multi-planet system in the Praesepe open cluster

Authors: John H. Livingston, Fei Dai, Teruyuki Hirano, Davide Gandolfi, Alessandro A. Trani, Grzegorz Nowak, William D. Cochran, Michael Endl, Simon Albrecht, Oscar Barragan, Juan Cabrera, Szilard Csizmadia, Jerome P. de Leon, Hans Deeg, Philipp Eigmüller, Anders Erikson, Malcolm Fridlund, Akihiko Fukui, Sascha Grziwa, Eike W. Guenther, Artie P. Hatzes, Judith Korth, Masayuki Kuzuhara, Pilar Montañes, Norio Narita , et al. (9 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 14 February, 2019; v1 submitted 6 September, 2018; originally announced September 2018.

Comments: 12 pages, 11 figures, 2 tables, published in MNRAS; see also secular dynamics at https://youtu.be/MSLodufHRes

Forming circumnuclear disks and rings in galactic nuclei: a competition between supermassive black hole and nuclear star cluster

Authors: Alessandro A. Trani, Michela Mapelli, Alessandro Ballone

Abstract: 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 an… ▽ More 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. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: Accepted for publication in ApJ. 12 pages, 6 figures, 3 tables. Comments welcome

Implementing Tidal and Gravitational Wave Energy Losses in Few-body Codes: A Fast and Easy Drag Force Model

Authors: Johan Samsing, Nathan W. C. Leigh, Alessandro A. Trani

Abstract: 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 en… ▽ More 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. △ Less

Submitted 22 August, 2018; v1 submitted 21 March, 2018; originally announced March 2018.

Comments: accepted version. MNRAS

Dynamics of tidally captured planets in the Galactic Center

Authors: Alessandro Alberto Trani, Michela Mapelli, Mario Spera, Alessandro Bressan

Abstract: Recent observations suggest ongoing planet formation in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. We investigate the chance of planet tidal captures by running three-body encounters of SMBH-star-planet systems with a high-accuracy regularized code. We s… ▽ More Recent observations suggest ongoing planet formation in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. We investigate the chance of planet tidal captures by running three-body encounters of SMBH-star-planet systems with a high-accuracy regularized code. We show that tidally captured planets have orbits close to those of their parent star. We conclude that the final periapsis distance of the captured planet from the SMBH will be much larger than 200 AU, unless its parent star was already on a highly eccentric orbit. △ Less

Submitted 25 July, 2016; originally announced July 2016.

Comments: 2 pages, 2 figures, to appear in: "Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics and Fundamental Physics", F.R. Ferraro & B. Lanzoni eds, Mem. SAIt, Vol 87

The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center

Authors: Alessandro Alberto Trani, Michela Mapelli, Alessandro Bressan, Federico Inti Pelupessy, Arjen van Elteren, Simon Portegies Zwart

Abstract: The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermo… ▽ More The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk. △ Less

Submitted 8 December, 2015; originally announced December 2015.

Comments: 10 pages, 10 figures, 1 table, accepted for publication in ApJ