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.

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.

NU Ori is a massive spectroscopic and visual binary in the Orion Nebula Cluster, with 4 components: Aa, Ab, B, and C. The B0.5 primary (Aa) is one of the most massive B-type stars reported to host a magnetic field. We report the detection of a spectroscopic contribution from the C component in high-resolution ESPaDOnS spectra, which is also detected in a Very Large Telescope Interferometer (VLTI) dataset. Radial velocity (RV) measurements of the inner binary (designated Aab) yield an orbital period of 14.3027(7) d. The orbit of the third component (designated C) was constrained using both RVs and interferometry. We find C to be on a mildly eccentric 476(1) d orbit. Thanks to spectral disentangling of mean line profiles obtained via least-squares deconvolution we show that the Zeeman Stokes $V$ signature is clearly associated with C, rather than Aa as previously assumed. The physical parameters of the stars were constrained using both orbital and evolutionary models, yielding $M_{\rm Aa} = 14.9 \pm 0.5 M_\odot$, $M_{\rm Ab} = 3.9 \pm 0.7 M_\odot$, and $M_{\rm C} = 7.8 \pm 0.7 M_\odot$. The rotational period obtained from longitudinal magnetic field $\langle B_z \rangle$ measurements is $P_{\rm rot} = 1.09468(7)$ d, consistent with previous results. Modeling of $\langle B_z \rangle$ indicates a surface dipole magnetic field strength of $\sim 8$ kG. NU Ori C has a magnetic field strength, rotational velocity, and luminosity similar to many other stars exhibiting magnetospheric H$\alpha$ emission, and we find marginal evidence of emission at the expected level ($\sim$1% of the continuum).

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

We investigate in the mid-IR the spatial properties of the PAHs emission in the disk of HD179218. We obtained mid-IR images in the PAH1, PAH2 and Si6 filters at 8.6, 11.3 and 12.5 mu, and N band low-resolution spectra using CanariCam on the GTC. We compared the PSFs measured in the PAH filters to the PSF derived in the Si6 filter, where the thermal continuum dominates. We performed radiative transfer modelling of the spectral energy distribution and produced synthetic images in the three filters to investigate different spatial scenarios. Our data show that the disk emission is spatially resolved in the PAHs filters, while unresolved in the Si6 filter. An average FHWM of 0.232", 0.280" and 0.293" is measured in the three filters. Gaussian disk fitting and quadratic subtraction of the science and calibrator suggest a lower-limit characteristic angular diameter of the emission of circa 100 mas (circa 40 au). The photometric and spectroscopic results are compatible with previous findings. Our radiative transfer (RT) modelling of the continuum suggests that the resolved emission results from PAH molecules on the disk atmosphere being UV-excited by the central star. Geometrical models of the PAH component compared to the underlying continuum point at a PAH emission uniformly extended out to the physical limits of the disk's model. Also, our RT best model of the continuum requires a negative exponent of the surface density power-law, in contrast to earlier modelling pointing at a positive exponent. Based on spatial and spectroscopic considerations as well as on qualitative comparison with IRS48 and HD97048, we favor a scenario in which PAHs extend out to large radii across the flared disk surface and are at the same time predominantly in an ionized charge state due to the strong UV radiation field of the 180 L_sun central star.

Modeling stellar atmospheres is a complex, intriguing task in modern astronomy. A systematic comparison of models with multi-technique observations is the only efficient way to constrain them. Aims: We performed a self-consistent modeling of the atmospheres of six C-rich AGB stars: R Lep, R Vol, Y Pav, AQ Sgr, U Hya and X TrA, with the aim of enlarging the knowledge of the dynamic processes occurring in their atmospheres. Methods. We used VLTI/MIDI interferometric observations, in combination with spectro-photometric data, and compared them with self-consistent dynamic models atmospheres. Results: The models can reproduce SED data well at wavelengths longwards of 1{\mu}m, and the interferometric observations between 8{\mu}m and 10{\mu}m. Discrepancies could be due to a combination of data- and model-related effects. The models best fitting the Miras are significantly extended, with a prominent shell-like structure, while the models best fitting the non-Miras are more compact, showing lower average mass-loss. The mass loss is of episodic or multi-periodic nature, but causes the visual amplitudes to be notably larger than the observed ones. Stellar parameters were derived from the model fitting: T_eff, L_bol, M, C/O, $\dot{\text{M}}$. Our findings agree well with literature values within the uncertainties. T_eff and L_bol are also in good agreement with the T derived from the angular diameter {\theta}_{(V-K)} and L_bol from the SED fitting, except for AQ Sgr. Finally, {\theta}_Ross and {\theta}_(V-K) agree with each other better for the Miras than for the non-Miras,which is probably connected to the episodic nature of the latter models. We also located the stars in the H-R diagram, comparing them with evolutionary tracks. We found that the main derived properties from the model fitting are in good agreement with TP-AGB evolutionary calculations for carbon stars (COLIBRI code).

Context. The mass loss from massive stars is not understood well. Eta Car is a unique object for studying the massive stellar wind during the LBV phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Methods. Observations of Eta Car were carried out with the ESO VLTI and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 different spectral channels distributed across the Br Gamma 2.166 micrometer emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is ~ 126 degree. The fan-shaped structure extends approximately 8.0 mas (~ 18.8 au) to the southeast and 5.8 mas (~ 13.6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Our long term aim is to derive model-independent stellar masses and distances for long period massive binaries by combining apparent astrometric orbit with double-lined radial velocity amplitudes (SB2). We follow-up ten O+O binaries with AMBER, PIONIER and GRAVITY at the VLTI. Here, we report about 130 astrometric observations over the last 7 years. We combine this dataset with distance estimates to compute the total mass of the systems. We also compute preliminary individual component masses for the five systems with available SB2 radial velocities. Nine over the ten binaries have their three dimensional orbit well constrained. Four of them are known colliding wind, non-thermal radio emitters, and thus constitute valuable targets for future high angular resolution radio imaging. Two binaries break the correlation between period and eccentricity tentatively observed in previous studies. It suggests either that massive star formation produce a wide range of systems, or that several binary formation mechanisms are at play. Finally, we found that the use of existing SB2 radial velocity amplitudes can lead to unrealistic masses and distances. If not understood, the biases in radial velocity amplitudes will represent an intrinsic limitation for estimating dynamical masses from SB2+interferometry or SB2+Gaia. Nevertheless, our results can be combined with future Gaia astrometry to measure the dynamical masses and distances of the individual components with an accuracy of 5 to 15\%, completely independently of the radial velocities.

In anticipation of the Gaia astrometric mission, a sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when combined with Gaia astrometric measurements. In order to validate the masses derived from Gaia, interferometric observations are obtained for three SB2s in our sample with F-K components: HIP 14157, HIP 20601 and HIP 117186. The masses of the six stellar components are derived. Due to its edge-on orientation, HIP 14157 is probably an eclipsing binary. We note that almost all the derived masses are a few percent larger than the expectations from the standard spectral-type-mass calibration and mass-luminosity relation. Our calculation also leads to accurate parallaxes for the three binaries, and the Hipparcos parallaxes are confirmed.

We analyse and report in detail new near-infrared (1.45 - 2.45 microns) observations of the Dusty S-cluster Object (DSO/G2) during its approach to the black hole at the center of the Galaxy that were carried out with ESO VLT/SINFONI between February and September 2014. Before May 2014 we detect spatially compact Br-gamma and Pa-alpha line emission from the DSO at about 40mas east of SgrA*. The velocity of the source, measured from the red-shifted emission, is 2700+-60 km/s. No blue-shifted emission above the noise level is detected at the position of SgrA* or upstream the presumed orbit. After May we find spatially compact Br-gamma blue-shifted line emission from the DSO at about 30mas west of SgrA* at a velocity of -3320+-60 km/s and no indication for significant red-shifted emission. We do not detect any significant extension of velocity gradient across the source. We find a Br-gamma-line full width at half maximum of 50+-10 Angstroem before and 15+-10 Angstroem after the peribothron transit, i.e. no significant line broadening with respect to last year is observed. Br-gamma line maps show that the bulk of the line emission originates from a region of less than 20mas diameter. This is consistent with a very compact source on an elliptical orbit with a peribothron time passage in 2014.39+-0.14. For the moment, the flaring activity of the black hole in the near-infrared regime has not shown any statistically significant increment. Increased accretion activity of SgrA* may still be upcoming. We discuss details of a source model according to which the DSO is rather a young accreting star than a coreless gas and dust cloud.

We present the latest update of the European Southern Observatory's Very Large Telescope interferometer (VLTI). The operations of VLTI have greatly improved in the past years: reduction of the execution time; better offering of telescopes configurations; improvements on AMBER limiting magnitudes; study of polarization effects and control for single mode fibres; fringe tracking real time data, etc. We present some of these improvements and also quantify the operational improvements using a performance metric. We take the opportunity of the first decade of operations to reflect on the VLTI community which is analyzed quantitatively and qualitatively. Finally, we present briefly the preparatory work for the arrival of the second generation instruments GRAVITY and MATISSE.

Context. The Gum 31 bubble containing the stellar cluster NGC 3324 is a poorly-studied young region close to the Carina Nebula. Aims. We are aiming to characterise the young stellar and protostellar population in and around Gum 31 and to investigate the star-formation process in this region. Methods. We identify candidate young stellar objects from Spitzer, WISE, and Herschel data. Combining these, we analyse the spectral energy distributions of the candidate young stellar objects. With density and temperature maps obtained from Herschel data and comparisons to a 'collect and collapse' scenario for the region we are able to further constrain the characteristics of the region as a whole. Results. 661 candidate young stellar objects are found from WISE data, 91 protostar candidates are detected through Herschel observations in a 1.0 deg x 1.1 deg area. Most of these objects are found in small clusters or are well aligned with the H II bubble. We also identify the sources of Herbig-Haro jets. The infrared morphology of the region suggests that it is part of the larger Carina Nebula complex. Conclusions. The location of the candidate young stellar objects in the rim of the H II bubble is suggestive of their being triggered by a 'collect and collapse' scenario, which agrees well with the observed parameters of the region. Some candidate young stellar objects are found in the heads of pillars, which points towards radiative triggering of star formation. Thus, we find evidence that in the region different mechanisms of triggered star formation are at work. Correcting the number of candidate young stellar objects for contamination we find ~ 600 young stellar objects in Gum 31 above our completeness limit of about 1 M_sol. Extrapolating the intital mass function down to 0.1 M_sol, we estimate a total population of ~ 5000 young stars for the region.

The characterization of multiple stellar systems is an important ingredient for testing current star formation models. Stars are more often found in multiple systems, the more massive they are. A complete knowledge of the multiplicity of high-mass stars over the full range of orbit separations is thus essential to understand their still debated formation process. Observations of the Orion Nebula Cluster can help to answer the question about the origin and evolution of multiple stars. Earlier studies provide a good knowledge about the multiplicity of the stars at very small (spectroscopic) and large separations (AO, speckle) and thus make the ONC a good target for such a project. We used the NIR interferometric instrument AMBER at VLTI to observe a sample of bright stars in the ONC. We complement our data set by archival NACO observations of \theta 1 Ori A to obtain more information about the orbit of the close visual companion. Our observations resolve the known multiple systems \theta 1 Ori C and \theta 1 Ori A and provide new orbit points, which confirm the predicted orbit and the determined stellar parameters for \theta 1 Ori C. Combining AMBER and NACO data for \theta 1 Ori A we were able to follow the motion of the companion from 2003 to 2011. We furthermore find hints for a companion around \theta 1 Ori D and a previously unknown companion to NU Ori. With a probability of ~90% we can exclude further companions with masses of > 3 Msun around our sample stars for separations between ~2 mas and ~110 mas. We conclude that the companion around \theta 1 Ori A is most likely physically related to the primary star. The newly discovered possible companions further increase the multiplicity in the ONC. For our sample of two O and three B-type stars we find on average 2.5 known companions per primary, which is around five times more than for low-mass stars.

Aims: Jets are excellent signposts for very young embedded protostars, so we want to identify jet-driving protostars as a tracer of the currently forming generation of stars in the Carina Nebula, which is one of the most massive galactic star-forming regions and which is characterised by particularly high levels of massive-star feedback on the surrounding clouds. Methods: We used archive data to construct large (> 2 deg x 2 deg) Spitzer IRAC mosaics of the Carina Nebula and performed a spatially complete search for objects with excesses in the 4.5 micron band, typical of shock-excited molecular hydrogen emission. We also identified the mid-infrared point sources that are the likely drivers of previously discovered Herbig-Haro jets and molecular hydrogen emission line objects. We combined the Spitzer photometry with our recent Herschel far-infrared data to construct the spectral energy distributions, and used the Robitaille radiative-transfer modelling tool to infer the properties of the objects. Results: The radiative-transfer modelling suggests that the jet sources are protostars with masses between ~1 M_sol and ~10 M_sol that are surrounded by circumstellar disks and embedded in circumstellar envelopes. Conclusions: The estimated protostar masses < 10 M_sol suggest that the current star-formation activity in the Carina Nebula is restricted to low- and intermediate-mass stars. More optical than infrared jets can be observed, indicating that star formation predominantly takes place close to the surfaces of clouds.

The optically invisible infrared-source NGC 2264 IRS 1 is thought to be a massive young stellar object (~10 Msun). Although strong infrared excess clearly shows that the central object is surrounded by large amounts of circumstellar material, no information about the spatial distribution of this circumstellar material has been available until now. We used the ESO Very Large Telescope Interferometer to perform long-baseline interferometric observations of NGC 2264 IRS 1 in the mid-infrared regime. Our observations resolve the circumstellar material around NGC 2264 IRS 1, provide the first direct measurement of the angular size of the mid-infrared emission, and yield direct constraints on the spatial distribution of the dust. We use different approaches (a geometrical model, a temperature-gradient model, and radiative transfer models) to jointly model the observed interferometric visibilities and the spectral energy distribution. The derived visibility values between ~0.02 and ~0.3 show that the mid-infrared emission is clearly resolved. The characteristic size of the MIR-emission region is ~30-60 AU; this value is typical for other YSOs with similar or somewhat lower luminosities. A comparison of the sizes for the two position angles shows a significant elongation of the dust distribution. Simple spherical envelope models are therefore inconsistent with the data. The radiative transfer modeling of our data suggests that we observe a geometrically thin and optically thick circumstellar disk with a mass of about 0.1 Msun. Our modeling indicates that NGC 2264 IRS 1 is surrounded by a flat circumstellar disk that has properties similar to disks typically found around lower-mass young stellar objects. This result supports the assumption that massive young stellar objects form via accretion from circumstellar disks.