Search

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−1 (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. [ABSTRACT FROM AUTHOR]

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 an |${\sim }100\hbox{-}\mathrm{M}_\odot$| BH to grow beyond |${\sim }1000 \, \mathrm{M}_\odot$| through successive mergers ranges from |${\sim }0.1~{{\ \rm per\ cent}}$| in low-density, low-mass clusters to nearly 90  per cent in high-mass, high-density clusters. However, in the most massive NSCs, the growth time-scale 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 behaviour. 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–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 \, \mathrm{M}_\odot }$| in an NSC. [ABSTRACT FROM AUTHOR]

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 Early Data Release 3 (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 (β ≈ −0.4 ± 0.1) in the outskirts. We also robustly measure a dark central mass of roughly |$800\pm 300 \, \rm M_\odot$| , but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), and a dark population of stellar remnants of extent |${\approx} 0.016\, {\rm pc} \simeq 3300\, {\rm 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 versus remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core. [ABSTRACT FROM AUTHOR]

We study close encounters of a 1 M ⊙ 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 BH = 100–104 M ⊙. 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 BH ≲ 103 M ⊙. 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 ⊙ star undergoes 16 (5) pericenter passages before ejection from a 10 M ⊙ (100 M ⊙) black hole. Observations of periodic flares from these repeated close passages could in principle indicate signatures of a partial tidal disruption event. [ABSTRACT FROM AUTHOR]

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 disc scales given sufficient 'fuelling' from larger scales, what remains unclear is whether or not BHs can actually capture sufficient gas from their surrounding interstellar medium (ISM). We explore this in a suite of multiphysics 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}\, \mathrm{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}\, \mathrm{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{{\ \rm per\ cent}}$| 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 intermediate-mass black hole (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. [ABSTRACT FROM AUTHOR]

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\, \rm M_\odot$| 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 reduces 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 Ă— 10âˆ'14 ≲ ÎĽ /eV ≲ 2 Ă— 10âˆ'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 ≈60% 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 ≈0.08 Gpcâˆ'3 yrâˆ'1 in the local universe. These results imply that a select range of ultralight boson masses 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. [ABSTRACT FROM AUTHOR]

There is increasing evidence that globular clusters retain sizeable black hole populations at present day. This is supported by dynamical simulations of cluster evolution, which have unveiled the spatial distribution and mass spectrum of black holes in clusters across cosmic age. However, black hole populations of young, high metallicity clusters remain unconstrained. Black holes hosted by these clusters mass segregate early in their evolutionary history, forming central subsystems of hundreds to thousands of black holes. We argue that after supernova feedback has subsided (≳50 Myr), the host cluster can accumulate gas from its dense surroundings, from which the black hole subsystem accretes at highly enhanced rates. The collective accretion luminosity can be substantial and provides a novel observational constraint for young massive clusters. We test this hypothesis by performing 3D hydrodynamic simulations where we embed discretized potentials, representing our black holes, within the potential of a massive cluster. This system moves supersonically with respect to a gaseous medium from which it accretes. We study the accretion of this black hole subsystem for different subsystem populations and determine the integrated accretion luminosity of the black hole subsystem. We apply our results to the young massive clusters of the Antennae Galaxies and find that a typical subsystem accretion luminosity should be in excess of ≈ 1040 erg s−1. We argue that no strong candidates of this luminous signal have been observed and constrain the subsystem population of a typical cluster in the Antennae Galaxies to ≲ 10–2 × 102 10 M⊙ black holes, given that feedback does not significantly impede accretion and that the gas remains optically thin. [ABSTRACT FROM AUTHOR]

As ancient, gravitationally bound stellar populations, globular clusters represent 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 CMC Cluster Catalog. Without performing 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 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. In addition, 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 python methods for examining the publicly accessible CMC Cluster Catalog , as well as any other models generated using CMC. [ABSTRACT FROM AUTHOR]

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 of 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 Gpc−3 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 1044 erg s−1, and bright optical luminosities (roughly 1041–1044 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. [ABSTRACT FROM AUTHOR]

Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above ∼40–50 M⊙ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae (PISNe). 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 ∼20% of all BH progenitors undergo one or more collisions prior to stellar collapse and up to ∼1% 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 the Laser Interferometer Gravitational-Wave Observatory/Virgo in GW190521. Finally, we describe how stellar collisions in clusters may provide a unique pathway to PISNe and briefly discuss the expected rate of these events and other electromagnetic transients. [ABSTRACT FROM AUTHOR]

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/Kamioka Gravitational-Wave Detector (LVK), and future missions such as LISA and the DECi-hertz Interferometer Gravitational wave Observatory. 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 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. [ABSTRACT FROM AUTHOR]

Depending on the stellar type, more than 15% of stars in the field have at least two stellar companions. Hierarchical triple 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 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. We show that a typical GC is expected to host tens of triples with at least one luminous component at present. We discuss how the Lidov–Kozai mechanism can drive the inner binary of these dynamically 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 merger products, such as blue stragglers, X-ray binaries, Type Ia supernovae, Thorne–Zytkow objects, and gravitational wave sources. [ABSTRACT FROM AUTHOR]

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 superelastic scattering, delaying the 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 not typically in directly detectable configurations. We demonstrated in Weatherford et al. that an anticorrelation 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 in 50 Milky Way GCs from the Advanced Camera for Surveys GC Survey. For each GC, Δ is measured between observed main-sequence populations and fed into correlations between Δ and BH retention found in our CMC Cluster Catalog's models. We demonstrate that the range in measured Δ from our models matches that for observed GCs to a remarkable degree. Our results constitute the largest sample of GCs for which BH populations have been predicted to date using a self-consistent and robust statistical approach. We identify NGC 2808, 5927, 5986, 6101, and 6205 to retain especially large BH populations, each with a total BH mass exceeding 103. [ABSTRACT FROM AUTHOR]

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, ultraluminous 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. 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 –. 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 , , , and , 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 ≳ 0.1), providing an indication of their dynamical origin. [ABSTRACT FROM AUTHOR]

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, leading to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward "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. Over the past decade, a large amount of work has suggested that stellar black holes (BHs) play a dynamically-significant role in clusters throughout their entire lifetimes. Here we review our latest understanding 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." [ABSTRACT FROM AUTHOR]

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 ~21 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 ~7 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. [ABSTRACT FROM AUTHOR]