Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C-and Q-bands along with spectral index maps to characterise their emission. We find that about 50% of the sources show evidence for non-thermal emission with 40% showing clear non-thermal lobes, especially sources of higher bolometric luminosity. The common or IRAS names of the sources that manifest non-thermal lobes are; V645Cyg, IRAS 22134+5834, NGC 7538 IRS 9, IRAS 23262+640, AFGL 402d and AFGL 490. All the central cores of the sources are thermal with corresponding mass-loss rates that lie in the range ∼ 3 × 10 −7 to 7 × 10 −6 M yr −1 . Given the presence of non-thermal lobes in some of the sources and the evidence of non-thermal emission from some spectral index maps, it seems that magnetic fields play a significant role in the jets of massive protostars. Also noted is that some of the sources show evidence of binarity and variability.

We report the results of the Jansky Very Large Array (JVLA) observation of five massive protostars at 6 and 22.2 GHz. The aim of the study was to compare their current fluxes and positions with previous observations to search for evidence of variability. Most of the observed sources present the morphologies of a thermal core, hosting the protostar and exhibiting no proper motion, and associated non-thermal radio lobes that are characterized by proper motions and located away from the thermal core. Some of the protostars drive jets whose lobes have dissimilar displacement vectors, implying precession of the jets or the presence of multiple jet drivers. The jets of the protostars were found to have proper motions that lie in the range of 170 ≤ v ≤ 650 km s−1, and precessions of periods of 40 ≤ p ≤ 50 yr and angles of 2 ≤ α ≤ 10°, assuming that their velocities v = 500 km s−1. The core of one of the sources, S255 NIRS3, which was in outburst at the time of our observations, showed a significant change in flux compared to the other sources. Its spectral index decreased during the outburst, consistent with the model of an expanding gas bubble. Modelling the emission of the outburst as that of a new non-thermal lobe that is emerging from a thermal core whose emission enshrouds that of the lobe also has the potential to account for the increase in flux and a decrease in the spectral index of the source’s outburst.