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.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

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.

HD50138 is a Herbig B[e] star with a circumstellar disc detected at IR and mm wavelength. Its brightness makes it a good candidate for NIR interferometry observations. We aim to resolve, spatially and spectrally, the continuum and hydrogen emission lines in the 2.12-2.47 micron region, to shed light on the immediate circumstellar environment of the star. VLTI/AMBER K-band observations provide spectra, visibilities, differential phases, and closure phases along three long baselines for the continuum, and HI emission in Br$\gamma$ and five high-n Pfund lines. By computing the pure-line visibilities, we derive the angular size of the different line-emitting regions. A simple LTE model was created to constrain the physical conditions of HI emitting region. The continuum region cannot be reproduced by a geometrical 2D elongated Gaussian fitting model. We estimate the size of the region to be 1 au. We find the Br$\gamma$ and Pfund lines come from a more compact region of size 0.4 au. The Br$\gamma$ line exhibits an S-shaped differential phase, indicative of rotation. The continuum and Br$\gamma$ line closure phase show offsets of $\sim$-25$\pm$5 $^o$ and 20$\pm$10$^o$, respectively. This is evidence of an asymmetry in their origin, but with opposing directions. We find that we cannot converge on constraints for the HI physical parameters without a more detailed model. Our analysis reveals that HD50138 hosts a complex circumstellar environment. Its continuum emission cannot be reproduced by a simple disc brightness distribution. Similarly, several components must be evoked to reproduce the interferometric observables within the Br$\gamma$, line. Combining the spectroscopic and interferometric data of the Br$\gamma$ and Pfund lines favours an origin in a wind region with a large opening angle. Finally, our results point to an evolved source.

In order to constrain the evolutionary state of the red supergiant Betelgeuse, we have produced a suite of models with ZAMS masses from 15 to 25 Msun in intervals of 1 Msun including the effects of rotation. The models were computed with the stellar evolutionary code MESA. For non-rotating models we find results that are similar to other work. It is somewhat difficult to find models that agree within 1 sigma of the observed values of R, Teff and L, but modestly easy within 3 sigma uncertainty. Incorporating the nominal observed rotational velocity, ~15 km/s, yields significantly different, and challenging, constraints. This velocity constraint is only matched when the models first approach the base of the red supergiant branch (RSB), having crossed the Hertzsprung gap, but not yet having ascended the RSB and most violate even generous error bars on R, Teff and L. Models at the tip of the RSB typically rotate at only ~0.1 km/s, independent of any reasonable choice of initial rotation. We discuss the possible uncertainties in our modeling and the observations, including the distance to Betelgeuse, the rotation velocity, and model parameters. We summarize various options to account for the rotational velocity and suggest that one possibility is that Betelgeuse merged with a companion star of about 1 Msun as it ascended the RSB, in the process producing the ring structure observed at about 7' away. A past coalescence would complicate attempts to understand the evolutionary history and future of Betelgeuse.