The onset of stellar multiplicity in massive star formation: A search for low-mass companions of massive young stellar objects with $L'$-band adaptive optics imaging
Given the high incidence of binaries among mature field massive stars, it is clear that multiplicity is an inevitable outcome of high-mass star formation. Understanding how massive multiples form requires the study of the birth environments of massive stars, covering the innermost to outermost regions. We aim to detect and characterise low-mass companions around massive young stellar objects (MYSOs) during and shortly after their formation phase. To investigate large spatial scales, we carried out an $L'$-band high-contrast direct imaging survey seeking low-mass companions (down to $L_{\text{bol}}\approx 10 L_{\odot}$, or late A-type) around thirteen previously identified MYSOs using the VLT/NACO instrument. From those images, we looked for the presence of companions on a wide orbit, covering scales from 300 to 56,000 au. Detection limits were determined for all targets and we tested the gravitational binding to the central object based on chance projection probabilities. We have discovered a total of thirty-nine potential companions around eight MYSOs, the large majority of which have never been reported to date. We derived a multiplicity frequency (MF) of $62\pm13$% and a companion fraction (CF) of $3.0\pm0.5$. The derived MF and CF are compared to other studies for similar separation ranges. The comparisons are effective for a fixed evolutionary stage spanning a wide range of masses and vice versa. We find an increased MF and CF compared to the previous studies targeting MYSOs, showing that the statement in which multiplicity scales with primary mass also extends to younger evolutionary stages. The separations at which the companions are found and their location with relation to the primary star allow us to discuss the implications for the massive star formation theories.