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.

We analyzed high angular resolution 45.5 GHz images of the W49 North massive star forming region obtained in 1998 and 2016 with the Very Large Array. Most of the ultracompact HII regions show no detectable changes over the time interval of the observations. However, subcomponents B1, B2, G2a and G2c have increased its peak flux densities by values in the range of 3.8 to 21.4 \%. Most interestingly, the cometary region C clearly shows proper motions that at the distance of the region are equivalent to a velocity of 76$\pm$6 km s$^{-1}$ in the plane of the sky. We interpret this region as the ionized bowshock produced by a runaway O6 ZAMS star that was ejected from the eastern edge of Welch's ring about 6,400 years ago.

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.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data.

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.

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.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

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.

The star V766 Cen (=HR 5171A) was originally classified as a yellow hypergiant but lately found to more likely be a 27-36 Msun red supergiant (RSG). Recent observations indicated a close eclipsing companion in the contact or common-envelope phase. Here, we aim at imaging observations of V766 Cen to confirm the presence of the close companion. We used near-infrared H -band aperture synthesis imaging at three epochs in 2014, 2016, and 2017, employing the PIONIER instrument at the Very Large Telescope Interferometer (VLTI). The visibility data indicate a mean Rosseland angular diameter of 4.1+/-0.8 mas, corresponding to a radius of 1575+/-400 Rsun. The data show an extended shell (MOLsphere) of about 2.5 times the Rosseland diameter, which contributes about 30% of the H-band flux. The reconstructed images at the 2014 epoch show a complex elongated structure within the photospheric disk with a contrast of about 10%. The second and third epochs show qualitatively and quantitatively different structures with a single very bright and narrow feature and high contrasts of 20-30%. This feature is located toward the south-western limb of the photospheric stellar disk. We estimate an angular size of the feature of 1.7+/-0.3 mas, corresponding to a radius of 650+/-150 Rsun, and giving a radius ratio of 0.42+0.35/-0.10} compared to the primary stellar disk. We interpret the images at the 2016 and 2017 epochs as showing the close companion, or a common envelope toward the companion, in front of the primary. At the 2014 epoch, the close companion is behind the primary and not visible. Instead, the structure and contrast at the 2014 epoch are typical of a single RSG harboring giant photospheric convection cells. The companion is most likely a cool giant or supergiant star with a mass of 5+15/-3 Msun.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

Characterization of the dynamics of massive star systems and the astrophysical properties of the interacting components are a prerequisite for understanding their formation and evolution. Optical interferometry at milliarcsecond resolution is a key observing technique for resolving high-mass multiple compact systems. Here we report on VLTI/GRAVITY, Magellan/FIRE, and MPG2.2m/FEROS observations of the late-O/early-B type system HD 93206 A, which is a member of the massive cluster Collinder 228 in the Carina nebula complex. With a total mass of about 90 M_sun, it is one of the most compact massive-quadruple systems known. In addition to measuring the separation and position angle of the outer binary Aa - Ac, we observe BrG and HeI variability in phase with the orbital motion of the two inner binaries. From the differential phases analysis, we conclude that the BrG emission arises from the interaction regions within the components of the individual binaries, which is consistent with previous models for the X-ray emission of the system based on wind-wind interaction. With an average 3-sigma deviation of ~15 deg, we establish an upper limit of p ~ 0.157 mas (0.35 AU) for the size of the BrG line-emitting region. Future interferometric observations with GRAVITY using the 8m UTs will allow us to constrain the line-emitting regions down to angular sizes of 20 uas (0.05 AU at the distance of the Carina nebula).

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

We present high spatial resolution mid-infrared images of the nebula around the late-type carbon-rich Wolf-Rayet (WC)-OB binary system WR~112 taken by the recently upgraded VLT spectrometer and imager for the mid-infrared (VISIR) with the PAH1, NeII\_2, and Q3 filters. The observations reveal a morphology resembling a series of arc-like filaments and broken shells. Dust temperatures and masses are derived for each of the identified filamentary structures, which exhibit temperatures ranging from $179_{-6}^{+8}$ K at the exterior W2 filament to $355_{-25}^{+37}$ K in the central 3". The total dust mass summed over the features is $2.6\pm0.4\times10^{-5}$ $\mathrm{M}_\odot$. A multi-epoch analysis of mid-IR photometry of WR~112 over the past $\sim20$ yr reveals no significant variability in the observed dust temperature and mass. The morphology of the mid-IR dust emission from WR~112 also exhibits no significant expansion from imaging data taken in 2001, 2007, and 2016, which disputes the current interpretation of the nebula as a high expansion velocity ($\sim1200$ km s$^{-1}$) "pinwheel"-shaped outflow driven by the central WC-OB colliding-wind binary. An upper limit of $\lesssim120$ km s$^{-1}$ is derived for the expansion velocity assuming a distance of $4.15$ kpc. The upper limit on the average total mass-loss rate from the central 3" of WR~112 is estimated to be $\lesssim8\times10^{-6}$ $\mathrm{M}_\odot$ yr$^{-1}$. We leave its true nature as an open question, but propose that the WR~112 nebula may have formed in the outflow during a previous red or yellow supergiant phase of the central Wolf-Rayet star.

Since their discovery, the Quintuplet proper members (QPMs) have been somewhat mysterious in nature. Originally dubbed the "cocoon stars" due to their cool featureless spectra, high-resolution near-infrared imaging observations have shown that at least two of the objects exhibit "pinwheel" nebulae consistent with binary systems with a carbon-rich Wolf-Rayet star and O/B companion. In this paper, we present 19.7, 25.2, 31.5, and 37.1 $\mu$m observations of the QPMs (with an angular resolution of 3.2-3.8") taken with the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST) in conjunction with high-resolution ($\sim$0.1-0.2") images at 8.8 and 11.7 $\mu$m from the Thermal-Region Camera Spectrograph (TReCS). DUSTY models of the thermal dust emission of two of the four detected QPMs, Q2 and Q3, are fitted by radial density profiles which are consistent with constant mass loss rates ($\rho_d \propto r^{-2}$). For the two remaining sources, Q1 and Q9, extended structures ($\sim$ 1") are detected around these objects in high-resolution imaging data. Based on the fitted dust masses, Q9 has an unusually large dust reservoir ($\mathrm{M_d}=1.3^{+0.8}_{-0.4}\times 10^{-3} \mathrm{M_{\odot}}$) compared to typical dusty Wolf-Rayet stars which suggests that it may have recently undergone an episode of enhanced mass loss.

Contemporary theory holds that massive stars gather mass during their initial phases via accreting disk-like structures. However, conclusive evidence for disks has remained elusive for the most massive young objects. This is mainly due to significant observational challenges. Incisive studies, even targeting individual objects, are therefore relevant to the progression of the field. NGC 3603 IRS 9A* is a young massive stellar object still surrounded by an envelope of molecular gas. Previous mid-infrared observations with long-baseline interferometry provided evidence for a disk of 50 mas diameter at its core. This work aims at a comprehensive study of the physics and morphology of IRS 9A at near-infrared wavelengths. New sparse aperture masking interferometry data taken with NACO/VLT at Ks and Lp filters were obtained and analysed together with archival CRIRES spectra of the H2 and BrG lines. The calibrated visibilities recorded at Ks and Lp bands suggest the presence of a partially resolved compact object of 30 mas at the core of IRS 9A, together with the presence of over-resolved flux. The spectroastrometric signal of the H2 line shows that this spectral feature proceeds from the large scale extended emission (300 mas) of IRS 9A, while the BrG line appears to be formed at the core of the object (20 mas). This scenario is consistent with the brightness distribution of the source for near- and mid-infrared wavelengths at various spatial scales. However, our model suffers from remaining inconsistencies between SED modelling and the interferometric data. Moreover, the BrG spectroastrometric signal indicates that the core of IRS 9A exhibits some form of complexity such as asymmetries in the disk. Future high-resolution observations are required to confirm the disk/envelope model and to flesh out the details of the physical form of the inner regions of IRS 9A.

We provide a complete characterization of the astrophysical properties of the $\sigma$ Ori Aa,Ab,B hierarchical triple system, and an improved set of orbital parameters for the highly eccentric $\sigma$ Ori Aa,Ab spectroscopic binary. We compiled a spectroscopic dataset comprising 90 high-resolution spectra covering a total time span of 1963 days. We applied the Lehman-Filh\'es method for a detailed orbital analysis of the radial velocity curves and performed a combined quantitative spectroscopic analysis of the {$\sigma$ Ori Aa,Ab,B} system by means of the stellar atmosphere code FASTWIND. We used our own plus other available information on photometry and distance to the system for measuring the radii, luminosities, and spectroscopic masses of the three components. We also inferred evolutionary masses and stellar ages using the Bayesian code BONNSAI. The orbital analysis of the new radial velocity curves led to a very accurate orbital solution of the $\sigma$ Ori Aa,Ab pair. We provided indirect arguments indicating that $\sigma$ Ori B is a fast rotating early-B dwarf. The FASTWIND+BONNSAI analysis showed that the Aa,Ab pair contains the hottest and most massive components of the triple system while $\sigma$ Ori B is a bit cooler and less massive. The derived stellar ages of the inner pair are intriguingly younger than the one widely accepted for the $\sigma$ Orionis cluster, at 3$\pm$1 Ma. The outcome of this study will be of key importance for a precise determination of the distance to the $\sigma$ Orionis cluster, the interpretation of the strong X-ray emission detected for $\sigma$ Ori Aa,Ab,B, and the investigation of the formation and evolution of multiple massive stellar systems and substellar objects.

Context: Massive stars are extremely important for the evolution of the galaxies; there are large gaps in our understanding of their properties and formation, however, mainly because they evolve rapidly, are rare, and distant. It may well be that almost all massive stars are born as triples or higher multiples, but their large distances require very high angular resolution to directly detect the companions at milliarcsecond scales. Aims: Herschel36 is a young massive system located at 1.3 kpc. It has a combined smallest predicted mass of 45 M_sun. Multi-epoch spectroscopic data suggest the existence of at least three gravitationally bound components. Two of them, system Ab, are tightly bound in a spectroscopic binary, and the third one, component Aa, orbits in a wider orbit. Our aim was to image and obtain astrometric and photometric measurements of components Aa and Ab using, for the first time, long-baseline optical interferometry to further constrain its nature. Methods: We observed Herschel 36 with the near-infrared instrument AMBER attached to the ESO VLT Interferometer, which provides an angular resolution of approx. 2 mas. We used the code BSMEM to perform the interferometric image reconstruction. We fitted the interferometric observables using proprietary IDL routines and the code LitPro. Results: We imaged the Aa+Ab components of Herschel 36 in H and K filters. Component Ab is located at a projected distance of 1.81 mas, at a position angle of approx. 222 deg. east of north, the flux ratio between components Aa and Ab is close to one. The small measured angular separation indicates that system Ab and Ab may be approaching the periastron of their orbits. These results, only achievable with long-baseline near-infrared interferometry, constitute the first step toward a thorough understanding of this massive triple system.

There are an enigmatic population of massive stars around the Galactic Center (GC) that were formed some Ma ago. A fraction of these stars has been found to orbit the supermassive black hole, SgrA*, in a projected clockwise disk, which suggests that they were formed in a formerly existing dense disk around SgrA*. We focus on the extended, near-infrared (NIR) sources IRS1W, IRS5, IRS10W, and IRS21 that have been suggested to be young, massive stars that form bow-shocks through their interaction with the ISM. Their nature has impeded accurate determination of their orbital parameters. We aim at establishing their nature and kinematics to test whether they form part of the clockwise disk. We performed NIR multi-wavelength imaging using adaptive optics (AO) and sparse aperture masking (SAM). We introduce a new method for self-calibration of the SAM PSF in dense stellar fields. The emission mechanism, morphology and kinematics of the targets were examined via 3D bow-shock models. We confirm previous findings that IRS21, IRS1W, and IRS5 are bow-shocks created by the interaction between mass-losing stars and the interstellar gas. The nature of IRS10W remains unclear. Our modeling shows that the bow-shock-emission is caused by thermal emission while the scattering of stellar light does not play any significant role. IRS 1W appears to be a bow-shock produced by an anisotropic stellar wind or by locally inhomogeneous ISM density. Our best-fit models provide an estimate of the local proper motion of the ISM in the NA in agreement with the published models. Assuming that all of the sources are tied to SgrA*, their orbital planes were obtained via a Monte-Carlo simulation. Our orbital analysis suggests that they are not part of any of the clockwise disk. We thus add more evidence to recent findings that a large part of the massive stars show apparently random orbital orientations.

Massive stars are of fundamental importance for almost all aspects of astrophysics, but there still exist large gaps in our understanding of their properties and formation because they are rare and therefore distant. It has been found that most O-stars are multiples. HD 150 136 is the nearest system to Earth with >100 M_sol, and provides a unique opportunity to study an extremely massive system. Recently, evidence for the existence of a third component in HD 150 136, in addition to the tight spectroscopic binary that forms the main component, was found in spectroscopic observations. Our aim was to image and obtain astrometric and photometric measurements of this component using long baseline optical interferometry to further constrain the nature of this component. We observed HD150136 with the near-infrared instrument AMBER attached to the ESO VLT Interferometer. The recovered closure phases are robust to systematic errors and provide unique information on the source asymmetry. Therefore, they are of crucial relevance for both image reconstruction and model fitting of the source structure. The third component in HD 150 136 is clearly detected in the high-quality data from AMBER. It is located at a projected angular distance of 7.3 mas, or about 13 AU at the line-of-sight distance of HD 150 136, at a position angle of 209 degrees East of North, and has a flux ratio of 0.25 with respect to the inner binary. We resolved the third component of HD 150 136 in J, H and K filters. The luminosity and color of the tertiary agrees with the predictions and shows that it is also an O main-sequence star. The small measured angular separation indicates that the tertiary may be approaching the periastron of its orbit. These results, only achievable with long baseline near infrared interferometry, constitute the first step towards the understanding of the massive star formation mechanisms.

Sparse aperture masking (SAM) interferometry combined with Adaptive Optics (AO) is a technique that is uniquely suited to investigate structures near the diffraction limit of large telescopes. The strengths of the technique are a robust calibration of the Point Spread Function (PSF) while maintaining a relatively high dynamic range. We used SAM+AO observations to investigate the circumstellar environment of several bright sources with infrared excess in the central parsec of the Galaxy. For our observations, unstable atmospheric conditions as well as significant residuals after the background subtraction presented serious problems for the standard approach of calibrating SAM data via interspersed observations of reference stars. We circumvented these difficulties by constructing a synthesized calibrator directly from sources within the field-of-view. When observing crowded fields, this novel method can boost the efficiency of SAM observations because it renders interspersed calibrator observations unnecessary. Here, we presented the first NaCo/SAM images reconstructed using this method.