Noninteracting Black Hole Binaries with Gaia and LAMOST

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in noninteracting BH binaries (nBHBs) through a precise analysis of the companion's motion. In this study, using an updated version of the StarTrack binary-star population modeling code and a detailed model of the Milky Way (MW) galaxy, we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty that is introduced by methods that Monte Carlo sample from binary population synthesis output to produce one potential specific realization of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realizations. Specifically, we predict ~41–340 nBHBs to be observed by Gaia, although the numbers may drop to ~10–70 if the recent ( 100 Myr) star formation is low (~1 M ⊙ yr⁻¹). For LAMOST we predict 14 detectable nBHBs, which is lower partially because its field of view covers just ~6% of the Galaxy