Multimass modelling of Milky Way globular clusters -- II. present-day black hole populations

Populations of stellar-mass black holes (BHs) in globular clusters (GCs) influence their dynamical evolution and have important implications on one of the main formation channels for gravitational wave sources. Inferring the size of these populations remains difficult, however. In this work, multimass models of 34 Milky Way GCs, first presented in Dickson et al., are used to explore the present-day BH populations. Direct constraints on both the total and visible mass components provided by several observables allow these models to accurately determine the distribution of the dark mass (including BHs) within clusters, as we demonstrate in a proof-of-concept fitting of the models to mock observations extracted from Monte Carlo cluster models. New constraints on the BH population retained to the present-day in each cluster are inferred from our models. We find that BH mass fractions ranging from 0 to 1 per cent of the total mass are typically required to explain the observations, except for Omega Cen, for which we infer a mass fraction above 5 per cent, in agreement with previous works. Relationships between the dark remnant populations and other cluster parameters are examined, demonstrating a clear anti-correlation between the amount of BHs and mass segregation between visible stars, as well as a correlation between remnant mass fractions and the dynamical age of clusters. Our inferred BH populations are in good agreement overall with other recent studies using different methodologies, but with notable discrepancies for individual clusters.