Context. The fate of a massive star during the latest stages of its evolution is highly dependent on its mass-loss rate/geometry and therefore knowing the geometry of the circumstellar material close to the star and its surroundings is crucial. Aims. We aim to study the nature (i.e. geometry, rates) of mass-loss episodes. In this context, yellow hypergiants are great targets. Methods. We analyse a large set of optical/near-infrared data, in spectroscopic and photometric (X-shooter/VLT), spectropolarimetric (ISIS/WHT), and interferometric GRAVITY-AMBER/VLTI) modes, toward the yellow hypergiant IRAS 17163-3907. We present the first model-independent reconstructed images of IRAS 17163-3907 at these wavelengths at milli-arcsecond scales. Lastly, we apply a 2D radiative transfer model to fit the dereddened photometry and the radial profiles of published VISIR images at 8.59 {\mu}m, 11.85 {\mu}m and 12.81 {\mu}m simultaneously, adopting the revised Gaia distance (DR2). Results. The interferometric observables around 2 {\mu}m show that the Br{\gamma} emission is more extended and asymmetric than the Na i and the continuum emission. In addition to the two known shells surrounding IRAS 17163-3907 we report on the existence of a third hot inner shell with a maximum dynamical age of only 30 yr. Conclusions. The interpretation of the presence of Na i emission at closer distances to the star compared to Br{\gamma} has been a challenge in various studies. We argue that the presence of a pseudophotosphere is not needed, but it is rather an optical depth effect. The three observed distinct mass-loss episodes are characterised by different mass-loss rates and can inform the theories on mass-loss mechanisms, which is a topic still under debate. We discuss these in the context of photospheric pulsations and wind bi-stability mechanisms.

Eta Car is one of the most luminous and massive stars in our Galaxy and is the brightest mid-infrared (mid-IR) source in the sky, outside our solar system. Since the late 1990s the central source has dramatically brightened at ultraviolet and optical wavelengths. This might be explained by a decrease in circumstellar dust extinction. We aim to establish the mid-IR flux evolution and further our understanding of the star's ultraviolet and optical brightening. Mid-IR images from $8-20~\mu$m were obtained in 2018 with VISIR at the Very Large Telescope. Archival data from 2003 and 2005 are retrieved from the ESO Science Archive Facility and historical records are collected from publications. We present the highest angular resolution mid-IR images of $\eta$ Car to date at the corresponding wavelengths ($\geq 0.22''$). We reconstruct the mid-IR evolution of the spectral energy distribution of the spatially integrated Homunculus nebula from 1968 to 2018 and find no long-term changes. Eta Car's bolometric luminosity has been stable over the past five decades. We do not observe a long-term decrease in the mid-IR flux densities that could be associated with the brightening at ultraviolet and optical wavelengths, but circumstellar dust must be declining in our line-of-sight only. Short-term flux variations within about 25% of the mean levels could be present.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

The formation of the most massive stars in the Universe remains an unsolved problem. Are they able to form in relative isolation in a manner similar to the formation of solar-type stars, or do they necessarily require a clustered environment? In order to shed light on this important question, we study the origin of two very massive stars (VMS): the O2.5If*/WN6 star RFS7 ($\sim$100 $M_{\odot}$), and the O3.5If* star RFS8 ($\sim$70 $M_{\odot}$), found within $\approx$ 53 and 58 pc respectively from the Galactic massive young cluster NGC 3603, using Gaia data. RFS7 is found to exhibit motions resembling a runaway star from NGC 3603. This is now the most massive runaway star candidate known in the Milky Way. Although RFS8 also appears to move away from the cluster core, it has proper-motion values that appear inconsistent with being a runaway from NGC 3603 at the $3\sigma$ level (but with substantial uncertainties due to distance and age). Furthermore, no evidence for a bow-shock or a cluster was found surrounding RFS8 from available near-infrared photometry. In summary, whilst RFS7 is likely a runaway star from NGC 3603, making it the first VMS runaway in the Milky Way, RFS8 is an extremely young ($\sim$2 Myr) VMS, which might also be a runaway, but this would need to be established from future spectroscopic and astrometric observations, as well as precise distances. If RFS8 were still not meeting the criteria for being a runaway from NGC 3603 from such future data, this would have important ramifications for current theories of massive star formation, as well as the way the stellar initial mass function (IMF) is sampled.

The rarity of young massive stars combined with the fact that they are often deeply embedded has limited the understanding of their formation. Ground based mid-infrared (IR) interferometry is one way of securing the spatial resolution required to study massive young stellar objects (MYSOs) and as the spatial-frequency coverage of such observations is often incomplete, direct-imaging can be supplementary to such a dataset. By consolidating these observations with modelling, the features of a massive protostellar environment can be constrained. This work simultaneously fits the aforementioned observations and a spectral energy distribution (SED) with a 2.5D radiative transfer model, providing an extensive view of the physical characteristics of the accreting regions of the MYSO G305.20+0.21. The high-resolution observations were obtained using the Very Large Telescope's MIDI and VISIR instruments, producing visibilities in the N-band and near-diffraction-limited imaging in the Q-band respectively. A model including a central protostar with a luminosity of $\sim$5$\times$10$^{4}$L$_{\odot}$ surrounded by a low-density bipolar cavity, a flared 1M$_{\odot}$ disk and an envelope provides a sufficient fit all three types of observation. The need to include a disk in the model implies that this MYSO follows a scaled-up version of the low-mass star formation process. The weak silicate absorption feature within the SED requires low-density envelope cavities to be successfully fit and is an atypical characteristic in comparison to previously studied MYSOs. Additionally, the inner radius of the disk must be three times the dust sublimation radius to satisfy the MIDI visibilities. The low density, low extinction environment implies the object is a more evolved MYSO and this combined with large inner radius of the disk suggests that it could be an example of a transitional disk around an MYSO.

Binarity and multiplicity appear to be a common outcome in star formation. In particular, the binary fraction of massive (OB-type) stars can be very high. In many cases, the further stellar evolution of these stars is affected by binary interactions at some stage during their lifetime. The origin of this high binarity and the binary parameters are poorly understood because observational constraints are scarce, which is predominantly due to a dearth of known young massive binary systems. We aim to identify and describe massive young binary systems in order to fill in the gaps of our knowledge of primordial binarity of massive stars, which is crucial for our understanding of massive star formation. We observed the two massive young stellar objects (MYSOs) PDS 27 and PDS 37 at the highest spatial resolution provided by VLTI/PIONIER in the H-band (1.3 mas). We applied geometrical models to fit the observed squared visibilities and closure phases. In addition, we performed a radial velocity analysis using published VLT/FORS2 spectropolarimetric and VLT/X-shooter spectroscopic observations. Our findings suggest binary companions for both objects at 12 mas (30 au) for PDS 27 and at 22-28 mas (42-54 au) for PDS 37. This means that they are among the closest MYSO binaries resolved to date. Our data spatially resolve PDS 27 and PDS 37 for the first time, revealing two of the closest and most massive ($>$8 M$_\odot$) YSO binary candidates to date. PDS 27 and PDS 37 are rare but great laboratories to quantitatively inform and test the theories on formation of such systems.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

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.

HH 50138 is one of the brightest B[e] stars at a distance of $\sim$ 380 pc with strong infrared excess. The star was observed in [O I] 63 $\mu$m and [C II] 158 $\mu$m with high velocity resolution with upGREAT on SOFIA. The velocity resolved [O I] emission provides evidence for a large gas-disk, $\sim$ 760 au in size, around HD 50138. Whereas previous interferometric observations give strong evidence for a hot gas and dust disk in Keplerian rotation, our bservations are the first to provide unambiguous evidence for a large warm disk around the star. Herschel/PACS observations showed that the [C II] emission is extended, therefore the [C II] emission most likely originates in an ionized gas shell created by a past outflow event. We confirm the isolated nature of HD 50138. It is far from any star forming region and has low proper motion. Neither is there any sign of a remnant cloud from which it could have formed. The extended disk around the star appears carbon poor. It shows OH and [O I] emission, but no CO. The CO abundance appears to be at least an order of magnitude lower than that of OH. Furthermore $^{13}$CO is enriched by more than a factor of five, confirming that the star is not a Herbig Be star. Finally we note that our high spectral resolution [O I] and [C II] observations provide a very accurate heliocentric velocity of the star, 40.8 $\pm$ 0.2 km~s$^{-1}$.

Using Gaia DR2 data, combined with OmegaCAM ground-based optical photometry from the AD-HOC survey, and detailed Radial Velocity measurements from ESO-Gaia, we analyse in detail a 10x5 deg region around the Wolf-Rayet star $\gamma^2$ Vel, including the previously known clusters Gamma Vel and NGC2547. Using clustering analysis that considers positions, proper motions and parallax, we discover 6 clusters or associations -- 4 of which appear new. Analysis of the colour-magnitude diagram for these clusters show that 4 of them formed coevally from the same molecular clouds 10 Myr ago, while NGC 2547 formed together with a newly discovered cluster 30 Myr ago. This study shows the incredible wealth of data provided by Gaia for the study of young stellar clusters.

Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Here we measure the variability of the oxygen-rich Mira star R Peg. We used K band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the VLTI. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.25 mu are steadily increasing with values of 8.7 mas, 9.4 mas, 9.8 mas, and 9.9 mas at visual phases of 0.15, 0.36, 0,45, 0.53, respectively. UD diameters at a bandpass around 2.05 mu, dominated by water vapour, follow the near-continuum variability at larger UD diameters. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

HD 100546 is one of the few known pre-main-sequence stars that may host a planetary system in its disk. We analyze new VLT/SPHERE/ZIMPOL polarimetric images of HD 100546 with filters in H$\alpha$ and the adjacent continuum. We have probed the disk gap and the surface layers of the outer disk, covering a region < 500 mas (< 55 au at 109 pc) from the star, at an angular resolution of ~ 20 mas. Our data show an asymmetry: the SE and NW regions of the outer disk are more polarized than the SW and NE. This can be explained from a preferential scattering angle close to 90$^o$, consistent with previous polarization images. The outer disk extends from 13 $\pm$ 2 to 45 $\pm$ 9 au, with a position angle and inclination of 137 $\pm$ 5$^o$ and 44 $\pm$ 8$^o$. The comparison with previous estimates suggests that the disk inclination could increase with the stellocentric distance, although the different measurements are still consistent within the error bars. In addition, no direct signature of the innermost candidate companion is detected from polarimetry, confirming recent results based on intensity imagery. We set an upper limit to its mass accretion rate < 10$^{-8}$ M$_{\odot}$/yr for a sub-stellar mass of 15M$_{Jup}$. Finally, we report the first detection (> 3$\sigma$) of a ~ 20 au bar-like structure that crosses the gap through the central region of HD 100546. It is tentatively suggested that the bar could be dust dragged by infalling gas that radially flows from the outer disk. This could represent an exceptional case in which a small-scale radial inflow is observed in a single system. If this scenario is confirmed, it could explain the presence of atomic gas in the inner disk that would otherwise accrete on to the central star on a timescale of a few months/years, as previously indicated from spectro-interferometric data, and could be related with additional (undetected) planets.

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.

As part of the Accretion Discs in H$\alpha$ with OmegaCAM (ADHOC) survey, we imaged in r, i and H-alpha a region of 12x8 square degrees around the Orion Nebula Cluster. Thanks to the high-quality photometry obtained, we discovered three well-separated pre-main sequences in the color-magnitude diagram. The populations are all concentrated towards the cluster's center. Although several explanations can be invoked to explain these sequences we are left with two competitive, but intriguing, scenarios: a population of unresolved binaries with an exotic mass ratio distribution or three populations with different ages. Independent high-resolution spectroscopy supports the presence of discrete episodes of star formation, each separated by about a million years. The stars from the two putative youngest populations rotate faster than the older ones, in agreement with the evolution of stellar rotation observed in pre-main sequence stars younger than 4 Myr in several star forming regions. Whatever the final explanation, our results prompt for a revised look at the formation mode and early evolution of stars in clusters.

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.

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41.

High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $\sim20$ and $\sim18 M_{\odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$\gamma$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L'$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $\sim3\times$ more massive and $\sim5\times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation.

We present a mid-IR interferometric survey of the circumstellar environment of a specific class of post-Asymptotic Giant Branch (post-AGB) binaries. For this class the presence of a compact dusty disk has been postulated on the basis of various spatially unresolved measurements. Our interferometric survey was performed with the MIDI instrument on the VLTI. In total 19 different systems were observed using variable baseline configurations. Combining all the visibilities at a single wavelength at 10.7 micron, we fitted two parametric models to the data: a uniform disk (UD) and a ring model mimicking a temperature gradient. We compared our observables of the whole sample, with synthetic data computed from a grid of radiative transfer models of passively irradiated disks in hydrostatic equilibrium. These models are computed with a Monte Carlo code that has been widely applied to describe the structure of protoplanetary disks around young stellar objects (YSO). The spatially resolved observations show that the majority of our targets cluster closely together in the distance-independent size-colour diagram, and have extremely compact N-band emission regions. The typical uniform disk diameter of the N-band emission region is about 40 mass which corresponds to a typical brightness temperature of 400-600~K. The resolved objects display very similar characteristics in the interferometric observables and in the spectral energy distributions. Therefore, the physical properties of the disks around our targets must be similar. The grid of protoplanetary disk models covers very well the observed objects. Much like for young stars, the spatially resolved N-band emission region is determined by the hot inner rim of the disk. Continued comparisons between post-AGB and protoplanetary disks will help to understand grain growth and disk evolution processes,

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.

Aims. The structural inhomogeneities and kinematics of massive star nebulae are tracers of their mass-loss history. We conduct a three-dimensional morpho-kinematic analysis of the ejecta of eta Car outside its famous Homunculus nebula. Methods. We carried out the first large-scale integral field unit observations of eta Car in the optical, covering a field of view of 1'x1' centered on the star. Observations with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) reveal the detailed three-dimensional structure of eta Car's outer ejecta. Morpho-kinematic modeling of these ejecta is conducted with the code SHAPE. Results. The largest coherent structure in eta Car's outer ejecta can be described as a bent cylinder with roughly the same symmetry axis as the Homunculus nebula. This large outer shell is interacting with the surrounding medium, creating soft X-ray emission. We establish the shape and extent of the ghost shell in front of the southern Homunculus lobe and confirm that the NN condensation can best be modeled as a bowshock in the orbital/equatorial plane. Conclusions. The SHAPE modeling of the MUSE observations indicates that the kinematics of the outer ejecta measured with MUSE can be described by a spatially coherent structure, and this structure also correlates with the extended soft X-ray emission associated with the outer debris field. The ghost shell just outside the southern Homunculus lobe hints at a sequence of eruptions within the time frame of the Great Eruption from 1837-1858 or possibly a later shock/reverse shock velocity separation. Our 3D morpho-kinematic modeling and the MUSE observations constitute an invaluable dataset to be confronted with future radiation-hydrodynamics simulations. Such a comparison may shed light on the yet elusive physical mechanism responsible for eta Car-like eruptions.

The upgrade of the VLTI infrastructure for the 2nd generation instruments is now complete with the transformation of the laboratory, and installation of star separators on both the 1.8-m Auxiliary Telescopes (ATs) and the 8-m Unit Telescopes (UTs). The Gravity fringe tracker has had a full semester of commissioning on the ATs, and a first look at the UTs. The CIAO infrared wavefront sensor is about to demonstrate its performance relative to the visible wavefront sensor MACAO. First astrometric measurements on the ATs and astrometric qualification of the UTs are on-going. Now is a good time to revisit the performance roadmap for VLTI that was initiated in 2014, which aimed at coherently driving the developments of the interferometer, and especially its performance, in support to the new generation of instruments: Gravity and MATISSE.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Luminous Blue Variables (LBVs) are massive stars caught in a post-main sequence phase, during which they are losing a significant amount of mass. As, on one hand, it is thought that the majority of massive stars are close binaries that will interact during their lifetime, and on the other, the most dramatic example of an LBV, Eta Car, is a binary, it would be useful to find other binary LBVs. We present here interferometric observations of the LBV HR Car done with the AMBER and PIONIER instruments attached to ESO's Very Large Telescope Interferometer (VLTI). Our observations, spanning two years, clearly reveal that HR Car is a binary star. It is not yet possible to constrain fully the orbit, and the orbital period may lie between a few years and several hundred years. We derive a radius for the primary in the system and possibly resolve as well the companion. The luminosity ratio in the H-band between the two components is changing with time, going from about 6 to 9. We also tentatively detect the presence of some background flux which remained at the 2% level until January 2016, but then increased to 6% in April 2016. Our AMBER results show that the emission line forming region of Br gamma is more extended than the continuum emitting region as seen by PIONIER and may indicate some wind-wind interaction. Most importantly, we constrain the total masses of both components, with the most likely range being 33.6 and 45 solar masses. Our results show that the LBV HR Car is possibly an Eta Car analog binary system with smaller masses, with variable components, and further monitoring of this object is definitively called for.

The operations model of the Paranal Observatory relies on the work of efficient staff to carry out all the daytime and nighttime tasks. This is highly dependent on adequate training. The Paranal Science Operations department (PSO) has a training group that devises a well-defined and continuously evolving training plan for new staff, in addition to broadening and reinforcing courses for the whole department. This paper presents the training activities for and by PSO, including recent astronomical and quality control training for operators, as well as adaptive optics and interferometry training of all staff. We also present some future plans.

To probe the circumstellar environment of IRAS 13481-6124, a 20 M_sun high-mass young stellar object (HMYSO) with a parsec-scale jet and accretion disc, we investigate the origin of its Br\gamma-emission line through NIR interferometry. We present the first AMBER/VLTI observations of the Br\gamma-emitting region in an HMYSO at R~1500. Our AMBER/VLTI observations reveal a spatially and spectrally resolved Br\gamma-line in emission with a strong P Cygni profile, indicating outflowing matter with a terminal velocity of ~500 km/s. Visibilities, differential phases, and closure phases are detected in our observations within the spectral line and in the adjacent continuum. Both total visibilities (continuum plus line emitting region) and pure-line visibilities indicate that the Br\gamma-emitting region is more compact (2-4 mas in diameter or ~6-13 au at 3.2 kpc) than the continuum-emitting region (~5.4 mas or ~17 au). The absorption feature is also spatially resolved at the longest baselines (81 and 85 m) and has a visibility that is slightly smaller than the continuum-emitting region. The differential phases at the four longest baselines display an \u2018S\u2019-shaped structure across the line, peaking in the blue- and red-shifted high-velocity components. The calibrated photocentre shifts are aligned with the known jet axis, i.e they are probably tracing an ionised jet. The high-velocity components (v_r~100-500 km/s) are located far from the source, whereas the low-velocity components (0-100 km/s) are observed to be closer, indicating a strong acceleration of the gas flow in the inner 10 au. Finally, a non-zero closure phase along the continuum is detected. By comparing our observations with the synthetic images of the continuum around 2.16 um, we confirm that this feature originates from the asymmetric brightness distribution of the continuum owing to the inclination of the inner disc.

We report on the OHANA interferometric snapshot survey, carried out by the VLTI group at the Paranal observatory. It makes use of observing time not useful for any other scheduled scientific or technical tasks in the sense of a backup programme, to characterize the mass-loss for early-type stars. The survey employs the combination of AMBER's high spectral and spatial resolution. The spatially unresolved central object provides a reference frame for the fringe properties observed in the light of the continuum.

This chapter discusses the use and possibilities of optical and infrared interferometry to study star formation. The chapter starts with a brief overview of the star formation process and highlights the open questions from an observational point of view. These are found at the smallest scales, as this is, inevitably, where all the action such as accretion and outflows, occurs. We then use basic astrophysical concepts to assess which scales and conditions can be probed with existing interferometric set-ups for which we use the ESO/VLTI instrument suite as example. We will concentrate on the more massive stars observed at high resolution with continuum interferometry. Throughout, some of the most recent interferometric results are used as examples of the various processes discussed.

We present new mid-infrared interferometric observations of the massive young stellar object IRAS 13481-6124, using VLTI/MIDI for spectrally-resolved, long-baseline measurements (projected baselines up to $\sim120$ m) and GSO/T-ReCS for aperture-masking interferometry in five narrow-band filters (projected baselines of $\sim1.8-6.4$ m) in the wavelength range of $7.5-13$ $\mu$m. We combine these measurements with previously-published interferometric observations in the $K$ and $N$ bands in order to assemble the largest collection of infrared interferometric observations for a massive YSO to date. Using a combination of geometric and radiative-transfer models, we confirm the detection at mid-infrared wavelengths of the disk previously inferred from near-infrared observations. We show that the outflow cavity is also detected at both near- and mid-infrared wavelengths, and in fact dominates the mid-infrared emission in terms of total flux. For the disk, we derive the inner radius ($\sim1.8$ mas or $\sim6.5$ AU at 3.6 kpc), temperature at the inner rim ($\sim1760$ K), inclination ($\sim48$ deg) and position angle ($\sim107$ deg). We determine that the mass of the disk cannot be constrained without high-resolution observations in the (sub-)millimeter regime or observations of the disk kinematics, and could be anywhere from $\sim10^{-3}$ to $20$ M$_\odot$. Finally, we discuss the prospects of interpreting the spectral energy distributions of deeply-embedded massive YSOs, and warn against attempting to infer disk properties from the SED.

Not all stars exhibiting the optical spectral characteristics of B[e] stars share the same evolutionary stage. The Galactic B[e] star MWC 137 is a prime example of an object with uncertain classification, with previous work suggesting pre- and post-main sequence classification. Our goal is to settle this debate and provide reliable evolutionary classification. Integral field spectrograph observations with VLT MUSE of the cluster SH 2-266 are used to analyze the nature of MWC 137. A collimated outflow is discovered that is geometrically centered on MWC 137. The central position of MWC 137 in the cluster SH 2-266 within the larger nebula suggests strongly that it is a member of this cluster and that it is both at the origin of the nebula and the newly discovered jet. Comparison of the color-magnitude diagram of the brightest cluster stars with stellar evolutionary models results in a distance of about 5.2$\pm$1.4 kpc. We estimate that the cluster is at least 3 Myr old. The jet extends over 66" (1.7 pc) projected on the plane of the sky, shows several knots, and projected velocities of up to $\pm$450 km s$^{-1}$. From the Balmer emission line decrement of the diffuse intracluster nebulosity we determine E(B-V)=1.4 mag for the inner 1' cluster region. The spectral energy distribution of the brightest cluster stars yield a slightly lower extinction of E(B-V)~1.2 mag. The extinction towards MWC 137 is estimated to be E(B-V)~1.8 mag (A$_V$~5.6 mag). Our analysis of the optical and near-infrared color-magnitude and color-color diagrams suggests a post-main sequence stage of MWC 137. The existence of a jet in this object implies the presence of an accretion disk.

An observational review is provided of the properties of accretion disks around young stars. It concerns the primordial disks of intermediate- and high-mass young stellar objects in embedded and optically revealed phases. The properties were derived from spatially resolved observations and therefore predominantly obtained with interferometric means, either in the radio/(sub)millimeter or in the optical/infrared wavelength regions. We make summaries and comparisons of the physical properties, kinematics, and dynamics of these circumstellar structures and delineate trends where possible. Amongst others, we report on a quadratic trend of mass accretion rates with mass from T Tauri stars to the highest mass young stellar objects and on the systematic difference in mass infall and accretion rates.

We present spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems.

Aims. Every 5.5 years eta Car's light curve and spectrum change remarkably across all observed wavelength bands. We compare the recent spectroscopic event in mid-2014 to the events in 2003 and 2009 and investigate long-term trends. Methods. Eta Car was observed with HST STIS, VLT UVES, and CTIO 1.5m CHIRON for a period of more than two years in 2012-2015. Archival observations with these instruments cover three orbital cycles. Results. Important spectroscopic diagnostics show significant changes in 2014 compared to previous events. While the timing of the first HeII 4686 flash was remarkably similar to previous events, the HeII equivalent widths were slightly larger and the line flux increased compared to 2003. The second HeII peak occurred at about the same phase as in 2009, but was stronger. The HeI line flux grew in 2009-2014 compared to 1998-2003. On the other hand, Halpha and FeII lines show the smallest emission strengths ever observed. Conclusions. The basic character of the spectroscopic events has changed in the past 2-3 cycles; ionizing UV radiation dramatically weakened during each pre-2014 event but not in 2014. The strengthening of HeI emission and the weakening of the lower-excitation wind features in our direct line of sight implies a substantial change in the physical parameters of the emitting regions. The polar spectrum at FOS4 shows less changes in the broad wind emission lines, which may be explained by the latitude-dependent wind structure of eta Car. The quick and strong recovery of the HeII emission in 2014 supports a scenario, in which the wind-wind shock may not have completely collapsed as was proposed for previous events. All this may be the consequence of just one elementary change, namely a strong decrease in the primary's mass-loss rate.

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.

Aims : The hydrogen-deficient supergiants known as R Coronae Borealis (RCB) stars might be the result of a double-degenerate merger of two white dwarfs (WDs), or a final helium shell flash in a planetary nebula central star. In this context, any information on the geometry of their circumstellar environment and, in particular, the potential detection of elongated structures, is of great importance. Methods : We obtained near-IR observations of V854 Cen with the AMBER recombiner located at the Very Large Telescope Interferometer (VLTI) array with the compact array (B$\leq$35m) in 2013 and the long array (B$\leq$140m) in 2014. At each time, V854 Cen was at maximum light. The $H$- and $K$-band continua were investigated by means of spectrally dependant geometric models. These data were supplemented with mid-IR VISIR/VLT images. Results : A dusty slightly elongated over density is discovered both in the $H$- and $K$-band images. With the compact array, the central star is unresolved ($\Theta\leq2.5$\,mas), but a flattened dusty environment of $8 \times 11$ mas is discovered whose flux increases from about $\sim$20% in the $H$ band to reach about $\sim$50% at 2.3$\micron$, which indicates hot (T$\sim$1500\,K) dust in the close vicinity of the star. The major axis is oriented at a position angle (P.A.) of 126$\pm$29$\deg$. Adding the long-array configuration dataset provides tighter constraints on the star diameter ($\Theta\leq1.0$ mas), a slight increase of the overdensity to $12 \times 15$ mas and a consistent P.A. of 133$\pm$49$\deg$. The closure phases, sensitive to asymmetries, are null and compatible with a centro-symmetric, unperturbed environment excluding point sources at the level of 3% of the total flux in 2013 and 2014. The VISIR images exhibit a flattened aspect ratio at the 15-20% level at larger distances ($\sim$1$\arcsec$) with a position angle of 92$\pm$19$\deg$, marginally consistent with the interferometric observations. Conclusions : This is the first time that a moderately elongated structure has been observed around an RCB star. These observations confirm the numerous suggestions for a bipolar structure proposed for this star in the literature, which were mainly based on polarimetric and spectroscopic observations.

VLTI/AMBER and VLTI/PIONIER observations of the LBV HR Car show an interferometric signature that could not possibly be explained by an extended wind, more or less symmetrically distributed around a single object. Instead, observations both in the Br$\gamma$ line and the H-band continuum are best explained by two point sources (or alternatively one point source and one slightly extended source) at about 2 mas separation and a contrast ratio of about 1:5. These observations establish that HR Car is a binary, but further interpretation will only be possible with future observations to constrain the orbit. Under the assumption that the current separation is close to the maximum one, the orbital period can be estimated to be of the order of 5 years, similar as in the $\eta$ Car system. This would make HR Car the second such LBV binary.

OHANA is an interferometric snapshot survey of the gaseous circumstellar environments of hot stars, carried out by the VLTI group at the Paranal observatory. It aims to characterize the mass-loss dynamics (winds/disks) at unexplored spatial scales for many stars. The survey employs the unique combination of AMBER's high spectral resolution with the unmatched spatial resolution provided by the VLTI. Because of the spatially unresolved central OBA-type star, with roughly neutral colour terms, their gaseous environments are among the easiest objects to be observed with AMBER, yet the extent and kinematics of the line emission regions are of high astrophysical interest.

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.

An overview is presented of the recent advances in understanding the B[e] phenomenon among blue supergiant stars in light of high-angular resolution observations and with an emphasis on the results obtained by means of long baseline optical stellar interferometry. The focus of the review is on the circumstellar material and evolutionary phase of B[e] supergiants, but recent results on dust production in regular blue supergiants are also highlighted.

Long-exposure stellar images recorded with large ground-based telescopes are blurred due to the turbulent nature of the atmosphere. The VLT employs active and adaptive optics (AO) systems to compensate for the deleterious effects of the atmosphere in real time. The speckle imaging technique provides an alternative way to achieve diffraction-limited imaging by post-processing a series of short-exposure images. The use of speckle imaging with the no-AO mode of NACO at the VLT is demonstrated. Application of this technique is particularly suited to the J-band and it provides versatile high angular resolution imaging under mediocre conditions and/or in imaging extended objects. The implementation of this mode underlines the continuing attractiveness of NACO at the VLT.

Aims. Eta Car's ultra-violet, optical, and X-ray light curves and its spectrum suggest a physical change in its stellar wind over the last decade. It was proposed that the mass-loss rate decreased by a factor of about 2 in the last 15 years. We complement these recent results by investigating the past evolution and the current state of eta Car in the near-infrared (IR). Methods. We present JHKL photometry of eta Car obtained at SAAO Sutherland from 2004-2013 with the Mk II photometer at the 0.75-m telescope and JHKs photometry with SIRIUS at the 1.4-m IRSF telescope from 2012-2013. The near-IR light curves since 1972 are analyzed. Results. The long-term brightening trends in eta Car's JHKL light curves were discontinuous around the 1998 periastron passage. After 1998, the star shows excess emission above the extrapolated trend from earlier dates, foremost in J and H, and the blueward, cyclical progression in its near-IR colors is accelerated. The near-IR color evolution is strongly correlated with the periastron passages. After correcting for the secular trend we find that the color evolution matches an apparent increase in blackbody temperature of an optically thick near-IR emitting plasma component from about 3500 to 6000 K over the last 20 years. Conclusions. We suggest that the changing near-IR emission may be caused by variability in optically thick bremsstrahlung emission. Periastron passages play a key role in the observed excess near-IR emission after 1998 and the long-term color evolution. We thus propose as a hypothesis that angular momentum transfer (via tidal acceleration) during periastron passages leads to sudden changes in eta Car's atmosphere resulting in a long-term decrease in the mass-loss rate.

We report the discovery of mid-infrared excess emission in the young object RZ Psc. The excess constitutes ~8% of its Lbol, and is well fit by a single 500K black-body implying a dust free region within 0.7AU for optically thick dust. The object displays dust obscuration events (UXOR behaviour) with a time-scale that suggests dusty material on orbits of 0.5AU. We also report a 12.4 year cyclical photometric variability which can be interpreted as due to perturbations in the dust distribution. The system is characterized by a high inclination, marginal extinction (during bright photometric states), a single temperature for the warm dust, and an age estimate which puts the star beyond the formation stage. We propose that the dust occultation events present a dynamical view of an active asteroid belt whose collisional products sporadically obscure the central star.

We are investigating the circumstellar material for a sample of B[e] stars using high spectral resolution data taken in the optical and near-infrared regions with ESO/FEROS and ESO/CRIRES spectrographs, respectively. B[e] stars are surrounded by dense disks of still unknown origin. While optical emission lines from [O I] and [Ca II] reflect the disk conditions close to the star (few stellar radii), the near-infrared data, especially the CO band emission, mirror the characteristics in the molecular part of the disk farther away from the star (several AU). Based on our high resolution spectroscopic data, we seek to derive the density and temperature structure of the disks, as well as their kinematics. This will allow us to obtain a better understanding of their structure, formation history and evolution. Here we present our preliminary results.

IRC +10420 is one of the few known massive stars in rapid transition from the Red Supergiant phase to the Wolf-Rayet or Luminous Blue Variable phase. The star has an ionised wind and using the Br gamma hydrogen recombination emission we assess the mass-loss on spatial scales of order 1 au. We present new VLT Interferometer AMBER data which are combined with all other AMBER data in the literature. The final dataset covers a position angle range of 180 degrees and baselines up to 110 meters. The spectrally dispersed visibilities, differential phases and line flux are conjointly analyzed and modelled. We also present AMBER/FINITO observations which cover a larger wavelength range and allow us to observe the Na I doublet at 2.2 micron. The data are complemented by X-Shooter data, which provide a higher spectral resolution view. The Brackett gamma line and the Na I doublet are both spatially resolved. After correcting the AMBER data for the fact that the lines are not spectrally resolved, we find that Br gamma traces a ring with a diameter of 4.18 milli-arcseconds. We consider a geometric model in which the Br gamma emission emerges from the top and bottom rings of an hour-glass shaped structure, viewed almost pole-on. It provides satisfactory fits to most visibilities and differential phases. The fact that we detect line emission from a neutral metal like Na I within the ionized region, a very unusual occurrence, suggests the presence of a dense pseudo-photosphere. The ionized wind can be reproduced with a polar wind, which could well have the shape of an hour-glass. The resolved Na I emission is found to occur on scales barely larger than the continuum. This fact and that many Yellow Hypergiants exhibit this comparatively rare emission hints at the presence of a "Yellow" or even "White Wall" in the Hertzsprung-Russell diagram, preventing them from visibly evolving to the blue.

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion.

The appearance of the B[e] phenomenon in evolved massive stars such as B[e] supergiants is still a mystery. While these stars are generally found to have disks that are cool and dense enough for efficient molecule and dust condensation, the origin of the disk material is still unclear. We aim at studying the kinematics and origin of the disk in the eccentric binary system GG Car, whose primary component is proposed to be a B[e] supergiant. Based on medium- and high-resolution near-infrared spectra we analyzed the CO-band emission detected from GG Car. The complete CO-band structure delivers information on the density and temperature of the emitting region, and the detectable 13CO bands allow us to constrain the evolutionary phase. In addition, the kinematics of the CO gas can be extracted from the shape of the first 12CO band head. We find that the CO gas is located in a ring surrounding the eccentric binary system, and its kinematics agrees with Keplerian rotation with a velocity, projected to the line of sight, of (80\pm 1) km/s. The CO ring has a column density of (5\pm 3)x10^21 cm^-2 and a temperature of 3200\pm 500 K. In addition, the material is chemically enriched in 13CO, which agrees with the primary component being slightly evolved off the main sequence. We discuss two possible scenarios for the origin of the circumbinary disk: (i) non-conservative Roche lobe overflow, and (ii) the possibility that the progenitor of the primary component could have been a classical Be star. Neither can be firmly excluded, but for Roche lobe overflow to occur, a combination of stellar and orbital parameter extremawould be required.

HD 327083 is a sgB[e] star that forms a binary system with an orbital semi-major axis of ~1.7 AU. Our previous observations using the VLTI and AMBER in the medium resolution K-band mode spatially resolved the environment of HD 327083. The continuum visibilities obtained indicate the presence of a circumbinary disc. CO bandhead emission was also observed. However, due to the limited spectral resolution of the previous observations, the kinematic structure of the emitting material was not constrained. In this paper, we address this and probe the source of the CO emission with high spectral resolution and spatial precision. We have observed HD 327083 with high spectral resolution (25 & 6 km/s) using AMBER and CRIRES. The observations are compared to kinematical models to constrain the source of the emission. It is shown that the CO bandhead emission can be reproduced using a model of a Keplerian disc with an inclination and size consistent with our previous VLTI observations. The model is compared to AMBER differential phase measurements, which have a precision as high as 30-micro-arcseconds. A differential phase signal corresponding to 0.15 milli-arcseconds (~5 sigma) is seen over the bandhead emission, which is in excellent agreement with the model that fits the CRIRES observations. In comparison, a model of an equatorial outflow, as envisaged in the standard sgB[e] scenario, does not reproduce the observations well. The excellent agreement between the disc model and observations in the spatial and spectral domains is compelling evidence that the CO bandhead emission of HD 327083 originates in a circumbinary Keplerian disc. In contrast, the model of an equatorial outflow cannot reproduce the observations well. This suggests that the standard sgB[e] scenario is not applicable to HD 327083, which supports the hypothesis that the B[e] behaviour of HD 327083 is due to binarity (ABRIDGED).

The closest examples of high-mass star birth occurs in deeply embedded environments at kiloparsec distances. Although much progress has been made, an observationally validated picture of the dominant processes which allows the central hydrostatic object to grow in mass has yet to be established. The observational technique of optical interferometry has demonstrated its potential in the field of high-mass star formation by delivering a milli-arcsecond infrared view on the complex accretion environment. We provide an overview of the scientific results obtained with multi-aperture telescope arrays and briefly discuss future instruments and their anticipated impact on our understanding of massive young stellar objects.

We present high angular resolution observations (0.5"x0.3") carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation.

HD 327083 is a luminous B type star which exhibits emission lines and an infrared excess and is therefore classified as a supergiant B[e] star. In addition, the star is the primary of a close binary system. It is not clear whether the B[e] behaviour of HD 327083 is related to its binarity or its evolutionary state. Here we address this issue by studying its circumstellar environment with high spatial resolution. To this end, we have observed HD 327083 with the VLTI and AMBER in the medium resolution K-band setting. 13CO bandhead emission is detected, confirming HD 327083 is a post-main sequence object. The observations spatially resolve the source of the NIR continuum and the Br-gamma and CO line emission. In addition, differential phase measurements allow us to probe the origin of the observed Br-gamma emission with sub-mas precision. Using geometrical models, we find that the visibilities and closure phases suggest that the close binary system is surrounded by a circum-binary disk. We also find that in the case of the binary HD 327083, the relative sizes of the continuum and Br-gamma emitting regions are different to those of a single supergiant B[e] star where the standard dual outflow scenario is thought to apply. These findings are consistent with the hypothesis that the mass loss of HD 327083 is related to its binary nature.

Knowledge of the mass function in open clusters constitutes one way to constrain the formation of low-mass stars and brown dwarfs as does the knowledge of the frequency of multiple systems and the properties of disks. The aim of the project is to determine the shape of the mass function in the low-mass and substellar regimes in the pre-main sequence (27 Myr) cluster IC4665, which is located at 350 pc from the Sun. We have cross-matched the near-infrared photometric data from the Eighth Data Release of the UKIRT Infrared Deep Sky Survey (UKIDSS) Galactic Clusters Survey with previous optical data obtained with the Canada-France-Hawaii wide-field camera to improve the determination of the luminosity and mass functions in the low-mass and substellar regimes. The availability of i and z photometry taken with the CFH12K camera on the Canada France Hawaii Telescope added strong constraints to the UKIDSS photometric selection in this cluster, which is located in a dense region of our Galaxy. We have derived the luminosity and mass functions of the cluster down to J=18.5 mag, corresponding to masses of ~0.025 Msun at the distance and age of IC4665 according to theoretical models. In addition, we have extracted new candidate members down to ~20 Jupiter masses in a previously unstudied region of the cluster. We have derived the mass function over the 0.6-0.04 Msun mass range and found that it is best represented by a log-normal function with a peak at 0.25-0.16 Msun, consistent with the determination in the Pleiades.

To investigate statistically whether magnetic fields in massive stars are ubiquitous or appear in stars with specific spectral classification, certain ages, or in a special environment, we acquired 41 new spectropolarimetric observations for 36 stars. Among the observed sample roughly half of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during three different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. To assess the membership in open clusters and associations, we used astrometric catalogues with the best currently available kinematic and photometric data. A field at a significance level of 3sigma was detected in ten O-type stars. Importantly, the largest longitudinal magnetic fields were measured in two Of?p stars: <B_z>=-381+-122G for CPD-282561 and <B_z>=-297+-62G for HD148937, previously detected by us as magnetic. The obtained observations of HD148937 on three different nights indicate that the magnetic field is slightly variable. Our new measurements support our previous conclusion that large-scale organized magnetic fields with polar field strengths in excess of 1kG are not widespread among O-type stars. Among the stars with a detected magnetic field, only one star, HD156154, belongs to an open cluster at high membership probability. According to previous kinematic studies, four magnetic O-type stars in the sample are well-known candidate runaway stars.

We apply mid-infrared spectro-interferometry to the massive young stellar object CRL2136. The observations were performed with the Very Large Telescope Interferometer and the MIDI instrument at a 42m baseline probing angular scales of 50 milli-arcseconds. We model the observed visibilities in parallel with diffraction-limited images at both 24.5micron and in the N-band (with resolutions of 0.6" and 0.3", respectively), as well as the spectral energy distribution. The arcsec-scale spatial information reveals the well-resolved emission from the dusty envelope. By simultaneously modelling the spatial and spectral data, we find that the bulk of the dust emission occurs at several dust sublimation radii (approximately 170 AU). This reproduces the high mid-infrared fluxes and at the same time the low visibilities observed in the MIDI data for wavelengths longward of 8.5micron. However, shortward of this wavelength the visibility data show a sharp up-turn indicative of compact emission. We discuss various potential sources of this emission. We exclude a dust disk being responsible for the observed spectral imprint on the visibilities. A cool supergiant star and an accretion disk are considered and both shown to be viable origins of the compact mid-infrared emission. We propose that CRL2136 is embedded in a dusty envelope, which truncates at several times the dust sublimation radius. A dust torus is manifest in the equatorial region. We find that the spectro-interferometric N-band signal can be reproduced by either a gaseous disk or a bloated central star. If the disk extends to the stellar surface, it accretes at a rate of 3.0 10^(-3) Msun/yr.

We present high resolution (R~50,000) spectroastrometry over the CO 1st overtone bandhead of a sample of seven intermediate/massive young stellar objects. These are primarily drawn from the red MSX source (RMS) survey, a systematic search for young massive stars which has returned a large, well selected sample of such objects. The mean luminosity of the sample is approximately 5 times 10^4 L_\odot, indicating the objects typically have a mass of ~15 solar masses. We fit the observed bandhead profiles with a model of a circumstellar disc, and find good agreement between the models and observations for all but one object. We compare the high angular precision (0.2-0.8 mas) spectroastrometric data to the spatial distribution of the emitting material in the best-fitting models. No spatial signatures of discs are detected, which is entirely consistent with the properties of the best-fitting models. Therefore, the observations suggest that the CO bandhead emission of massive young stellar objects originates in small-scale disks, in agreement with previous work. This provides further evidence that massive stars form via disc accretion, as suggested by recent simulations.

The circumstellar structure on 100 AU scales of the massive young stellar object W33A is probed using the VLTI and the MIDI instrument. N-band visibilities on 4 baselines are presented which are inconsistent with a spherically symmetric geometry. The visibility spectra and SED are simultaneously compared to 2D axi-symmetric dust radiative transfer models with a geometry including a rotationally flattened envelope and outflow cavities. We assume an O7.5 ZAMS star as the central source, consistent with the observed bolometric luminosity. The observations are also compared to models with and without (dusty and gaseous) accretion disks. A satisfactory model is constructed which reproduces the visibility spectra for each (u,v) point. It fits the silicate absorption, the mid-IR slope, the far-infrared peak, and the (sub)mm of the SED. It produces a 350 micron morphology consistent with observations. The 10 micron emission on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the outflow. The visibilities rule out the presence of dust disks with total (gas and dust) masses more than 0.01 Msun. However, optically thick accretion disks, interior to the dust sublimation radius, are allowed to accrete at rates equalling the envelope's mass infall rate (up to 10^(-3) Msun/yr) without substantially affecting the visibilities due to the extinction by the extremely massive envelope of W33A.

The Young Stellar Object (YSO) W33A is one of the best known examples of a massive star still in the process of forming. Here we present Gemini North ALTAIR/NIFS laser-guide star adaptive-optics assisted K-band integral-field spectroscopy of W33A and its inner reflection nebula. In our data we make the first detections of a rotationally-flattened outer envelope and fast bi-polar jet of a massive YSO at near-infrared wavelengths. The predominant spectral features observed are Br-gamma, H_2, and a combination of emission and absorption from CO gas. We perform a 3-D spectro-astrometric analysis of the line emission, the first study of its kind. We find that the object's Br-gamma emission reveals evidence for a fast bi-polar jet on sub-milliarcsecond scales, which is aligned with the larger-scale outflow. The hybrid CO features can be explained as a combination of hot CO emission arising in a disk close to the central star, while cold CO absorption originates in the cooler outer envelope. Kinematic analysis of these features reveals that both structures are rotating, and consistent with being aligned perpendicularly to both the ionised jet and the large-scale outflow. Assuming Keplerian rotation, we find that the circumstellar disk orbits a central mass of >10Msun, while the outer envelope encloses a mass of ~15Msun. Our results suggest a scenario of a central star accreting material from a circumstellar disk at the centre of a cool extended rotating torus, while driving a fast bi-polar wind. These results therefore provide strong supporting evidence for the hypothesis that the formation mechanism for high-mass stars is qualitatively similar to that of low-mass stars.

The imminent availability of large arrays of large light collectors deployed to exploit atmospheric Cherenkov radiation for gamma-ray astronomy at more than 100GeV, motivates the growing interest in application of intensity interferometry in astronomy. Indeed, planned arrays numbering up to one hundred telescopes will offer close to 5,000 baselines, ranging from less than 50m to more than 1000m. Recent and continuing signal processing technology developments reinforce this interest. Revisiting Stellar Intensity Interferometry for imaging is well motivated scientifically. It will fill the short wavelength (B/V bands) and high angular resolution (< 0.1mas) gap left open by amplitude interferometers. It would also constitute a first and important step toward exploiting quantum optics for astronomical observations, thus leading the way for future observatories. In this paper we outline science cases, technical approaches and schedule for an intensity interferometer to be constructed and operated in the visible using gamma-ray astronomy Air Cherenkov Telescopes as receivers.

We discuss VLTI AMBER and MIDI interferometry in addition to single-dish Subaru observations of massive young stellar objects. The observations probe linear size scales between 10 to 1000 AU for the average distance of our sources.

The space-time correlations of streams of photons can provide fundamentally new channels of information about the Universe. Today's astronomical observations essentially measure certain amplitude coherence functions produced by a source. The spatial correlations of wave fields has traditionally been exploited in Michelson-style amplitude interferometry. However the technology of the past was largely incapable of fine timing resolution and recording multiple beams. When time and space correlations are combined it is possible to achieve spectacular measurements that are impossible by any other means. Stellar intensity interferometry is ripe for development and is one of the few unexploited mechanisms to obtain potentially revolutionary new information in astronomy. As we discuss below, the modern use of stellar intensity interferometry can yield unprecedented measures of stellar diameters, binary stars, distance measures including Cepheids, rapidly rotating stars, pulsating stars, and short-time scale fluctuations that have never been measured before.

Massive young stellar objects (MYSO) are surrounded by massive dusty envelopes. Our aim is to establish their density structure on scales of ~1000 AU, i.e. a factor 10 increase in angular resolution compared to similar studies performed in the (sub)mm. We have obtained diffraction-limited (0.6") 24.5 micron images of 14 well-known massive star formation regions with Subaru/COMICS. The images reveal the presence of discrete MYSO sources which are resolved on arcsecond scales. For many sources, radiative transfer models are capable of satisfactorily reproducing the observations. They are described by density powerlaw distributions (n(r) ~ r^(-p)) with p = 1.0 +/-0.25. Such distributions are shallower than those found on larger scales probed with single-dish (sub)mm studies. Other sources have density laws that are shallower/steeper than p = 1.0 and there is evidence that these MYSOs are viewed near edge-on or near face-on, respectively. The images also reveal a diffuse component tracing somewhat larger scale structures, particularly visible in the regions S140, AFGL 2136, IRAS 20126+4104, Mon R2, and Cep A. We thus find a flattening of the MYSO envelope density law going from ~10 000 AU down to scales of ~1000 AU. We propose that this may be evidence of rotational support of the envelope (abridged).

In this poster contribution we highlight the equivalence between an Imaging Air Cherenkov Telescope (IACT) array and an Intensity Interferometer for a range of technical requirements. We touch on the differences between a Michelson and an Intensity Interferometer and give a brief overview of the current IACT arrays, their upgrades and next generation concepts (CTA, AGIS, completion 2015). The latter are foreseen to include 30-90 telescopes that will provide 400-4000 different baselines that range in length between 50m and a kilometre. Intensity interferometry with such arrays of telescopes attains 50 micro-arcseconds resolution for a limiting V magnitude of ~8.5. This technique opens the possibility of a wide range of studies, amongst others, probing the stellar surface activity and the dynamic AU scale circumstellar environment of stars in various crucial evolutionary stages. Here we discuss possibilities for using IACT arrays as optical Intensity Interferometers.

We present diffraction limited (0.6") 24.5micron Subaru/COMICS images of the red supergiant mu Cep. We report the detection of a circumstellar nebula, that was not detected at shorter wavelengths. It extends to a radius of at least 6" in the thermal infrared. On these angular scales, the nebula is roughly spherical, in contrast, it displays a pronounced asymmetric morphology closer in. We simultaneously model the azimuthally averaged intensity profile of the nebula and the observed spectral energy distribution using spherical dust radiative transfer models. The models indicate a constant mass-loss process over the past 1000 years, for mass-loss rates a few times 10^(-7) Msun/yr. This work supports the idea that at least part of the asymmetries in shells of evolved massive stars and supernovae may be due to the mass-loss process in the red supergiant phase.

We present interferometric and single-dish mid-infrared observations of a sample of massive young stellar objects (BN-type objects), using VLTI-MIDI (10 micron) and Subaru-COMICS (24.5 micron). We discuss the regions S140, Mon R2, M8E-IR, and W33A. The observations probe the inner regions of the dusty envelope at scales of 50 milli arcsecond and 0.6 arcsec (100-1000 AU), respectively. Simultaneous model fits to spectral energy distributions and spatial data are achieved using self-consistent spherical envelope modelling. We conclude that those MYSO envelopes that are best described by a spherical geometry, the commensurate density distribution is a powerlaw with index -1.0. Such a powerlaw is predicted if the envelope is supported by turbulence on the 100-1000AU scales probed with MIDI and COMICS, but the role of rotation at these spatial scales need testing.

The yellow hypergiants are found in a stage between the massive Red Supergiants and the Wolf-Rayet stars. This review addresses current issues concerning the evolution of massive stars, concentrating on the transitional post-Red Supergiant phase. Few yellow hypergiants are known and even fewer show direct evidence for having evolved off the Red Supergiant branch. Indeed, only two such rare objects with clear evidence for having gone through of a previous mass losing phase are known, IRC +10420 and HD 179821. We will review their properties and present recent results employing near-infrared interferometry, integral field spectroscopy and polarimetry. Finally, their real-time evolution is discussed.

The so-called lithium depletion boundary (LDB) provides a secure and independent tool for deriving the ages of young open clusters.In this context, our goal is to determine membership for a sample of 147 photometrically selected candidates of the young open cluster IC 4665 and to use confirmed members to establish an age based on the LDB. Employing the FLAMES multi-object spectrograph on VLT/UT2, we have obtained intermediate-resolution spectra of the cluster candidates. The spectra were used to measure radial velocities and to infer the presence of the Li I 670.8 nm doublet and Halpha emission. We have identified 39 bona fide cluster members based on radial velocity, Halpha emission, and Li absorption. The mean radial velocity of IC 4665 is found to be vrad=-15.95 +/- 1.13 km/s. Confirmed cluster members display a sharp transition in magnitude between stars with and without lithium, both in the Im vs. Im-z and in the Ks vs. Im-Ks diagrams.From this boundary, we deduce a cluster age of 27.7^(+4.2)_(-3.5) +/- 1.1 +/- 2 Myr. IC 4665 is the fifth cluster for which an LDB age has been determined, and it is the youngest cluster among these five. Thus, the LDB is established from relatively bright stars still in the contracting pre-main sequence phase. The mass of the boundary is M*=0.24 +/- 0.04 Msun. The LDB age agrees well with the ages derived from isochrone fitting of both low and high mass, turn-off stars, a result similar to what is found in the slightly older NGC 2547.

The paper investigates the milli-arcsecond scale structure of the present-day mass-loss of the post-Red Supergiant IRC+10420. We use three telescopes of the VLT Interferometer in combination with the AMBER near-infrared beam combiner to measure spectrally dispersed correlated fluxes in the K-band around the Br gamma transition. The resulting visibilities are compared to the predicted visibilities of emission structures with various simple models in order to infer the size of the observed emission region. The Br gamma line is resolved by VLTI+AMBER on all three baselines, with the maximum projected baseline extending 69 meter and a P.A. ranging between 10 and 30 degrees. A differential phase between line and continuum is detected on the longest baseline. The Br gamma emission region is found to have a diameter of 3.3 milli-arcseconds (FWHM), when compared to a Gaussian intensity distribution. A uniform disk and a ring-like intensity distribution do not fit the line visibilities. Comparing the AMBER equivalent width of Br gamma with measurements from various epochs, we find that the stellar photosphere contributes about 60% of the total continuum light at 2.2 micron. The remaining 40% continuum emission is found on scales larger than the 66mas AMBER field of view. Using simple arguments, and assuming optically thick line emission, we find that the line emitting region is elongated. We briefly discuss the possibilities whether such a structure is due to a bi-polar flow or a circumstellar disk. (Abridged).

We report on resolved interferometric observations with VLTI/MIDI of the massive young stellar object (MYSO) W33A. The MIDI observations deliver spectrally dispersed visibilities with values between 0.03 and 0.06, for a baseline of 45m over the wavelength range 8-13 micron. The visibilities indicate that W33A has a FWHM size of approximately 120AU (0.030'') at 8 micron which increases to 240AU at 13 micron, scales previously unexplored among MYSOs. This observed trend is consistent with the temperature falling off with distance. 1D dust radiative transfer models are simultaneously fit to the visibility spectrum, the strong silicate feature and the shape of the mid infrared spectral energy distribution (SED). For any powerlaw density distribution, we find that the sizes (as implied by the visibilities) and the stellar luminosity are incompatible. A reduction to a third of W33A's previously adopted luminosity is required to match the visibilities; such a reduction is consistent with new high resolution 70 micron data from Spitzer's MIPSGAL survey. We obtain best fits for models with shallow dust density distributions of r^(-0.5) and r^(-1.0) and for increased optical depth in the silicate feature produced by decreasing the ISM ratio of graphite to silicates and using optical grain properties by Ossenkopf et al. (1992).

The present generation of ground-based Very High Energy (VHE) gamma-ray observatories consist of arrays of up to four large (> 12m diameter) light collectors quite similar to those used by R. Hanbury Brown to measure stellar diameters by Intensity Interferometry in the late 60's. VHE gamma-ray observatories to be constructed over the coming decade will involve several tens of telescopes of similar or greater sizes. Used as intensity interferometers, they will provide hundreds of independent baselines. Now is the right time to re-assess the potential of intensity interferometry so that it can be taken into consideration in the design of these large facilities.

Intensity interferometry exploits a quantum optical effect in order to measure objects with extremely small angular scales. The first experiment to use this technique was the Narrabri intensity interferometer, which was successfully used in the 1970s to measure 32 stellar diameters at optical wavelengths; some as small as 0.4 milli-arcseconds. The advantage of this technique, in comparison with Michelson interferometers, is that it requires only relatively crude, but large, light collectors equipped with fast (nanosecond) photon detectors. Ground-based gamma-ray telescope arrays have similar specifications, and a number of these observatories are now operating worldwide, with more extensive installations planned for the future. These future instruments (CTA, AGIS, completion 2015) with 30-90 telescopes will provide 400-4000 different baselines that range in length between 50m and a kilometre. Intensity interferometry with such arrays of telescopes attains $50 \mu$-arcsecond resolution for a limiting visual magnitude ~8.5. Phase information can be extracted from the interferometric measurement with phase closure, allowing image reconstruction. This technique opens the possibility of a wide range of studies amongst others, probing the stellar surface activity and the dynamic AU scale circumstellar environment of stars in various crucial evolutionary stages. Here we focuse on the astrophysical potential of an intensity interferometer utilising planned new gamma-ray instrumentation.

We present an analysis of optical lightcurves of Small Magellanic Cloud Be-type stars. Observations show that (1) optical excess flux is correlated with near-IR excess flux indicating a similar mechanism and (2) the lightcurves can trace out ``loops'' in a colour-magnitude diagram. A simple model for the time dependent bound-free and free-free emission produced by an outflowing circumstellar disk gives reasonable fits to the observations.

(Abridged) Some 240 blue stars in the Small Magellanic Cloud are investigated on their fantastic irregular continuum variability. We report here two results regarding these stars. First, their optical flux excess is correlated to their near-IR flux excess, as determined from optical EROSII light curves and 2MASS measurements. Second, the relation between optical colour and magnitude is observed to be bi-valued in 40% of the cases, resulting in a "loop" when the light curve is presented in a colour-magnitude diagram. We argue that optical variability for a large fraction of the variable stars is due to variations in the amount of bound-free and free-free radiation. We do simple model calculations that allow us to interpret the observed colour-magnitude variability as due to an outflowing ionized circumstellar disk. The mass loss of the central star is variable, i.e. on or off. Once the star stops losing mass, the disk evolves naturally into a ring. The observed bi-valued colour-magnitude relation is the transition of a partially optically thick to a fully optically thin disk. Significantly, the loop is traversed clockwise by outflowing matter, but anti-clockwise by infalling matter. It is observed that the material is generally outflowing, but few cases of inflow are also observed.

We present the first deep, optical, wide-field imaging survey of the young open cluster Collinder 359, complemented by near-infrared follow-up observations. This study is part of a large programme aimed at examining the dependence of the mass function on environment and time. We have surveyed 1.6 square degrees in the cluster, in the I and z filters, with the CFH12K camera on the Canada-France-Hawaii 3.6-m telescope down to completeness and detection limits in both filters of 22.0 and 24.0 mag, respectively. Based on their location in the optical (I-z,I) colour-magnitude diagram, we have extracted new cluster member candidates in Collinder 359 spanning 1.3-0.03 Msun assuming an age of 60 Myr and a distance of 450 pc for the cluster. We have used the 2MASS database as well as our own near-infrared photometry to examine the membership status of the optically-selected cluster candidates. Comparison of the location of the most massive members in Collinder 359 in a (B-V,V) diagram with theoretical isochrones suggests that Collinder 359 is older than Alpha Per but younger than the Pleiades. We discuss the possible relationship between Collinder 359 and IC 4665 as both clusters harbour similar parameters, including proper motion, distance, and age.

We present a study of the young (30-100Myr) open cluster IC 4665 with the aim to determine the shape of the mass function well into the brown dwarf regime. We photometrically select 691 low-mass stellar and 94 brown dwarf candidate members over an area of 3.82 square degrees centred on the cluster. K-band follow-up photometry and Two-Micron All-Sky Survey data allow a first filtering of contaminant objects from our catalogues. A second filtering is performed for the brightest stars using proper motion data provided by the Tycho-2 and UCAC2 public catalogues. Contamination by the field population for the lowest mass objects is estimated using same latitude control fields. We fit observed surface densities of various cluster populations with King profiles and find a consistent tidal radius of 1.0 degree. The presence of possible mass segregation is discussed. In most respects investigated, IC 4665 is similar to other young open clusters at this age: (1) a power law fit to the mass function between 1 and 0.04Msun results in best fit for a slope of -0.6; (2) a cusp in the mass function is noticed at about the substellar boundary with respect to the power law description, the interpretation of which is discussed; (3) a fraction between 10-19% for BDs with M > 0.03Msun to total members; (4) a best-fit lognormal function to the full mass distribution shows an average member mass of 0.32Msun, if IC 4665 has an age of 50Myr.

The young stellar object MWC 297 has been observed with the VLT interferometer equipped with the AMBER instrument. MWC 297 has been spatially resolved in the continuum with a visibility of 0.50 as well as in the Br gamma emission line where the visibility decrease to a lower value of 0.33. This change in the visibility with the wavelength can be interpreted by the presence of an optically thick disk responsible for the visibility in the continuum and of a stellar wind traced by Br gamma and whose apparent size is 40% larger. We validate this interpretation by building a model of the stellar environment that combines a geometrically thin, optically thick accretion disk model consisting of gas and dust, and a latitude-dependent stellar wind outflowing above the disk surface. The continuum emission and visibilities obtained from this model are fully consistent with the interferometric AMBER data. They agree also with existing optical, near-infrared spectra and other broad-band near-infrared interferometric visibilities. A picture emerges in which MWC 297 is surrounded by an equatorial flat disk that is possibly still accreting and an outflowing wind which has a much higher velocity in the polar region than at the equator. [abridged]

The young stellar object MWC297 is a B1.5Ve star exhibiting strong hydrogen emission lines. This object has been observed by the AMBER/VLTI instrument in 2-telescope mode in a sub-region of the K spectral band centered around the Br gamma line at 2.1656 microns. The object has not only been resolved in the continuum with a visibility of 0.50+/-0.10, but also in the Br gamma line, where the flux is about twice larger, with a visibility about twice smaller (0.33+/-0.06). The continuum emission is consistent with the expectation of an optically thick thermal emission from dust in a circumstellar disk. The hydrogen emission can be understood by the emission of a halo above the disk surface. It can be modelled as a latitudinal-dependant wind model and it explains the width, the strength and the visibibility through the emission lines. The AMBER data associated with a high resolution ISAAC spectrum constrains the apparent size of the wind but also its kinematics.

In two papers we try to confirm that all Galactic high-mass stars are formed in a cluster environment, by excluding that O-type stars found in the Galactic field actually formed there. In de Wit et al. (2004) we presented deep K-band imaging of 5 arcmin fields centred on 43 massive O-type field stars that revealed that the large majority of these objects are single objects. In this contribution we explore the possibility that the field O stars are dynamically ejected from young clusters, by investigating their peculiar space velocity distribution, their distance from the Galactic plane, and their spatial vicinity to known young stellar clusters. We (re-)identify 22 field O-type stars as candidate runaway OB-stars. The statistics show that ~4% of all O-type stars with V < 8 can be considered as formed outside a cluster environment. Most are spectroscopically single objects, some are visual binaries. The derived percentage for O-type stars that form isolated in the field based on our statistical analyses is in agreement with what is expected from calculations adopting a universal cluster richness distribution with power index of beta = 1.7, assuming that the cluster richness distribution is continuous down to the smallest clusters containing one single star.

We present a study aimed at clarifying the birthplace for 43 massive O-type field stars. In this first paper we present the observational part: a search for stellar clusters near the target stars. We derive stellar density maps at two different resolving scales, viz. 0.25pc and 1.0pc from NTT and TNG imaging and the 2MASS catalogue. These scales are typical for cluster sizes. The main result is that the large majority of the O-type field population are isolated stars: only 12% (5 out of 43) of the O-type field stars is found to harbour a small-scale stellar cluster. We review the literature and aim at characterizing the stellar field of each O-type field star with the emphasis on star formation and the presence of known young stellar clusters. An analysis of the result of this paper and a discussion of the O-type field population as products of a dynamical ejection event is presented in an accompanying paper.

We report the discovery of 5 new Herbig Ae/Be candidate stars in the Small Magellanic Cloud in addition to the 2 reported in Beaulieu et al. (2001). We discuss these 7 HAeBe candidate stars in terms of (1) their irregular photometric variability, (2) their near infrared emission, (3) their Halpha emission and (4) their spectral type. One star has the typical photometric behaviour that is observed only among Pre-Main Sequence UX Orionis type stars. The objects are more luminous than Galactic HAeBe stars and Large Magellanic Cloud HAeBe candidates of the same spectral type. The stars were discovered in a systematic search for variable stars in a subset of the EROS2 database consisting of 115,612 stars in a field of 24x24 arcmin in the Small Magellanic Cloud. In total we discovered 504 variable stars. After classifying the different objects according to their type of variability, we concentrate on 7 blue objects with irregular photometric behaviour. We cross-identified these objects with emission line catalogues from Simbad and JHK photometry from 2MASS. The analysis is supplemented with obtained narrow and broad band imaging. We discuss their variability in terms of dust obscuration and bound-free and free-free emission. We estimate the influence of metallicity on the circumstellar dust emission from pre-main sequence stars.

We have performed intermediate resolution (R = 5000), high signal-to-noise H-band spectroscopy of a small, initial sample of three massive embedded young stellar objects (YSOs), using VLT/ISAAC. The sample has been selected from sources characterised in previous literature as being likely of OB spectral type, to be unambiguously associated with bright (H < 14) single point sources in the 2MASS database, and to have no optical counterparts. Of the targets observed, one object shows a ~B3 spectrum, similar to a main sequence object of the same spectral type. A second object exhibits weak HeI and H emission, indicating an early-type source: we detect HeII absorption, which supports a previous indirect derivation of the spectral type as mid-O. The third object does not show absorption lines, so no spectral type can de derived. It does, however, exhibit a rich spectrum of strong, broad emission lines and is likely to be surrounded by dense circumstellar material and at a very early evolutionary stage. Our results from this very small sample are in agreement with those of Kaper et al. (2002), who also find spectra similar to optically visible main sequence stars, together with emission line objects representing a very early evolutionary phase, in their much larger sample of K-band spectra.

We studied the HST-WFPC2 observations of the inner kpc of the interacting galaxy M51 in six bands from 2550 to 8140 Angstrom. The images show an oval shaped area ("bulge") of 11x16 arcsec or 450x650 pc around the nucleus, dominated by a smooth population of old stars with overimposed dust lanes. We found 30 bright point-like sources in the bulge of M51; many of these are bright in the UV. They are located in elongated "strings" which follow the general pattern of the dust lanes. The spectral energy distributions of the pointsources are compared with those of clusters or single stars. We present arguments that the pointsources are isolated massive stars or small groups of a few massive stars. The energy distributions are best fitted with those of single stars with 4000<Teff<50000 K, 4.2<log L/Lsun<7.2, and 12<Mstar<200 Msun. We have derived upper limits to the total mass of low mass stars that could be "hiding" within the point sources. For the bluest sources the upper limit is only a few hundred Msun. We conclude that massive stars are formed outside clusters (or in very low mass clusters) in the bulge of M51. The estimated star formation rate in the bulge is 1 to 2x10^(-3) Msun/yr. This star formation rate could be sustained for about 200 to 400 Myrs, suggesting that the ongoing massive star formation in the bulge of M51 is fed/triggered by the interaction with its companion about 400 Mrs ago. Theoretical cloud-models show that CO is destroyed in low extinction clouds near the bright nucleus of M51, so that H_2, [OI] and [CII] are the dominants coolants. This condition favors the formation of massive stars, in agreement with the observations. This mode of massive-star formation in the bulge of M51 may resemble the star formation in the early Universe, when the CO and dust contents were low due to the low metallicity.

Using data from the EROS2 microlensing survey, we report the discovery of two blue objects with irregular photometric behaviour of $\Delta V \sim 0.1$-0.4 mag on time scales of 20 to 200 days. They show a bluer when fainter behaviour. Subsequent spectra taken with the ESO 3.6m telescope show spectral type B4eIII and B2eIV-V with strong $H \alpha$ emission. These objects resemble the Herbig AeBe but also classical Be stars. At this stage, it is not possible to distinguish unambiguously between pre-main sequence and classical Be nature. If we favour the pre-main sequence interpretation, they are more luminous than the luminosity upper limit for Galactic HAeBe stars. The same was found for the HAeBe candidates in the LMC. This might be due to a shorter accretion time scale ($\tau = M_*/\dot{M}$), or the smaller dust content during the pre-main sequence evolution of SMC and LMC stars.

We will review some results based on the EROS and the on going EROS-2 microlensing surveys to search for dark matter in the Galactic halo via microlensing effects on LMC/SMC stars. Microlensing surveys provide systematic observations for millions of stars over a long period in both Magellanic Clouds and give birth to very powerful database. We will quickly review the results obtained on pulsating stars and we will detail the similarity and differences of the Herbig Ae/Be stars in the Magellanic Clouds compared to their Galactic counterpart.

In the early nineties several teams started large scale systematic surveys of the Magellanic Clouds and the Galactic Bulge to search for microlensing effects. As a by product, these groups have created enormous time-series databases of photometric measurements of stars with a temporal sampling duration and accuracy which are unprecedented. They provide the opportunity to test the accuracy of primary distance indicators, such as Cepheids, RRLyrae stars, the detached eclipsing binaries, or the luminosity of the red clump. We will review the contribution of the microlensing surveys to the understanding of the physics of the primary distance indicators, recent differential studies and direct distance determinations to the Magellanic Clouds and the Galactic Bulge.