The First Interferometric Survey of Massive YSOs in the K-band Hot Dust, Ionised Gas, and Binarity at Au Scales

Circumstellar discs are essential for high mass star formation, while multiplicity, in particular binarity, appears to be an inevitable outcome since the vast majority of massive stars (> 8 Msun) are found in binaries (up to 100%). We constrain the sizes of the dust and ionised gas (Brgamma) emission of the innermost regions towards a sample of six MYSOs, and provide high-mass binary statistics of young stars at 2-300 au scales using VLTI (GRAVITY, AMBER) observations. We determine the inner radius of the dust emission and place MYSOs with K-band measurements in a size-luminosity diagram for the first time, and compare our findings to T Tauris and Herbig AeBes. We also compare the observed K-band sizes to the sublimation radius predicted by three different disc scenarios. Lastly, we apply binary geometries to trace close binarity among MYSOs. The inner sizes of MYSOs, Herbig AeBe and T Tauri stars appear to follow a universal trend at which the sizes scale with the square-root of the stellar luminosity. The Brgamma emission originates from somewhat smaller and co-planar area compared to the 2.2 {\mu}m continuum emission. We discuss this new finding with respect to disc-wind or jet origin. Finally, we report an MYSO binary fraction of 17-25% at milli-arcsecond separations (2-300 au). The size-luminosity diagram indicates that the inner regions of discs around young stars scale with luminosity independently of the stellar mass. At the targeted scales (2-300 au), the MYSO binary fraction is lower than what was previously reported for the more evolved main sequence massive stars, which, if further confirmed, could implicate the predictions from massive binary formation theories. Lastly, we spatially resolve the crucial star/disc interface in a sample of MYSOs, showing that au-scale discs are prominent in high-mass star formation and similar to their low-mass equivalents.