Search

Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $ii_{\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., 9 pages, 7 figures, 1 table

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., Comment: 16 pages, 8 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Accepted for publication in ApJL; 13 pages, 5 figures, 1 table

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$ 103) 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 1 standard deviation.

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., Comment: 15 pages, 7 figures; accepted for publication in MNRAS

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, Accepted for publication in MNRAS

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

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., 13 pages, 7 figures. Accepted by MNRAS

More than 200 supermassive black holes (SMBHs) of masses $\gtrsim 10^9\, \mathrm{M_{\odot }}$ have been discovered at z ≳ 6. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses $\sim 10^{3}{-}10^{5} \, \mathrm{M_{\odot }}$ into seed black holes which could grow upto few times $10^9\, \mathrm{M_{\odot }}$ SMBHs observed at z ∼ 7. In this paper, we explore how SMSs with masses $\sim 10^{3}{-}10^{5} \, \mathrm{M_{\odot }}$ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealized 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 time-scale remains considerably shorter than the time-scale for collisions with the most massive object (MMO). However, overall the time-scale for collisions between any two stars in the cluster can become comparable or shorter than the accretion time-scale, 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 $\sim 10^{3}{-}10^{5} \, \mathrm{M_{\odot }}$ can form for all three accretion scenarios considered in this work.

We investigate the old open cluster M67 using ultraviolet photometric data of Ultra-Violet Imaging Telescope in multi-filter far-UV bands. M67, well known for the presence of several blue straggler stars (BSS), has been put to detailed tests to understand their formation pathways. Currently, there are three accepted formation channels: mass transfer due to Roche-lobe overflow in binary systems, stellar mergers either due to dynamical collisions or through coalescence of close binaries. So far, there had not been any confirmed detection of a white dwarf (WD) companion to any of the BSSs in this cluster. Here, we present the detection of WD companions to 5 bright BSSs in M67. The multiwavelength spectral energy distributions covering 0.12 -11.5 $\mu$m range, were found to require binary spectral fits for 5 BSSs, consisting of a cool (BSS) and a hot companion. The parameters (Luminosity, Temperature, Radius and Mass) of the hot companions suggest them to be WDs with mass in the range 0.2 - 0.35 M$_{\odot}$ with T$_{eff}$ $\sim$ 11000 - 24000 K., Comment: 4 pages, 2 figures, to appear in IAU Symposium No. 351 proceedings - Star Clusters: From the Milky Way to the Early Universe, 2019

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., 10 pages, 7 figures

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

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., Comment: submitted to MNRAS

We present the open source few-body gravity integration toolkit {\tt SpaceHub}. {\tt SpaceHub} offers a variety of algorithmic methods, including the unique algorithms AR-Radau, AR-Sym6, AR-ABITS and AR-chain$^+$ which we show out-perform other methods in the literature and allow for fast, precise and accurate computations to deal with few-body problems ranging from interacting black holes to planetary dynamics. We show that AR-Sym6 and AR-chain$^+$, with algorithmic regularization, chain algorithm, active round-off error compensation and a symplectic kernel implementation, are the fastest and most accurate algorithms to treat black hole dynamics with extreme mass ratios, extreme eccentricities and very close encounters. AR-Radau, the first regularized Radau integrator with round off error control down to 64 bits floating point machine precision, has the ability to handle extremely eccentric orbits and close approaches in long-term integrations. AR-ABITS, a bit efficient arbitrary precision method, achieves any precision with the least CPU cost compared to other open source arbitrary precision few-body codes. With the implementation of deep numerical and code optimization, these new algorithms in {\tt SpaceHub} prove superior to other popular high precision few-body codes in terms of performance, accuracy and speed., Comment: Submitted to MNRAS. Comments are welcome

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.

Black hole (BH) mergers detectable with the Laser Interferometer Gravitational-wave Observatory (LIGO) can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum, and the spin spectrum of BHs in AGN disks. The predicted merger rate spans ~10−3–104 Gpc−1 yr−1, so upper limits from LIGO (

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., 16 pages. Accepted for publication in Astronomy and Astrophysics

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

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 $\Gamma\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (\sigma/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $\sigma$ 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., Comment: Accepted to ApJ Letters

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., Submitted to MNRAS. Comments are welcome

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, 6 pages, 5 figures; accepted for publication in MNRAS

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., 17 pages, 18 figures, 2 tables; accepted for publication in MNRAS

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 $\chi_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $\chi_{\rm eff}$ BBH mergers are most likely for massive binaries., Comment: submitted to ApJL. Comments are welcome

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., 28 pages, 23 figures. MNRAS, in press

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., Comment: 13 main pages, 4 figures, 1 table plus Appendices. Accepted for publication in ApJ Letters

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, \chi_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 )., Comment: 16 pages, 9 figures, accepted to ApJ

The old open cluster M67, populated with blue straggler stars (BSSs), is a well known test bed to study the BSS formation pathways. Here, we report the first direct detection of a white dwarf (WD) companion to a BSS in M67, using far-UV images from the \textit{Ultra Violet Imaging telescope} ({\it UVIT}) on {\it ASTROSAT}. Near-simultaneous observations in three far-UV bands combined with {\it GALEX}, {\it IUE}, ground and space based photometric data covering 0.14 -11.5 $\mu$m range for WOCS1007 were found to require a binary fit to its spectral energy distribution (SED), consisting of a BSS and a hot companion. On the other hand, a single spectral fit was found to be satisfactory for the SEDs of two other BSSs, WOCS1006 and WOCS2011, with the latter showing a deficient far-UV flux. The hot companion of WOCS1007 is found to have a T$_{eff}$ $\sim$ 13250-13750K and a radius of 0.09$\pm$0.01 R$_\odot$. A comparison with WD models suggests it to be a low mass WD ($\sim$ 0.18$M_\odot$), in agreement with the kinematic mass from the literature. As a low mass WD ($, Comment: 11 pages, 4 figures, Accepted for publication in Astrophysical Journal

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., 13 pages, 11 figures, Accepted for publication in MNRAS

We combine sophisticated high-precision scattering experiments, together with results from the Millenium-II simulation, to compute the cosmic merger rate of bound compact object (CO) binaries dynamically interacting with supermassive black hole binaries (SMBHBs). We consider binaries composed of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). The overall merger rates for WD–WD, NS–NS, BH–BH, BH–NS binaries, and EBBH (eccentric binaries of BHs) from redshift ∼5 are found to be $4.32\times 10^3\, \mathrm{yr}^{-1}(5.93\times 10^2\, \mathrm{yr}^{-1}$ for Type Ia SNe), $82.7\, \mathrm{yr}^{-1}$, $96.3\, \mathrm{yr}^{-1}$, $13.1\, \mathrm{yr}^{-1}$, and $148\, \mathrm{yr}^{-1}$, respectively, for a nominal CO binary fraction in the Galactic centre of 0.1. We calculate the distance (R) distribution of the merger sites with respect to the host galaxies of the binaries. The distribution shows a wide range of distances up to ∼Mpc; this tail is produced by escaped hypervelocity CO binaries. Due to the differences in the matter density of the surrounding environment, merger events with different R are expected to display significantly different signatures in their EM counterparts. In particular, merger events (and especially NS–NS) producing a relativistic jet but occurring in the intergalactic medium will have very weak afterglow radiation relative to their prompt emission. These events, which we call ‘off-centre’, can only be produced from a close encounter between CO binaries and SMBHBs; hence the detection of such merger events would indicate the existence of nearby SMBHBs, and in particular with high mass ratio, produced in the aftermath of a major galaxy merger.

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., 16 pages, 10 figures, 2 tables, submitted to Astrophysical Journal

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., Comment: 23 pages, 17 figures. Accepted for publication in the Astrophysical Journal

Active galactic nucleus (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, as 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.

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}({\sigma}/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $\sigma$ 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., Comment: Accepted for publication in ApJS

We report the discovery of a hot white dwarf (WD) companion to a blue straggler star (BSS) in the globular cluster (GC) NGC 5466, based on observations from the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat. The Spectral Energy Distribution (SED) of the Far-UV detected BSS NH 84 was constructed by combining the flux measurements from 4 filters of UVIT, with GALEX, GAIA and other ground-based observations. The SED of NH 84 reveals the presence of a hot companion to the BSS. The temperature and radius of the BSS (T$_{\mathrm{eff}} = 8000^{+1000}_{-250}$ K, R/R$_\odot = 1.44 \pm 0.05$) derived from Gemini spectra and SED fitting using Kurucz atmospheric models are consistent with each other. The temperature and radius of the hotter companion of NH 84 (T$_{\mathrm{eff}} = 32,000 \pm 2000$ K, R/R$_\odot = 0.021 \pm 0.007$) derived by fitting Koester WD models to the SED suggest that it is likely to be a hot WD. The radial velocity derived from the spectra along with the proper motion from GAIA DR2 confirms NH 84 to be a kinematic member of the cluster. This is the second detection of a BSS-WD candidate in a GC, and the first in the outskirts of a low density GC. The location of this BSS in NGC 5466 along with its dynamical age supports the mass-transfer pathway for BSS formation in low density environments., Comment: 12 pages, 20 figures. accepted for publication in ApJ

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., Comment: 54 pages, 8 figures

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., Comment: 13 pages, 8 figures, 1 table, accepted by MNRAS

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

We present a novel method to constrain the past collisional evolution of observed globular cluster (GC) systems, in particular their mass functions. We apply our method to a pair of galaxies hypothesized to have recently undergone an episode of violent relaxation due to a strong galaxy-galaxy interaction, namely NGC 1052-DF2 and NGC 1052-DF4. We begin by exploring the observational evidence for a collisional origin for these two recently discovered ultra-diffuse galaxies observed in the NGC 1052 group, posited in the literature to be dark matter (DM)-free. We compute the timescales for infall to the central nucleus due to dynamical friction (DF) for the GCs in these galaxies, using the shortest of these times to constrain how long ago a galaxy-galaxy interaction could have occurred. We go on to quantify the initial GC numbers and densities needed for significant collisional evolution to occur within the allotted times, and show that, if the hypothesis of a previous galaxy-galaxy interaction is correct, a paucity of low-mass GCs should be revealed by deeper observational surveys. If any are found, they should be more spatially extended than the currently observed GC population. Finally, we apply our method to these galaxies, in order to illustrate its efficacy in constraining their dynamical evolution. Our results motivate more complete observations of the GC luminosity functions in these galaxies, in addition to future studies aimed at combining the method presented here with a suite of numerical simulations in order to further constrain the origins of the curious GC populations in these (and other) galaxies., Comment: 12 pages, 5 figures, 2 table, accepted by ApJ

Nuclear star clusters that surround supermassive black holes (SMBHs) in galactic nuclei are thought to contain large numbers of black holes (BHs) and neutron stars (NSs), a fraction of which form binaries and could merge by Kozai-Lidov oscillations (KL). Triple compact objects are likely to be present, given what is known about the multiplicity of massive stars, whose life ends either as a NS or a BH. In this paper, we present a new possible scenario for merging BHs and NSs in galactic nuclei. We study the evolution of a triple black hole (BH) or neutron star (NS) system orbiting an SMBH in a galactic nucleus by means of direct high-precision $N$-body simulations, including Post-Newtonian terms. We find that the four-body dynamical interactions can increase the KL angle window for mergers compared to the binary case and make BH and NS binaries merge on shorter timescales. We show that the merger fraction can be up to $\sim 5$--$8$ times higher for triples than for binaries. Therefore, even if the triple fraction is only $\sim 10\%$--$20\%$ of the binary fraction, they could contribute to the merger events observed by LIGO/VIRGO in comparable numbers., Comment: 13 pages, 7 figures, 2 tables, accepted by MNRAS

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., 10 pages, 8 figures, 1 table; accepted for publication in MNRAS

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$)., Comment: 12 pages, 6 figures; submitted to ApJ

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.)., 13 pages, 5 figures, 3 tables, accepted for publication at ApJ

Hypervelocity binaries (HVBs) can be produced by a close encounter between stellar binaries and supermassive black hole (SMBH)-intermediate massive black hole (IMBH) binaries. We perform scattering experiments between tight stellar binaries and an SMBH-IMBH binary (within the observationally allowed parameters), using high-precision N-body integrations. Simulation results show that only a few HVBs can be produced in the Galactic center if it harbors a $4\times10^3M_\odot$ IMBH. However, we find that tight HVBs can be efficiently produced in prograde coplanar orbits with respect to the SMBH-IMBH binary. Therefore, the trajectories of the HVBs provide a promising way to constrain the IMBH orbital plane, if one exists. More importantly, we find that the ejected tight binaries are eccentric. For other galaxies hosting equal-mass SMBH binaries for which the HVB ejection rate is much higher relative to SMBH-IMBH binaries, these ejected compact binaries constitute promising extragalactic merger sources that could be detected by gravitational wave detectors.

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., 10 pages, 5 figures, 1 table; accepted for publication in MNRAS

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., accepted version. MNRAS

We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function is strongly correlated with the fraction of mass that the cluster has lost, independent of the cluster's initial mass, and nearly independent of its orbit and initial size. Thus, the mass function - initial mass relation can be used to determine a Galactic cluster's initial total stellar mass, if the initial stellar mass function is known. We apply the mass function - initial mass relation presented here to determine the initial stellar masses of 33 Galactic globular clusters, assuming an universal Kroupa initial mass function. Our study suggests that globular clusters had initial masses that were on average a factor of 4.5 times larger than their present day mass, with three clusters showing evidence for being 10 times more massive at birth., Comment: Revised Version: 11 pages, 4 figures, Accepted for publication in MNRAS