An Estimate of the Binary Star Fraction among Young Stars at the Galactic Center: Possible Evidence of a Radial Dependence

Abstract We present the first estimate of the intrinsic binary fraction of young stars across the central ≈0.4 pc surrounding the supermassive black hole (SMBH) at the Milky Way Galactic center (GC). This experiment searched for photometric variability in 102 spectroscopically confirmed young stars, using 119 nights of 10″ wide adaptive optics imaging observations taken at W. M. Keck Observatory over 16 yr in the K ′ -2.1 μ m and H -1.6 μ m bands. We photometrically detected three binary stars, all of which are situated more than 1″ (0.04 pc) from the SMBH and one of which, S2-36, is newly reported here with spectroscopic confirmation. All are contact binaries or have photometric variability originating from stellar irradiation. To convert the observed binary fraction into an estimate of the underlying binary fraction, we determined the experimental sensitivity through detailed light-curve simulations, incorporating photometric effects of eclipses, irradiation, and tidal distortion in binaries. The simulations assumed a population of young binaries, with stellar ages (4 Myr) and masses matched to the most probable values measured for the GC young star population, and underlying binary system parameters (periods, mass ratios, and eccentricities) similar to those of local massive stars. As might be expected, our experimental sensitivity decreases for eclipses narrower in phase. The detections and simulations imply that the young, massive stars in the GC have a stellar binary fraction ≥71% (68% confidence), or ≥42% (95% confidence). This inferred GC young star binary fraction is consistent with that typically seen in young stellar populations in the solar neighborhood. Furthermore, our measured binary fraction is significantly higher than that recently reported by Chu et al. based on radial velocity measurements for stars ≲1″ of the SMBH. Constrained with these two studies, the probability that the same underlying young star binary fraction extends across the entire region is <1.4%. This tension provides support for a radial dependence of the binary star fraction, and therefore, for the dynamical predictions of binary merger and evaporation events close to the SMBH.