Search

The existence of compact stellar remnants in the mass range $2-5\,M_{\odot}$ has long been debated. This so-called lower mass gap was initially suggested by the lack of low-mass X-ray binary observations with accretors in this mass range, but it has recently been called into question following newer observations, including a lower-mass-gap candidate with a millisecond pulsar companion in the dense globular cluster NGC~1851. Here we model NGC~1851 with a grid of similar dense star clusters utilizing the state-of-the-art Monte Carlo $N$-body code \texttt{CMC}, and we specifically study the formation of lower-mass-gap black holes. We demonstrate that both massive star evolution and dynamical interactions can contribute to forming lower-mass-gap black holes. In general, the collapse of massive remnants formed through mergers of neutron stars or massive white dwarfs produces the largest number of lower-mass-gap black holes among all formation channels. However, in more massive clusters, supernova core collapse can contribute comparable numbers. Our NGC~1851-like models can reproduce millisecond pulsar -- lower-mass-gap black hole binaries similar to the observed system. Additionally, the lower-mass-gap black holes can also become components of dynamically assembled binaries, and some will be in merging black hole--neutron star systems similar to the recently detected gravitational wave source GW230529. However, the corresponding merger rate is probably $\lesssim 1~{\rm Gpc^{-3}\,yr^{-1}}$., Comment: 18 pages, 7 figures, 3 tables. Submitted to ApJ. Comments welcome

We report the discovery of a high velocity, very low-mass star or brown dwarf whose kinematics suggest it is unbound to the Milky Way. CWISE J124909.08+362116.0 was identified by citizen scientists in the Backyard Worlds: Planet 9 program as a high proper motion ($\mu$ $=$ 0''9/yr) faint red source. Moderate resolution spectroscopy with Keck/NIRES reveals it to be a metal-poor early L subdwarf with a large radial velocity ($-$103$\pm$10 km/s), and its estimated distance of 125$\pm$8 pc yields a speed of 456$\pm$27 km/s in the Galactic rest frame, near the local escape velocity for the Milky Way. We explore several potential scenarios for the origin of this source, including ejection from the Galactic center $\gtrsim$3 Gyr in the past, survival as the mass donor companion to an exploded white dwarf. acceleration through a three-body interaction with a black hole binary in a globular cluster, and accretion from a Milky Way satellite system. CWISE J1249+3621 is the first hypervelocity very low mass star or brown dwarf to be found, and the nearest of all such systems. It may represent a broader population of very high velocity, low-mass objects that have undergone extreme accelerations., Comment: 17 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

The origins and mergers of supermassive black holes (BHs) remain a mystery. We describe a scenario from a novel multi-physics simulation featuring rapid ($\lesssim 1\,$Myr) hyper-Eddington gas capture by a $\sim 1000\,{\rm M}_{\odot}$ ``seed'' BH up to supermassive ($\gtrsim 10^{6}\,M_{\odot}$) masses, in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient and most of the gas turns into stars in star clusters which rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be ``captured'' by merging sub-clusters and migrate to the center in $\sim1$ free-fall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with $\sim 0.1$\,pc separation. The centrally-concentrated stellar density profile (akin to a ``proto-bulge'') allows the cluster as a whole to capture and retain gas and build up a large (pc-scale) circum-binary accretion disk with gas coherently funnelled to the central BH (even when the BH radius of influence is small). The disk is ``hyper-magnetized'' and ``flux-frozen'': dominated by a toroidal magnetic field with plasma $\beta \sim 10^{-3}$, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates $\gtrsim 1\,\rm M_\odot yr^{-1}$, which also drive the two BHs to nearly-equal masses. The late-stage system appears remarkably similar to recently-observed high-redshift ``little red dots.'' This scenario can provide an explanation for rapid SMBH formation, growth and mergers in high-redshift galaxies., Comment: 9 pages, 5 figures, final version matching the publication

The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kilo-parsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density ($\sim 10^4\,{\rm {\rm M_\odot}\, pc}^{-2}$), facilitating ``hyper-Eddington'' accretion via efficient feeding by dense clumps which are driven by turbulence and stellar feedback. This article investigates the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. We run a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the sub-grid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Utilizing radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities as high as $10^{41}$ -- $10^{44}\,\rm erg/s$, depending on the GMC properties and specific feedback model assumed. We find the maximum possible mass growth of seed BHs ($\Delta M_{\rm BH}^{\rm max}$) is regulated by the momentum deposition rate from BH feedback, $\dot{p}_{\rm feedback}/(\dot{M}_{\rm BH} c)$, which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of $\sim 10^4\,\rm M_\odot$ intermediate-massive BHs (IMBHs) from stellar-mass BHs within $\sim \rm Myr$. Furthermore, we examine the impacts of sub-grid accretion models and how BH feedback may influence star formation within these cloud complexes., Comment: Comments welcome

Context. The star-black hole (S-BH) binary known as Gaia BH3, discovered by the Gaia Collaboration is chemically and kinematically associated with the metal-poor ED-2 stream in the Milky Way halo. Aims. We explore the possibility that Gaia BH3 was assembled dynamically in the progenitor globular cluster (GC) of the ED-2 stream. Methods. We used a public suite of star-by-star dynamical Monte Carlo models to identify S-BH binaries in GCs with different initial masses and (half-mass) radii. Results. We show that a likely progenitor of the ED-2 stream was a relatively low-mass ($\lesssim10^5M_\odot$) GC with an initial half-mass radius of ~4 pc. Such a GC can dynamically retain a large fraction of its BH population and dissolve on the orbit of ED-2. From the suite of models we find that GCs produce ~ 3 - 30 S-BH binaries, approximately independent of initial GC mass and inversely correlated with initial cluster radius. Scaling the results to the Milky Way GC population, we find that ~75% of the S-BH binaries formed in GCs are ejected from their host GC, all in the early phases of evolution ($\lesssim1$ Gyr); these are expected to no longer be close to streams. The ~25% of S-BH binaries retained until dissolution are expected to form part of streams, such that for an initial mass of the progenitor of ED-2 of a few $10^4M_\odot$, we expect ~2-3 S-BH to end up in the stream. GC models with metallicities similar to Gaia BH3 ($\lesssim1\%$ solar) include S-BH binaries with similar BH masses ($\gtrsim30M_\odot$), orbital periods, and eccentricities. Conclusions. We predict the Galactic halo contains of order $10^5$ S-BH binaries that formed dynamically in GCs, a fraction of which may readily be detected in Gaia DR4. The detection of these sources provides valuable tests of BH dynamics in clusters and the possible role in formation of gravitational wave sources., Comment: 5 pages, 3 figures, accepted in A&A

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., Comment: Accepted for ApJ. Comments welcome. 19 pages, 10 figures

Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and wander, ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are swept up into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short ($\ll$ yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at $\lesssim 10$ au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models., Comment: 39 pages, 23 figures, replaced with version accepted to the Open Journal of Astrophysics. Additional images and movies from the simulations available at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories -- but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black hole in gas-free dynamical environments, the presence of anti-aligned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger's dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin ($\chi_{\rm eff}$) and in-plane spin ($\chi_{\rm p}$) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers' $\chi_{\rm eff}$ and $\chi_{rm p}$. The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network's third observing period and projections for design sensitivity. Using a "likelihood-ratio"-based statistic, we find that $\sim2\%$ of the recovered population by the current gravitational-wave detector network has a statistically significant $\chi_{\rm p}$ measurement, whereas no $\chi_{\rm eff}$ measurement was capable of confidently determining a system to be anti-aligned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession through $\chi_{\rm p}$ is a more detectable signature of a hierarchical mergers and dynamical formation than anti-aligned spins., Comment: 10 pages, 2 figures. If you would like to give us comments, just go ahead now

Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $\alpha$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $\beta \ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates., Comment: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $\beta \sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk., Comment: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

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 $\Delta t$ of PEs, we find new behavior related to chaotic scattering in the three-body problem; the $\Delta 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 $\Delta t$ in an evolving cluster, finding that $\Delta 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 $\Delta 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., Comment: 31 pages, 19 Figures (+ linked video), 1 Table, accepted to ApJ

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., Comment: 19 pages, 12 figures, 4 tables. Published at ApJ

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe., Comment: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core-collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper-motion data from Gaia EDR3 and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions ($\beta\approx -0.4$$\pm$$0.1$) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm300 \,$M$_{\odot}$, but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), or a dark population of stellar remnants of extent $\approx 0.016\,\rm pc \simeq 3300\,AU$. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended ($\sim 0.034\, \rm{pc} \approx 7000\,\rm AU$). We use Monte Carlo $N$-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e., IMBH vs. remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core., Comment: 19 page, 15 figures, 3 tables. Accepted for publication in MNRAS

The detections of four apparently young radio pulsars in the Milky Way globular clusters are difficult to reconcile with standard neutron star formation scenarios associated with massive star evolution. Here we discuss formation of these young pulsars through white dwarf mergers in dynamically-old clusters that have undergone core collapse. Based on observed properties of magnetic white dwarfs, we argue neutron stars formed via white dwarf merger are born with spin periods of roughly $10-100\,$ms and magnetic fields of roughly $10^{11}-10^{13}\,$G. As these neutron stars spin down via magnetic dipole radiation, they naturally reproduce the four observed young pulsars in the Milky Way clusters. Rates inferred from $N$-body cluster simulations as well as the binarity, host cluster properties, and cluster offsets observed for these young pulsars hint further at a white dwarf merger origin. These young pulsars may be descendants of neutron stars capable of powering fast radio bursts analogous to the bursts observed recently in a globular cluster in M81., Comment: 5 pages, 2 figures, 1 table. Accepted for publication in MNRAS

Tidal disruptions of stars by stellar-mass black holes are expected to occur frequently in dense star clusters. Building upon previous studies that performed hydrodynamic simulations of these encounters, we explore the formation and long-term evolution of the thick, super-Eddington accretion disks formed. We build a disk model that includes fallback of material from the tidal disruption, accretion onto the black hole, and disk mass losses through winds launched in association with the super-Eddington flow. We demonstrate that bright transients are expected when radiation from the central engine powered by accretion onto the black hole is reprocessed at large radii by the optically-thick disk wind. By combining hydrodynamic simulations of these disruption events with our disk+wind model, we compute light curves of these wind-reprocessed transients for a wide range of stellar masses and encounter penetration depths. We find typical peak bolometric luminosities of roughly $10^{41}-10^{44}\,$erg/s (depending mostly on accretion physics parameters) and temperatures of roughly $10^5-10^6\,$K, suggesting peak emission in the ultraviolet/blue bands. We predict all-sky surveys such as the Vera Rubin Observatory and ULTRASAT will detect up to thousands of these events per year in dense star clusters out to distances of several Gpc., Comment: 16 Pages, 13 figures, 2 tables. Accepted for publication in MNRAS

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., Comment: 28 pages, 6 figures, 4 tables, accepted to ApJ

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

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., Comment: 19 pages, 8 figures, 2 tables. Submitted to ApJ

The recent detection of a repeating fast radio burst (FRB) in an old globular cluster in M81 challenges traditional FRB formation mechanisms based on magnetic activity in young neutron stars formed recently in core-collapse supernovae. Furthermore, the detection of this repeater in such a nearby galaxy implies a high local universe rate of similar events in globular clusters. Building off the properties inferred from the M81 FRB, we predict the number of FRB sources in nearby ($d \lesssim 20\,$Mpc) galaxies with large globular cluster systems known. Incorporating the uncertain burst energy distribution, we estimate the rate of bursts detectable in these galaxies by radio instruments such as FAST and MeerKat. Of all local galaxies, we find M87 is the best candidate for FRB detections. We predict M87's globular cluster system contains $\mathcal{O}(10)$ FRB sources at present and that a dedicated radio survey (by FAST or MeerKat) of $\mathcal{O}(10)\,$hr has a $90\%$ probability of detecting a globular cluster FRB in M87. The detection of even a handful of additional globular cluster FRBs would provide invaluable constraints on FRB mechanisms and population properties. Previous studies have demonstrated young neutron stars formed following collapse of dynamically-formed massive white dwarf binary mergers may provide the most natural mechanism for these bursts. We explore the white dwarf merger scenario using a suite of $N$-body cluster models, focusing in particular on such mergers in M87 clusters. We describe a number of outstanding features of this scenario that in principle may be testable with an ensemble of observed FRBs in nearby globular clusters., Comment: 22 pages, 9 figures, 2 tables. Accepted for publication in ApJ

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$., Comment: Accepted for ApJ. Comments welcome. 16 pages and 10 figures

Formation of supermassive black holes (BHs) remains a theoretical challenge. In many models, especially beginning from stellar relic "seeds," this requires sustained super-Eddington accretion. While studies have shown BHs can violate the Eddington limit on accretion disk scales given sufficient "fueling" from larger scales, what remains unclear is whether or not BHs can actually capture sufficient gas from their surrounding ISM. We explore this in a suite of multi-physics high-resolution simulations of BH growth in magnetized, star-forming dense gas complexes including dynamical stellar feedback from radiation, stellar mass-loss, and supernovae, exploring populations of seeds with masses $\sim 1-10^{4}\,M_{\odot}$. In this initial study, we neglect feedback from the BHs: so this sets a strong upper limit to the accretion rates seeds can sustain. We show that stellar feedback plays a key role. Complexes with gravitational pressure/surface density below $\sim 10^{3}\,M_{\odot}\,{\rm pc^{-2}}$ are disrupted with low star formation efficiencies so provide poor environments for BH growth. But in denser cloud complexes, early stellar feedback does not rapidly destroy the clouds but does generate strong shocks and dense clumps, allowing $\sim 1\%$ of randomly-initialized seeds to encounter a dense clump with low relative velocity and produce runaway, hyper-Eddington accretion (growing by orders of magnitude). Remarkably, mass growth under these conditions is almost independent of initial BH mass, allowing rapid IMBH formation even for stellar-mass seeds. This defines a necessary (but perhaps not sufficient) set of criteria for runaway BH growth: we provide analytic estimates for the probability of runaway growth under different ISM conditions., Comment: Final version matching the publication

We report the discovery of a new 5.78 ms-period millisecond pulsar (MSP), PSR J1740-5340B (NGC 6397B), in an eclipsing binary system discovered with the Parkes radio telescope (now also known as Murriyang), Australia, and confirmed with the MeerKAT radio telescope in South Africa. The measured orbital period, 1.97 days, is the longest among all eclipsing binaries in globular clusters (GCs) and consistent with that of the coincident X-ray source U18, previously suggested to be a 'hidden MSP'. Our XMM-Newton observations during NGC 6397B's radio quiescent epochs detected no X-ray flares. NGC 6397B is either a transitional MSP or an eclipsing binary in its initial stage of mass transfer after the companion star left the main sequence. The discovery of NGC 6397B potentially reveals a subgroup of extremely faint and heavily obscured binary pulsars, thus providing a plausible explanation to the apparent dearth of binary neutron stars in core-collapsed GCs as well as a critical constraint on the evolution of GCs.

Compact-object binaries including a white dwarf component are unique among gravitational-wave sources because their evolution is governed not just by general relativity and tides, but also by mass transfer. While the black hole and neutron star binaries observed with ground-based gravitational-wave detectors are driven to inspiral due to the emission of gravitational radiation--manifesting as a "chirp-like" gravitational-wave signal--the astrophysical processes at work in double white dwarf (DWD) systems can cause the inspiral to stall and even reverse into an outspiral. The dynamics of the DWD outspiral thus encode information about tides, which tell us about the behaviour of electron-degenerate matter. We carry out a population study to determine the effect of the strength of tides on the distributions of the DWD binary parameters that the Laser Interferometer Space Antenna (LISA) will be able to constrain. We find that the strength of tidal coupling parameterized via the tidal synchronization timescale at the onset of mass transfer affects the distribution of gravitational-wave frequencies and frequency derivatives for detectably mass-transferring DWD systems. Using a hierarchical Bayesian framework informed by binary population synthesis simulations, we demonstrate how this parameter can be inferred using LISA observations. By measuring the population properties of DWDs, LISA will be able to probe the behavior of electron-degenerate matter., Comment: version accepted in ApJ

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., Comment: 14 pages, 4 figures, 2 tables. Accepted for publication in ApJL

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

We analyse Gaia EDR3 and re-calibrated HST proper motion data from the core-collapsed and non core-collapsed globular clusters NGC 6397 and NGC 3201, respectively, with the Bayesian mass-orbit modelling code MAMPOSSt-PM. We use Bayesian evidence and realistic mock data sets constructed with AGAMA to select between different mass models. In both clusters, the velocities are consistent with isotropy within the extent of our data. We robustly detect a dark central mass (DCM) of roughly 1000 solar masses in both clusters. Our MAMPOSSt-PM fits strongly prefer an extended DCM in NGC 6397, while only presenting a mild preference for it in NGC 3201, with respective sizes of a roughly one and a few per cent of the cluster effective radius. We explore the astrophysics behind our results with the CMC Monte Carlo N-body code, whose snapshots best matching the phase space observations lead to similar values for the mass and size of the DCM. The internal kinematics are thus consistent with a population of hundreds of massive white dwarfs in NGC 6397, and roughly 100 segregated stellar-mass black holes in NGC 3201, as previously found with CMC. Such analyses confirm the accuracy of both mass-orbit modelling and Monte Carlo N-body techniques, which together provide more robust predictions on the DCM of globular clusters (core-collapsed or not). This opens possibilities to understand a vast range of interesting astrophysical phenomena in clusters, such as fast radio bursts, compact object mergers, and gravitational waves., Comment: 28 pages, 16 figures, 4 tables, accepted for publication in MNRAS

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., Comment: 23 pages, 11 figures, 3 tables. Accepted for publication in ApJ

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., Comment: v2: 17 pages, 5 figures

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., Comment: 23 pages, 15 figures, 4 tables, published at ApJ

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduce the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass $ 2\times10^{-14}\lesssim \mu/\textrm{eV} \lesssim2\times10^{-13}$ would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to $\approx60\%$ more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate mass black hole mergers to $\approx0.08$\,Gpc$^{-3}$\,yr$^{-1}$ in the local Universe. These results imply that a select range of ultralight boson mass can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario., Comment: 20 pages, 7 figures

The repeating fast radio burst (FRB) localized to a globular cluster in M81 challenges our understanding of FRB models. In this Letter, we explore dynamical formation scenarios for objects in old globular clusters that may plausibly power FRBs. Using N-body simulations, we demonstrate that young neutron stars may form in globular clusters at a rate of up to $\sim50\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ through a combination of binary white dwarf mergers, white dwarf--neutron star mergers, binary neutron star mergers, and accretion induced collapse of massive white dwarfs in binary systems. We consider two FRB emission mechanisms: First, we show that a magnetically-powered source (e.g., a magnetar with field strength $\gtrsim10^{14}\,$G) is viable for radio emission efficiencies $\gtrsim10^{-4}$. This would require magnetic activity lifetimes longer than the associated spin-down timescales and longer than empirically-constrained lifetimes of Galactic magnetars. Alternatively, if these dynamical formation channels produce young rotation-powered neutron stars with spin periods of $\sim10\,$ms and magnetic fields of $\sim10^{11}\,$G (corresponding to spin-down lifetimes of $\gtrsim10^5\,$yr), the inferred event rate and energetics can be reasonably reproduced for order unity duty cycles. Additionally, we show that recycled millisecond pulsars or low-mass X-ray binaries similar to those well-observed in Galactic globular clusters may also be plausible channels, but only if their duty cycle for producing bursts similar to the M81 FRB is small., Comment: 11 pages, 2 tables, 1 figure. Accepted for publication in ApJ Letters. Comments welcome

Orbital eccentricity is one of the most robust discriminators for distinguishing between dynamical and isolated formation scenarios of binary black hole mergers using gravitational-wave observatories such as LIGO and Virgo. Using state-of-the-art cluster models, we show how selection effects impact the detectable distribution of eccentric mergers from clusters. We show that the observation (or lack thereof) of eccentric binary black hole mergers can significantly constrain the fraction of detectable systems that originate from dynamical environments, such as dense star clusters. After roughly 150 observations, observing no eccentric binary signals would indicate that clusters cannot make up the majority of the merging binary black hole population in the local universe (95% credibility). However, if dense star clusters dominate the rate of eccentric mergers and a single system is confirmed to be measurably eccentric in the first and second gravitational-wave transient catalogs, clusters must account for at least 14% of detectable binary black hole mergers. The constraints on the fraction of detectable systems from dense star clusters become significantly tighter as the number of eccentric observations grows and will be constrained to within 0.5 dex once 10 eccentric binary black holes are observed., Comment: 15 pages (9 main text + 6 appendices/references), 4 figures, ApJ Letters in press

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\'enon'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)., Comment: Code is available at https://clustermontecarlo.github.io/ 25 pages, 8 Figures, Matches version accepted by ApJS

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., Comment: Submitted to ApJ, 26 pages, 12 figures, 3 tables. Comments welcome

Supersonically Induced Gas Objects (SIGOs), are structures with little to no dark matter component predicted to exist in regions of the Universe with large relative velocities between baryons and dark matter at the time of recombination. They have been suggested to be the progenitors of present-day globular clusters. Using simulations, SIGOs have been studied on small scales (around 2 Mpc), where these relative velocities are coherent. However, it is challenging to study SIGOs using simulations on large scales due to the varying relative velocities at scales larger than a few Mpc. Here, we study SIGO abundances semi-analytically: using perturbation theory, we predict the number density of SIGOs analytically, and compare these results to small-box numerical simulations. We use the agreement between the numerical and analytic calculations to extrapolate the large-scale variation of SIGO abundances over different stream velocities. As a result, we predict similar large-scale variations of objects with high gas densities before reionization that could possibly be observed by JWST. If indeed SIGOs are progenitors of globular clusters, then we expect a similar variation of globular cluster abundances over large scales. Significantly, we find that the expected number density of SIGOs is consistent with observed globular cluster number densities. As a proof-of-concept, and because globular clusters were proposed to be natural formation sites for gravitational wave sources from binary black hole (BBH) mergers, we show that SIGOs should imprint an anisotropy on the gravitational wave signal on the sky, consistent with SIGOs' distribution., Comment: 17 pages, 8 figures, Accepted to ApJ

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., Comment: 3 pages, 1 figure, to submit to AAS journals

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}., Comment: 22 pages, 11 figures, 5 tables; accepted to ApJ

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., Comment: 4 pages, one figure, matches version in RNAAS

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., Comment: 8 pages, 4 figures, 1 table; accepted for publication to ApJL

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., Comment: 10 pages, 2 figures, 1 table. Accepted for publication in ApJ Letters. Comments welcome

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., Comment: 17 pages, 5 figures, 3 tables. Accepted for publication in ApJ

Globular clusters are considered to be likely breeding grounds for compact binary mergers. In this paper, we demonstrate how the gravitational-wave signals produced by compact object mergers can act as tracers of globular cluster formation and evolution. Globular cluster formation is a long-standing mystery in astrophysics, with multiple competing theories describing when and how globular clusters formed. The limited sensitivity of electromagnetic telescopes inhibits our ability to directly observe globular cluster formation. However, with future audio-band detectors sensitive out to redshifts of $z \approx 50$ for GW150914-like signals, gravitational-wave astronomy will enable us to probe the Universe when the first globular clusters formed. We simulate a population of binary black hole mergers from theoretically-motivated globular cluster formation models, and construct redshift measurements consistent with the predicted accuracy of third-generation detectors. We show that we can locate the peak time of a cluster formation epoch during reionisation to within 0.05Gyr after one year of observations. The peak of a formation epoch that coincides with the Universal star formation rate can be measured to within 0.4Gyr-10.5Gyr after one year of observations, depending on the relative weighting of the model components.

The recent release of the second Gravitational-Wave Transient Catalog (GWTC-2) has increased significantly the number of known GW events, enabling unprecedented constraints on formation models of compact binaries. One pressing question is to understand the fraction of binaries originating from different formation channels, such as isolated field formation versus dynamical formation in dense stellar clusters. In this paper, we combine the $\texttt{COSMIC}$ binary population synthesis suite and the $\texttt{CMC}$ code for globular cluster evolution to create a mixture model for black hole binary formation under both formation scenarios. For the first time, these code bodies are combined self-consistently, with $\texttt{CMC}$ itself employing $\texttt{COSMIC}$ to track stellar evolution. We then use a deep-learning enhanced hierarchical Bayesian analysis to constrain the mixture fraction $f$ between formation models, while simultaneously constraining the common envelope efficiency $\alpha$ assumed in $\texttt{COSMIC}$ and the initial cluster virial radius $r_v$ assumed in $\texttt{CMC}$. Under specific assumptions about other uncertain aspects of isolated binary and globular cluster evolution, we report the median and $90\%$ confidence interval of three physical parameters $(f,\alpha,r_v)=(0.20^{+0.32}_{-0.18},2.26^{+2.65}_{-1.84},2.71^{+0.83}_{-1.17})$. This simultaneous constraint agrees with observed properties of globular clusters in the Milky Way and is an important first step in the pathway toward learning astrophysics of compact binary formation from GW observations., Comment: 8 pages, 5 figures

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., Comment: 21 Pages, 11 Figures, 2 Tables, Accepted for publication in ApJ

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., Comment: 28 pages, 12 figures, 2 tables, accepted by ApJ

LIGO's third observing run (O3) has reported several neutron star-black hole (NSBH) merger candidates. From a theoretical point of view, NSBH mergers have received less attention in the community than either binary black holes (BBHs), or binary neutron stars (BNSs). Here we examine single-single (sin-sin) gravitational wave (GW) captures in different types of star clusters -- galactic nuclei (GN), globular clusters (GC), and young stellar clusters (YSC) -- and compare the merger rates from this channel to other proposed merger channels in the literature. There are currently large uncertainties associated with every merger channel, making a definitive conclusion about the origin of NSBH mergers impossible. However, keeping these uncertainties in mind, we find that sin-sin GW capture is unlikely to significantly contribute to the overall NSBH merger rate. In general, it appears that isolated binary evolution in the field or in clusters, and dynamically interacting binaries in triple configurations, may result in a higher merger rate., Comment: Accepted to ApJ

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., Comment: 25 pages, 6 figures, accepted for publication in ApJ. Comments welcome

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., Comment: 9 pages, 4 figures, matches version accepted by ApJL