We present multi-wavelength images observed with SOFIA-FORCAST from ~10 to 40 $\mu$m of 14 protostars, selected as intermediate-mass protostar candidates, as part of the SOFIA Massive (SOMA) Star Formation Survey. We build protostellar spectral energy distributions (SEDs) with the SOFIA observations, together with archival data from Spitzer, Herschel and IRAS. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key properties of the protostars. The SEDs generally indicate the validity of these RT models down to intermediate-mass and/or early-stage protostars. The protostars analyzed so far in the SOMA survey span a range of luminosities from ~$10^{2}$ to ~$10^{6} L_{\odot}$, current protostellar masses from ~0.5 to ~30 $M_{\odot}$, and ambient clump mass surface densities, $\Sigma_{\rm cl}$ of 0.1 - 3 g cm$^{-2}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. The 19 to 37 $\mu$m spectral index of the sources correlates with outflow cavity opening angle, ratio of this angle to viewing angle, and evolutionary stage. We have also added a sample of ~50 protostellar sources identified from within Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models of the protostars, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to ~25 $M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm cl}$ needs to be at least 1 g cm$^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Aims: The inner regions of high-mass protostars are often invisible in the near-infrared. We aim to investigate the inner gaseous disc of IRAS11101-5829 through scattered light from the outflow cavity walls. Methods: We observed the environment of the high-mass young stellar object IRAS11101-5829 and the closest knots of its jet, HH135-136, with the VLT/SINFONI. We also retrieved archival data from the high-resolution long-slit spectrograph VLT/X-shooter. Results: We detect the first three bandheads of the $\upsilon=2-0$ CO vibrational emission for the first time in this object. It is coincident with continuum and Br$\gamma$ emission and extends up to $\sim10000$ au to the north-east and $\sim10 000$ au to the south-west. The line profiles have been modelled as a Keplerian rotating disc assuming a single ring in LTE. The model output gives a temperature of $\sim3000$ K, a CO column density of $\sim1\times10^{22}\mathrm{ cm^{-2}}$, and a projected Keplerian velocity $v_\mathrm{K}\sin i_\mathrm{disc} \sim 25\mathrm{ km s^{-1}}$, which is consistent with previous modelling in other high-mass protostars. In particular, the low value of $v_\mathrm{K}\sin i_\mathrm{disc}$ suggests that the disc is observed almost face-on, whereas the well-constrained geometry of the jet imposes that the disc must be close to edge-on. This apparent discrepancy is interpreted as the CO seen reflected in the mirror of the outflow cavity wall. Conclusions: From both jet geometry and disc modelling, we conclude that all the CO emission is seen through reflection by the cavity walls and not directly. This result implies that in the case of highly embedded objects, as for many high-mass protostars, line profile modelling alone might be deceptive and the observed emission could affect the derived physical and geometrical properties; in particular the inclination of the system can be incorrectly interpreted.

It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.

RW Aur A is a CTTS that has suddenly undergone three major dimming events since 2010. We aim to understand the dimming properties, examine accretion variability, and derive the physical properties of the inner disc traced by the CO ro-vibrational emission at NIR wavelengths (2.3 mic). We compared two epochs of X-Shooter observations, during and after the dimming. We modelled the rarely detected CO bandhead emission in both epochs to examine whether the inner disc properties had changed. The SED was used to derive the extinction properties of the dimmed spectrum and compare the infrared excess between the two epochs. Lines tracing accretion were used to derive the mass accretion rate in both states. The CO originates from a region with physical properties of T=3000 K, N$_{CO}$=1x10$^{21}$ cm$^{-2}$ and vsini=113 km/s. The extinction properties of the dimming layer were derived with the effective optical depth ranging from teff 2.5-1.5 from the UV to the NIR. The inferred mass accretion rate Macc is $1.5x 10^{-8}$ Msun/yr and $\sim 2x 10^{-8}$ Msun/yr after and during the dimming respectively. By fitting the SED, additional emission is observed in the IR during the dimming event from dust grains with temperatures of 500-700K. The physical conditions traced by the CO are similar for both epochs, indicating that the inner gaseous disc properties do not change during the dimming events. The extinction curve is flatter than that of the ISM, and large grains of a few hundred microns are thus required. When we correct for the observed extinction, Macc is constant in the two epochs, suggesting that the accretion is stable and therefore does not cause the dimming. The additional hot emission in the NIR is located at about 0.5 au from the star. The dimming events could be due to a dust-laden wind, a severe puffing-up of the inner rim, or a perturbation caused by the recent star-disc encounter.

Context: Protostellar jets in high-mass young stellar objects (HMYSOs) play a key role in the understanding of star formation and provide us with an excellent tool to study fundamental properties of HMYSOs. Aims: We aim at studying the physical and kinematic properties of the near-IR (NIR) jet of IRAS\,13481-6124 from au to parsec scales. Methods: Our study includes NIR data from the Very Large Telescope instruments SINFONI, CRIRES, and ISAAC. Information about the source and its immediate environment is retrieved with SINFONI. The technique of spectro-astrometry is performed with CRIRES to study the jet on au scales. The parsec-scale jet and its kinematic and dynamic properties are investigated using ISAAC. Results: The SINFONI spectra in the $H$ and $K$ band are rich in emission lines that are mainly associated with ejection and accretion processes. Spectro-astrometry is applied to the Br$\gamma$ line, and for the first time, to the Br$\alpha$ line, revealing their jet origin with milliarcsecond-scale photocentre displacements ($11-15$\,au). This allows us to constrain the kinematics of the au-scale jet and to derive its position angle ($\sim216^{\circ}$). ISAAC spectroscopy reveals H$_2$ emission along the parsec-scale jet, which allows us to infer kinematic and dynamic properties of the NIR parsec-scale jet. The mass-loss rate inferred for the NIR jet is $\dot{M}_\mathrm{ejec}\sim10^{-4}\mathrm{\,M_\odot\,yr^{-1}}$ and the thrust is $\dot{P}\sim10^{-2}\mathrm{\,M_\odot\,yr^{-1}\,km\,s^{-1}}$ , which is roughly constant for the formation history of the young star. A tentative estimate of the ionisation fraction is derived for the massive jet by comparing the radio and NIR mass-loss rates. An ionisation fraction $\lesssim8\%$ is obtained, which means that the bulk of the ejecta is traced by the NIR jet and that the radio jet only delineates a small portion of it.