The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

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

Theory surrounding the origin of the dust-laden winds from evolved stars remains mired in controversy. Characterizing the formation loci and the dust distribution within approximately the first stellar radius above the surface is crucial for understanding the physics that underlie the mass-loss phenomenon. By exploiting interferometric polarimetry, we derive the fundamental parameters that govern the dust structure at the wind base of a red supergiant. We present near-infrared aperture-masking observations of Betelgeuse in polarimetric mode obtained with the NACO/SAMPol instrument. We used both parametric models and radiative transfer simulations to predict polarimetric differential visibility data and compared them to SPHERE/ZIMPOL measurements. Using a thin dust shell model, we report the discovery of a dust halo that is located at only 0.5 R$_{\star}$ above the photosphere (i.e. an inner radius of the dust halo of 1.5 R$_{\star}$). By fitting the data under the assumption of Mie scattering, we estimate the grain size and density for various dust species. By extrapolating to the visible wavelengths using radiative transfer simulations, we compare our model with SPHERE/ZIMPOL data and find that models based on dust mixtures that are dominated by forsterite are most favored. Such a close dusty atmosphere has profound implications for the dust formation mechanisms around red supergiants.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

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

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

The GRAVITY instrument has been commissioned on the VLTI during 2016 and is now available to the astronomical community. It is the first optical interferometer capable of observing sources as faint as magnitude 19 in K-band. This is possible thanks to the fringe tracker which compensates the differential piston based on measurements of a brighter off-axis astronomical reference source. The goal of this paper is to consign the main developments made in the context of the GRAVITY fringe tracker. This could serve as basis for future fringe tracking systems. The paper therefore covers all aspects of the fringe tracker, from hardware, to control software and on-sky observations. Special emphasis is placed on the interaction between the group delay controller and the phase delay controller. The group delay control loop is a simple but robust integrator. The phase delay controller is a state-space control loop based on an auto-regressive representation of the atmospheric and vibrational perturbations. A Kalman filter provides optimal determination of the state of the system. The fringe tracker shows good tracking performance on sources with coherent K magnitudes of 11 on the UTs and 9.5 on the ATs. It can track fringes with an SNR level of 1.5 per DIT, limited by photon and background noises. On the ATs, during good seeing conditions, the optical path delay residuals can be as low as 75 nm root mean square. On the UTs, the performance is limited to around 250 nm because of structural vibrations.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

On 6 Sept. 2017, an exceptional coronal mass ejection departed from the Sun toward the Earth and Uranus, whose magnetospheres are sensitive to the solar wind. The resulting interplanetary shock triggered geomagnetic storm and intense aurora at Earth the next day and was predicted by MHD models to reach Uranus around 10-11 Nov. This event provided a unique opportunity to investigate the auroral response of the asymmetric Uranian magnetosphere in its intermediate equinox-to-solstice configuration. Coordinated multi-spectral observations were acquired with the Hubble Space Telescope (HST) in the far-UV (FUV), with the Very Large Telescope (VLT) and Gemini North in the near-IR (NIR) and with Chandra in the X-ray domain. In this study, we focus on the analysis of NIR images obtained between 9 and 17 Nov. 2017 which are compared to one FUV image acquired on 11 Nov. The latter reveals a bright southern auroral spot in the H2 bands, which we use as a reference to locate auroral precipitations. The NIR images were aimed at mapping H3+ emission from the Uranian ionosphere and at updating the results built from a couple of pioneer images taken 25 years ago. These new high resolution images reveal H3+ from the whole disc although brighter near the southern pole, but show no evidence of localized auroral emission.

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

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

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

Knowledge of the Earth's atmospheric optical turbulence is critical for astronomical instrumentation. Not only does it enable performance verification and optimisation of existing systems but it is required for the design of future instruments. As a minimum this includes integrated astro-atmospheric parameters such as seeing, coherence time and isoplanatic angle, but for more sophisticated systems such as wide field adaptive optics enabled instrumentation the vertical structure of the turbulence is also required. Stereo-SCIDAR is a technique specifically designed to characterise the Earth's atmospheric turbulence with high altitude resolution and high sensitivity. Together with ESO, Durham University has commissioned a Stereo-SCIDAR instrument at Cerro Paranal, Chile, the site of the Very Large Telescope (VLT), and only 20~km from the site of the future Extremely Large Telescope (ELT). Here we provide results from the first 18 months of operation at ESO Paranal including instrument comparisons and atmospheric statistics. Based on a sample of 83 nights spread over 22 months covering all seasons, we find the median seeing to be 0.64" with 50% of the turbulence confined to an altitude below 2 km and 40% below 600 m. The median coherence time and isoplanatic angle are found as 4.18 ms and 1.75" respectively. A substantial campaign of inter-instrument comparison was also undertaken to assure the validity of the data. The Stereo-SCIDAR profiles (optical turbulence strength and velocity as a function of altitude) have been compared with the Surface-Layer SLODAR, MASS-DIMM and the ECMWF weather forecast model. The correlation coefficients are between 0.61 (isoplanatic angle) and 0.84 (seeing).

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.

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.

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.

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.

Context. The mass-loss mechanism of cool massive evolved stars is poorly understood. The proximity of Betelgeuse makes it an appealing target to study its atmosphere, map the shape of its envelope, and follow the structure of its wind from the photosphere out to the interstellar medium. Aims. A link is suspected between the powerful convective motions in Betelgeuse and its mass loss. We aim to constrain the spatial structure and temporal evolution of the convective pattern on the photosphere and to search for evidence of this link. Methods. We report new interferometric observations in the infrared H band using the VLTI/PIONIER instrument. We monitored the photosphere of Betelgeuse between 2012 January and 2014 November to look for evolutions that may trigger the outflow. Results. Our interferometric observations at low spatial frequencies are compatible with the presence of a hot spot on the photosphere that has a characteristic width of one stellar radius. It appears to be superposed on the smaller scale convective pattern. In the higher spatial frequency domain, we observe a significant difference between the observations and the predictions of 3D hydrodynamical simulations. Conclusions. We bring new evidence for the presence of a convective pattern in the photosphere of red supergiants. The inferred hot spot is probably the top of a giant convection cell although an asymmetric extension of the star cannot be excluded by these interferometric observations alone. The properties of the observed surface features show a stronger contrast and inhomogeneity as predicted by 3D radiative hydrodynamical simulations. We propose that the large observed feature is modifying the signature of the convective pattern at the surface of the star in a way that simulations cannot reproduce.

The bright southern star L2 Pup is a particularly prominent asymptotic giant branch (AGB) star, as its distance of 64 pc makes it the nearest of its type. We report new adaptive optics observations of L2 Pup at visible wavelengths with the SPHERE/ZIMPOL instrument of the VLT that confirm the presence of the circumstellar dust disk at high inclination discovered recently by Kervella et al. (2014b). The signature of the three-dimensional structure of the disk is clearly observed in the map of the degree of linear polarization pL. We identify the inner rim of the disk through its polarimetric signature at a radius of 6 AU from the AGB star. The ZIMPOL intensity images in the V and R bands also reveal a close-in secondary source at a projected separation of 2 AU from the primary. The identification of the spectral type of this companion is uncertain due to the strong reddening from the disk, but its photometry suggests that it is a late K giant, of comparable mass to the AGB star. We present refined physical parameters for the dust disk derived using the RADMC-3D radiative transfer code. We also interpret the pL map using a simple polarization model to infer the three-dimensional structure of the envelope. Interactions between the inner binary system and the disk apparently form spiral structures that propagate along the orthogonal axis to the disk to form streamers. Two dust plumes propagating orthogonally to the disk are also detected. They originate in the inner stellar system, and are possibly related to the interaction of the wind of the two stars with the material in the disk. Based on the morphology of the envelope of L2 Pup, we propose that this star is at an early stage of the formation of a bipolar planetary nebula.

The physical mechanism through which the outgoing material of massive red supergiants is accelerated above the escape velocity is unclear. Thanks to the transparency of its circumstellar envelope, the nearby red supergiant Betelgeuse gives the opportunity to probe the innermost layers of the envelope of a typical red supergiant down to the photosphere, i.e. where the acceleration of the wind is expected to occur. We took advantage of the SPHERE/ZIMPOL adaptive optics imaging polarimeter to resolve the visible photosphere and close envelope of Betelgeuse. We detect an asymmetric gaseous envelope inside a radius of 2 to 3 times the near-infrared photospheric radius of the star (R*), and a significant Halpha emission mostly contained within 3 R*. From the polarimetric signal, we also identify the signature of dust scattering in an asymmetric and incomplete dust shell located at a similar radius. The presence of dust so close to the star may have a significant impact on the wind acceleration through radiative pressure on the grains. The 3 R* radius emerges as a major interface between the hot gaseous and dusty envelopes. The detected asymmetries strengthen previous indications that the mass loss of Betelgeuse is likely tied to the vigorous convective motions in its atmosphere.

The mass-loss process in Mira stars probably occurs in an asymmetric way where dust can form in inhomogeneous circumstellar molecular clumps. Following asymmetries along the pulsation cycle can give us clues about these mass-loss processes. We imaged the Mira star X Hya and its environnement at different epochs to follow the evolution of the morphology in the continuum and in the molecular bands. We observed X Hya with AMBER in J-H-K at low resolution at two epochs. We modelled squared visibilities with geometrical and physical models. We also present imaging reconstruction results obtained with MiRA and based on the physical a priori images. We report on the angular scale change of X Hya between the two epochs. 1D CODEX profiles allowed us to understand and model the spectral variation of squared visibilities and constrain the stellar parameters. Reconstructed model-dependent images enabled us to reproduce closure phase signals and the azimuthal dependence of squared visibilities. They show evidence for material inhomogeneities located in the immediate environment of the star.

We report the detection of seven low mass companions to intermediate-mass stars (SpT B/A/F; $M$$\approx$1.5-4.5 solar masses) in the Scorpius-Centaurus Association using nonredundant aperture masking interferometry. Our newly detected objects have contrasts $\Delta L'$$\approx$4-6, corresponding to masses as low as $\sim$20 Jupiter masses and mass ratios of $q$$\approx$0.01-0.08, depending on the assumed age of the target stars. With projected separations $\rho$$\approx$10-30 AU, our aperture masking detections sample an orbital region previously unprobed by conventional adaptive optics imaging of intermediate mass Scorpius-Centaurus stars covering much larger orbital radii ($\approx$30-3000 AU). At such orbital separations, these objects resemble higher mass versions of the directly imaged planetary mass companions to the 10-30 Myr, intermediate-mass stars HR 8799, $\beta$ Pictoris, and HD95086. These newly discovered companions span the brown dwarf desert, and their masses and orbital radii provide a new constraint on models of the formation of low-mass stellar and substellar companions to intermediate-mass stars.

Optical interferometry is a powerful tool to investigate the close environment of AGB stars. With a spatial resolution of a few milli-arcseconds, it is even possible to image directly the surface of angularly large objects. This is of special interest forMira stars and red supergiants for which the dust-wind is initiated from or very close to the photosphere by an interplay between pulsation and convection. Based on two-epoch interferometric observations of the Mira star X Hya, we present how the variation of the angular size with wavelength challenges pulsation models and how reconstructed images can reveal the evolution of the object shape and of its asymmetric structures.

The idea behind FIRST (Fibered Imager foR a Single Telescope) is to use single-mode fibers to combine multiple apertures in a pupil plane as such as to synthesize a bigger aperture. The advantages with respect to a pure imager are i) relaxed tolerance on the pointing and cophasing, ii) higher accuracy in phase measurement, and iii) availability of compact, precise, and active single-mode optics like Lithium Niobate. The latter point being a huge asset in the context of a space mission. One of the problems of DARWIN or SIM-like projects was the difficulty to find low cost pathfinders missions. But the fact that Lithium Niobate optic is small and compact makes it easy to test through small nanosats missions. Moreover, they are commonly used in the telecom industry, and have already been tested on communication satellites. The idea of the FIRST-S demonstrator is to spatialize a 3U CubeSat with a Lithium Niobate nulling interferometer. The technical challenges of the project are: star tracking, beam combination, and nulling capabilities. The optical baseline of the interferometer would be 30 cm, giving a 2.2 AU spatial resolution at distance of 10 pc. The scientific objective of this mission would be to study the visible emission of exozodiacal light in the habitable zone around the closest stars.

The stellar population in the central parsec of the Galaxy is dominated by an old (several Gyr) population, but young, massive stars dominate the luminosity function. We have studied the most luminous of these stars, GCIRS 7, in order to constrain the age of the recent star formation event in the Galactic Centre and to characterise it as an interferometric reference for observations of the Galactic Centre with the instrument GRAVITY, which will equip the Very Large Telescope Interferometer in the near future. We present the first H-band interferometric observations of GCIRS 7, obtained using the PIONIER visitor instrument on the VLTI using the four 8.2-m unit telescopes. In addition, we present unpublished K-band VLTI/AMBER data, build JHKL light-curves based on data spanning 4 decades, and measured the star's effective temperature using SINFONI spectroscopy. GCIRS 7 is marginally resolved at H-band (in 2013: uniform-disk diameter=1.076+/-0.093mas, R=960+/-92Rsun at 8.33+/-0.35kpc). We detect a significant circumstellar contribution at K-band. The star and its environment are variable in brightness and in size. The photospheric H-band variations are well modelled with two periods: P0~470+/-10 days (amplitude ~0.64mag) and long secondary period LSP~2700-2850 days (~1.1mag). As measured from CO equivalent width, <Teff>=3600+/-195K. The size, periods, luminosity (<Mbol>=-8.44+/-0.22) and effective temperature are consistent with an M1 supergiant with an initial mass of 22.5+/-2.5Msun and an age of 6.5-10Myr (depending on rotation). This age is in remarkable agreement with most estimates for the recent star formation event in the central parsec. Caution should be taken when using this star as an interferometric reference as it is variable in size, is surrounded by a variable circumstellar environment and large convection cells may form on its photosphere.

As the nearest known AGB star (d=64pc) and one of the brightest (mK-2), L2 Pup is a particularly interesting benchmark object to monitor the final stages of stellar evolution. We report new lucky imaging observations of this star with the VLT/NACO adaptive optics system in twelve narrow band filters covering the 1.0-4.0 microns wavelength range. These diffraction limited images reveal an extended circumstellar dust lane in front of the star, that exhibits a high opacity in the J band and becomes translucent in the H and K bands. In the L band, extended thermal emission from the dust is detected. We reproduce these observations using Monte-Carlo radiative transfer modeling of a dust disk with the RADMC-3D code. We also present new interferometric observations with the VLTI/VINCI and MIDI instruments. We measure in the K band an upper limit to the limb-darkened angular diameter of theta_LD = 17.9 +/- 1.6 mas, converting to a maximum linear radius of R = 123 +/- 14 Rsun. Considering the geometry of the extended K band emission in the NACO images, this upper limit is probably close to the actual angular diameter of the star. The position of L2 Pup in the Herzsprung-Russell diagram indicates that this star has a mass around 2 Msun and is probably experiencing an early stage of the asymptotic giant branch. We do not detect any stellar companion of L2 Pup in our adaptive optics and interferometric observations, and we attribute its apparent astrometric wobble in the Hipparcos data to variable lighting effects on its circumstellar material. We however do not exclude the presence of a binary companion, as the large loop structure extending to more than 10 AU to the North-East of the disk in our L band images may be the result of interaction between the stellar wind of L2 Pup and a hidden secondary object. The geometric configuration that we propose, with a large dust disk seen almost edge-on, appears particularly favorable to test and develop our understanding of the formation of bipolar nebulae.

Be stars exhibit variability for a great number of observables. Putting the pieces together of the disk dynamics is not an easy task and requires arduous modeling before achieving a good fit of the observational data. In order to guide the modeling process and make it more efficient, it is very instructive to investigate reference dynamical cases. This paper focuses on continuum polarimetric quantities and is the second of a series that aims to demonstrate the capacity of deriving the dynamical history and fundamental parameters of a classical Be star through the follow-up of various observables. After a detailed study of the different opacities at play in the formation of polarized spectra, we investigate predictions of polarimetric observables in the continuum for different dynamical scenarios. Our models are based on a coupling of a hydrodynamic viscous decretion simulations in a disk and a 3-D non-LTE radiative transfer code. Introducing the polarization color diagram (PCD), we show that certain combinations of polarimetric observables exhibit features that are characteristic of a mass loss history. This diagram also enables to estimate fundamental parameters such as the inclination angle, the disk density scale and the $\alpha$ viscous diffusion parameter. We present the PCD as a powerful diagnosis tool to track the dynamical phases of a Be star such as disk building-up, dissipation, periodic and episodic outbursts. Finally we confront our models with observations of 4 Be stars that exhibited long-term polarimetric activity.

Recent observational and theoretical studies of classical Be stars have established the utility of polarization color diagrams (PCD) in helping to constrain the time-dependent mass decretion rates of these systems. We expand on our pilot observational study of this phenomenon, and report the detailed analysis of a long-term (1989-2004) spectropolarimetric survey of 9 additional classical Be stars, including systems exhibiting evidence of partial disk-loss/disk-growth episodes as well as systems exhibiting long-term stable disks. After carefully characterizing and removing the interstellar polarization along the line of sight to each of these targets, we analyze their intrinsic polarization behavior. We find that many steady-state Be disks pause at the top of the PCD, as predicted by theory. We also observe sharp declines in the Balmer jump polarization for later spectral type, near edge-on steady-state disks, again as recently predicted by theory, likely caused when the base density of the disk is very high, and the outer region of the edge-on disk starts to self absorb a significant number of Balmer jump photons. The intrinsic $V$-band polarization and polarization position angle of $\gamma$ Cas exhibits variations that seem to phase with the orbital period of a known one-armed density structure in this disk, similar to the theoretical predictions of Halonen & Jones. We also observe stochastic jumps in the intrinsic polarization across the Balmer jump of several known Be+sdO systems, and speculate that the thermal inflation of part of the outer region of these disks could be responsible for producing this observational phenomenon. Finally, we estimate the base densities of this sample of stars to be between $\approx 8\times 10^{-11}$ to $\approx 4 \times 10^{-12}\,\rm g cm^{-3}$ during quasi steady state periods given there maximum observed polarization.

Recent improvements on the sensitivity and spectral resolution of X-ray observations have led to a better understanding of the properties of matter in the vicinity of High Mass X-ray Binaries hosting a supergiant star and a compact object. However the geometry and physical properties of their environment at larger scales are currently only predicted by simulations. We aim at exploring the environment of Vela X-1 at a few stellar radii of the supergiant using spatially resolved observations in the near-infrared and at studying its dynamical evolution along the 9-day orbital period of the system. We observed Vela X-1 in 2010 and 2012 using long baseline interferometry at VLTI, respectively with the AMBER instrument in the K band and the PIONIER instrument in the H band. The PIONIER observations span through one orbital period to monitor possible evolutions in the geometry of the system. We resolved a structure of $8\pm3~R_\star$ from the AMBER data and $2.0\,_{-1.2}^{+0.7}~R_\star$ from the PIONIER data. From the closure phase we found that the environment of Vela X-1 is symmetrical. We observed comparable measurements between the continuum and the spectral lines in the K band, meaning that both emissions originate from the same forming region. From the monitoring of the system over one period in 2012, we found the signal to be constant with the orbital phase within the error bars. We propose three scenarios for the discrepancy between the two measurements: either there is a strong temperature gradient in the supergiant wind leading to a hot component much more compact than the cool part of the wind observed in the K band, or we observed a diffuse shell in 2010 possibly triggered by an off-state in the accretion rate of the pulsar that was dissolved in the interstellar medium in 2012, or the structure observed in the H band was the stellar photosphere instead of the supergiant wind.

We present the first multiplicity-dedicated long baseline optical interferometric survey of the Scorpius-Centaurus-Lupus-Crux association. We used the Sydney University Stellar Interferometer to undertake a survey for new companions to 58 Sco-Cen B- type stars and have detected 24 companions at separations ranging from 7-130mas, 14 of which are new detections. Furthermore, we use a Bayesian analysis and all available information in the literature to determine the multiplicity distribution of the 58 stars in our sample, showing that the companion frequency is F = 1.35 and the mass ratio distribution is best described as a power law with exponent equal to -0.46, agreeing with previous Sco-Cen high mass work and differing significantly from lower-mass stars in Tau-Aur. Based on our analysis, we estimate that among young B-type stars in moving groups, up to 23% are apparently single stars. This has strong implications for the understanding of high-mass star formation, which requires angular momentum dispersal through some mechanism such as formation of multiple systems.

Aims. We report on simultaneous observations and modeling of mid-infrared (MIR), near-infrared (NIR), and submillimeter (submm) emission of the source Sgr A* associated with the supermassive black hole at the center of our Galaxy. Our goal was to monitor the activity of Sgr A* at different wavelengths in order to constrain the emitting processes and gain insight into the nature of the close environment of Sgr A*. Methods. We used the MIR instrument VISIR in the BURST imaging mode, the adaptive optics assisted NIR camera NACO, and the sub-mm antenna APEX to monitor Sgr A* over several nights in July 2007. Results. The observations reveal remarkable variability in the NIR and sub-mm during the five nights of observation. No source was detected in the MIR, but we derived the lowest upper limit for a flare at 8.59 microns (22.4 mJy with A_8.59mu = 1.6+/- 0.5). This observational constraint makes us discard the observed NIR emission as coming from a thermal component emitting at sub-mm frequencies. Moreover, comparison of the sub-mm and NIR variability shows that the highest NIR fluxes (flares) are coincident with the lowest sub-mm levels of our five-night campaign involving three flares. We explain this behavior by a loss of electrons to the system and/or by a decrease in the magnetic field, as might conceivably occur in scenarios involving fast outflows and/or magnetic reconnection.

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard $\alpha$-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk build-up as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a 3-D non-LTE radiative transfer code. We show that lightcurves for different wavelengths are specific of a mass loss history, inclination angle and $\alpha$ viscosity parameter. The diagnostic potential of those lightcurves is also discussed.

Be stars possess gaseous circumstellar decretion disks, which are well described using standard $\alpha$-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter $\alpha=1.0\pm0.2$, corresponding to a mass injection rate $\dot{M}=(3.5\pm 1.3) \times 10^{-8}\ M_\sun\,\mathrm{yr}^{-1}$. Such a large value of $\alpha$ suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

Stellar limb-darkening is an important tool for constraining the properties of a stellar atmosphere. We present a novel method for relating the fundamental stellar parameters mass and radius to limb-darkening laws using grids of spherical model stellar atmospheres. This method is applied to interferometric observations of the red supergiant Betelgeuse, where an unique measure of the stellar mass is determined.

Context: Mass-loss occurring in red supergiants (RSGs) is a major contributor to the enrichment of the interstellar medium in dust and molecules. The physical mechanism of this mass loss is however relatively poorly known. Betelgeuse is the nearest RSG, and as such a prime object for high angular resolution observations of its surface (by interferometry) and close circumstellar environment. Aims: The goal of our program is to understand how the material expelled from Betelgeuse is transported from its surface to the interstellar medium, and how it evolves chemically in this process. Methods: We obtained diffraction-limited images of Betelgeuse and a PSF calibrator (Aldebaran) in six filters in the N band (7.76 to 12.81 mic) and two filters in the Q band (17.65 and 19.50 mic), using the VLT/VISIR instrument. Results: Our images show a bright, extended and complex circumstellar envelope at all wavelengths. It is particularly prominent longwards of 9-10 mic, pointing at the presence of O-rich dust, such as silicates or alumina. A partial circular shell is observed between 0.5 and 1.0" from the star, and could correspond to the inner radius of the dust envelope. Several knots and filamentary structures are identified in the nebula, one of the knots, located at a distance of 0.9" west of the star, being particularly bright and compact. Conclusions: The circumstellar envelope around Betelgeuse extends at least up to several tens of stellar radii. Its relatively high degree of clumpiness indicates an inhomogeneous spatial distribution of the material lost by the star. Its extension corresponds to an important intermediate scale, where most of the dust is probably formed, between the hot and compact gaseous envelope observed previously in the near infrared and the interstellar medium.

We analyze the intrinsic polarization of two classical Be stars in the process of losing their circumstellar disks via a Be to normal B star transition originally reported by Wisniewski et al. During each of five polarimetric outbursts which interrupt these disk-loss events, we find that the ratio of the polarization across the Balmer jump (BJ+/BJ-) versus the V-band polarization traces a distinct loop structure as a function of time. Since the polarization change across the Balmer jump is a tracer of the innermost disk density whereas the V-band polarization is a tracer of the total scattering mass of the disk, we suggest such correlated loop structures in Balmer jump-V band polarization diagrams (BJV diagrams) provide a unique diagnostic of the radial distribution of mass within Be disks. We use the 3-D Monte Carlo radiation transfer code HDUST to reproduce the observed clockwise loops simply by turning "on/off" the mass decretion from the disk. We speculate that counter-clockwise loop structures we observe in BJV diagrams might be caused by the mass decretion rate changing between subsequent "on/off" sequences. Applying this new diagnostic to a larger sample of Be disk systems will provide insight into the time-dependent nature of each system's stellar decretion rate.

In Wisniewski et al. 2010, paper I, we analyzed 15 years of spectroscopic and spectropolarimetric data from the Ritter and Pine Bluff Observatories of 2 Be stars, 60 Cygni and {\pi} Aquarii, when a transition from Be to B star occurred. Here we anaylize the intrinsic polarization, where we observe loop-like structures caused by the rise and fall of the polarization Balmer Jump and continuum V-band polarization being mismatched temporally with polarimetric outbursts. We also see polarization angle deviations from the mean, reported in paper I, which may be indicative of warps in the disk, blobs injected at an inclined orbit, or spiral density waves. We show our ongoing efforts to model time dependent behavior of the disk to constrain the phenomena, using 3D Monte Carlo radiative transfer codes.

We present a novel theoretical tool to analyze the dynamical behaviour of a Be disk fed by non-constant decretion rates. It is mainly based on the computer code HDUST, a fully three-dimensional radiative transfer code that has been successfully applied to study several Be systems so far, and the SINGLEBE code that solves the 1D viscous diffusion problem. We have computed models of the temporal evolution of different types of Be star disks for different dynamical scenarios. By showing the behaviour of a large number of observables (interferometry, polarization, photometry and spectral line profiles), we show how it is possible to infer from observations some key dynamical parameters of the disk.

In the course of an imaging survey we have detected a visual companion to the intermediate-mass star HR 6037. In this letter, we present two epoch observations of the binary with NACO/VLT, and near-IR spectroscopy of the secondary with ISAAC/VLT. The NACO observations allow us to confirm HR 6037B as a co-moving companion. Its J and H band ISAAC spectra suggest the object has an spectral type of M9+-1, with a surface gravity intermediate between that of 10 Myr dwarfs and field dwarfs with identical spectral type. The comparison of its Ks-band photometry with evolutionary tracks allows us to derive a mass, effective temperature, and surface gravity of 62+-20 MJup, Teff = 2330+-200 K, and log g = 5.1+-0.2, respectively. The small mass ratio of the binary, -0.03, and its long orbital period, -5000 yr, makes HR 6037 a rare and uncommon binary system.

GRAVITY is an adaptive optics assisted Beam Combiner for the second generation VLTI instrumentation. The instrument will provide high-precision narrow-angle astrometry and phase-referenced interferometric imaging in the astronomical K-band for faint objects. We describe the wide range of science that will be tackled with this instrument, highlighting the unique capabilities of the VLTI in combination with GRAVITY. The most prominent goal is to observe highly relativistic motions of matter close to the event horizon of Sgr A*, the massive black hole at center of the Milky Way. We present the preliminary design that fulfils the requirements that follow from the key science drivers: It includes an integrated optics, 4-telescope, dual feed beam combiner operated in a cryogenic vessel; near-infrared wavefrontsensing adaptive optics; fringe-tracking on secondary sources within the field of view of the VLTI and a novel metrology concept. Simulations show that 10 {\mu}as astrometry within few minutes is feasible for a source with a magnitude of mK = 15 like Sgr A*, given the availability of suitable phase reference sources (mK = 10). Using the same setup, imaging of mK = 18 stellar sources in the interferometric field of view is possible, assuming a full night of observations and the corresponding UV coverage of the VLTI.

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.

This paper reports on H-band interferometric observations of Betelgeuse made at the three-telescope interferometer IOTA. We image Betelgeuse and its asymmetries to understand the spatial variation of the photosphere, including its diameter, limb darkening, effective temperature, surrounding brightness, and bright (or dark) star spots. We used different theoretical simulations of the photosphere and dusty environment to model the visibility data. We made images with parametric modeling and two image reconstruction algorithms: MIRA and WISARD. We measure an average limb-darkened diameter of 44.28 +/- 0.15 mas with linear and quadratic models and a Rosseland diameter of 45.03 +/- 0.12 mas with a MARCS model. These measurements lead us to derive an updated effective temperature of 3600 +/- 66 K. We detect a fully-resolved environment to which the silicate dust shell is likely to contribute. By using two imaging reconstruction algorithms, we unveiled two bright spots on the surface of Betelgeuse. One spot has a diameter of about 11 mas and accounts for about 8.5% of the total flux. The second one is unresolved (diameter < 9 mas) with 4.5% of the total flux. Resolved images of Betelgeuse in the H band are asymmetric at the level of a few percent. The MOLsphere is not detected in this wavelength range. The amount of measured limb-darkening is in good agreement with model predictions. The two spots imaged at the surface of the star are potential signatures of convective cells.

We present infrared interferometric imaging of the S-type Mira star Chi Cygni. The object was observed at four different epochs in 2005-2006 with the IOTA optical interferometer (H band). Images show up to 40% variation in the stellar diameter, as well as significant changes in the limb darkening and stellar inhomogeneities. Model fitting gave precise time-dependent values of the stellar diameter, and reveals presence and displacement of a warm molecular layer. The star radius, corrected for limb darkening, has a mean value of 12.1 mas and shows a 5.1mas amplitude pulsation. Minimum diameter was observed at phase 0.94+/-0.01. Maximum temperature was observed several days later at phase 1.02+/-0.02. We also show that combining the angular acceleration of the molecular layer with CO (Delta v = 3) radial velocity measurements yields a 5.9+/-1.5 mas parallax. The constant acceleration of the CO molecules -- during 80% of the pulsation cycle -- lead us to argument for a free-falling layer. The acceleration is compatible with a gravitational field produced by a 2.1(+1.5/-0.7) solar mass star. This last value is in agreement with fundamental mode pulsator models. We foresee increased development of techniques consisting in combining radial velocity with interferometric angular measurements, ultimately allowing total mapping of the speed, density, and position of the diverse species in pulsation driven atmospheres.

Context: Betelgeuse is one the largest stars in the sky in terms of angular diameter. Structures on the stellar photosphere have been detected in the visible and near-infrared as well as a compact molecular environment called the MOLsphere. Mid-infrared observations have revealed the nature of some of the molecules in the MOLsphere, some being the precursor of dust. Aims: Betelgeuse is an excellent candidate to understand the process of mass loss in red supergiants. Using diffraction-limited adaptive optics (AO) in the near-infrared, we probe the photosphere and close environment of Betelgeuse to study the wavelength dependence of its extension, and to search for asymmetries. Methods: We obtained AO images with the VLT/NACO instrument, taking advantage of the "cube" mode of the CONICA camera to record separately a large number of short-exposure frames. This allowed us to adopt a "lucky imaging" approach for the data reduction, and obtain diffraction-limited images over the spectral range 1.04-2.17 $\mu$m in 10 narrow-band filters. Results: In all filters, the photosphere of Betelgeuse appears partly resolved. We identify an asymmetric envelope around the star, with in particular a relatively bright "plume" extending in the southwestern quadrant up to a radius of approximately six times the photosphere. The CN molecule provides an excellent match to the 1.09 mic bandhead in absorption in front of the stellar photosphere, but the emission spectrum of the plume is more difficult to interpret. Conclusions: Our AO images show that the envelope surrounding Betelgeuse has a complex and irregular structure. We propose that the southwestern plume is linked either to the presence of a convective hot spot on the photosphere, or to the rotation of the star.

This paper reports measurements of Sgr A* made with NACO in L' -band (3.80 um), Ks-band (2.12 um) and H-band (1.66 um) and with VISIR in N-band (11.88 um) at the ESO VLT, as well as with XMM-Newton at X-ray (2-10 keV) wavelengths. On 4 April, 2007, a very bright flare was observed from Sgr A* simultaneously at L'-band and X-ray wavelengths. No emission was detected using VISIR. The resulting SED has a blue slope (beta > 0 for nuL_nu ~ nu^beta, consistent with nuL_nu ~ nu^0.4) between 12 micron and 3.8 micron. For the first time our high quality data allow a detailed comparison of infrared and X-ray light curves with a resolution of a few minutes. The IR and X-ray flares are simultaneous to within 3 minutes. However the IR flare lasts significantly longer than the X-ray flare (both before and after the X-ray peak) and prominent substructures in the 3.8 micron light curve are clearly not seen in the X-ray data. From the shortest timescale variations in the L'-band lightcurve we find that the flaring region must be no more than 1.2 R_S in size. The high X-ray to infrared flux ratio, blue nuL_nu slope MIR to L' -band, and the soft nuL_nu spectral index of the X-ray flare together place strong constraints on possible flare emission mechanisms. We find that it is quantitatively difficult to explain this bright X-ray flare with inverse Compton processes. A synchrotron emission scenario from an electron distribution with a cooling break is a more viable scenario.

Two of the three instruments proposed to ESO for the second generation instrumentation of the VLTI would use integrated optics for beam combination. Several design are studied, including co-axial and multi-axial recombination. An extensive quantity of combiners are therefore under test in our laboratories. We will present the various components, and the method used to validate and compare the different combiners. Finally, we will discuss the performances and their implication for both VSI and Gravity VLTI instruments.

We present the second-generation VLTI instrument GRAVITY, which currently is in the preliminary design phase. GRAVITY is specifically designed to observe highly relativistic motions of matter close to the event horizon of Sgr A*, the massive black hole at center of the Milky Way. We have identified the key design features needed to achieve this goal and present the resulting instrument concept. It includes an integrated optics, 4-telescope, dual feed beam combiner operated in a cryogenic vessel; near infrared wavefront sensing adaptive optics; fringe tracking on secondary sources within the field of view of the VLTI and a novel metrology concept. Simulations show that the planned design matches the scientific needs; in particular that 10 microarcsecond astrometry is feasible for a source with a magnitude of K=15 like Sgr A*, given the availability of suitable phase reference sources.

This paper is an H band interferometric examination of Arcturus, a star frequently used as a spatial and spectral calibrator. Using the IOTA 3 telescope interferometer, we performed spectro-interferometric observation (R~35) of Arcturus. Atmospheric models and prescriptions were fitted to the data to derive the brightness distribution of the photosphere. Image reconstruction was also obtained using two software algorithms: Wisard and Mira. An achromatic power law proved to be a good model of the brightness distribution, with a limb darkening compatible with the one derived from atmospheric model simulations using our Marcs model. A Rosseland diameter of 21.05+/-0.21 was derived, corresponding to an effective temperature of T_eff=4295+/-26 K. No companion was detected from the closure phases, with an upper limit on the brightness ratio of 8e-4 at 1AU. Dynamic range at such distance from the photosphere was established at 1.5e-4 (1sigma rms). An upper limit of 1.7e-3 was also derived for the level of brightness asymmetries present on the photosphere.

The multiplicity of early-type stars is still not well established. The studies of individual star forming regions suggest a connection with the age and the environment. To fill in this gap, we started the first detailed adaptive-optic-assisted imaging survey of 308 BA-type field stars within 300 pc from the Sun, to derive their multiplicity in a homogeneous way. Our first epoch observations yield 195 companion candidates around 117 sample targets. The second epoch observations are underway.