2025ApJ...985L..42X

The dynamical formation of binary black holes (BBHs) in globular clusters (GCs) may contribute significantly to the observed gravitational-wave (GW) merger rate. Furthermore, the Laser Interferometer Space Antenna (LISA) may detect many BBH sources from GCs at mHz frequencies, enabling the characterization of such systems within the Milky Way and nearby Universe. In this work, we use Monte Carlo N-body simulations to construct a realistic sample of Galactic clusters, thus estimating the population, detectability, and parameter measurement accuracy of BBHs hosted within them. In particular, we show that the GW signal from 0.7 ± 0.7, 2.0 ± 1.7, 3.6 ± 2.3, and 13.4 ± 4.7 BBHs in Milky Way GCs can exceed the signal-to-noise ratio (SNR) threshold of SNR = 30, 5, 3, and 1 for a 10 yr LISA observation, with ∼50% of detectable sources exhibiting high eccentricities (e ≳ 0.9). Moreover, the Fisher matrix and Bayesian analyses of the GW signals indicate that these systems typically feature highly resolved orbital frequencies (δf_orb/f_orb ∼ 10^–7 to 10^–5) and eccentricities (δe/e ∼ 10^–3 to 0.1), as well as a measurable total mass when SNR exceeds ∼20. Notably, we show that high-SNR BBHs can be confidently localized to specific Milky Way GCs with a sky localization accuracy of δΩ ∼ 1 deg², and we address the large uncertainties in their distance measurement (δR ∼ 0.3-20 kpc). The detection and localization of even a single BBH in a Galactic GC would allow accurate tracking of its long-term orbital evolution, enable a direct test of the role of GCs in BBH formation, and provide a unique probe into the evolutionary history of Galactic clusters.