Showing 1–50 of 219 results for author: Rasio, F A

IMBH Progenitors from Stellar Collisions in Dense Star Clusters

Authors: Elena González Prieto, Newlin C. Weatherford, Giacomo Fragione, Kyle Kremer, Frederic A. Rasio

Abstract: Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes (IMBHs), which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies hav… ▽ More Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes (IMBHs), which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo (CMC) star cluster simulations -- the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 solar masses. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: Accepted for ApJ. Comments welcome. 19 pages, 10 figures

Binary black hole mergers and intermediate-mass black holes in dense star clusters with collisional runaways

Authors: Rujuta A. Purohit, Giacomo Fragione, Frederic A. Rasio, Grayson C. Petter, Ryan C. Hickox

Abstract: Intermediate-mass black holes (IMBHs) are believed to be the missing link between the supermassive black holes (BHs) found at the centers of massive galaxies and BHs formed through stellar core collapse. One of the proposed mechanisms for their formation is a collisional runaway process in high-density young star clusters, where an unusually massive object forms through repeated stellar collisions… ▽ More Intermediate-mass black holes (IMBHs) are believed to be the missing link between the supermassive black holes (BHs) found at the centers of massive galaxies and BHs formed through stellar core collapse. One of the proposed mechanisms for their formation is a collisional runaway process in high-density young star clusters, where an unusually massive object forms through repeated stellar collisions and mergers, eventually collapsing to form an IMBH. This seed IMBH could then grow further through binary mergers with other stellar-mass BHs. Here we investigate the gravitational-wave (GW) signals produced during these later IMBH--BH mergers. We use a state-of-the-art semi-analytic approach to study the stellar dynamics and to characterize the rates and properties of IMBH--BH mergers. We also study the prospects for detection of these mergers by current and future GW observatories, both space-based (LISA) and ground-based (LIGO Voyager, Einstein Telescope, and Cosmic Explorer). We find that most of the merger signals could be detected, with some of them being multi-band sources. Therefore, GWs represent a unique tool to test the collisional runaway scenario and to constrain the population of dynamically assembled IMBHs. △ Less

Submitted 13 March, 2024; v1 submitted 12 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astronomical Journal

Stellar Escape from Globular Clusters. II. Clusters May Eat Their Own Tails

Authors: Newlin C. Weatherford, Frederic A. Rasio, Sourav Chatterjee, Giacomo Fragione, Fulya Kıroğlu, Kyle Kremer

Abstract: We apply for the first time the Monte Carlo star cluster modeling method to study tidal tail and stellar stream formation from globular clusters, assuming a circular orbit in a smooth Galactic potential. Approximating energetically unbound bodies (potential escapers; PEs) as collisionless enables this fast but spherically symmetric method to capture asymmetric tidal phenomena with unprecedented de… ▽ More We apply for the first time the Monte Carlo star cluster modeling method to study tidal tail and stellar stream formation from globular clusters, assuming a circular orbit in a smooth Galactic potential. Approximating energetically unbound bodies (potential escapers; PEs) as collisionless enables this fast but spherically symmetric method to capture asymmetric tidal phenomena with unprecedented detail. Beyond reproducing known stream features, including epicyclic overdensities, we show how 'returning tidal tails' may form after the stream fully circumnavigates the Galaxy back to the cluster, enhancing the stream's velocity dispersion. While a realistically clumpy, time-dependent Galactic potential may disrupt such tails, they warrant scrutiny as potentially excellent constraints on the Galactic potential's history and substructure. Re-examining the escape timescale $Δt$ of PEs, we find new behavior related to chaotic scattering in the three-body problem; the $Δt$ distribution features sharp plateaus corresponding to distinct locally smooth patches of the chaotic saddle separating the phase space basins of escape. We study for the first time $Δt$ in an evolving cluster, finding that $Δt\sim(E_{\rm J}^{-0.1},E_{\rm J}^{-0.4})$ for PEs with (low, high) Jacobi energy $E_{\rm J}$, flatter than for a static cluster ($E_{\rm J}^{-2}$). Accounting for cluster mass loss and internal evolution -- and (roughly) for ongoing relaxation among PEs -- lowers the median $Δt$ from ${\sim}10\,$Gyr to ${\lesssim}100\,$Myr. We finally outline future improvements to escape physics in the Monte Carlo method intended to enable both the first large-parameter-space studies of tidal tail/stellar stream formation from full globular cluster simulations and detailed comparisons to stream observations. △ Less

Submitted 8 April, 2024; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 31 pages, 19 Figures (+ linked video), 1 Table, accepted to ApJ

Single Millisecond Pulsars from Dynamical Interaction Processes in Dense Star Clusters

Authors: Claire S. Ye, Kyle Kremer, Scott M. Ransom, Frederic A. Rasio

Abstract: Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we cr… ▽ More Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we create N-body models that reproduce the ratio of single-to-binary pulsars in Milky-Way-like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play the dominant role: (1) tidal disruption of main-sequence stars by neutron stars; and (2) gravitational collapse of heavy white-dwarf-binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC. △ Less

Submitted 19 January, 2024; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 19 pages, 12 figures, 4 tables. Published at ApJ

Demographics of Hierarchical Black Hole Mergers in Dense Star Clusters

Authors: Giacomo Fragione, Frederic A. Rasio

Abstract: With about a hundred mergers of binary black holes (BBHs) detected via gravitational waves by the LIGO-Virgo-KAGRA (LVK) Collaboration, our understanding of the darkest objects in the Universe has taken unparalleled steps forward. While most of the events are expected to consist of BHs directly formed from the collapse of massive stars, some may contain the remnants of previous BBH mergers. In the… ▽ More With about a hundred mergers of binary black holes (BBHs) detected via gravitational waves by the LIGO-Virgo-KAGRA (LVK) Collaboration, our understanding of the darkest objects in the Universe has taken unparalleled steps forward. While most of the events are expected to consist of BHs directly formed from the collapse of massive stars, some may contain the remnants of previous BBH mergers. In the most massive globular clusters and in nuclear star clusters, successive mergers can produce second- (2G) or higher-generation BHs, and even form intermediate-mass BHs. Overall, we predict that up to $\sim 10\%$, $\sim 1\%$ or $\sim 0.1\%$ of the BBH mergers have one component being a 2G, 3G, or 4G BH, respectively. Assuming that $\sim 500$ BBH mergers will be detected in O4 by LVK, this means that $\sim 50$, $\sim 5$, or $\sim 0.5$ events, respectively, will involve a 2G, 3G, or 4G BH, if most sources are produced dynamically in dense star clusters. With their distinctive signatures of higher masses and spins, such hierarchical mergers offer an unprecedented opportunity to learn about the BH populations in the densest stellar systems and to shed light on the elusive intermediate-mass BHs that may form therein. △ Less

Submitted 30 May, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: 19 pages, 13 figures, accepted by ApJ

Stellar Escape from Globular Clusters. I. Escape Mechanisms and Properties at Ejection

Authors: Newlin C. Weatherford, Fulya Kıroğlu, Giacomo Fragione, Sourav Chatterjee, Kyle Kremer, Frederic A. Rasio

Abstract: The theory of stellar escape from globular clusters (GCs) dates back nearly a century, especially the gradual evaporation of GCs via two-body relaxation coupled with external tides. More violent ejection can also occur via strong gravitational scattering, supernovae, gravitational wave-driven mergers, tidal disruption events, and physical collisions, but comprehensive study of the many escape mech… ▽ More The theory of stellar escape from globular clusters (GCs) dates back nearly a century, especially the gradual evaporation of GCs via two-body relaxation coupled with external tides. More violent ejection can also occur via strong gravitational scattering, supernovae, gravitational wave-driven mergers, tidal disruption events, and physical collisions, but comprehensive study of the many escape mechanisms has been limited. Recent exquisite kinematic data from the Gaia space telescope has revealed numerous stellar streams in the Milky Way (MW) and traced the origin of many to specific MWGCs, highlighting the need for further examination of stellar escape from these clusters. In this study, the first of a series, we lay the groundwork for detailed follow-up comparisons between Cluster Monte Carlo (CMC) GC models and the latest Gaia data on the outskirts of MWGCs, their tidal tails, and associated streams. We thoroughly review escape mechanisms from GCs and examine their relative contributions to the escape rate, ejection velocities, and escaper demographics. We show for the first time that three-body binary formation may dominate high-speed ejection from typical MWGCs, potentially explaining some of the hypervelocity stars in the MW. Due to their mass, black holes strongly catalyze this process, and their loss at the onset of observable core collapse, characterized by a steep central brightness profile, dramatically curtails three-body binary formation, despite the increased post-collapse density. We also demonstrate that even when born from a thermal eccentricity distribution, escaping binaries have significantly nonthermal eccentricities consistent with the roughly uniform distribution observed in the Galactic field. △ Less

Submitted 27 February, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 28 pages, 6 figures, 4 tables, accepted to ApJ

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

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

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

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

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

Partial Tidal Disruption of Main-Sequence Stars by Intermediate-Mass Black Holes

Authors: Fulya Kıroğlu, James C. Lombardi Jr., Kyle Kremer, Giacomo Fragione, Shane Fogarty, Frederic A. Rasio

Abstract: We study close encounters of a $1\,M_{\odot}$ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, $M_{\rm BH}= 100-10^4\,M_{\odot}$. Possible outcomes vary from a small tidal perturbation for weak enc… ▽ More We study close encounters of a $1\,M_{\odot}$ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, $M_{\rm BH}= 100-10^4\,M_{\odot}$. Possible outcomes vary from a small tidal perturbation for weak encounters all the way to partial or full disruption for stronger encounters. We find that stronger encounters lead to increased mass loss at the first pericenter passage, in many cases ejecting the partially disrupted star on an unbound orbit. For encounters that initially produce a bound system, with only partial stripping of the star, the fraction of mass stripped from the star increases with each subsequent pericenter passage and a stellar remnant of finite mass is ultimately ejected in all cases. The critical penetration depth that separates bound and unbound remnants has a dependence on the black hole mass when $M_{\rm BH} \lesssim 10^3\,M_{\odot}$. We also find that the number of successive close passages before ejection decreases as we go from the stellar-mass black hole to the intermediate-mass black hole regime. For instance, after an initial encounter right at the classical tidal disruption limit, a $1\,M_{\odot}$ star undergoes 16 (5) pericenter passages before ejection from a $10\,M_{\odot}$ ($100\,M_{\odot}$) black hole. Observations of periodic flares from these repeated close passages could in principle indicate signatures of a partial tidal disruption event. △ Less

Submitted 12 May, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 19 pages, 8 figures, 2 tables. Submitted to ApJ

Intermediate-mass Black Holes on the Run from Young Star Clusters

Authors: Elena González Prieto, Kyle Kremer, Giacomo Fragione, Miguel A. S. Martinez, Newlin C. Weatherford, Michael Zevin, Frederic A. Rasio

Abstract: The existence of black holes (BHs) with masses in the range between stellar remnants and supermassive BHs has only recently become unambiguously established. GW190521, a gravitational wave signal detected by the LIGO/Virgo Collaboration, provides the first direct evidence for the existence of such intermediate-mass BHs (IMBHs). This event sparked and continues to fuel discussion on the possible fo… ▽ More The existence of black holes (BHs) with masses in the range between stellar remnants and supermassive BHs has only recently become unambiguously established. GW190521, a gravitational wave signal detected by the LIGO/Virgo Collaboration, provides the first direct evidence for the existence of such intermediate-mass BHs (IMBHs). This event sparked and continues to fuel discussion on the possible formation channels for such massive BHs. As the detection revealed, IMBHs can form via binary mergers of BHs in the "upper mass gap" ($\approx40 -120\,M_{\odot}$). Alternatively, IMBHs may form via the collapse of a very massive star formed through stellar collisions and mergers in dense star clusters. In this study, we explore the formation of IMBHs with masses between $120$ and $500\,M_{\odot}$ in young, massive star clusters using state-of-the-art Cluster Monte Carlo ($\texttt{CMC}$) models. We examine the evolution of IMBHs throughout their dynamical lifetimes, ending with their ejection from the parent cluster due to gravitational radiation recoil from BH mergers, or dynamical recoil kicks from few-body scattering encounters. We find that $ \textit{all}$ of the IMBHs in our models are ejected from the host cluster within the first $\sim 500$ Myr, indicating a low retention probability of IMBHs in this mass range for globular clusters today. We estimate the peak IMBH merger rate to be $\mathcal{R} \approx 2 \, \rm{Gpc}^{-3}\,\rm{yr}^{-1}$ at redshift $z \approx 2$. △ Less

Submitted 17 October, 2022; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for ApJ. Comments welcome. 16 pages and 10 figures

Formation of Low-mass Black Holes and Single Millisecond Pulsars in Globular Clusters

Authors: Kyle Kremer, Claire S. Ye, Fulya Kıroğlu, James C. Lombardi Jr., Scott M. Ransom, Frederic A. Rasio

Abstract: Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in tidal disruption of the star. Using $N$-body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on host cluster properties. We find that tidal disruption events occur most frequently in… ▽ More Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in tidal disruption of the star. Using $N$-body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on host cluster properties. We find that tidal disruption events occur most frequently in core-collapsed globular clusters and that roughly $25\%$ of the disrupted stars are merger products (i.e., blue straggler stars). Using hydrodynamic simulations, we model the tidal disruptions themselves (over timescales of days) to determine the mass bound to the neutron star and the properties of the accretion disks formed. In general, we find that roughly $80-90\%$ of the initial stellar mass becomes bound to the neutron star following disruption. Additionally, we find that neutron stars receive impulsive kicks of up to about $20\,$km/s as a result of the asymmetry of unbound ejecta; these kicks place these neutron stars on elongated orbits within their host cluster, with apocenter distances well outside the cluster core. Finally, we model the evolution of the (hypercritical) accretion disks on longer timescales (days to years after disruption) to estimate the accretion rate onto the neutron stars and accompanying spin-up. As long as $\gtrsim1\%$ of the bound mass accretes onto the neutron star, millisecond spin periods can be attained. We argue the growing numbers of isolated millisecond pulsars observed in globular clusters may have formed, at least in part, through this mechanism. In the case of significant mass growth, some of these neutron stars may collapse to form low-mass ($\lesssim3\,M_{\odot}$) black holes. △ Less

Submitted 27 June, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 14 pages, 4 figures, 2 tables. Accepted for publication in ApJL

Eccentric Mergers of Intermediate-Mass Black Holes from Evection Resonances in AGN Disks

Authors: Diego J. Muñoz, Nicholas C. Stone, Cristobal Petrovich, Frederic A. Rasio

Abstract: We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) si… ▽ More We apply the theory of nonlinear resonance capture to the problem of a black hole binary (BHB) orbiting a supermassive black hole (SMBH) while embedded in the accretion disk of an active galactic nucleus (AGN). If successful, resonance capture can trigger dramatic growth in the BHB eccentricity, with important consequences for the BHB merger timescale, as well as for the gravitational wave (GW) signature of such an eccentric merger. This resonance capture may occur when the orbital period around the SMBH (the "outer binary") and the apsidal precession of the BHB (the "inner binary") are in a 1:1 commensurability. This effect is analogous to the phenomenon of lunar evection resonance in the early Sun-Earth-Moon system, with the distinction that in the present case, the BHB apsidal precession is due to general relativity, rather than rotationally-induced distortion. In contrast to the case of lunar evection, however, the inner binary undergoes orbital decay driven by GW emission, rather than orbital expansion driven by tidal dissipation. This distinction fundamentally alters the three-body dynamics, forbidding resonance capture, and limiting eccentricity growth. However, if the BHB migrates through of a gaseous AGN disk, the change in the outer binary can counterbalance the suppressing effect of BHB decay, permitting evection resonance capture and the production of eccentric BHB mergers. We compute the likelihood of resonance capture assuming an agnostic distribution of parameters for the three bodies involved and for the properties of the AGN disk. We find that intermediate-mass ratio BHBs (involving an intermediate-mass black hole and a stellar-mass black hole) are the most likely to be captured into evection resonance and thus undergo an eccentric merger. We also compute the GW signature of these mergers, showing that they can enter the LISA band while eccentric. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 15 pages (before references and appendices), 11 figures, comments welcome

Merger rates of intermediate-mass black hole binaries in nuclear star clusters

Authors: Giacomo Fragione, Abraham Loeb, Bence Kocsis, Frederic A. Rasio

Abstract: Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravita… ▽ More Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravitational waves (GWs), which represent an unprecedented opportunity to reveal IMBHs. In this paper, we analyze the statistical results of a wide range of numerical simulations, which encompass different cluster metallicities, initial BH seed masses, and initial BH spins, and we compute the merger rate of IMBH binaries. We find that merger rates are in the range $0.01$-$10$\,Gpc$^{-3}$\,yr$^{-1}$ depending on IMBH masses. We also compute the number of multi-band detections in ground-based and space-based observatories. Our model predicts that a few merger events per year should be detectable with LISA, DECIGO, ET, and LIGO for IMBHs with masses $\lesssim 1000\msun$, and a few tens of merger events per year with DECIGO, ET, and LIGO only. △ Less

Submitted 7 June, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 12 pages, 6 figures, accepted by ApJ

Hydrodynamics of Collisions and Close Encounters between Stellar Black Holes and Main-sequence Stars

Authors: Kyle Kremer, James C. Lombardi Jr., Wenbin Lu, Anthony L. Piro, Frederic A. Rasio

Abstract: Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a… ▽ More Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a wide range of hydrodynamic simulations of close encounters between black holes and main-sequence stars that collectively cover the parameter space of interest, and we identify and classify the various possible outcomes. In the case of nearly head-on collisions, the star is completely disrupted with roughly half of the stellar material becoming bound to the black hole. For more distant encounters near the classical tidal-disruption radius, the star is only partially disrupted on the first pericenter passage. Depending upon the interaction details, the partially disrupted stellar remnant may be tidally captured by the black hole or become unbound (in some cases, receiving a sufficiently large impulsive kick from asymmetric mass loss to be ejected from its host cluster). In the former case, the star will undergo additional pericenter passages before ultimately being disrupted fully. Based on the properties of the material bound to the black hole at the end of our simulations (in particular, the total bound mass and angular momentum), we comment upon the expected accretion process and associated electromagnetic signatures that are likely to result. △ Less

Submitted 24 June, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 23 pages, 11 figures, 3 tables. Accepted for publication in ApJ

Gravitational Microlensing Rates in Milky Way Globular Clusters

Authors: Fulya Kıroğlu, Newlin C. Weatherford, Kyle Kremer, Claire S. Ye, Giacomo Fragione, Frederic A. Rasio

Abstract: Many recent observational and theoretical studies suggest that globular clusters (GCs) host compact object populations large enough to play dominant roles in their overall dynamical evolution. Yet direct detection, particularly of black holes and neutron stars, remains rare and limited to special cases, such as when these objects reside in close binaries with bright companions. Here we examine the… ▽ More Many recent observational and theoretical studies suggest that globular clusters (GCs) host compact object populations large enough to play dominant roles in their overall dynamical evolution. Yet direct detection, particularly of black holes and neutron stars, remains rare and limited to special cases, such as when these objects reside in close binaries with bright companions. Here we examine the potential of microlensing detections to further constrain these dark populations. Based on state-of-the-art GC models from the CMC Cluster Catalog, we estimate the microlensing event rates for black holes, neutron stars, white dwarfs, and, for comparison, also for M dwarfs in Milky Way GCs, as well as the effects of different initial conditions on these rates. Among compact objects, we find that white dwarfs dominate the microlensing rates, simply because they largely dominate by numbers. We show that microlensing detections are in general more likely in GCs with higher initial densities, especially in clusters that undergo core collapse. We also estimate microlensing rates in the specific cases of M22 and 47 Tuc using our best-fitting models for these GCs. Because their positions on the sky lie near the rich stellar backgrounds of the Galactic bulge and the Small Magellanic Cloud, respectively, these clusters are among the Galactic GCs best-suited for dedicated microlensing surveys. The upcoming 10-year survey with the Rubin Observatory may be ideal for detecting lensing events in GCs. △ Less

Submitted 27 February, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: v2: 17 pages, 5 figures

Compact Object Modeling in the Globular Cluster 47 Tucanae

Authors: Claire S. Ye, Kyle Kremer, Carl L. Rodriguez, Nicholas Z. Rui, Newlin C. Weatherford, Sourav Chatterjee, Giacomo Fragione, Frederic A. Rasio

Abstract: The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc… ▽ More The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47~Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47~Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $70\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47~Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables. △ Less

Submitted 6 June, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 23 pages, 15 figures, 4 tables, published at ApJ

Impact of natal kicks on merger rates and spin-orbit misalignments of black hole -- neutron star mergers

Authors: Giacomo Fragione, Abraham Loeb, Frederic A. Rasio

Abstract: The long wait for the detection of merging black hole -- neutron star (BH--NS) binaries is finally over with the announcement by the LIGO/Virgo/Kagra collaboration of GW200105 and GW200115. Remarkably, the primary of GW200115 has a negative spin projection onto the orbital angular momentum, with about $90\%$ probability. Merging BH--NS binaries are expected to form mainly through the evolution of… ▽ More The long wait for the detection of merging black hole -- neutron star (BH--NS) binaries is finally over with the announcement by the LIGO/Virgo/Kagra collaboration of GW200105 and GW200115. Remarkably, the primary of GW200115 has a negative spin projection onto the orbital angular momentum, with about $90\%$ probability. Merging BH--NS binaries are expected to form mainly through the evolution of massive binary stars in the field, since their dynamical formation in dense star clusters is strongly suppressed by mass segregation. In this paper, we carry out a systematic statistical study of the binary stars that evolve to form a BH--NS binary, considering different metallicities and taking into account the uncertainties on the natal kick distributions for BHs and NSs and on the common envelope phase of binary evolution. Under the assumption that the initial stellar spins are aligned with the binary angular momentum, we show that both large natal kicks ($\gtrsim 150$ km s$^{-1}$) and high efficiencies for common envelope ejection are required to simultaneously explain the inferred high merger rates and the large spin-orbit misalignment of GW200115. △ Less

Submitted 31 August, 2021; v1 submitted 14 August, 2021; originally announced August 2021.

Comments: 5 pages, 4 figures, accepted by ApJL

Repeated mergers, mass-gap black holes, and formation of intermediate-mass black holes in nuclear star clusters

Authors: Giacomo Fragione, Bence Kocsis, Frederic A. Rasio, Joseph Silk

Abstract: Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly $50\,M_\odot$ and $100\,M_\odot$, while, above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational wave emission with the current LIGO/Virgo/Kagra interferometers and future detecto… ▽ More Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly $50\,M_\odot$ and $100\,M_\odot$, while, above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. We show that a massive stellar BH seed can easily grow to $\sim 10^3 - 10^4\,M_\odot$ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to larger final BH masses. We also show that the growing BH spin tends to decrease in magnitude with the number of mergers, so that a negative correlation exists between final mass and spin of the resulting IMBHs. Assumptions about the birth spins of stellar BHs affect our results significantly, with low birth spins leading to the production of a larger population of massive BHs. △ Less

Submitted 30 January, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 13 pages, 8 figures, 2 tables, accepted by ApJ

Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code

Authors: Carl L. Rodriguez, Newlin C. Weatherford, Scott C. Coughlin, Pau Amaro Seoane, Katelyn Breivik, Sourav Chatterjee, Giacomo Fragione, Fulya Kıroğlu, Kyle Kremer, Nicholas Z. Rui, Claire S. Ye, Michael Zevin, Frederic A. Rasio

Abstract: We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using Hénon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star a… ▽ More We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using Hénon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest ($N = 10^8$) star-by-star $N$-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap). △ Less

Submitted 11 October, 2021; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: Code is available at https://clustermontecarlo.github.io/ 25 pages, 8 Figures, Matches version accepted by ApJS

White Dwarf Subsystems in Core-Collapsed Globular Clusters

Authors: Kyle Kremer, Nicholas Z. Rui, Newlin C. Weatherford, Sourav Chatterjee, Giacomo Fragione, Frederic A. Rasio, Carl L. Rodriguez, Claire S. Ye

Abstract: Numerical and observational evidence suggests that massive white dwarfs dominate the innermost regions of core-collapsed globular clusters by both number and total mass. Using NGC 6397 as a test case, we constrain the features of white dwarf populations in core-collapsed clusters, both at present day and throughout their lifetimes. The dynamics of these white dwarf subsystems have a number of astr… ▽ More Numerical and observational evidence suggests that massive white dwarfs dominate the innermost regions of core-collapsed globular clusters by both number and total mass. Using NGC 6397 as a test case, we constrain the features of white dwarf populations in core-collapsed clusters, both at present day and throughout their lifetimes. The dynamics of these white dwarf subsystems have a number of astrophysical implications. We demonstrate that the collapse of globular cluster cores is ultimately halted by the dynamical burning of white dwarf binaries. We predict core-collapsed clusters in the local universe yield a white dwarf merger rate of $\mathcal{O}(10\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$, roughly $0.1-1\%$ of the observed Type Ia supernova rate. We show that prior to merger, inspiraling white dwarf binaries will be observable as gravitational wave sources at milli- and decihertz frequencies. Over $90\%$ of these mergers have a total mass greater than the Chandrasekhar limit. If the merger/collision remnants are not destroyed completely in an explosive transient, we argue the remnants may be observed in core-collapsed clusters as either young neutron stars/pulsars/magnetars (in the event of accretion-induced collapse) or as young massive white dwarfs offset from the standard white dwarf cooling sequence. Finally, we show collisions between white dwarfs and main sequence stars, which may be detectable as bright transients, occur at a rate of $\mathcal{O}(100\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$ in the local universe. We find that these collisions lead to depletion of blue straggler stars and main sequence star binaries in the centers of core-collapsed clusters. △ Less

Submitted 21 May, 2021; v1 submitted 23 April, 2021; originally announced April 2021.

Comments: Submitted to ApJ, 26 pages, 12 figures, 3 tables. Comments welcome

No Black Holes in NGC 6397

Authors: Nicholas Z. Rui, Newlin Weatherford, Kyle Kremer, Sourav Chatterjee, Giacomo Fragione, Frederic A. Rasio, Carl L. Rodriguez, Claire S. Ye

Abstract: Recently, \citet{vitral2021does} detected a central concentration of dark objects in the core-collapsed globular cluster NGC 6397, which could be interpreted as a subcluster of stellar-mass black holes. However, it is well established theoretically that any significant number of black holes in the cluster would provide strong dynamical heating and is fundamentally inconsistent with this cluster's… ▽ More Recently, \citet{vitral2021does} detected a central concentration of dark objects in the core-collapsed globular cluster NGC 6397, which could be interpreted as a subcluster of stellar-mass black holes. However, it is well established theoretically that any significant number of black holes in the cluster would provide strong dynamical heating and is fundamentally inconsistent with this cluster's core-collapsed profile. Claims of intermediate-mass black holes in core-collapsed clusters should similarly be treated with suspicion, for reasons that have been understood theoretically for many decades. Instead, the central dark population in NGC 6397 is exactly accounted for by a compact subsystem of white dwarfs, as we demonstrate here by inspection of a previously published model that provides a good fit to this cluster. These central subclusters of heavy white dwarfs are in fact a generic feature of core-collapsed clusters, while central black hole subclusters are present in all {\em non\/}-collapsed clusters. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 3 pages, 1 figure, to submit to AAS journals

Matching Globular Cluster Models to Observations

Authors: Nicholas Z. Rui, Kyle Kremer, Newlin C. Weatherford, Sourav Chatterjee, Frederic A. Rasio, Carl L. Rodriguez, Claire S. Ye

Abstract: As ancient, gravitationally bound stellar populations, globular clusters are abundant, vibrant laboratories characterized by high frequencies of dynamical interactions coupled to complex stellar evolution. Using surface brightness and velocity dispersion profiles from the literature, we fit $59$ Milky Way globular clusters to dynamical models from the \texttt{CMC Cluster Catalog}. Without doing an… ▽ More As ancient, gravitationally bound stellar populations, globular clusters are abundant, vibrant laboratories characterized by high frequencies of dynamical interactions coupled to complex stellar evolution. Using surface brightness and velocity dispersion profiles from the literature, we fit $59$ Milky Way globular clusters to dynamical models from the \texttt{CMC Cluster Catalog}. Without doing any interpolation, and without any directed effort to fit any particular cluster, $26$ globular clusters are well-matched by at least one of our models. We discuss in particular the core-collapsed clusters NGC 6293, NGC 6397, NGC 6681, and NGC 6624, and the non-core-collapsed clusters NGC 288, NGC 4372, and NGC 5897. As NGC 6624 lacks well-fitting snapshots on the main \texttt{CMC Cluster Catalog}, we run six additional models in order to refine the fit. We calculate metrics for mass segregation, explore the production of compact object sources such as millisecond pulsars, cataclysmic variables, low-mass X-ray binaries, and stellar-mass black holes, finding reasonable agreement with observations. Additionally, closely mimicking observational cuts, we extract the binary fraction from our models, finding good agreement except in the dense core regions of core-collapsed clusters. Accompanying this paper are a number of \textsf{python} methods for examining the publicly accessible \texttt{CMC Cluster Catalog}, as well as any other models generated using \texttt{CMC}. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 22 pages, 11 figures, 5 tables; accepted to ApJ

The Observed Rate of Binary Black Hole Mergers can be Entirely Explained by Globular Clusters

Authors: Carl L. Rodriguez, Kyle Kremer, Sourav Chatterjee, Giacomo Fragione, Abraham Loeb, Frederic A. Rasio, Newlin C. Weatherford, Claire S. Ye

Abstract: Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshif… ▽ More Since the first signal in 2015, the gravitational-wave detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical work on the astrophysical origin of these objects. We show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters are the dominant formation channel, we emphasize that many formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior. △ Less

Submitted 27 January, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: 4 pages, one figure, matches version in RNAAS

Journal ref: Res. Notes AAS 5 19 (2021)

Black Hole Mergers from Star Clusters with Top-Heavy Initial Mass Functions

Authors: Newlin C. Weatherford, Giacomo Fragione, Kyle Kremer, Sourav Chatterjee, Claire S. Ye, Carl L. Rodriguez, Frederic A. Rasio

Abstract: Recent observations of globular clusters (GCs) provide evidence that the stellar initial mass function (IMF) may not be universal, suggesting specifically that the IMF grows increasingly top-heavy with decreasing metallicity and increasing gas density. Non-canonical IMFs can greatly affect the evolution of GCs, mainly because the high end determines how many black holes (BHs) form. Here we compute… ▽ More Recent observations of globular clusters (GCs) provide evidence that the stellar initial mass function (IMF) may not be universal, suggesting specifically that the IMF grows increasingly top-heavy with decreasing metallicity and increasing gas density. Non-canonical IMFs can greatly affect the evolution of GCs, mainly because the high end determines how many black holes (BHs) form. Here we compute a new set of GC models, varying the IMF within observational uncertainties. We find that GCs with top-heavy IMFs lose most of their mass within a few Gyr through stellar winds and tidal stripping. Heating of the cluster through BH mass segregation greatly enhances this process. We show that, as they approach complete dissolution, GCs with top-heavy IMFs can evolve into 'dark clusters' consisting of mostly BHs by mass. In addition to producing more BHs, GCs with top-heavy IMFs also produce many more binary BH (BBH) mergers. Even though these clusters are short-lived, mergers of ejected BBHs continue at a rate comparable to, or greater than, what is found for long-lived GCs with canonical IMFs. Therefore these clusters, although they are no longer visible today, could still contribute significantly to the local BBH merger rate detectable by LIGO/Virgo, especially for sources with higher component masses well into the BH mass gap. We also report that one of our GC models with a top-heavy IMF produces dozens of intermediate-mass black holes (IMBHs) with masses $M>100\,{\rm M_\odot}$, including one with $M>500\,{\rm M_\odot}$. Ultimately, additional gravitational wave observations will provide strong constraints on the stellar IMF in old GCs and the formation of IMBHs at high redshift. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 8 pages, 4 figures, 1 table; accepted for publication to ApJL

Intermediate-mass Black Holes from High Massive-star Binary Fractions in Young Star Clusters

Authors: Elena González, Kyle Kremer, Sourav Chatterjee, Giacomo Fragione, Carl L. Rodriguez, Newlin C. Weatherford, Claire S. Ye, Frederic A. Rasio

Abstract: Black holes formed in dense star clusters, where dynamical interactions are frequent, may have fundamentally different properties than those formed through isolated stellar evolution. Theoretical models for single star evolution predict a gap in the black hole mass spectrum from roughly $40-120\,M_{\odot}$ caused by (pulsational) pair-instability supernovae. Motivated by the recent LIGO/Virgo even… ▽ More Black holes formed in dense star clusters, where dynamical interactions are frequent, may have fundamentally different properties than those formed through isolated stellar evolution. Theoretical models for single star evolution predict a gap in the black hole mass spectrum from roughly $40-120\,M_{\odot}$ caused by (pulsational) pair-instability supernovae. Motivated by the recent LIGO/Virgo event GW190521, we investigate whether black holes with masses within or in excess of this "upper-mass gap" can be formed dynamically in young star clusters through strong interactions of massive stars in binaries. We perform a set of $N$-body simulations using the CMC cluster-dynamics code to study the effects of the high-mass binary fraction on the formation and collision histories of the most massive stars and their remnants. We find that typical young star clusters with low metallicities and high binary fractions in massive stars can form several black holes in the upper-mass gap and often form at least one intermediate-mass black hole. These results provide strong evidence that dynamical interactions in young star clusters naturally lead to the formation of more massive black hole remnants. △ Less

Submitted 21 January, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 10 pages, 2 figures, 1 table. Accepted for publication in ApJ Letters. Comments welcome

Fast Optical Transients from Stellar-Mass Black Hole Tidal Disruption Events in Young Star Clusters

Authors: Kyle Kremer, Wenbin Lu, Anthony L. Piro, Sourav Chatterjee, Frederic A. Rasio, Claire S. Ye

Abstract: Observational evidence suggests that the majority of stars may have been born in stellar clusters or associations. Within these dense environments, dynamical interactions lead to high rates of close stellar encounters. A variety of recent observational and theoretical indications suggest stellar-mass black holes may be present and play an active dynamical role in stellar clusters of all masses. In… ▽ More Observational evidence suggests that the majority of stars may have been born in stellar clusters or associations. Within these dense environments, dynamical interactions lead to high rates of close stellar encounters. A variety of recent observational and theoretical indications suggest stellar-mass black holes may be present and play an active dynamical role in stellar clusters of all masses. In this study, we explore the tidal disruption of main sequence stars by stellar-mass black holes in young star clusters. We compute a suite of over 3000 independent $N$-body simulations that cover a range in cluster mass, metallicity, and half-mass radii. We find stellar-mass black hole tidal disruption events (TDEs) occur at an overall rate of up to roughly $200\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ in young stellar clusters in the local universe. These TDEs are expected to have several characteristic features, namely fast rise times of order a day, peak X-ray luminosities of at least $10^{44}\,\rm{erg\,s}^{-1}$, and bright optical luminosities (roughly $10^{41}-10^{44}\,\rm{erg\,s}^{-1}$) associated with reprocessing by a disk wind. In particular, we show these events share many features in common with the emerging class of Fast Blue Optical Transients. △ Less

Submitted 23 March, 2021; v1 submitted 4 December, 2020; originally announced December 2020.

Comments: 17 pages, 5 figures, 3 tables. Accepted for publication in ApJ

Black Hole Mergers from Hierarchical Triples in Dense Star Clusters

Authors: Miguel A. S. Martinez, Giacomo Fragione, Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Johan Samsing, Claire S. Ye, Newlin C. Weatherford, Michael Zevin, Smadar Naoz, Frederic A. Rasio

Abstract: Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical trip… ▽ More Hierarchical triples are expected to be produced by the frequent binary-mediated interactions in the cores of globular clusters. In some of these triples, the tertiary companion can drive the inner binary to merger following large eccentricity oscillations, as a result of the eccentric Kozai-Lidov mechanism. In this paper, we study the dynamics and merger rates of black hole (BH) hierarchical triples, formed via binary--binary encounters in the CMC Cluster Catalog, a suite of cluster simulations with present-day properties representative of the Milky Way's globular clusters. We compare the properties of the mergers from triples to the other merger channels in dense star clusters, and show that triple systems do not produce significant differences in terms of mass and effective spin distribution. However, they represent an important pathway for forming eccentric mergers, which could be detected by LIGO--Virgo/KAGRA (LVK), and future missions such as LISA and DECIGO. We derive a conservative lower limit for the merger rate from this channel of $0.35$ Gpc$^{-3}$yr$^{-1}$ in the local Universe and up to $\sim9\%$ of these events may have a detectable eccentricity at LVK design sensitivity. Additionally, we find that triple systems could play an important role in retaining second-generation BHs, which can later merge again in the core of the host cluster. △ Less

Submitted 21 September, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 21 Pages, 11 Figures, 2 Tables, Accepted for publication in ApJ

On the Origin of GW190521-like events from repeated black hole mergers in star clusters

Authors: Giacomo Fragione, Abraham Loeb, Frederic A. Rasio

Abstract: LIGO and Virgo have reported the detection of GW190521, from the merger of a binary black hole (BBH) with a total mass around $150$ M$_\odot$. While current stellar models limit the mass of any black hole (BH) remnant to about $40 - 50$ M$_\odot$, more massive BHs can be produced dynamically through repeated mergers in the core of a dense star cluster. The process is limited by the recoil kick (du… ▽ More LIGO and Virgo have reported the detection of GW190521, from the merger of a binary black hole (BBH) with a total mass around $150$ M$_\odot$. While current stellar models limit the mass of any black hole (BH) remnant to about $40 - 50$ M$_\odot$, more massive BHs can be produced dynamically through repeated mergers in the core of a dense star cluster. The process is limited by the recoil kick (due to anisotropic emission of gravitational radiation) imparted to merger remnants, which can escape the parent cluster, thereby terminating growth. We study the role of the host cluster metallicity and escape speed in the buildup of massive BHs through repeated mergers. Almost independent of host metallicity, we find that a BBH of about $150$ M$_\odot$ could be formed dynamically in any star cluster with escape speed $\gtrsim 200$ km s$^{-1}$, as found in galactic nuclear star clusters as well as the most massive globular clusters and super star clusters. Using an inspiral-only waveform, we compute the detection probability for different primary masses ($\ge 60$ M$_\odot$) as a function of secondary mass and find that the detection probability increases with secondary mass and decreases for larger primary mass and redshift. Future additional detections of massive BBH mergers will be of fundamental importance for understanding the growth of massive BHs through dynamics and the formation of intermediate-mass BHs. △ Less

Submitted 25 September, 2020; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 8 pages, 5 figures, accepted by ApJ Lett

Demographics of triple systems in dense star clusters

Authors: Giacomo Fragione, Miguel A. S. Martinez, Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Claire S. Ye, Newlin C. Weatherford, Smadar Naoz, Frederic A. Rasio

Abstract: Depending on the stellar type, more than $\sim 50$\% and $\sim 15\%$ of stars in the field have at least one and two stellar companions, respectively. Hierarchical systems can be assembled dynamically in dense star clusters, as a result of few-body encounters among stars and/or compact remnants in the cluster core. In this paper, we present the demographics of stellar and compact-object triples fo… ▽ More Depending on the stellar type, more than $\sim 50$\% and $\sim 15\%$ of stars in the field have at least one and two stellar companions, respectively. Hierarchical systems can be assembled dynamically in dense star clusters, as a result of few-body encounters among stars and/or compact remnants in the cluster core. In this paper, we present the demographics of stellar and compact-object triples formed via binary--binary encounters in the \texttt{CMC Cluster Catalog}, a suite of cluster simulations with present-day properties representative of the globular clusters (GCs) observed in the Milky Way. We show how the initial properties of the host cluster set the typical orbital parameters and formation times of the formed triples. We find that a cluster typically assembles hundreds of triples with at least one black hole (BH) in the inner binary, while only clusters with sufficiently small virial radii are efficient in producing triples with no BHs, as a result of the BH-burning process. We show that a typical GC is expected to host tens of triples with at least one luminous component at present day. We discuss how the Lidov-Kozai mechanism can drive the inner binary of the formed triples to high eccentricities, whenever it takes place before the triple is dynamically reprocessed by encountering another cluster member. Some of these systems can reach sufficiently large eccentricities to form a variety of transients and sources, such as blue stragglers, X-ray binaries, Type Ia Supernovae, Thorne-Zytkow objects, and LIGO/Virgo sources. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 28 pages, 12 figures, 2 tables, accepted by ApJ

Populating the upper black hole mass gap through stellar collisions in young star clusters

Authors: Kyle Kremer, Mario Spera, Devin Becker, Sourav Chatterjee, Ugo N. Di Carlo, Giacomo Fragione, Carl L. Rodriguez, Claire S. Ye, Frederic A. Rasio

Abstract: Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of B… ▽ More Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through collisions of young massive stars in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. This has clear relevance for the formation of very massive BH binaries as recently detected by LIGO/Virgo in GW190521. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients. △ Less

Submitted 15 September, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 25 pages, 6 figures, accepted for publication in ApJ. Comments welcome

GW190412 as a Third-Generation Black Hole Merger from a Super Star Cluster

Authors: Carl L. Rodriguez, Kyle Kremer, Michael Y. Grudić, Zachary Hafen, Sourav Chatterjee, Giacomo Fragione, Astrid Lamberts, Miguel A. S. Martinez, Frederic A. Rasio, Newlin Weatherford, Claire S. Ye

Abstract: We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary… ▽ More We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics. △ Less

Submitted 27 May, 2020; v1 submitted 8 May, 2020; originally announced May 2020.

Comments: 9 pages, 4 figures, matches version accepted by ApJL

Merging Black Holes in the Low-mass and High-mass Gaps from 2+2 Quadruple Systems

Authors: Giacomo Fragione, Abraham Loeb, Frederic A. Rasio

Abstract: The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron star… ▽ More The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function. △ Less

Submitted 6 May, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 8 pages, 1 Table, 4 Figures, accepted by ApJ Lett

Gravitational-wave captures by intermediate-mass black holes in galactic nuclei

Authors: Giacomo Fragione, Abraham Loeb, Kyle Kremer, Frederic A. Rasio

Abstract: Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO… ▽ More Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin. △ Less

Submitted 18 May, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: 10 pages, 5 figures, accepted by ApJ

A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

Authors: Newlin C. Weatherford, Sourav Chatterjee, Kyle Kremer, Frederic A. Rasio

Abstract: Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars.… ▽ More Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($Δ$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $Δ$ between populations of main sequence stars, using correlations found between $Δ$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $Δ$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$. △ Less

Submitted 23 June, 2020; v1 submitted 20 November, 2019; originally announced November 2019.

Comments: 25 pages, 8 figures, 6 tables; accepted for publication to ApJ

COSMIC Variance in Binary Population Synthesis

Authors: Katelyn Breivik, Scott Coughlin, Michael Zevin, Carl L. Rodriguez, Kyle Kremer, Claire S. Ye, Jeff J. Andrews, Michael Kurkowski, Matthew C. Digman, Shane L. Larson, Frederic A. Rasio

Abstract: The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observat… ▽ More The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA. △ Less

Submitted 10 September, 2020; v1 submitted 3 November, 2019; originally announced November 2019.

Comments: reflects the version submitted to ApJ; 17 pages, 2 Tables, 4 Figures; corrects typo in SNR calculation

Modeling Dense Star Clusters in the Milky Way and Beyond with the $\texttt{CMC}$ Cluster Catalog

Authors: Kyle Kremer, Claire S. Ye, Nicholas Z. Rui, Newlin C. Weatherford, Sourav Chatterjee, Giacomo Fragione, Carl L. Rodriguez, Mario Spera, Frederic A. Rasio

Abstract: We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models… ▽ More We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters. △ Less

Submitted 8 March, 2020; v1 submitted 31 October, 2019; originally announced November 2019.

Comments: 50 pages, 15 figures, 10 tables. Accepted for publication in ApJS

On the Rate of Neutron Star Binary Mergers from Globular Clusters

Authors: Claire S. Ye, Wen-fai Fong, Kyle Kremer, Carl L. Rodriguez, Sourav Chatterjee, Giacomo Fragione, Frederic A. Rasio

Abstract: The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations… ▽ More The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates. △ Less

Submitted 6 December, 2019; v1 submitted 23 October, 2019; originally announced October 2019.

Comments: 17 pages, 4 figures, 2 tables, accepted by ApJL

Journal ref: The Astrophysical Journal Letters, Volume 888, Issue 1, article id. L10, 13 pp. (2020)

Probing the Survival of Planetary Systems in Globular Clusters with Tidal Disruption Events

Authors: Kyle Kremer, Daniel J. D'Orazio, Johan Samsing, Sourav Chatterjee, Frederic A. Rasio

Abstract: Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clu… ▽ More Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clusters that may no longer be bound to a host star. Here, we use the cluster Monte Carlo code CMC to explore the evolution of planetary systems in dense star clusters. Depending on a number of initial conditions, we show that $10-50\%$ of primordial planetary systems are broken through dynamical encounters over a cluster's full lifetime, populating clusters with "free-floating" planets. Furthermore, a large number ($30-80\%$) of planets are ejected from their host cluster through strong dynamical encounters and/or tidal loss. Additionally, we show that planets naturally mix with stellar-mass black holes (BHs) in the central regions of their host cluster. As a consequence, up to a few hundreds of planets will be tidally disrupted through close passages of BHs. We show these BH-planet tidal disruption events (TDEs) occur in clusters at a rate of up to $10^{-5}\,\rm{yr}^{-1}$ in a Milky Way-type galaxy. We predict BH-planet TDEs should be detected by upcoming transient surveys such as LSST at a rate of roughly a few events per year. The observed rate of BH-planet TDEs would place new constraints upon the formation and survival of both planetary systems and BHs in dense star clusters. Additionally, depending on various assumptions including the initial number of planets and their orbital properties, we predict that typical globular clusters may contain a few dynamically-formed NS-planet systems at present as well as up to roughly 100 dynamically-formed WD-planet systems. △ Less

Submitted 22 August, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: 18 pages, 8 figures. Submitted to ApJ. Comments welcome

The Role of "Black Hole Burning" in the Evolution of Dense Star Clusters

Authors: Kyle Kremer, Claire S. Ye, Sourav Chatterjee, Carl L. Rodriguez, Frederic A. Rasio

Abstract: As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in… ▽ More As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning." △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: To appear in the Star Clusters: From the Milky Way to the Early Universe IAU Proceedings. Comments welcome

Black Holes: The Next Generation -- Repeated Mergers in Dense Star Clusters and their Gravitational-Wave Properties

Authors: Carl L. Rodriguez, Michael Zevin, Pau Amaro-Seoane, Sourav Chatterjee, Kyle Kremer, Frederic A. Rasio, Claire S. Ye

Abstract: When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the p… ▽ More When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the production of binary black hole mergers where at least one component of the binary was forged in a previous merger. We create four hypothetical universes where every black hole born in the collapse of a massive star has a dimensionless Kerr spin parameter, $χ_{\rm birth}$, of 0.0, 0.1, 0.2, or 0.5. We show that if all stellar-born black holes are non-spinning ($χ_{\rm birth} = 0.0$), then more than 10% of merging binary black holes from clusters have components formed from previous mergers, accounting for more than 20% of the mergers from globular clusters detectable by LIGO/Virgo. Furthermore, nearly 7% of detectable mergers would have a component with a mass $\gtrsim 55M_{\odot}$, placing them clearly in the mass "gap" region where black holes cannot form from isolated collapsing stars due to the pulsational-pair instability mechanism. On the other hand, if black holes are born spinning, then the contribution from these second-generation mergers decreases, making up as little as 1% of all detections from globular clusters when $χ_{\rm birth} = 0.5$. We make quantitative predictions for the detected masses, mass ratios, and spin properties of first- and second-generation mergers from dense star clusters, and show how these distributions are highly sensitive to the birth spins of black holes. △ Less

Submitted 4 August, 2019; v1 submitted 24 June, 2019; originally announced June 2019.

Comments: 15 pages, 8 figures, matches version accepted by PRD

Journal ref: Phys. Rev. D 100, 043027 (2019)

Millisecond Pulsars and Black Holes in Globular Clusters

Authors: Claire S. Ye, Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Frederic A. Rasio

Abstract: Over a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic $N$-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (i… ▽ More Over a hundred millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic $N$-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses $\approx 2 \times 10^5\,M_\odot$ can produce up to $10-20$ MSPs, while a very massive GC model with mass $\approx 10^6\,M_\odot$ can produce close to $100$. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated towards the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS--NS and NS--BH binaries in GCs, although the predicted numbers are very small. △ Less

Submitted 18 April, 2019; v1 submitted 15 February, 2019; originally announced February 2019.

Comments: 14 pages, 6 figures, Accepted for publication in ApJ

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

Post-Newtonian Dynamics in Dense Star Clusters: Binary Black Holes in the LISA Band

Authors: Kyle Kremer, Carl L. Rodriguez, Pau Amaro-Seoane, Katelyn Breivik, Sourav Chatterjee, Michael L. Katz, Shane L. Larson, Frederic A. Rasio, Johan Samsing, Claire S. Ye, Michael Zevin

Abstract: The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all B… ▽ More The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission. △ Less

Submitted 26 February, 2019; v1 submitted 28 November, 2018; originally announced November 2018.

Comments: Accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 99, 063003 (2019)

Post-Newtonian Dynamics in Dense Star Clusters: Formation, Masses, and Merger Rates of Highly-Eccentric Black Hole Binaries

Authors: Carl L. Rodriguez, Pau Amaro-Seoane, Sourav Chatterjee, Kyle Kremer, Frederic A. Rasio, Johan Samsing, Claire S. Ye, Michael Zevin

Abstract: Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission betw… ▽ More Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc^-3 yr^-1 to the binary merger rate in the local universe, increasing to ~10 Gpc^-3 yr^-1 at z~3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture. △ Less

Submitted 15 November, 2018; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: 16 Pages, 10 Figures. Accepted by Phys Rev D

How initial size governs core collapse in globular clusters

Authors: Kyle Kremer, Sourav Chatterjee, Claire S. Ye, Carl L. Rodriguez, Frederic A. Rasio

Abstract: Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called "core-collapsed" and "non-core-collapsed" clusters. Here, we use our Hénon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions i… ▽ More Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called "core-collapsed" and "non-core-collapsed" clusters. Here, we use our Hénon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, $r_v$) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe ($r_v \approx 0.5-5$ pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify "best-fit" models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of $\sim 50-100$ stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main sequence stars, white dwarfs, and sub-subgiants. △ Less

Submitted 3 December, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: v2. Accepted for publication in ApJ. Comments welcome

Survival of non-coplanar, closely-packed planetary systems after a close encounter

Authors: David R. Rice, Frederic A. Rasio, Jason H. Steffen

Abstract: Planetary systems with more than two bodies will experience orbital crossings at a time related to the initial orbital separations of the planets. After a crossing, the system enters a period of chaotic evolution ending in the reshaping of the system's architecture via planetary collisions or ejections. We carry out N-body integrations on a large number of systems with equally-spaced planets (in u… ▽ More Planetary systems with more than two bodies will experience orbital crossings at a time related to the initial orbital separations of the planets. After a crossing, the system enters a period of chaotic evolution ending in the reshaping of the system's architecture via planetary collisions or ejections. We carry out N-body integrations on a large number of systems with equally-spaced planets (in units of the Hill radius) to determine the distribution of instability times for a given planet separation. We investigate both the time to the initiation of instability through a close encounter and the time to a planet-planet collision. We find that a significant portion of systems with non-zero mutual inclinations survive after a close encounter and do not promptly experience a planet-planet collision. Systems with significant inclinations can continue to evolve for over 1,000 times longer than the encounter time. The fraction of long lived systems is dependent on the absolute system scale and the initial inclination of the planets. These results have implications to the assumed stability of observed planetary systems. △ Less

Submitted 31 August, 2018; v1 submitted 19 July, 2018; originally announced July 2018.

Comments: 9 pages, 10 figures, 1 table, Accepted for publication in MNRAS

Journal ref: MNRAS, 481, 2205-2212 (2018)

MeerTime - the MeerKAT Key Science Program on Pulsar Timing

Authors: M. Bailes, E. Barr, N. D. R. Bhat, J. Brink, S. Buchner, M. Burgay, F. Camilo, D. J. Champion, J. Hessels, G. H. Janssen, A. Jameson, S. Johnston, A. Karastergiou, R. Karuppusamy, V. Kaspi, M. J. Keith, M. Kramer, M. A. McLaughlin, K. Moodley, S. Oslowski, A. Possenti, S. M. Ransom, F. A. Rasio, J. Sievers, M. Serylak , et al. (6 additional authors not shown)

Abstract: The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to fou… ▽ More The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results. △ Less

Submitted 18 March, 2018; originally announced March 2018.

Comments: 15 pages, MeerKAT Science: On the Pathway to the SKA, 25-27 May, 2016, Stellenbosch, South Africa, available at: https://pos.sissa.it/277/011/pdf

Journal ref: 2018, Proceedings of Science, Volume 277, 011, https://pos.sissa.it/277/011/pdf

How black holes shape globular clusters: Modeling NGC 3201

Authors: Kyle Kremer, Claire S. Ye, Sourav Chatterjee, Carl L. Rodriguez, Frederic A. Rasio

Abstract: Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains… ▽ More Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains $\gtrsim 200$ stellar-mass BHs. Furthermore, we explore the dynamical formation of main sequence-BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue-straggler-BH binaries unique to core-collapsed clusters, which otherwise retain few BHs. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: Accepted for publication in ApJL. Comments welcome

LISA sources in Milky Way globular clusters

Authors: Kyle Kremer, Sourav Chatterjee, Katelyn Breivik, Carl L. Rodriguez, Shane L. Larson, Frederic A. Rasio

Abstract: We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These… ▽ More We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These detectable sources contain all combinations of black hole (BH), neutron star, and white dwarf components. We predict $\sim7$ of these sources will be BH-BH binaries. Furthermore, we show that some of these BH-BH binaries can have signal-to-noise ratios large enough to be detectable at the distance of the Andromeda galaxy or even the Virgo cluster. △ Less

Submitted 29 May, 2018; v1 submitted 15 February, 2018; originally announced February 2018.

Comments: 7 pages, 4 figures, 2 appendices. As published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 120, 191103 (2018)

Low-mass X-ray binaries ejected from globular clusters

Authors: Kyle Kremer, Sourav Chatterjee, Carl L. Rodriguez, Frederic A. Rasio

Abstract: We explore the population of mass-transferring binaries ejected from globular clusters (GCs) with both black hole (BH) and neutron star (NS) accretors. We use a set of 137 fully evolved globular cluster models which span a large range in cluster properties and, overall, match very well the properties of old GCs observed in the Milky Way. We identify all binaries ejected from our set of models that… ▽ More We explore the population of mass-transferring binaries ejected from globular clusters (GCs) with both black hole (BH) and neutron star (NS) accretors. We use a set of 137 fully evolved globular cluster models which span a large range in cluster properties and, overall, match very well the properties of old GCs observed in the Milky Way. We identify all binaries ejected from our set of models that eventually undergo mass-transfer. These binaries are ejected from their host clusters over a wide range of ejection times and include white dwarf, giant, and main sequence donors. We calculate the orbits of these ejected systems in the Galactic potential to determine their present-day positions in the Galaxy and compare to the distribution of observed low-mass X-ray binaries (XRBs) in the Milky Way. We estimate $\sim 300$ mass-transferring NS binaries and $\sim 180$ mass-transferring BH binaries may currently be present in the Milky Way that originated from within GCs. Of these, we estimate, based on mass-transfer rates and duty cycles at the present time, at most a few would be observable as BH--XRBs and NS--XRBs at the present day. Based on our results, XRBs that originated from GCs are unlikely to contribute significantly to the total population of low-mass XRBs in the Galactic field. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: Submitted to ApJ. Comments welcome

Post-Newtonian Dynamics in Dense Star Clusters: Highly-Eccentric, Highly-Spinning, and Repeated Binary Black Hole Mergers

Authors: Carl L. Rodriguez, Pau Amaro-Seoane, Sourav Chatterjee, Frederic A. Rasio

Abstract: We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities… ▽ More We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae. △ Less

Submitted 14 March, 2018; v1 submitted 13 December, 2017; originally announced December 2017.

Comments: 9 pages, 3 figures, 2 appendices. To appear in Physical Review Letters

Journal ref: Phys. Rev. D 98, 123005 (2018)