Showing 1–50 of 76 results for author: Leigh, N W C

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

New candidate hypervelocity red clump stars in the inner Galactic bulge

Authors: A. Luna, T. Marchetti, M. Rejkuba, N. W. C. Leigh, J. Alonso-García, A. Valenzuela Navarro, D. Minniti, L. C. Smith

Abstract: We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the de… ▽ More We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the detection of HVSs, ejection mechanism models can be constrained by exploring the stellar dynamics in the Galactic centre through a recent stellar interaction with the SMBH. Based on a previously developed methodology by our group, we searched with a sample of preliminary data from version 2 of the Vista Variables in the Via Lactea (VVV) Infrared Astrometric Catalogue (VIRAC2) and Gaia DR3 data, including accurate optical and NIR proper motions. This search resulted in a sample of 46 stars with transverse velocities larger than the local escape velocity within the Galactic bulge, of which 4 are prime candidate HVSs with high-proper motions consistent with being ejections from the Galactic centre. Adding to that, we studied a sample of reddened stars without a Gaia DR3 counterpart and found 481 stars with transverse velocities larger than the local escape velocity, from which 65 stars have proper motions pointing out of the Galactic centre and are candidate HVSs. In total, we found 69 candidate HVSs pointing away from the Galactic centre with transverse velocities larger than the local escape velocity. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. 18 pages (5 of Appendix), 15 figures

A catalogue of Galactic GEMS: Globular cluster Extra-tidal Mock Stars

Authors: Steffani M. Grondin, Jeremy J. Webb, James M. M. Lane, Joshua S. Speagle, Nathan W. C. Leigh

Abstract: This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are p… ▽ More This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample N=50,000 extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way's dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way. △ Less

Submitted 16 January, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 23 pages, 12 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society. For the associated GEMS catalogue, please visit https://zenodo.org/records/8436703

In-plane Tidal Disruption of Stars in Disks of Active Galactic Nuclei

Authors: Taeho Ryu, Barry McKernan, Saavik Ford, Matteo Cantiello, Matthew Graham, Daniel Stern, Nathan W. C Leigh

Abstract: Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk g… ▽ More Stars embedded in active galactic nucleus (AGN) disks or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with {\small AREPO}, we investigate the dependence of debris properties in in-plane TDEs in AGN disks on the disk density and the orientation of stellar orbits relative to the disk gas (pro- and retro-grade). Key findings are: 1) Debris experiences continuous perturbations from the disk gas, which can result in significant and continuous changes in debris energy and angular momentum compared to `naked' TDEs. 2) Above a critical density of a disk around a SMBH with mass $M_{\bullet}$ ($ρ_{\rm crit} \sim 10^{-8}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$) for retrograde stars, both bound and unbound debris is fully mixed into the disk. The density threshold for no bound debris return, inhibiting the accretion component of TDEs, is $ρ_{\rm crit,bound} \sim 10^{-9}{\rm g~cm^{-3}}(M_{\bullet}/10^{6}{\rm M}_{\odot})^{-2.5}$. 3) Observationally, AGN-TDEs transition from resembling naked TDEs in the limit of $ρ_{\rm disk}\lesssim 10^{-2}ρ_{\rm crit,bound}$ to fully muffled TDEs with associated inner disk state changes at $ρ_{\rm disk}\gtrsimρ_{\rm crit,bound}$, with a superposition of AGN+TDE in between. Stellar or remnant passages themselves can significantly perturb the inner disk. This can lead to an immediate X-ray signature and optically detectable inner disk state changes, potentially contributing to the changing-look AGN phenomenon. 4) Debris mixing can enriches the average disk metallicity over time if the star's metallicity exceeds that of the disk gas. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Comments: 16 pages, 12 figures, 2 tables. Submitted to MNRAS. Movies: https://youtube.com/playlist?list=PLxLK3qI02cQfhzVB3JAk7ThAnv7oLMUQL&si=DO6AHjXhkgpBgN30

The dominant mechanism(s) for populating the outskirts of star clusters with neutron star binaries

Authors: Nathan W. C. Leigh, Claire S. Ye, Steffani M. Grondin, Giacomo Fragione, Jeremy J. Webb, Craig O. Heinke

Abstract: It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of or… ▽ More It has been argued that heavy binaries composed of neutron stars (NSs) and millisecond pulsars (MSPs) can end up in the outskirts of star clusters via an interaction with a massive black hole (BH) binary expelling them from the core. We argue here, however, that this mechanism will rarely account for such observed objects. Only for primary masses $\lesssim$ 100 M$_{\odot}$ and a narrow range of orbital separations should a BH-BH binary be both dynamically hard and produce a sufficiently low recoil velocity to retain the NS binary in the cluster. Hence, BH binaries are in general likely to eject NSs from clusters. We explore several alternative mechanisms that would cause NS/MSP binaries to be observed in the outskirts of their host clusters after a Hubble time. The most likely mechanism is a three-body interaction involving the NS/MSP binary and a normal star. We compare to Monte Carlo simulations of cluster evolution for the globular clusters NGC 6752 and 47 Tuc, and show that the models not only confirm that normal three-body interactions involving all stellar-mass objects are the dominant mechanism for putting NS/MSP binaries into the cluster outskirts, they also reproduce the observed NS/MSP binary radial distributions without needing to invoke the presence of a massive BH binary. Higher central densities and an episode of core-collapse can broaden the radial distributions of NSs/MSPs and NS/MSP binaries due to three-body interactions, making these clusters more likely to host NSs in the cluster outskirts. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures, 2 tables, submitted to MNRAS

Globular Cluster UVIT legacy Survey (GlobUleS) III. Omega Centauri in Far-Ultraviolet

Authors: Deepthi S. Prabhu, Annapurni Subramaniam, Snehalata Sahu, Chul Chung, Nathan W. C. Leigh, Emanuele Dalessandro, Sourav Chatterjee, N. Kameswara Rao, Michael Shara, Patrick Cote, Samyaday Choudhury, Gajendra Pandey, Aldo A. R. Valcarce, Gaurav Singh, Joesph E. Postma, Sharmila Rani, Avrajit Bandyopadhyay, Aaron M. Geller, John Hutchings, Thomas Puzia, Mirko Simunovic, Young-Jong Sohn, Sivarani Thirupathi, Ramakant Singh Yadav

Abstract: We present the first comprehensive study of the most massive globular cluster Omega Centauri in the far-ultraviolet (FUV) extending from the center to ~ 28% of the tidal radius using the Ultraviolet Imaging Telescope aboard AstroSat. A comparison of the FUV-optical color-magnitude diagrams with available canonical models reveals that the horizontal branch (HB) stars bluer than the knee (hHBs) and… ▽ More We present the first comprehensive study of the most massive globular cluster Omega Centauri in the far-ultraviolet (FUV) extending from the center to ~ 28% of the tidal radius using the Ultraviolet Imaging Telescope aboard AstroSat. A comparison of the FUV-optical color-magnitude diagrams with available canonical models reveals that the horizontal branch (HB) stars bluer than the knee (hHBs) and the white dwarfs (WDs) are fainter in the FUV by ~ 0.5 mag than model predictions. They are also fainter than their counterparts in M13, another massive cluster. We simulated HB with at least five subpopulations, including three He-rich populations with a substantial He enrichment of Y up to 0.43 dex, to reproduce the observed FUV distribution. We find the He-rich younger subpopulations to be radially more segregated than the He-normal older ones, suggesting an in-situ enrichment from older generations. The Omega Cen hHBs span the same effective temperature range as their M13 counterparts, but some have smaller radii and lower luminosities. This may suggest that a fraction of Omega Cen hHBs are less massive than those of M13, similar to the result derived from earlier spectroscopic studies of outer extreme HB stars. The WDs in Omega Cen and M13 have similar luminosity-radius-effective temperature parameters, and 0.44 - 0.46 M$_\odot$ He-core WD model tracks evolving from progenitors with Y = 0.4 dex are found to fit the majority of these. This study provides constraints on the formation models of Omega Cen based on the estimated range in age, [Fe/H] and Y (in particular), for the HB stars. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted for publication in ApJL; 13 pages, 5 figures, 1 table

Searching for the extra-tidal stars of globular clusters using high-dimensional analysis and a core particle spray code

Authors: Steffani M. Grondin, Jeremy J. Webb, Nathan W. C. Leigh, Joshua S. Speagle, Reem J. Khalifeh

Abstract: Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data ana… ▽ More Three-body interactions can eject stars from the core of a globular cluster, causing them to enter the Galactic halo as extra-tidal stars. While finding extra-tidal stars is imperative for understanding cluster evolution, connecting isolated extra-tidal field stars back to their birth cluster is extremely difficult. In this work, we present a new methodology consisting of high-dimensional data analysis and a particle spray code to identify extra-tidal stars of any Galactic globular cluster using M3 as a case study. Using the t-Stochastic Neighbour Embedding (t-SNE) and Uniform Manifold Approximation and Projection (UMAP) machine learning dimensionality reduction algorithms, we first identify a set of 103 extra-tidal candidates in the APOGEE DR17 data catalogue with chemical abundances similar to M3 stars. To confirm each candidate's extra-tidal nature, we introduce corespray - a new Python-based three-body particle spray code that simulates extra-tidal stars for any Galactic globular cluster. Using Gaia EDR3 proper motions and APOGEE DR17 radial velocities, we apply multivariate Gaussian modelling and an extreme deconvolution to identify the extra-tidal candidates that are more likely to be associated with a distribution of corespray-simulated M3 extra-tidal stars than the field. Through these methods, we identify 10 new high-probability extra-tidal stars produced via three-body interactions in M3. We also explore whether any of our extra-tidal candidates are consistent with being ejected from M3 through different dynamical processes. Future applications of corespray will yield better understandings of core dynamics, star formation histories and binary fractions in globular clusters. △ Less

Submitted 16 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 16 pages, 8 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

Authors: Syeda S. Nasim, Gaia Fabj, Freddy Caban, Amy Secunda, K. E. Saavik Ford, Barry McKernan, Jillian M. Bellovary, Nathan W. C. Leigh, Wladimir Lyra

Abstract: Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population w… ▽ More Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population which drives sBH-sBH mergers detectable in gravitational waves using LIGO-Virgo-KAGRA (LVK) or sBH-SMBH mergers detectable with LISA (Laser Interferometer Space Antenna). Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive 'immortal' stars. Here, we investigate the behaviour of polar and retrograde orbiters $(i \geq 90^{\circ})$ interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation $(i<90^{\circ})$ during capture, and decrease $a$ dramatically towards the SMBH. For sBH, we find a critical angle $i_{\rm ret} \sim 113^{\circ}$, below which retrograde sBH decay towards embedded prograde orbits $(i \to 0^{\circ})$, while for $i_{\rm o}>i_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution. △ Less

Submitted 17 May, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures, 1 table

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, July 2023, Pages 5393-5401

The thermodynamics of stellar multiplicity: dynamical evolution of binary star populations in dense stellar environments

Authors: N. W. C. Leigh, N. C. Stone, J. J. Webb, W. Lyra

Abstract: We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant… ▽ More We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single-binary interactions, that drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step toward the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 12 pages, 8 figures, 2 appendices

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

The mean free path approximation and stellar collisions in star clusters: Numerical exploration of the analytic rates and the role of perturbations on binary star mergers

Authors: Bastián Reinoso, Nathan W. C. Leigh, Carlos M. Barrera-Retamal, Dominik Schleicher, Ralf S. Klessen, Amelia M. Stutz

Abstract: In this paper we compute predictions for the number of stellar collisions derived from analytic models based on the mean free path (MFP) approximation and compare them to the results of $N$-body simulations. Our goal is to identify the cluster conditions under which the MFP approximation remains valid. Adopting a range of particle numbers ($100\leq N\leq5000$) and different combinations of particl… ▽ More In this paper we compute predictions for the number of stellar collisions derived from analytic models based on the mean free path (MFP) approximation and compare them to the results of $N$-body simulations. Our goal is to identify the cluster conditions under which the MFP approximation remains valid. Adopting a range of particle numbers ($100\leq N\leq5000$) and different combinations of particle masses and radii, we explore three different channels leading to stellar collisions, all of which are expected to occur in realistic stellar environments. At high densities, binaries form from isolated three-body interactions of single stars. Hence, we consider collisions between single stars and collisions involving binary stars, after they form in our simulations. For the latter, we consider two channels for mergers, namely direct stellar collisions during chaotic single-binary interactions and perturbation-driven mergers of binaries due to random walks in eccentricity approaching unity. In the densest systems considered here, a very massive object is formed at the cluster centre, causing local stellar orbits to become increasingly Keplerian and the assumptions going into our analytic model to break down. Before reaching this limit, we obtain excellent agreement between our theoretical predictions and the simulations: the analytic rates are typically accurate to within one standard deviation for the entire parameter space considered here, but the agreement is best for short integration times. Our results have direct implications for blue straggler formation in dense star clusters, and stellar mergers in galactic nuclei hosting massive black holes. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 13 pages, 7 figures. Accepted by MNRAS

Starfall: A heavy rain of stars in 'turning on' AGN

Authors: B. McKernan, K. E. S. Ford, M. Cantiello, M. J. Graham, A. S. Jermyn, N. W. C. Leigh, T. Ryu, D. Stern

Abstract: As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and… ▽ More As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and onto the SMBH. A large AGN loss cone ($θ_{\rm AGN,lc}$) can result from binary scatterings in the inner disk and yield tidal disruption events (TDEs). Signatures of starfall include optical/UV flares that rise in luminosity over time, particularly in the inner disk. If the SMBH mass is $M_{\rm SMBH} \ge 10^{8}M_{\odot}$, flares truncate abruptly and the star is swallowed. If $M_{\rm SMBH}<10^{8}M_{\odot}$, and if the infalling orbit lies within $θ_{\rm AGN,lc}$, the flare is followed by a TDE which can be prograde or retrograde relative to the AGN inner disk. Retrograde AGN TDEs are over-luminous and short-lived as in-plane ejecta collide with the inner disk and a lower AGN state follows. Prograde AGN TDEs add angular momentum to inner disk gas and so start off looking like regular TDEs but are followed by an AGN high state. Searches for such flare signatures test models of AGN `turn on', SMBH mass, as well as disk properties and the embedded population. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 10 pages, 1 figure, MNRAS submitted

Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei

Authors: Yihan Wang, Barry McKernan, Saavik Ford, Rosalba Perna, Nathan W. C. Leigh, Mordecai-Mark Mac Low

Abstract: Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-maj… ▽ More Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: submitted to ApJL. Comments are welcome

Hot Jupiter formation in dense clusters: secular chaos in multi-planetary systems

Authors: Yihan Wang, Rosalba Perna, Nathan W. C. Leigh, Michael M. Shara

Abstract: Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity ex… ▽ More Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close fly-by stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys in dense clusters have a significant impact on hot Jupiter formation via activation of high eccentricity excitation mechanisms such as the Zeipel-Lidov-Kozai (ZLK) effect and planet-planet scattering. Previous studies have shown that if there are two giant planets in the planetary system, close flybys can efficiently activate the ZLK mechanism, thus triggering high eccentricity tidal migration and ultimately form hot Jupiters in star clusters. Here we extend our previous study with a multi-planet (triple) system. We perform high precision, high-accuracy few-body simulations of stellar flybys and subsequent planetary migration within the perturbed planetary systems using the code {\tt SpaceHub}. Our simulations demonstrate that a single close flyby on a multi-planet system in a cluster can activate secular chaos and ultimately lead to hot Jupiter formation via high eccentricity migration. We find that the hot Jupiter formation rate per system increases with both the size of the planetary system as well as with the mass of the outer planet, and we quantify the relative formation fractions for a range of parameters. Hot Jupiters formed via secular chaos are expected to be accompanied by massive companions with very long periods. Our study further shows that this flyby-induced secular chaos is preferred in low-density clusters where multi-planet systems are more likely to survive, and that it contributes a significant fraction of the hot Jupiter formation in star clusters compared to the flyby-induced ZLK mechanism. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: submitted to MNRAS

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

Stellar collisions in flattened and rotating Pop. III star clusters

Authors: M. Z. C. Vergara, D. R. G. Schleicher, T. C. N. Boekholt, B. Reinoso, M. Fellhauer, R. S. Klessen, N. W. C. Leigh

Abstract: Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive c… ▽ More Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive cluster, in which the most massive star moves to the center of the cluster due to dynamical friction. This increases the rate of collisions and mergers since massive stars have larger collisional cross sections. This can lead to runaway growth of a very massive star which may collapse to become an intermediate-mass black hole. Here we investigate the dynamical evolution of Miyamoto-Nagai models that allow us to describe dense stellar clusters, including flattening and different degrees of rotation. We find that the collisions in these clusters depend mostly on the number of stars and the initial stellar radii for a given radial size of the cluster. By comparison, rotation seems to affect the collision rate by at most $20\%$. For flatness, we compared spherical models with systems that have a scale height of about $10\%$ of their radial extent, in this case finding a change in the collision rate of less than $25\%$. Overall, we conclude that the parameters only have a minor effect on the number of collisions. Our results also suggest that rotation helps to retain more stars in the system, reducing the number of escapers by a factor of $2-3$ depending on the model and the specific realization. After two million years, a typical lifetime of a very massive star, we find that about $630$ collisions occur in typical models with $N=10^4$, $R=100$ $\rm~R_\odot$ and a half-mass radius of $0.1$ $\rm~pc$, leading to a mass of about $6.3\times10^3$ $\rm~M_\odot$ for the most massive object. △ Less

Submitted 3 April, 2021; originally announced April 2021.

Comments: 10 pages, 7 figures

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

Formation of supermassive black hole seeds in nuclear star clusters via gas accretion and runaway collisions

Authors: Arpan Das, Dominik R. G. Schleicher, Nathan W. C. Leigh, Tjarda C. N. Boekholt

Abstract: More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. I… ▽ More More than two hundred supermassive black holes (SMBHs) of masses $\gtrsim 10^9\,\mathrm{M_{\odot}}$ have been discovered at $z \gtrsim 6$. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ into seed black holes which could grow upto few times $10^9\,\mathrm{M_{\odot}}$ SMBHs observed at $z\sim 7$. In this paper, we explore how SMSs with masses $\sim 10^{3-5}\,\mathrm{M_{\odot}}$ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealised N-body simulations. We explore physically motivated accretion scenarios, e.g. Bondi-Hoyle-Lyttleton accretion and Eddington accretion, as well as simplified scenarios such as constant accretions. While gas is present, the accretion timescale remains considerably shorter than the timescale for collisions with the most massive object (MMO). However, overall the timescale for collisions between any two stars in the cluster can become comparable or shorter than the accretion timescale, hence collisions still play a crucial role in determining the final mass of the SMSs. We find that the problem is highly sensitive to the initial conditions and our assumed recipe for the accretion, due to the highly chaotic nature of the problem. The key variables that determine the mass growth mechanism are the mass of the MMO and the gas reservoir that is available for the accretion. Depending on different conditions, SMSs of masses $\sim10^{3-5} \,\mathrm{M_{\odot}}$ can form for all three accretion scenarios considered in this work. △ Less

Submitted 14 May, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

A phylogenetic analysis of galaxies in the Coma Cluster and the field: a new approach to galaxy evolution

Authors: M. Martínez-Marín, R. Demarco, G. Cabrera-Vives, P. Cerulo, N. W. C. Leigh, R. Herrera-Camus

Abstract: We propose a phylogenetic approach (PA) as a novel and robust tool to detect galaxy populations (GPs) based on their chemical composition. The branches of the tree are interpreted as different GPs and the length between nodes as the internal chemical variation along a branch. We apply the PA using 30 abundance indices from the Sloan Digital Sky Survey to 475 galaxies in the Coma Cluster and 438 ga… ▽ More We propose a phylogenetic approach (PA) as a novel and robust tool to detect galaxy populations (GPs) based on their chemical composition. The branches of the tree are interpreted as different GPs and the length between nodes as the internal chemical variation along a branch. We apply the PA using 30 abundance indices from the Sloan Digital Sky Survey to 475 galaxies in the Coma Cluster and 438 galaxies in the field. We find that a dense environment, such as Coma, shows several GPs, which indicates that the environment is promoting galaxy evolution. Each population shares common properties that can be identified in colour magnitude space, in addition to minor structures inside the red sequence. The field is more homogeneous, presenting one main GP. We also apply a principal component analysis (PCA) to both samples, and find that the PCA does not have the same power in identifying GPs. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, 2020, Issue 4, pp.5607-5622

Hot Jupiter and ultra-cold Saturn formation in dense star clusters

Authors: Yi-Han Wang, Nathan W. C. Leigh, Rosalba Perna, Michael M. Shara

Abstract: The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequen… ▽ More The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJS

Small-N Collisional Dynamics V: Beyond the Realm of Not-So-Small-N

Authors: Carlos Barrera, Nathan W. C. Leigh, Bastián Reinoso, Amelia M. Stutz, Dominik Schleicher

Abstract: Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important application is stellar collisions, which occur commonly in dense star clusters, and their relevance for the formation of various types of stellar exotica. In this paper,… ▽ More Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important application is stellar collisions, which occur commonly in dense star clusters, and their relevance for the formation of various types of stellar exotica. In this paper, we return to our study of the collision rates and probabilities during small-number chaotic gravitational interactions ($N$ $\lesssim$ 10), moving beyond the small-number particle limit and into the realm of larger particle numbers ($N$ $\gtrsim$ 10$^3$) to test the extent of validity of our analytic model as a function of the particle properties and the number of interacting particles. This is done using direct $N$-body simulations of stellar collisions in dense star clusters, by varying the relative numbers of particles with different particle masses and radii. We compute the predicted rate of collisions using the mean free path approximation, adopting the point-particle limit and using the sticky-star approximation as our collision criterion. We evaluate its efficacy in the regime where gravitational-focusing is important by comparing the theoretical rates to numerical simulations. Using the tools developed in previous papers in this series, in particular Collision Rate Diagrams, we illustrate that our predicted and simulated rates are in excellent agreement, typically consistent with each other to within one standard deviation. △ Less

Submitted 30 October, 2020; originally announced November 2020.

Comments: 11 pages, 6 figures and 3 tables; submitted for publication; comments welcome

Massive Stellar Triples Leading to Sequential Binary Black-Hole Mergers in the Field

Authors: Alejandro Vigna-Gómez, Silvia Toonen, Enrico Ramirez-Ruiz, Nathan W. C. Leigh, Jeff Riley, Carl-Johan Haster

Abstract: Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are… ▽ More Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well. △ Less

Submitted 22 December, 2020; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 13 main pages, 4 figures, 1 table plus Appendices. Accepted for publication in ApJ Letters

Filament rotation in the California L1482 cloud

Authors: Rodrigo H. Álvarez-Gutiérrez, Amelia M. Stutz, Chi Y. Law, Stefan Reissl, Ralf S. Klessen, Nathan W. C. Leigh, Hong-Li Liu, Rodrigo A. Reeves

Abstract: We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C… ▽ More We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation. △ Less

Submitted 17 December, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 23 pages, 17 figures. Accepted for publication in the Astrophysical Journal

AGN as Potential Factories for Eccentric Black Hole Mergers

Authors: J. Samsing, I. Bartos, D. J. D'Orazio, Z. Haiman, B. Kocsis, N. W. C. Leigh, B. Liu, M. E. Pessah, H. Tagawa

Abstract: There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to th… ▽ More There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to the merger, as basic physics would argue for circularization. Here we show that AGN-disk environments can lead to an excess of eccentric mergers, if the interactions between single and binary black holes are frequent, and occur with mutual inclinations of less than a few degrees. We further illustrate that this eccentric population has a different distribution of the inclination between the spin vectors of the black holes and their orbital angular momentum at merger, referred to as the spin-orbit tilt, compared to the remaining circular mergers. △ Less

Submitted 21 April, 2022; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted

Journal ref: Nature, Volume 603, Issue 7900, p.237-240, 2022

Evolution of Retrograde Orbiters in an AGN Disk

Authors: Amy Secunda, Betsy Hernandez, Jeremy Goodman, Nathan W. C. Leigh, Barry McKernan, K. E. Saavik Ford, Jose I. Adorno

Abstract: AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifyin… ▽ More AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifying the retrograde contribution is important, since roughly half of all orbiters should initially be on retrograde orbits when the disk forms. We perform an analytic calculation of the evolution of ROs in an AGN disk. Because this evolution could cause the orbits of ROs to cross those of prograde BBHs, we derive the collision rate between a given RO and a given BBH orbiting in the prograde direction. In the examples given here, ROs in the inner region of the disk experience a rapid decrease in the semimajor axis of their orbits while also becoming highly eccentric in less than a million years. This rapid orbital evolution could lead to extreme mass ratio inspirals detectable by the Laser Interferometer Space Antenna. The collision rates of our example ROs with prograde BBHs in the migration trap depend strongly on the volume of the inner radiation-pressure-dominated region which depends on the mass of the supermassive black hole (SMBH). Rates are lowest for larger mass SMBHs, which dominate the AGN merger channel, suggesting that merger rates for this channel may not be significantly altered by ROs. △ Less

Submitted 4 February, 2021; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: 12 pages, 3 figures, accepted to ApJL

The effects of a background potential in star cluster evolution: a delay in the relaxation time-scale and runaway collision processes

Authors: B. Reinoso, D. R. G. Schleicher, M. Fellhauer, N. W. C. Leigh, R. S. Klessen

Abstract: Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explic… ▽ More Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: 16 pages. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 639, A92 (2020)

Mergers of Equal-Mass Binaries with Compact Object Companions from Mass Transfer in Triple Star Systems

Authors: N. W. C. Leigh, S. Toonen, S. F. Portegies Zwart, R. Perna

Abstract: In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analyti… ▽ More In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion. △ Less

Submitted 23 June, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: 15 pages, 7 figures; accepted for publication in MNRAS

Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes II. Spins and Incoming Objects

Authors: Amy Secunda, Jillian Bellovary, Mordecai-Mark Mac Low, K. E. Saavik Ford, Barry McKernan, Nathan W. C. Leigh, Wladimir Lyra, Zsolt Sandor, Jose I. Adorno

Abstract: The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previou… ▽ More The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. Orbiters entering the inner disk form BBHs with sBHs on resonant orbits near the migration trap. When these sBHs reach ~80 Msun, they form BBHs with sBHs in the migration trap, which over 10 Myr reach ~1000 Msun. We find 68% of the BBHs in our simulation orbit in the retrograde direction, which implies BBHs in our merger channel will have small dimensionless aligned spins, χ_eff. Overall, our models produce BBHs that resemble both the majority of BBH mergers detected thus far (0.66 to 120 Gpc^-3 yr^-1 ) and two recent unusual detections, GW190412 (~0.3 Gpc^-3 yr^-1 ) and GW190521 (~0.1 Gpc^-3 yr^-1 ). △ Less

Submitted 25 September, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures, accepted to ApJ

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

Giant planet swaps during close stellar encounters

Authors: Yi-Han Wang, Rosalba Perna, Nathan W. C. Leigh

Abstract: The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not… ▽ More The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $Γ\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (σ/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $σ$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ Letters

Planetary Architectures in Interacting Stellar Environments

Authors: Yi-Han Wang, Rosalba Perna, Nathan W. C. Leigh

Abstract: The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers… ▽ More The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers (stars and binaries), we explore the reconfiguration to the planetary system due to the encounters. In particular, via an exploration focused on the strong scattering regime, we derive the fraction of encounters which result in planet ejections, planet transfers and collisions by the interloper star/binary, as a function of the characteristics of the environment (density, velocity dispersion), and for different masses of the flyby star/binary. We find that binary interlopers can significantly increase the cross section of planet ejections and collisions, while they only slightly change the cross section for planet transfers. Therefore, in environments with high binary fractions, floating planets are expected to be relatively common, while in environments with low binary fractions, where the cross sections of planet ejection and transfer are comparable, the rate of planet exchanges between two stars will be comparable to the rate of production of free-floating planets. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: Submitted to MNRAS. Comments are welcome

The X-ray Emissivity of Low-Density Stellar Populations

Authors: C. O. Heinke, M. G. Ivanov, E. W. Koch, R. Andrews, L. Chomiuk, H. N. Cohn, S. Crothers, T. de Boer, N. Ivanova, A. K. H. Kong, N. W. C. Leigh, P. M. Lugger, L. Nelson, C. J. Parr, E. W. Rosolowsky, A. J. Ruiter, C. L. Sarazin, A. W. Shaw, G. R. Sivakoff, M. van den Berg

Abstract: The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations,… ▽ More The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: 28 pages, 23 figures. MNRAS, in press

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

Authors: Matteo Bonetti, Alexander Rasskazov, Alberto Sesana, Massimo Dotti, Francesco Haardt, Nathan W. C. Leigh, Manuel Arca Sedda, Giacomo Fragione, Elena Rossi

Abstract: We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrap… ▽ More We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$). △ Less

Submitted 4 February, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 6 pages, 5 figures; accepted for publication in MNRAS

A Statistical Solution to the Chaotic, Non-Hierarchical Three-Body Problem

Authors: Nicholas C. Stone, Nathan W. C. Leigh

Abstract: The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost a… ▽ More The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost always disintegrate into a single escaping star and a stable, bound binary, but the chaotic nature of the three-body problem prevents the derivation of tractable analytic formulae deterministically mapping initial conditions to final outcomes. However, chaos also motivates the assumption of ergodicity, suggesting that the distribution of outcomes is uniform across the accessible phase volume. Here, we use the ergodic hypothesis to derive a complete statistical solution to the non-hierarchical three-body problem, one which provides closed-form distributions of outcomes (e.g. binary orbital elements) given the conserved integrals of motion. We compare our outcome distributions to large ensembles of numerical three-body integrations, and find good agreement, so long as we restrict ourselves to "resonant" encounters (the ~50% of scatterings that undergo chaotic evolution). In analyzing our scattering experiments, we identify "scrambles" (periods in time where no pairwise binaries exist) as the key dynamical state that ergodicizes a non-hierarchical triple. The generally super-thermal distributions of survivor binary eccentricity that we predict have notable applications to many astrophysical scenarios. For example, non-hierarchical triples produced dynamically in globular clusters are a primary formation channel for black hole mergers, but the rates and properties of the resulting gravitational waves depend on the distribution of post-disintegration eccentricities. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: 54 pages, 8 figures

Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

Authors: Giacomo Fragione, Brian D. Metzger, Rosalba Perna, Nathan W. C. Leigh, Bence Kocsis

Abstract: Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effect… ▽ More Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies. △ Less

Submitted 27 April, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 13 pages, 8 figures, 1 table, accepted by MNRAS

Kinematic study of the association Cyg OB3 with Gaia DR2

Authors: Anjali Rao, Poshak Gandhi, Christian Knigge, John A. Paice, Nathan W. C. Leigh, Douglas Boubert

Abstract: We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean pec… ▽ More We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is ~22 km/s, dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds <20 km/s with respect to the association for a majority of the sample stars. △ Less

Submitted 1 May, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: Accepted for publication 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

Constraining Mass Transfer Histories of Blue Straggler Stars with COS Spectroscopy of White Dwarf Companions

Authors: Natalie M. Gosnell, Emily M. Leiner, Robert D. Mathieu, Aaron M. Geller, Christian Knigge, Alison Sills, Nathan W. C. Leigh

Abstract: Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with… ▽ More Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with HST COS far-ultraviolet spectroscopy. Adopting a Gaia-based cluster distance of $1847\pm107$ pc, we determine that one system, WOCS 4540, formed through Case C mass transfer resulting in a CO-core white dwarf with $T_{\text{eff}}=17000^{+140}_{-200}$ K and a $\log g=7.80^{+0.06}_{-0.06}$, corresponding to a mass of $0.53^{+0.03}_{-0.03}$ $M_{\odot}$ and a cooling age of $105^{+6}_{-5}$ Myr. The other system, WOCS 5379, formed through Case B mass transfer resulting in a He-core white dwarf with $T_{\text{eff}}=15500^{+170}_{-150}$ K and a $\log g=7.50^{+0.06}_{-0.05}$, corresponding to a mass of $0.42^{+0.02}_{-0.02}$ $M_{\odot}$ and an age of $250^{+20}_{-20}$ Myr. The WD parameters are consistent across four different cluster distance assumptions. We determine possible progenitor binary systems with a grid of accretion models using MESA, and investigate whether these systems would lead to stable or unstable mass transfer. WOCS 4540 likely resulted from stable mass transfer during periastron passage in an eccentric binary system, while WOCS 5379 challenges our current understanding of the expected outcomes for mass transfer from red giant branch stars. Both systems are examples of the value in using detailed analyses to fine-tune our physical understanding of binary evolutionary processes. △ Less

Submitted 19 August, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 16 pages, 10 figures, 2 tables, submitted to Astrophysical Journal

In Search of the Thermal Eccentricity Distribution

Authors: Aaron M. Geller, Nathan W. C. Leigh, Mirek Giersz, Kyle Kremer, Frederic A. Rasio

Abstract: About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary… ▽ More About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two. △ Less

Submitted 31 January, 2019; originally announced February 2019.

Comments: 14 pages, 7 figures, accepted for publication in The Astrophysical Journal

A Triple Origin for Twin Blue Stragglers in Close Binaries

Authors: Simon Portegies Zwart, Nathan W. C. Leigh

Abstract: We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed… ▽ More We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed model. We predict the following properties for twin BSs: (1) For the outer tertiary orbit, the initial orbital period should lie between 220 days $\lesssim$ P$_{\rm out}$ $\lesssim$ 1100 days, assuming initial masses for the inner binary components of $m_{\rm 1} = 1.1$ M$_{\odot}$ and $m_{\rm 2} =$ 0.9 M$_{\odot}$ and an outer tertiary mass of $m_{\rm 3} = 1.4$ M$_{\odot}$. After Roche-lobe overflow, the outer star turns into a white dwarf (WD) of mass 0.43 to 0.54\,\MSun. There is a correlation between the mass of this WD and the outer orbital period: more massive WDs will be on wider orbits. (3) The rotational axes of both BSs will be aligned with each other and the orbital plane of the outer tertiary WD. (4) The BSs will have roughly equal masses, independent of their initial masses (since the lower mass star accretes the most). The dominant accretor should, therefore, be more enriched by the accreted material. Hence, one of the BSs will appear to be more enriched by either He, C and O or by s-process elements, if the donor started Roche lobe overflow on, respectively, the red giant or asymptotic giant branch. (5) Relative to old clusters, twin BSs in close binaries formed from the proposed mechanism should be more frequent in the Galactic field and open clusters with ages $\lesssim$ 4-6 Gyr, since then the donor will have a radiative envelope. (6) The orbit of the binary BS will have a small semi-major axis (typically $\aplt 0.3$\,au) and be close to circular ($e \aplt 0.2$). △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: 12 pages, 6 figures; submitted to ApJ

Black hole and neutron star mergers in galactic nuclei

Authors: Giacomo Fragione, Evgeni Grishin, Nathan W. C. Leigh, Hagai. B. Perets, Rosalba Perna

Abstract: Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies,… ▽ More Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes. △ Less

Submitted 12 June, 2019; v1 submitted 26 November, 2018; originally announced November 2018.

Comments: 19 pages, 8 figures, 5 Tables, accepted by MNRAS

Hypervelocity Stars from a Supermassive Black Hole-Intermediate-mass Black Hole binary

Authors: Alexander Rasskazov, Giacomo Fragione, Nathan W. C. Leigh, Hiromichi Tagawa, Alberto Sesana, Adrian Price-Whelan, Elena Maria Rossi

Abstract: In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various… ▽ More In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.). △ Less

Submitted 23 May, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 13 pages, 5 figures, 3 tables, accepted for publication at ApJ

Journal ref: The Astrophysical Journal, Volume 878, Issue 1, article id. 17, 12 pp. (2019)

Constraints on blue straggler formation mechanisms in Galactic globular clusters from proper motion velocity distributions

Authors: N. W. C. Leigh, T. Panurach, M. Simunovic, A. M. Geller, D. Zurek, M. M. Shara, A. Sills, C. Knigge, N. M. Gosnell, R. Mathieu, T. H. Puzia, J. Ventura, Q. Minor

Abstract: For a sample of 38 Galactic globular clusters (GCs), we confront the observed distributions of blue straggler (BS) proper motions and masses (derived from isochrone fitting) from the BS catalog of Simunovic & Puzia with theoretical predictions for each of the two main competing BS formation mechanisms. These are mass transfer from an evolved donor on to a main-sequence (MS) star in a close binary… ▽ More For a sample of 38 Galactic globular clusters (GCs), we confront the observed distributions of blue straggler (BS) proper motions and masses (derived from isochrone fitting) from the BS catalog of Simunovic & Puzia with theoretical predictions for each of the two main competing BS formation mechanisms. These are mass transfer from an evolved donor on to a main-sequence (MS) star in a close binary system, and direct collisions involving MS stars during binary encounters. We use the \texttt{FEWBODY} code to perform simulations of single-binary and binary-binary interactions. This provides collisional velocity and mass distributions for comparison to the observed distributions. Most clusters are consistent with BSs derived from a dynamically relaxed population, supportive of the binary mass-transfer scenario. In a few clusters, including all the post-core collapse clusters in our sample, the collisional velocities provide the best fit. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 10 pages, 8 figures, 1 table; accepted for publication in MNRAS

Small-N collisional dynamics IV: Order in the realm of not-so-small-N

Authors: Nathan W. C. Leigh, Aaron M. Geller, Michael M. Shara, Lukas Baugher, Vianny Hierro, De'Andre Ferreira, Elizabeth Teperino

Abstract: In this paper, the fourth in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We focus once again on the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and interactions that produce direct collisions between any two particles. Our long-term goal is to predict the probability of a given collision event occurring over the course of an… ▽ More In this paper, the fourth in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We focus once again on the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and interactions that produce direct collisions between any two particles. Our long-term goal is to predict the probability of a given collision event occurring over the course of an interaction, as a function of the numbers and properties of the particles. In previous papers, we varied the number of interacting particles, as well as the distributions of particle radii and masses. Here, we refine the methods developed in these preceding studies, and arrive at a final and robust methodology that can be used to study collisional dynamics in a variety of astrophysical contexts, ranging from stars in star clusters, galaxies in galaxy groups and clusters and even the collisional growth of planetesimals in protoplanetary disks. We further present and refine the concept of a Collision Rate Diagram (CRD), the primary tool we use to quantify the relative rates for different collision scenarios to occur. The agreement between our final theoretical model and the results of numerical scattering simulations is excellent. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: 10 pages, 5 figures, 1 table; accepted for publication in MNRAS

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

When do star clusters become multiple star systems? II. Toward a half-life formalism with four bodies

Authors: Timur Ibragimov, Nathan W. C. Leigh, Taeho Ryu, Teresa Panurach, Rosalba Perna

Abstract: We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-… ▽ More We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-life for the distribution of disruption times. As in radioactive decay, the half-lives should provide a direct prediction for the relative fractions of each decay product. We test this prediction with simulated data and find good agreement with our hypothesis. We briefly discuss applications of this feature of the gravitational four-body problem to populations of black holes in globular clusters. This paper, the second in the series, builds on extending the remarkable similarity between gravitational chaos at the macroscopic scale and radioactive decay at the microscopic scale to larger-N systems. △ Less

Submitted 14 March, 2018; originally announced March 2018.

Comments: 17 pages, 18 figures, 2 tables; accepted for publication in MNRAS

Illustrating chaos: A schematic discretization of the general three-body problem in Newtonian gravity

Authors: Nathan W. C. Leigh, Shalma Wegsman

Abstract: We present a formalism for constructing schematic diagrams to depict chaotic three-body interactions in Newtonian gravity. This is done by decomposing each interaction in to a series of discrete transformations in energy- and angular momentum-space. Each time a transformation is applied, the system changes state as the particles re-distribute their energy and angular momenta. These diagrams have t… ▽ More We present a formalism for constructing schematic diagrams to depict chaotic three-body interactions in Newtonian gravity. This is done by decomposing each interaction in to a series of discrete transformations in energy- and angular momentum-space. Each time a transformation is applied, the system changes state as the particles re-distribute their energy and angular momenta. These diagrams have the virtue of containing all of the quantitative information needed to fully characterize most bound or unbound interactions through time and space, including the total duration of the interaction, the initial and final stable states in addition to every intervening temporary meta-stable state. As shown via an illustrative example for the bound case, prolonged excursions of one of the particles, which by far dominates the computational cost of the simulations, are reduced to a single discrete transformation in energy- and angular momentum-space, thereby potentially mitigating any computational expense. We further generalize our formalism to sequences of (unbound) three-body interactions, as occur in dense stellar environments during binary hardening. Finally, we provide a method for dynamically evolving entire populations of binaries via three-body scattering interactions, using a purely analytic formalism. In principle, the techniques presented here are adaptable to other three-body problems that conserve energy and angular momentum. △ Less

Submitted 22 January, 2018; originally announced January 2018.

Comments: 9 pages, 4 figures; accepted for publication in MNRAS

On the rate of black hole binary mergers in galactic nuclei due to dynamical hardening

Authors: Nathan W. C. Leigh, Aaron M. Geller, B. McKernan, K. E. S. Ford, M. -M. Mac Low, J. Bellovary, Z. Haiman, W. Lyra, J. Samsing, M. O'Dowd, B. Kocsis, S. Endlich

Abstract: We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to ha… ▽ More We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS