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

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.

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

Recent discoveries of black hole (BH) candidates in Galactic and extragalactic globular clusters (GCs) have ignited interest in understanding how BHs dynamically evolve in a GC and the number of BHs ($N_{\rm{BH}}$) that may still be retained by today's GCs. Numerical models show that even if stellar-mass BHs are retained in today's GCs, they are typically in configurations that are not directly detectable. We show that a suitably defined measure of mass segregation ($\Delta$) between, e.g., giants and low-mass main-sequence stars, can be an effective probe to indirectly estimate $N_{\rm{BH}}$ in a GC aided by calibrations from numerical models. Using numerical models including all relevant physics we first show that $N_{\rm{BH}}$ is strongly anticorrelated with $\Delta$ between giant stars and low-mass main-sequence stars. We apply the distributions of $\Delta$ vs $N_{\rm{BH}}$ obtained from models to three Milky Way GCs to predict the $N_{\rm{BH}}$ retained by them at present. We calculate $\Delta$ using the publicly available ACS survey data for 47 Tuc, M 10, and M 22, all with identified stellar-mass BH candidates. Using these measured $\Delta$ and distributions of $\Delta$ vs $N_{\rm{BH}}$ from models as calibration we predict distributions for $N_{\rm{BH}}$ expected to be retained in these GCs. For 47 Tuc, M 10, and M 22 our predicted distributions peak at $N_{\rm{BH}}\approx20$, $24$, and $50$, whereas, within the $2\sigma$ confidence level, $N_{\rm{BH}}$ can be up to $\sim150$, $50$, and $200$, respectively.