Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

Double-lined spectroscopic binaries (SB2s) are one of the main sources of stellar masses, as additional observations are only needed to give the inclinations of the orbital planes in order to obtain the individual masses of the components. For this reason, we are observing a selection of SB2s using the SOPHIE spectrograph at the Haute-Provence observatory in order to precisely determine their orbital elements. Our objective is to finally obtain masses with an accuracy of the order of one percent by combining our radial velocity (RV) measurements and the astrometric measurements that will come from the Gaia satellite. We present here the RVs and the re-determined orbits of 10 SB2s. In order to verify the masses we will derive from Gaia, we obtained interferometric measurements of the ESO VLTI for one of these SB2s. Adding the interferometric or speckle measurements already published by us or by others for 4 other stars, we finally obtain the masses of the components of 5 binary stars, with masses ranging from 0.51 to 2.2 solar masses, including main-sequence dwarfs and some more evolved stars whose location in the HR diagram has been estimated.

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.

The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale.

Young stellar objects are thought to accrete material from their circumstellar disks through their strong stellar magnetospheres. We aim to directly probe the magnetospheric accretion region on a scale of a few 0.01 au in a young stellar system using long-baseline optical interferometry. We observed the pre-transitional disk system DoAr 44 with VLTI/GRAVITY on two consecutive nights in the K-band. We computed interferometric visibilities and phases in the continuum and in the BrG line in order to constrain the extent and geometry of the emitting regions. We resolve the continuum emission of the inner dusty disk and measure a half-flux radius of 0.14 au. We derive the inclination and position angle of the inner disk, which provides direct evidence that the inner and outer disks are misaligned in this pre-transitional system. This may account for the shadows previously detected in the outer disk. We show that BrG emission arises from an even more compact region than the inner disk, with an upper limit of 0.047 au (5 Rstar). Differential phase measurements between the BrG line and the continuum allow us to measure the astrometric displacement of the BrG line-emitting region relative to the continuum on a scale of a few tens of microarcsec, corresponding to a fraction of the stellar radius. Our results can be accounted for by a simple geometric model where the BrG line emission arises from a compact region interior to the inner disk edge, on a scale of a few stellar radii, fully consistent with the concept of magnetospheric accretion process in low-mass young stellar systems.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

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.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

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.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

Post-AGB binaries are surrounded by circumbinary disks of gas and dust that are similar to protoplanetary disks found around young stars. We aim to understand the structure of these disks and identify the physical phenomena at play in their very inner regions. We want to understand the disk-binary interaction and to further investigate the comparison with protoplanetary disks. We have conducted an interferometric snapshot survey of 23 post-AGB binaries in the near-infrared (H-band) using VLTI/PIONIER. We have fitted the multiwavelength visibilities and closure phases with purely geometrical models with an increasing complexity in order to retrieve the sizes, temperatures and flux ratios of the different components All sources are resolved and the different components contributing to the H-band flux are dissected. The environment of these targets is very complex: 13/23 targets need models with thirteen or more parameters to fit the data. We find that the inner disk rims follow and extend the size-luminosity relation established for disks around young stars with an offset toward larger sizes. The measured temperature of the near-infrared circumstellar emission of post-AGB binaries is lower (Tsub~1200K) than for young stars, probably due to a different dust mineralogy and/or gas density in the dust sublimation region. The dusty inner rims of the circumbinary disks around post-AGB binaries are ruled by dust sublimation physics. Additionally, a significant amount of the circumstellar H-band flux is over-resolved (14 targets have more than 10% of their non-stellar flux over-resolved) hinting for more structure from a yet unknown origin (disk structure or outflow). The amount of over-resolved flux is larger than around young stars. Due to the complexity of these targets, interferometric imaging is a necessary tool to reveal the interacting inner regions in a model-independent way.

$\varepsilon$ Lupi A is a binary system consisting of two main sequence early B-type stars Aa and Ab in a short period, moderately eccentric orbit. The close binary pair is the only doubly-magnetic massive binary currently known. Using photometric data from the BRITE-Constellation we identify a modest heartbeat variation. Combining the photometry with radial velocities of both components we determine a full orbital solution including empirical masses and radii. These results are compared with stellar evolution models as well as interferometry and the differences discussed. We also find additional photometric variability at several frequencies, finding it unlikely these frequencies can be caused by tidally excited oscillations. We do, however, determine that these signals are consistent with gravity mode pulsations typical for slowly pulsating B stars. Finally we discuss how the evolution of this system will be affected by magnetism, determining that tidal interactions will still be dominant.

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.

We implement an electron avalanche photodiode (e-APD) in the MIRC-X instrument, upgrade of the 6-telescope near-infrared imager MIRC, at the CHARA array. This technology should improve the sensitivity of near-infrared interferometry. We first used the classical Mean-Variance analysis to measure the system gain and the amplification gain. We then developed a physical model of the statistical distribution of the camera output signal. This model is based on multiple convolutions of the Poisson statistic, the intrinsic avalanche gain distribution, and the observed distribution of the background signal. At low flux level, this model constraints independently the incident illumination level, the total gain, and the excess noise factor of the amplification. We measure a total transmission of $48\pm3\%$ including the cold filter and the Quantum Efficiency. We measure a system gain of 0.49 ADU/e, a readout noise of $10$ ADU, and amplification gains as high as 200. These results are consistent between the two methods and therefore validate our modeling approach. The measured excess noise factor based on the modeling is $1.47\pm0.03$, with no obvious dependency with flux level or amplification gain. The presented model allows measuring the characteristics of the e-APD array at low flux level independently of preexisting calibration. With $<0.3$\,electron equivalent readout noise at kilohertz frame rates, we confirm the revolutionary performances of the camera with respect to the PICNIC or HAWAII technologies. However, the measured excess noise factor is significantly higher than the one claimed in the literature ($<$1.25), and explains why counting multiple photons remains challenging with this camera.

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 report on 13 new high-precision measurements of stellar diameters for low-mass dwarfs obtained by means of near-infrared long-baseline interferometry with PIONIER at the Very Large Telescope Interferometer. Together with accurate parallaxes from Gaia DR2, these measurements provide precise estimates for their linear radii, effective temperatures, masses, and luminosities. This allows us to refine the effective temperature scale, in particular towards the coolest M-dwarfs. We measure for late-type stars with enhanced metallicity slightly inflated radii, whereas for stars with decreased metallicity we measure smaller radii. We further show that Gaia DR2 effective temperatures for M-dwarfs are underestimated by $\sim$ 8.2 % and give an empirical $M_{G}$-$T_{\rm eff}$ relation which is better suited for M-dwarfs with $T_{\rm eff}$ between 2600 and 4000 K. Most importantly, we are able to observationally identify a discontinuity in the $T_{\rm eff}$-radius plane, which is likely due to the transition from partially convective M-dwarfs to the fully convective regime. We found this transition to happen between 3200 K and 3340 K, or equivalently for stars with masses $\approx 0.23 M_{\odot}$. We find that in this transition region the stellar radii are in the range from 0.18 to 0.42$R_{\odot}$ for similar stellar effective temperatures.

The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the stars HD 163296 and HD 190073, including 30 previously unpublished nights of observations of the former and 45 nights of the latter with the CHARA long-baseline interferometer, in addition to archival VLTI data. We confirm previous observations suggesting significant near-infrared emission originates within the putative dust evaporation front of HD 163296 and show this is the case for HD 190073 as well. The H- and K-band sizes are the same within $(3 \pm 3)\%$ for HD 163296 and within $(6 \pm 10)\%$ for HD 190073. The radial surface brightness profiles for both disks are remarkably Gaussian-like with little or no sign of the sharp edge expected for a dust evaporation front. Coupled with spectral energy distribution analysis, our direct measurements of the stellar flux component at H and K bands suggest that HD 190073 is much younger (<400 kyr) and more massive (~5.6 M$_\odot$) than previously thought, mainly as a consequence of the new Gaia distance (891 pc).

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

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.

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.

We present the design for MYSTIC, the Michigan Young STar Imager at CHARA. MYSTIC will be a K-band, cryogenic, 6-beam combiner for the Georgia State University CHARA telescope array. The design follows the image-plane combination scheme of the MIRC instrument where single-mode fibers bring starlight into a non-redundant fringe pattern to feed a spectrograph. Beams will be injected in polarization-maintaining fibers outside the cryogenic dewar and then be transported through a vacuum feedthrough into the ~220K cold volume where combination is achieved and the light is dispersed. We will use a C-RED One camera (First Light Imaging) based on the eAPD SAPHIRA detector to allow for near-photon-counting performance. We also intend to support a 4-telescope mode using a leftover integrated optics component designed for the VLTI-GRAVITY experiment, allowing better sensitivity for the faintest targets. Our primary science driver motivation is to image disks around young stars in order to better understand planet formation and how forming planets might influence disk structures.

We characterize a near-infrared C-RED ONE camera from First Light Imaging (FLI). This camera uses a SAPHIRA electron avalanche photo-diode array (e-APD) from Leonardo (previously Selex). To do so, we developed a model of the signal distribution. This model allows a measurement of the gain and the Excess Noise Factor (ENF) independently of preexisting calibration such as the system gain. The results of this study show a gain which is 0.53 +/- 0.04 times the gain reported by the manufacturer. The measured ENF is 1.47 +/- 0.03 when we expected 1.25. For an avalanche gain of 60 and a frame rate larger than 100 Hz, the total noise can be lower than 1 e-/frame/pixel. The lowest dark current level is 90e-/s/pixel, in agreement with the expected H-band background passing through the camera window. These performance values provide a significant improvement compared to earlier-generation PICNIC camera and allowed us to improve the performance of the Michigan infrared combiner (MIRC) instrument at the Center for High Angular Resolution Astronomy (CHARA), as part of our MIRC-X instrumentation project.

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.

Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5~$\mu$m) and advanced fringe tracking strategies. We also address the strong exoplanet science case (young exoplanets, planet formation, and exozodiacal disks) offered by this wavelength regime as well as other possible science cases such as stellar physics (fundamental parameters and multiplicity) and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.

MIRC-X is an upgrade of the six-telescope infrared beam combiner at the CHARA telescope array, the world's largest baseline interferometer in the optical/infrared, located at the Mount Wilson Observatory in Los Angeles. The upgraded instrument features an ultra-low noise and fast frame rate infrared camera (SAPHIRA detector) based on e-APD technology. We report the MIRC-X sensitivity upgrade work and first light results in detail focusing on the detector characteristics and software architecture.

MIRC-X is a new beam combination instrument at the CHARA array that enables 6-telescope interferometric imaging on object classes that until now have been out of reach for milliarcsecond-resolution imaging. As part of an instrumentation effort lead by the University of Exeter and University of Michigan, we equipped the MIRC instrument with an ultra-low read-noise detector system and extended the wavelength range to the J and H-band. The first phase of the MIRC-X commissioning was successfully completed in June 2017. In 2018 we will commission polarisation control to improve the visibility calibration and implement a 'cross-talk resiliant' mode that will minimise visibility cross-talk and enable exoplanet searches using precision closure phases. Here we outline our key science drivers and give an overview about our commissioning timeline. We comment on operational aspects, such as remote observing, and the prospects of co-phased parallel operations with the upcoming MYSTIC combiner.

The Very Large Telescope Interferometer Auxiliary Telescopes will soon be equipped with an adaptive optics system called NAOMI. The corrective optics deformable mirror is the commercial DM241 from ALPAO. Being part of an interferometer operating from visible to mid-infrared, the DMs of NAOMI face several challenges (high level of reliability, open-loop chopping, piston-free control, WFS/DM pupil rotation, high desired bandwidth and stroke). We here describe our extensive characterization of the DMs through measurements and simulations. We summarize the operational scenario we have defined to handle the specific mirror properties. We conclude that the ALPAO DMs have overall excellent properties that fulfill most of the stringent requirements and that deviations from specifications are easily handled. To our knowledge, NAOMI will be the first astronomical system with a command in true Zernike modes (allowing software rotation), and the first astronomical system in which a chopping is performed with the deformable mirror (5'' sky, at 5~Hz).

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10$^{-3}$ in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at least a few 10$^{-4}$ are within reach using specific beam combination and data acquisition techniques. In this paper, we explore the possibility to reach similar or higher contrasts on the VLTI. After reviewing the state-of-the-art in high-contrast infrared interferometry, we discuss key features that made the success of other high-contrast interferometric instruments (e.g., integrated optics, nulling, closure phase, and statistical data reduction) and address possible avenues to improve the contrast of the VLTI by at least one order of magnitude. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5 $\mu$m), a sweet spot to image and characterize young extra-solar planetary systems. Finally, we address the science cases of a high-contrast VLTI imaging instrument and focus particularly on exoplanet science (young exoplanets, planet formation, and exozodiacal disks), stellar physics (fundamental parameters and multiplicity), and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.

Atmospheric turbulence and precise measurement of the astrometric baseline vector between any two telescopes are two major challenges in implementing phase referenced interferometric astrometry and imaging. They limit the performance of a fibre-fed interferometer by degrading the instrument sensitivity and astrometric measurements precision and by introducing image reconstruction errors due to inaccurate phases. A multiple beam acquisition and guiding camera was built to meet these challenges for a recently commissioned four beam combiner instrument, GRAVITY, at the ESO Very Large Telescope Interferometer. For each telescope beam it measures: a) field tip-tilts by imaging stars in the sky; b) telescope pupil shifts by imaging pupil reference laser beacons installed on each telescope using a $2 \times 2$ lenslet; c) higher order aberrations using a $9 \times 9$ Shack-Hartmann. The telescope pupils are imaged for a visual monitoring while observing. These measurements enable active field and pupil guiding by actuating a train of tip-tilt mirrors placed in the pupil and field planes, respectively. The Shack-Hartmann measured quasi-static aberrations are used to focus the Auxiliary Telescopes and allow the possibility of correcting the non-common path errors between the Unit Telescopes adaptive optics systems and GRAVITY. The guiding stabilizes light injection into single-mode fibres, increasing sensitivity and reducing the astrometric and image reconstruction errors. The beam guiding enables to achieve astrometric error less than $50\,\mu$as. Here, we report on the data reduction methods and laboratory tests of the multiple beam acquisition and guiding camera and its performance on-sky.

Accurate stellar masses are requested in order to improve our understanding of stellar interiors, but they are still rather rare. Fortunately, the forthcoming Gaia Mission will provide astrometric measurements permitting the derivation of the orbital inclinations of nearby binaries which are also observed as double-lined spectroscopic binaries (SB2s) with ground-based telescopes. A programme of radial velocity (RV) measurements was initiated in 2010 with the Sophie spectrograph of the Haute-Provence observatory in order to derive accurate SB2 orbits for a large set of stars. Therefore, combined SB2+astrometric orbits will be derived thanks to Gaia, and masses with errors around 1~\% are expected for both components. The programme includes 70 SB2s, and the accurate SB2 orbits of 24 of them were already derived. In addition, two complementary programmes devoted to southern stars or to late-type dwarf stars were also initiated with the HERMES and the CARMENES spectrographs, respectively. Interferometric measurements were obained with the VLTI/PIONIER for 7 SB2s, and were taken from other sources for 4 others. Currently, combined "visual binary" (VB) +SB2 solutions were derived for 7 binaries, leading to the masses of the components and to the parallaxes. The parallaxes from the Hipparcos 2 catalogue were corrected for orbital motion and compared to our solution, confirming the high quality of Hipparcos 2.

Context. Convection is a candidate to explain the trigger of red supergiant star mass loss. Owing to the small size of the convective cells on the photosphere, few of their characteristics are known. Aims. Using near infrared interferometry, we intend to resolve the photosphere of red supergiant stars and to bring new constraints on their modeling. Methods. We observed the nearby red supergiant Antares using the 4-telescope instrument VLTI/PIONIER. We collected data on the three available configurations of the 1.8m telescopes in the H band. Results. We obtained unprecedented angular resolution on the disk of a star (6% of the star angular diameter) that limits the mean size of convective cells, brings new constraints on numerical simulations. Using an analytical model with a distribution of bright spots we determine their effect on the visibility signal. Conclusion. We determine that the interferometric signal on Antares is compatible with convective cells with various sizes from 45% to 5% of the angular diameter. We also conclude that convective cells can strongly affect the angular diameter and limb-darkening measurements. In particular, the apparent angular diameter becomes dependent on the sampled position angles.

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

We present 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.

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.

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.

We present near-infrared interferometry of the carbon-rich asymptotic giant branch (AGB) star R Sculptoris. The visibility data indicate a broadly circular resolved stellar disk with a complex substructure. The observed AMBER squared visibility values show drops at the positions of CO and CN bands, indicating that these lines form in extended layers above the photosphere. The AMBER visibility values are best fit by a model without a wind. The PIONIER data are consistent with the same model. We obtain a Rosseland angular diameter of 8.9+-0.3 mas, corresponding to a Rosseland radius of 355+-55 Rsun, an effective temperature of 2640+-80 K, and a luminosity of log L/Lsun=3.74+-0.18. These parameters match evolutionary tracks of initial mass 1.5+-0.5 Msun and current mass 1.3+-0.7 Msun. The reconstructed PIONIER images exhibit a complex structure within the stellar disk including a dominant bright spot located at the western part of the stellar disk. The spot has an H-band peak intensity of 40% to 60% above the average intensity of the limb-darkening-corrected stellar disk. The contrast between the minimum and maximum intensity on the stellar disk is about 1:2.5. Our observations are broadly consistent with predictions by dynamic atmosphere and wind models, although models with wind appear to have a circumstellar envelope that is too extended compared to our observations. The detected complex structure within the stellar disk is most likely caused by giant convection cells, resulting in large-scale shock fronts, and their effects on clumpy molecule and dust formation seen against the photosphere at distances of 2-3 stellar radii.

To understand the chemical composition of planets, it is important to know the chemical composition of the region where they form in protoplanetary disks. Due to its fundamental role in chemical and biological processes, carbon is a key element to trace. We aim to identify the carriers and processes behind the extended NIR flux observed around several Herbig stars. We compare the extended NIR flux from objects in the PIONIER Herbig Ae/Be survey with their flux in the PAH features. HD 100453 is used as a benchmark case to investigate the influence of quantum heated particles, like PAHs or very small carbonaceous grains, in more detail. We use the Monte Carlo radiative transfer code MCMax to do a parameter study of the QHP size and scale- height and examine the influence of quantum heating on the amount of extended flux in the NIR visibilities. There is a correlation between the PAH feature flux of a disk and the amount of its extended NIR flux. We find that very small carbonaceous grains create the observed extended NIR flux around HD 100453 and still lead to a realistic SED. These results can not be achieved without using quantum heating effects, e.g. only with scattered light and grains in thermal equilibrium. It is possible to explain the extended NIR emission around Herbig stars with the presence of carbonaceous, quantum heated particles. Interferometric observations can be used to constrain the spatial distribution and typical size of carbonaceous material in the terrestrial planet forming region.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

We have used (a) HST ACS imaging and STIS spectroscopy, (b) ground-based PIONIER/VLT long-baseline interferometry, and (c) ground-based spectroscopy from different instruments to study the orbit of the extreme multiple system HD 93 129 Aa,Ab, which is composed of (at least) two very massive stars in a long-period orbit with e>0.92 that will pass through periastron in 2017/2018. In several ways, the system is an eta Car precursor. Around the time of periastron passage the two very strong winds will collide and generate an outburst of non-thermal hard X-ray emission without precedent in an O+O binary since astronomers have been able to observe above Earth's atmosphere. A coordinated multiwavelength monitoring in the next two years will enable a breakthrough understanding of the wind interactions in such extreme close encounters. Furthermore, we have found evidence that HD 93 129 Aa may be a binary system itself. In that case, we could witness a three-body interaction that may yield a runaway star or a stellar collision close to or shortly after the periastron passage. Either of those outcomes would be unprecedented, as they are predicted to be low-frequency events in the Milky Way.

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

Context: Extended circumstellar emission has been detected within a few 100 milli-arcsec around > 10% of nearby main sequence stars using near-infrared interferometry. Follow-up observations using other techniques, should they yield similar results or non-detections, can provide strong constraints on the origin of the emission. They can also reveal the variability of the phenomenon. Aims: We aim to demonstrate the persistence of the phenomenon over time scales of a few years and to search for variability of our previously detected excesses. Methods: Using VLTI/PIONIER in H band we have carried out multi-epoch observations of the stars for which a near-infrared excess was previously detected with the same observing technique and instrument. The detection rates and distribution of the excesses from our original survey and the follow-up observations are compared statistically. A search for variability of the excesses in our time series is carried out based on the level of the broadband excesses. Results: In 12 of 16 follow-up observations, an excess is re-detected with a significance of > 2 sigma, and in 7 of 16 follow-up observations significant excess (> 3 sigma) has been re-detected. We statistically demonstrate with very high confidence that the phenomenon persists for the majority of the systems. We also present the first detection of potential variability in two sources. Conclusions: We conclude that the phenomenon responsible for the excesses persists over time scales of a few years for the majority of the systems. However, we also find that variability intrinsic to a target can cause it to have no significant excess at the time of a specific observation.

The nearby red dwarf binary GJ65AB (UB+BL Ceti, M5.5Ve+M6Ve) is a cornerstone system to probe the physics of very low mass stars. We present new interferometric measurements of the angular diameters of the two components of GJ65AB with the PIONIER instrument in the H band: thetaUD (A) = 0.558 +/- 0.008 +/- 0.020 mas and thetaUD (B) = 0.539 +/- 0.009 +/- 0.020 mas. They translate into limb-darkened angular diameters of LD (A) = 0.573 +/- 0.021 mas and LD (B) = 0.554 +/- 0.022 mas. From the known parallax, the radii are R(A) = 0.165 +/- 0.006 Rsun and R(B) = 0.159 +/- 0.006 Rsun (sigma(R)/R = 4%). We also observed GJ65AB with the VLT/NACO adaptive optics and refine the orbital parameters and infrared magnitudes of the system. We derive masses for the two components of m(A) = 0.1225 +/- 0.0043 Msun and m(B) = 0.1195 +/- 0.0043 Msun (sigma(m)/m = 4%). To derive the radial and rotational velocities of the two stars, as well as their relative metallicity with respect to Proxima, we also present new individual UVES high resolution spectra. Placing GJ65 A&B in the mass-radius diagram shows that their radii are larger respectively by 14 +/- 4% and 12 +/- 4% compared to expectations from the models by Baraffe et al. (2015). Following Chabrier et al. (2007), we propose that this discrepancy is caused by the inhibition of convective energy transport by a strong internal magnetic field generated by dynamo effect in these two fast rotating stars. A comparison with the almost identical twin but slowly-rotating Proxima strengthens this hypothesis, as the radius of Proxima does not appear inflated compared to models.

MWC158 is a star with the B[e] phenomenon that shows strong spectrophotometric variability (in lines and in UV and visible continuum) attributed to phases of shell ejection. The evolutionary stage of this star was never clearly determined. Previous interferometric, spectropolarimetric and spectro-interferometric studies suggest a disk morphology for its environment. We investigate the origin of the variability within the inner astronomical unit of the central star using near-infrared interferometric observations with PIONIER at the VLTI over a two-year period. We performed an image reconstruction of the circumstellar environment using the SPARCO method. We discovered that the morphology of the circumstellar environment could vary on timescales of weeks or days. We carried out a parametric fit of the data with a model consisting of a star, a disk and a bright spot that represents a brighter emission in the disk. We detect strong morphological changes in the first astronomical unit around the star, that happen on a timescale of few months. We cannot account for such variability well with a binary model. Our parametric model fits the data well and allows us to extract the location of the asymmetry for different epochs. For the first time, we detect a morphological variability in the environment of MWC158. This variability is reproduced by a model of a disk and a bright spot. The locations of the bright spot suggest that it is located in the disk, but its precise motion is not determined. The origin of the asymmetry in the disk is complex and may be related to asymmetric shell ejections.

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.

We present the first near-IR milli-arcsecond-scale image of a post-AGB binary that is surrounded by hot circumbinary dust. A very rich interferometric data set in six spectral channels was acquired of IRAS08544-4431 with the new RAPID camera on the PIONIER beam combiner at the Very Large Telescope Interferometer (VLTI). A broadband image in the \textit{H} band was reconstructed by combining the data of all spectral channels using the SPARCO method. We spatially separate all the building blocks of the IRAS08544-4431 system in our milliarcsecond-resolution image. Our dissection reveals a dust sublimation front that is strikingly similar to that expected in early-stage protoplanetary disks, as well as an unexpected flux signal of $\sim$4\% from the secondary star. The energy output from this companion indicates the presence of a compact circum-companion accretion disk, which is likely the origin of the fast outflow detected in H$\alpha$. Our image provides the most detailed view into the heart of a dusty circumstellar disk to date. Our results demonstrate that binary evolution processes and circumstellar disk evolution can be studied in detail in space and over time.

The distance to pulsating stars is classically estimated using the parallax-of-pulsation (PoP) method, which combines spectroscopic radial velocity measurements and angular diameter estimates to derive the distance of the star. An important application of this method is the determination of Cepheid distances, in view of the calibration of their distance scale. However, the conversion of radial to pulsational velocities in the PoP method relies on a poorly calibrated parameter, the projection factor (p-factor). We aim to measure empirically the value of the p-factors of a homogeneous sample of nine Galactic Cepheids for which trigonometric parallaxes were measured with the Hubble Space Telescope Fine Guidance Sensor. We use the SPIPS algorithm, a robust implementation of the PoP method that combines photometry, interferometry, and radial velocity measurements in a global modeling of the pulsation. We obtained new interferometric angular diameters using the PIONIER instrument at the Very Large Telescope Interferometer, completed by data from the literature. Using the known distance as an input, we derive the value of the p-factor and study its dependence with the pulsation period. We find the following p-factors: 1.20 $\pm$ 0.12 for RT Aur, 1.48 $\pm$ 0.18 for T Vul, 1.14 $\pm$ 0.10 for FF Aql, 1.31 $\pm$ 0.19 for Y Sgr, 1.39 $\pm$ 0.09 for X Sgr, 1.35 $\pm$ 0.13 for W Sgr, 1.36 $\pm$ 0.08 for $\beta$ Dor, 1.41 $\pm$ 0.10 for $\zeta$ Gem, and 1.23 $\pm$ 0.12 for $\ell$ Car. These values are consistently close to p = 1.324 $\pm$ 0.024. We observe some dispersion around this average value, but the observed distribution is statistically consistent with a constant value of the p-factor as a function of the pulsation period. The error budget of our determination of the p-factor values is presently dominated by the uncertainty on the parallax, a limitation that will soon be waived by Gaia.

Context. Luminous blue variables (LBVs) are rare massive stars with very high luminosity. They are characterized by strong photo-metric and spectroscopic variability related to transient eruptions. The mechanisms at the origin of these eruptions is not well known. In addition, their formation is still problematic and the presence of a companion could help to explain how they form. Aims. This article presents a study of seven LBVs (about 20% of the known Galactic population), some Wolf-Rayet stars, and massive binaries. We probe the environments that surround these massive stars with near-, mid-, and far-infrared images, investigating potential nebula/shells and the companion stars. Methods. To investigate large spatial scales, we used seeing-limited and near diffraction-limited adaptive optics images to obtain a differential diagnostic on the presence of circumstellar matter and to determine their extent. From those images, we also looked for the presence of binary companions on a wide orbit. Once a companion was detected, its gravitational binding to the central star was tested. Tests include the chance projection probability, the proper motion estimates with multi-epoch observations, flux ratio, and star separations. Results. We find that two out of seven of LBVs may have a wide orbit companion. Most of the LBVs display a large circumstellar envelope or several shells. In particular, HD168625, known for its rings, possesses several shells with possibly a large cold shell at the edge of which the rings are formed. For the first time, we have directly imaged the companion of LBV stars.

Spatially resolving the inner dust cavity of the transitional disks is a key to understanding the connection between planetary formation and disk dispersal. The disk around the Herbig star HD 139614 is of particular interest since it presents a pretransitional nature with an au-sized gap, in the dust, that was spatially resolved by mid-IR interferometry. Using new NIR interferometric observations, we aim to characterize the 0.1-10~au region of the HD~139614 disk further and identify viable mechanisms for the inner disk clearing. We report the first multiwavelength radiative transfer modeling of the interferometric data acquired on HD~139614 with PIONIER, AMBER, and MIDI, complemented by Herschel/PACS photometries. We confirm a gap structure in the um-sized dust, extending from about 2.5 au to 6 au, and constrained the properties of the inner dust component: e.g., a radially increasing surface density profile, and a depletion of 10^3 relative to the outer disk. Since self-shadowing and photoevaporation appears unlikely to be responsible for the au-sized gap of HD~139614, we thus tested if dynamical clearing could be a viable mechanism using hydrodynamical simulations to predict the gaseous disk structure. Indeed, a narrow au-sized gap is expected when a single giant planet interacts with the disk. Assuming that small dust grains are well coupled to the gas, we found that a ~ 3~Mjup planet located at 4.5 au from the star could, in less than 1 Myr, reproduce most of the aspects of the dust surface density profile, while no significant depletion in gas occurred in the inner disk, in contrast to the dust. However, the dust-depleted inner disk could be explained by the expected dust filtration by the gap and the efficient dust growth/fragmentation in the inner disk regions. Our results support the hypothesis of a giant planet opening a gap and shaping the inner region of the HD~139614 disk.

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

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

The distances of pulsating stars, in particular Cepheids, are commonly measured using the parallax of pulsation technique. The different versions of this technique combine measurements of the linear diameter variation (from spectroscopy) and the angular diameter variation (from photometry or interferometry) amplitudes, to retrieve the distance in a quasi-geometrical way. However, the linear diameter amplitude is directly proportional to the projection factor (hereafter p-factor), which is used to convert spectroscopic radial velocities (i.e., disk integrated) into pulsating (i.e., photospheric) velocities. The value of the p-factor and its possible dependence on the pulsation period are still widely debated. Our goal is to measure an observational value of the p-factor of the type-II Cepheid kappa Pavonis, whose parallax was measured with an accuracy of 5% using HST/FGS. We used this parallax as a starting point to derive the p-factor of kappa Pav, using the SPIPS technique, which is a robust version of the parallax-of-pulsation method that employs radial velocity, interferometric and photometric data. We applied this technique to a combination of new VLTI/PIONIER optical interferometric angular diameters, new CORALIE and HARPS radial velocities, as well as multi-colour photometry and radial velocities from the literature. We obtain a value of p = 1.26 +/- 0.07 for the p-factor of kappa Pav. This result agrees with several of the recently derived Period-p-factor relationships from the literature, as well as previous observational determinations for Cepheids. Individual estimates of the p-factor are fundamental to calibrating the parallax of pulsation distances of Cepheids. Together with previous observational estimates, the projection factor we obtain points to a weak dependence of the p-factor on period.

Context : The properties of the inner disks of bright Herbig AeBe stars have been studied with near infrared (NIR) interferometry and high resolution spectroscopy. The continuum and a few molecular gas species have been studied close to the central star; however, sensitivity problems limit direct information about the inner disks of the fainter T Tauri stars. Aims : Our aim is to measure some of the properties of the inner regions of disks surrounding southern T Tauri stars. Methods : We performed a survey with the PIONIER recombiner instrument at H-band of 21 T Tauri stars. The baselines used ranged from 11 m to 129 m, corresponding to a maximum resolution of 3mas (0.45 au at 150 pc). Results : Thirteen disks are resolved well and the visibility curves are fully sampled as a function of baseline in the range 45-130 m for these 13 objects. A simple qualitative examination of visibility profiles allows us to identify a rapid drop-off in the visibilities at short baselines in 8 resolved disks. This is indicative of a significant contribution from an extended contribution of light from the disk. We demonstrate that this component is compatible with scattered light, providing strong support to a prediction made by Pinte et al. (2008). The amplitude of the drop-off and the amount of dust thermal emission changes from source to source suggesting that each disk is different. A by-product of the survey is the identification of a new milli-arcsec separation binary: WW Cha. Spectroscopic and interferometric data of AK Sco have also been fitted with a binary and disk model. Conclusions : Visibility data are reproduced well when thermal emission and scattering form dust are fully considered. The inner radii measured are consistent with the expected dust sublimation radii. Modelling of AK Sco suggests a likely coplanarity between the disk and the binary's orbital plane

The main goal of the EXOZODI survey is to detect and characterize circumstellar dust and to propose the first statistical study of exozodiacal disks in the near-infrared using telescopes in both hemispheres. For this purpose, Ertel et al. have conducted in 2012 a survey of nearby main sequence stars with VLTI/PIONIER to search for the presence of circumstellar dust. This survey, carried out during 12 nights, comprises about 100 stars. For each star, we obtained typically three OBs and we searched for circumstellar emission based on the measurement of squared visibilities at short baselines. A drop in the measured visibilities with respect to the expected photospheric visibility indicates the presence of resolved emission around the target star. It is however generally not possible to conclude on the morphology of the detected emission based solely on the squared visibilities. Here, we focus on closure phases to search for faint companions around the whole sample. Indeed, to derive robust statistics on the occurrence rate of bright exozodiacal disks, we need to discriminate between companions and disks. For this reason, the main goal of this paper is to discriminate between circumstellar disks (which show no closure phase provided that they are point-symmetric) and faint companions (point-like sources, creating non-zero closure phases). We also aim to reveal new companions that do not necessarily produce a significant signature in the squared visibilities, as the signature of the companion may show up more prominently in the closure phases. In this process, we reveal four new stellar companions with contrasts ranging from 2% to 95% (i.e., up to equal flux binaries). We also tentatively detect faint companions around one other target that will require follow-up observations to be confirmed or infirmed. We discuss the implications of these discoveries on the results of the exozodi survey.

HD 50138 is a B[e] star surrounded by a large amount of circumstellar gas and dust. Its spectrum shows characteristics which may indicate either a pre- or a post-main-sequence system. Mapping the kinematics of the gas in the inner few au of the system contributes to a better understanding of its physical nature. We present the first high spatial and spectral resolution interferometric observations of the Br-gamma line of HD~50138, obtained with VLTI/AMBER. The line emission originates from a region more compact (up to 3 au) than the continuum-emitting region. Blue- and red-shifted emission originates from the two different hemispheres of an elongated structure perpendicular to the polarization angle. The velocity of the emitting medium decreases radially. An overall offset along the NW direction between the line- and continuum-emitting regions is observed. We compare the data with a geometric model of a thin Keplerian disk and a spherical halo on top of a Gaussian continuum. Most of the data are well reproduced by this model, except for the variability, the global offset and the visibility at the systemic velocity. The evolutionary state of the system is discussed; most diagnostics are ambiguous and may point either to a post-main-sequence or a pre-main-sequence nature.

Context. Detecting and characterizing circumstellar dust is a way to study the architecture and evolution of planetary systems. Cold dust in debris disks only traces the outer regions. Warm and hot exozodiacal dust needs to be studied in order to trace regions close to the habitable zone. Aims. We aim to determine the prevalence and to constrain the properties of hot exozodiacal dust around nearby main-sequence stars. Methods. We search a magnitude limited (H < 5) sample of 92 stars for bright exozodiacal dust using our VLTI visitor instrument PIONIER in the H-band. We derive statistics of the detection rate with respect to parameters such as the stellar spectral type and age or the presence of a debris disk in the outer regions of the systems. We derive more robust statistics by combining our sample with the results from our CHARA/FLUOR survey in the K-band. In addition, our spectrally dispersed data allows us to put constraints on the emission mechanism and the dust properties in the detected systems. Results. We find an over-all detection rate of bright exozodiacal dust in the H-band of 11% (9 out of 85 targets) and three tentative detections. The detection rate decreases from early type to late type stars and increases with the age of the host star. We do not confirm the tentative correlation between the presence of cold and hot dust found in our earlier analysis of the FLUOR sample alone. Our spectrally dispersed data suggest that either the dust is extremely hot or the emission is dominated by the scattered light in most cases. The implications of our results for the target selection of future terrestrial planet finding missions using direct imaging are discussed.

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 mas remain mostly unknown due to intrinsic observational limitations. [...] The Southern MAssive Stars at High angular resolution survey (SMASH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250mas and brightness contrasts of Delta H < 4 and Delta H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8''-radius down to Delta H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1mas to 8'' and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (DEC < 0 deg; H<7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation < 8'' increases to f_m = 0.91 +/- 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30mas while their average number of physically connected companions within 8'' is f_c = 2.2 +/- 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. Additionally, the nine non-thermal (NT) radio emitters observed by SMASH+ are all resolved [...]

The Exozodi survey aims to determine the occurrence rate of bright exozodiacal discs around nearby main sequence stars using infrared interferometry. Although the Exozodi survey targets have been carefully selected to avoid the presence of binary stars, the results of this survey can still be biased by the presence of unidentified stellar companions. Using the PIONIER data set collected within the Exozodi survey, we aim to search for the signature of point-like companions around the Exozodi target stars. We use both the closure phases and squared visibilities collected by PIONIER to search for companions within the ~100 mas interferometric field of view. The presence of a companion is assessed by computing the goodness of fit to the data for a series of binary models with various separations and contrasts. Five stellar companions are resolved for the first time around five A-type stars: HD 4150, HD 16555, HD 29388, HD 202730, and HD 224392 (although the companion to HD 16555 was independently resolved by speckle interferometry while we were carrying out the survey). In the most likely case of main sequence companions, their spectral types range from A5V to K4V. Three of these stars were already suspected to be binaries from Hipparcos astrometric measurements, although no information was available on the companions themselves so far. In addition to debiasing the statistics of the Exozodi survey, these results can also be used to revise the fraction of visual binaries among A-type stars, suggesting that an extra ~13% A-type stars are visual binaries in addition to the ones detected in previous direct imaging surveys. We estimate that about half the population of nearby A-type stars could be resolved as visual binaries using a combination of state-of-the-art interferometry and single-aperture imaging, and we suggest that a significant fraction of these binaries remains undetected to date.

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.

I report the detection of a companion around HD142384 from Optical Long Baseline Interferometry in the near-infrared (flux difference DeltaH~4 and separation ~20mas). HD142384 has been used as PSF reference for published observations of the well-known young star HD142527. This is now strongly discouraged.

Delta Cir is known as an O7.5 III eclipsing and spectroscopic binary with an eccentric orbit. Penny et al. discovered the presence of a third component in the IUE spectra. The eclipsing binary and the third body revolve around a common centre of gravity with a period of 1644 days in an eccentric orbit with the semimajor axis of 10 AU. We demonstrate the presence of the apsidal-line rotation with a period of about 141 years, which is considerably longer than its theoretically predicted value, based on the published radii of the binary components derived from the Hipparchos Hp light curve. However, our new solution of the same light curve resulted in smaller radii and a better agreement between the observed and predicted period of the apsidal- line advance. There are indications that the third body is a binary. The object was resolved by VLTI with the PIONIER combiner; in June 2012 the separation was 3.78 mas, magnitude difference in the H region 1.75 mag. This result means that (assuming the distance 770 pc) the inclination of the long orbit is 87.7 degrees.

Analyses of multi-epoch, high-resolution (R ~ 50.000) optical spectra of the O-type star HD 152246 (O9 IV according to the most recent classification), complemented by a limited number of earlier published radial velocities, led to the finding that the object is a hierarchical triple system, where a close inner pair (Ba-Bb) with a slightly eccentric orbit (e = 0.11) and a period of 6.0049 days revolves in a 470-day highly eccentric orbit (e = 0.865) with another massive and brighter component A. The mass ratio of the inner system must be low since we were unable to find any traces of the secondary spectrum. The mass ratio A/(Ba+Bb) is 0.89. The outer system has recently been resolved using long-baseline interferometry on three occasions. The interferometry confirms the spectroscopic results and specifies elements of the system. Our orbital solutions, including the combined radial-velocity and interferometric solution indicate an orbital inclination of the outer orbit of 112{\deg} and stellar masses of 20.4 and 22.8 solar masses. We also disentangled the spectra of components A and Ba and compare them to synthetic spectra from two independent programmes, TLUSTY and FASTWIND. In either case, the fit was not satisfactory and we postpone a better determination of the system properties for a future study, after obtaining observations during the periastron passage of the outer orbit (the nearest chance being March 2015). For the moment, we can only conclude that component A is an O9 IV star with v*sin(i) = 210 +\- 10 km/s and effective temperature of 33000 +\- 500 K, while component Ba is an O9 V object with v*sin(i) = 65 +/- 3 km/s and T_eff = 33600 +\- 600 K.

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.

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.

Symbiotic stars display absorption lines of a cool red giant together with emission lines of a nebula ionized by a hotter star, indicative of an active binary star system in which mass transfer is occurring. PIONIER at the VLT has been used to combine the light of four telescopes at a time to study in unprecedented detail how mass is transferred in symbiotic stars. The results of a mini-survey of symbiotic stars with PIONIER are summarised and some tentative general results about the role of Roche lobe overflow are presented.

HD 135344B is an accreting (pre-) transition disk that displays the emission of warm CO extending tens of AU inside its 30 AU dust cavity. We used the dust radiative transfer code MCFOST and the thermochemical code ProDiMo to derive the disk structure from the simultaneous modeling of the spectral energy distribution (SED), VLT/CRIRES CO P(10) 4.75 micron, Herschel/PACS [O I] 63 micron, Spitzer-IRS, and JCMT 12CO J=3-2 spectra, VLTI/PIONIER H-band visibilities, and constraints from (sub-)mm continuum interferometry and near-IR imaging. We found a disk model able to describe the current observations simultaneously. This disk has the following structure. (1) To reproduce the SED, the near-IR interferometry data, and the CO ro-vibrational emission, refractory grains (we suggest carbon) are present inside the silicate sublimation radius (0.08<R<0.2 AU). (2) The dust cavity (R<30 AU) is filled with gas, the surface density of this gas must increase with radius to fit the CO P(10) line profile, a small gap of a few AU in the gas is compatible with current data, and a large gap in the gas is not likely. (4) The gas/dust ratio inside the cavity is > 100 to account for the 870 micron continuum upper limit and the CO P(10) line flux. (5) The gas/dust ratio at 30<R<200 AU is < 10 to simultaneously describe the [O I] 63 micron line flux and the CO P(10) line profile. (6) In the outer disk, most of the mass should be located in the midplane, and a significant fraction of the dust is in large grains. Conclusions: Simultaneous modeling of the gas and dust is required to break the model degeneracies and constrain the disk structure. An increasing gas surface density with radius in the inner dust cavity echoes the effect of a migrating Jovian planet. The low gas mass (a few MJupiter) in the HD 135344B's disk suggests that it is an evolved disk that has already lost a large portion of its mass.

We aim at unveiling the observational imprint of physical mechanisms that govern planetary formation in young, multiple systems. In particular, we investigate the impact of tidal truncation on the inner circumstellar disks. We observed the emblematic system GG Tau at high-angular resolution: a hierarchical quadruple system composed of low-mass T Tauri binary stars surrounded by a well-studied, massive circumbinary disk in Keplerian rotation. We used the near-IR 4-telescope combiner PIONIER on the VLTI and sparse-aperture-masking techniques on VLT/NaCo to probe this proto-planetary system at sub-AU scales. We report the discovery of a significant closure-phase signal in H and Ks bands that can be reproduced with an additional low-mass companion orbiting GG Tau Ab, at a (projected) separation rho = 31.7 +/- 0.2mas (4.4 au) and PA = 219.6 +/- 0.3deg. This finding offers a simple explanation for several key questions in this system, including the missing-stellar-mass problem and the asymmetry of continuum emission from the inner dust disks observed at millimeter wavelengths. Composed of now five co-eval stars with 0.02 <= Mstar <= 0.7 Msun, the quintuple system GG Tau has become an ideal test case to constrain stellar evolution models at young ages (few 10^6yr).

The emergence of optical interferometers with three and more telescopes allows image reconstruction of astronomical objects at the milliarcsecond scale. However, some objects contain components with very different spectral energy distributions (SED; i.e. different temperatures), which produces strong chromatic effects on the interferograms that have to be managed with care by image reconstruction algorithms. For example, the gray approximation for the image reconstruction process results in a degraded image if the total (u, v)-coverage given by the spectral supersynthesis is used. The relative flux contribution of the central object and an extended structure changes with wavelength for different temperatures. For young stellar objects, the known characteristics of the central object (i.e., stellar SED), or even the fit of the spectral index and the relative flux ratio, can be used to model the central star while reconstructing the image of the extended structure separately. Methods. We present a new method, called SPARCO (semi-parametric algorithm for the image reconstruction of chromatic objects), which describes the spectral characteristics of both the central object and the extended structure to consider them properly when reconstructing the image of the surrounding environment. We adapted two image-reconstruction codes (Macim, Squeeze, and MiRA) to implement this new prescription. SPARCO is applied using Macim, Squeeze and MiRA on a young stellar object model and also on literature data on HR5999 in the near-infrared with the VLTI. This method paves the way to improved aperture-synthesis imaging of several young stellar objects with existing datasets. More generally, the approach can be used on astrophysical sources with similar features such as active galactic nuclei, planetary nebulae, and asymptotic giant branch stars.

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

Using the PIONIER visitor instrument that combines the light of the four Auxiliary Telescopes of ESO's Very Large Telescope Interferometer, we measure precisely the diameters of several symbiotic and related stars: HD 352, HD 190658, V1261 Ori, ER Del, FG Ser, and AG Peg. These diameters - in the range of 0.6 to 2.3 milli-arcseconds - are used to assess the filling factor of the Roche lobe of the mass-losing giants and provide indications on the nature of the ongoing mass transfer. We also provide the first spectroscopic orbit of ER Del, based on CORAVEL and HERMES/Mercator observations. The system is found to have an eccentric orbit with a period of 5.7 years. In the case of the symbiotic star FG Ser, we find that the diameter is changing by 13% over the course of 41 days, while the observations of HD 352 are indicative of an elongation. Both these stars are found to have a Roche filling factor close to 1, as is most likely the case for HD 190658 as well, while the three other stars have factors below 0.5-0.6. Our observations reveal the power of interferometry for the study of interacting binary stars - the main limitation in our conclusions being the poorly known distances of the objects.

We aim to refine the dynamical masses of the individual component of the low-mass pre-main sequence binary Haro 1-14 C. We combine the data of the preliminary orbit presented previously with new interferometric observations obtained with the four 8m telescopes of the Very Large Telescope Interferometer. The derived masses are $M_a=0.905\pm0.043\,\Msun$ and $M_b=0.308\pm0.011\,\Msun$ for the primary and secondary components, respectively. This is about five times better than the uncertainties of the preliminary orbit. Moreover, the possibility of larger masses is now securely discarded. The new dynamical distance, $d=96\pm\,9\,$pc, is smaller than the distance to the Ophiuchus core with a significance of $2.6\,\sigma$. Fitting the spectral energy distribution yields apparent diameters of $\phi_a=0.13\pm0.01\mas$ and $\phi_b=0.10\pm0.01\mas$ (corresponding to $\Ra=1.50\,\Rsun$ and $\Rb=1.13\,\Rsun$) and a visual extinction of $A_v\approx1.75$. Although the revised orbit has a nearly edge-on geometry, the system is unlikely to be a long-period eclipsing binary. The secondary in Haro~1-14C is one of the few low-mass, pre-main sequence stars with an accurately determined dynamical mass and distance.

Aims: We aim at detecting and characterizing the main-sequence companion of the Cepheid AX Cir ($P_\mathrm{orb} \sim $ 18 yrs). The long-term objective is to estimate the mass of both components and the distance to the system. Methods: We used the PIONIER combiner at the VLT Interferometer to obtain the first interferometric measurements of the short-period Cepheid AX Cir and its orbiting component. Results: The companion is resolved by PIONIER at a projected separation $\rho = 29.2 \pm 0.2$ mas and projection angle $PA = 167.6 \pm 0.3^{\circ}$. We measured $H$-band flux ratios between the companion and the Cepheid of $0.90 \pm 0.10$ % and $0.75 \pm 0.17$ %, respectively at a pulsation phase for the Cepheid $\phi = 0.24$ and 0.48. The lower contrast at $\phi = 0.48$ is due to increased brightness of the Cepheid compared to the $\phi = 0.24$. This gives an average apparent magnitude $m\mathrm{_H (comp)} = 9.06 \pm 0.24$ mag. The limb-darkened angular diameter of the Cepheid at the two pulsation phases was measured to be $\theta_\mathrm{LD} = 0.839 \pm 0.023$ mas and $\theta_\mathrm{LD} = 0.742 \pm 0.020$ mas, respectively at $\phi = 0.24$ and 0.48. A lower limit on the total mass of the system was also derived based on our measured separation, we found $M_\mathrm{T} \geq 9.7 \pm 0.6 M_\odot$.

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.

Young close binaries clear central cavities in their surrounding circumbinary disk from which the stars can still accrete material. This process takes place within the very first astronomical units, and is still not well constrained as the observational evidence has been gathered, until now, only by means of spectroscopy. The young object HD200775 (MWC361) is a massive spectroscopic binary (separation of ~15.9mas, ~5.0~AU), with uncertain classification (early/late Be), that shows a strong and variable Halpha emission. We aim to study the mechanisms that produce the Halpha line at the AU-scale. Combining the radial velocity measurements and astrometric data available in the literature, we determined new orbital parameters. With the VEGA instrument on the CHARA array, we spatially and spectrally resolved the Halpha emission of HD200775, at low and medium spectral resolutions (R~1600 and 5000) over a full orbital period (~3.6 years). We observe that the Halpha equivalent width varies with the orbital phase, and increases close to periastron, as expected from theoretical models that predict an increase of the mass transfer from the circumbinary disk to the primary disk. In addition, using spectral visibilities and differential phases, we find marginal variations of the typical extent of the Halpha emission (at 1 to 2-sigma level) and location (at 1 to 5-sigma level). The spatial extent of the Halpha emission, as probed by a Gaussian FWHM, is minimum at the ascending node (0.67+/-0.20 mas, i.e., 0.22+/-0.06 AU), and more than doubles at periastron. In addition, the Gaussian photocenter is slightly displaced in the direction opposite to the secondary, ruling out the scenario in which all or most of the Halpha emission is due to accretion onto the secondary. These findings, together with the wide Halpha line profile, may be due to a non-spherical wind enhanced at periastron.

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

We aim at measuring the stellar parameters of the two Chemically Peculiar components of the B9.5Vp HgMn + A2 Vm double-lined spectroscopic binary HD141556, whose period is 15.25 days. We combined historical radial velocity measurements with new spatially resolved astrometric observations from PIONIER/VLTI to reconstruct the three-dimensional orbit of the binary, and thus obtained the individual masses. We fit the available photometric points together with the flux ratios provided by interferometry to constrain the individual sizes, which we compared to predictions from evolutionary models.The individual masses of the components are $\Ma = 2.84 \pm 0.12\ \Msun$ and $\Mb = 1.94 \pm 0.09\ \Msun$. The dynamical distance is compatible with the Hipparcos parallax. We find linear stellar radii of $\Ra=2.85 \pm 0.15\ \Rsun$ and $\Rb=1.75 \pm 0.18\ \Rsun$. This result validates a posteriori the flux ratio used in previous detailed abundance studies. We determine a sub-solar initial metallicity $Z=0.012\pm0.003$ and an age of $(2.8\pm0.3)\times10^8\ $years. Our results imply that the primary rotates more slowly than its synchronous velocity, while the secondary is probably synchronous. We show that strong tidal coupling during the pre-main sequence evolution followed by a full decoupling at zero-age main sequence provides a plausible explanation for these very low rotation rates.

(Abridged) Circumstellar disks are believed to be the birthplace of planets and are expected to dissipate on a timescale of a few Myr. The processes responsible for the removal of the dust and gas will strongly modify the radial distribution of the dust and consequently the SED. In particular, a young planet will open a gap, resulting in an inner disk dominating the near-IR emission and an outer disk emitting mostly in the far-IR. We analyze a full set of data (including VLTI/Pionier, VLTI/Midi, and VLT/NaCo/Sam) to constrain the structure of the transition disk around TCha. We used the Mcfost radiative transfer code to simultaneously model the SED and the interferometric observations. We find that the dust responsible for the emission in excess in the near-IR must have a narrow temperature distribution with a maximum close to the silicate sublimation temperature. This translates into a narrow inner dusty disk (0.07-0.11 AU). We find that the outer disk starts at about 12 AU and is partially resolved by the Pionier, Sam, and Midi instruments. We show that the Sam closure phases, interpreted as the signature of a candidate companion, may actually trace the asymmetry generated by forward scattering by dust grains in the upper layers of the outer disk. These observations help constrain the inclination and position angle of the outer disk. The presence of matter inside the gap is difficult to assess with present-day observations. Our model suggests the outer disk contaminates the interferometric signature of any potential companion that could be responsible for the gap opening, and such a companion still has to be unambiguously detected. We stress the difficulty to observe point sources in bright massive disks, and the consequent need to account for disk asymmetries (e.g. anisotropic scattering) in model-dependent search for companions.

Nearby companions alter the evolution of massive stars in binary systems. Using a sample of Galactic massive stars in nearby young clusters, we simultaneously measure all intrinsic binary characteristics relevant to quantify the frequency and nature of binary interactions. We find a large intrinsic binary fraction, a strong preference for short orbital periods and a flat distribution for the mass-ratios. Our results do not support the presence of a significant peak of equal-mass `twin' binaries. As a result of the measured distributions, we find that over seventy per cent of all massive stars exchange mass with a companion. Such a rate greatly exceeds previous estimates and implies that the majority of massive stars have their evolution strongly affected by interaction with a nearby companion.

We aim at resolving the circumstellar environment around beta Pic in the near-infrared in order to study the inner planetary system (< 200 mas, i.e., ~4 AU). Precise interferometric fringe visibility measurements were obtained over seven spectral channels dispersed across the H band with the four-telescope VLTI/PIONIER interferometer. Thorough analysis of interferometric data was performed to measure the stellar angular diameter and to search for circumstellar material. We detected near-infrared circumstellar emission around beta Pic that accounts for 1.37% +/- 0.16% of the near-infrared stellar flux and that is located within the field-of-view of PIONIER (i.e., ~200 mas in radius). The flux ratio between this excess and the photosphere emission is shown to be stable over a period of 1 year and to vary only weakly across the H band, suggesting that the source is either very hot (> 1500 K) or dominated by the scattering of the stellar flux. In addition, we derived the limb-darkened angular diameter of beta Pic with an unprecedented accuracy (theta_LD= 0.736 +/- 0.019 mas). The presence of a small H-band excess originating in the vicinity of beta Pic is revealed for the first time thanks to the high-precision visibilities enabled by VLTI/PIONIER. This excess emission is likely due to the scattering of stellar light by circumstellar dust and/or the thermal emission from a yet unknown population of hot dust, although hot gas emitting in the continuum cannot be firmly excluded.

In the course of our VLTI young stellar object PIONIER imaging program, we have identified a strong visibility chromatic dependency that appeared in certain sources. This effect, rising value of visibilities with decreasing wavelengths over one base, is also present in previous published and archival AMBER data. For Herbig AeBe stars, the H band is generally located at the transition between the star and the disk predominance in flux for Herbig AeBe stars. We believe that this phenomenon is responsible for the visibility rise effect. We present a method to correct the visibilities from this effect in order to allow "gray" image reconstruction software, like Mira, to be used. In parallel we probe the interest of carrying an image reconstruction in each spectral channel and then combine them to obtain the final broadband one. As an illustration we apply these imaging methods to MWC158, a (possibly Herbig) B[e] star intensively observed with PIONIER. Finally, we compare our result with a parametric model fitted onto the data.

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, X-ray binaries and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. Over seventy per cent of all massive stars will exchange mass with a companion, leading to a binary merger in one third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

The visitor instrument PIONIER provides VLTI with improved imaging capabilities and sensitivity. The instrument started routinely delivering scientific data in November 2010, that is less than 12 months after being approved by the ESO Science and Technical Committee. We recall the challenges that had to be tackled to design, built and commission PIONIER. We summarize the typical performances and some astrophysical results obtained so far. We conclude this paper by summarizing lessons learned.

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

We explore the sensitivity and completeness of long baseline interferometric observations for detecting unknown, faint companions around bright unresolved stars. We derive a linear expression for the closure phase signature of a faint companion in the high contrast regime (<0.1), and provide a quantitative estimation of the detection efficiency for the currently offered four-telescope configurations at the Very Large Telescope Interferometer. The results are compared to the performances provided by linear and Y-shaped interferometric configurations in order to identify the ideal array. We find that all configurations have a similar efficiency in discovering companions wider than 10mas. Assuming a closure phase accuracy of 0.25deg, that is typical of state-of-the-art instruments, we predict a median dynamic range of up to six magnitudes when stacking observations obtained at five different hour angles. Surveying bright stars to search for faint companions can be considered as an ideal filler programme for modern interferometric facilities because that places few constraints on the choice of the interferometric configuration.

Using optical long baseline interferometry, we resolved for the first time the two wide components of HD167971, a candidate hierarchical triple system known to efficiently accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a gravitational link between the O8 supergiant and the close eclipsing O + O binary. The separation varies from 8 to 15 mas over the three-year baseline of our observations, suggesting that the components evolve on a wide and very eccentric orbit (most probably e>0.5). These results provide evidence that the wide orbit revealed by our study is not coplanar with the orbit of the inner eclipsing binary. From our measurements of the near-infrared luminosity ratio, we constrain the spectral classification of the components in the close binary to be O6-O7, and confirm that these stars are likely main-sequence objects. Our results are discussed in the context of the bright non-thermal radio emission already reported for this system, and we provide arguments in favour of a maximum radio emission coincident with periastron passage. HD167971 turns out to be an efficient O-type particle accelerator that constitutes a valuable target for future high angular resolution radio imaging using VLBI facilities.

We observed the microlensing event MOA-2007-BLG-192 at high angular resolution in JHKs with the NACO adaptive optics system on the VLT while the object was still amplified by a factor 1.23 and then at baseline 18 months later. We analyzed and calibrated the NACO photometry in the standard 2MASS system in order to accurately constrain the source and the lens star fluxes. We detect light from the host star of MOA-2007-BLG-192Lb, which significantly reduces the uncertainties in its char- acteristics as compared to earlier analyses. We find that MOA-2007-BLG-192L is most likely a very low mass late type M-dwarf (0.084 [+0.015] [-0.012] M\odot) at a distance of 660 [+100] [-70] pc orbited by a 3.2 [+5.2] [-1.8] M\oplus super-Earth at 0.66 [+0.51] [-0.22] AU. We then discuss the properties of this cold planetary system.

Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, multi-band BVRI polarimetry, broad-band and high-resolution interferometric visibilities and phases. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years. The broad H-band interferometry gives a disk diameter of (1.72+-0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50+-9) deg and a relatively low disk flattening of 1.66+-0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii).

Context. Determining the mass transfer in a close binary system is of prime importance for understanding its evolution. SS Leporis, a symbiotic star showing the Algol paradox and presenting clear evidence of ongoing mass transfer, in which the donor has been thought to fill its Roche lobe, is a target particularly suited to this kind of study. Aims. Since previous spectroscopic and interferometric observations have not been able to fully constrain the system morphology and characteristics, we go one step further to determine its orbital parameters, for which we need new interferometric observations directly probing the inner parts of the system with a much higher number of spatial frequencies. Methods. We use data obtained at eight different epochs with the VLTI instruments AMBER and PIONIER in the H- and K-bands. We performed aperture synthesis imaging to obtain the first model-independent view of this system. We then modelled it as a binary (whose giant is spatially resolved) that is surrounded by a circumbinary disc. Results. Combining these interferometric measurements with previous radial velocities, we fully constrain the orbit of the system. We then determine the mass of each star and significantly revise the mass ratio. The M giant also appears to be almost twice smaller than previously thought. Additionally, the low spectral resolution of the data allows the flux of both stars and of the dusty disc to be determined along the H and K bands, and thereby extracting their temperatures. Conclusions. We find that the M giant actually does not stricto sensus fill its Roche lobe. The mass transfer is more likely to occur through the accretion of an important part of the giant wind. We finally rise the possibility for an enhanced mass loss from the giant, and we show that an accretion disc should have formed around the A star.

Context. A new four-telescope interferometric instrument called PIONIER has recently been installed at VLTI. It provides improved imaging capabilities together with high precision. Aims. We search for low-mass companions around a few bright stars using different strategies, and determine the dynamic range currently reachable with PIONIER. Methods. Our method is based on the closure phase, which is the most robust interferometric quantity when searching for faint companions. We computed the chi^2 goodness of fit for a series of binary star models at different positions and with various flux ratios. The resulting chi^2 cube was used to identify the best-fit binary model and evaluate its significance, or to determine upper limits on the companion flux in case of non detections. Results. No companion is found around Fomalhaut, tau Cet and Regulus. The median upper limits at 3 sigma on the companion flux ratio are respectively of 2.3e-3 (in 4 h), 3.5e-3 (in 3 h) and 5.4e-3 (in 1.5 h) on the search region extending from 5 to 100 mas. Our observations confirm that the previously detected near-infrared excess emissions around Fomalhaut and tau Cet are not related to a low-mass companion, and instead come from an extended source such as an exozodiacal disk. In the case of del Aqr, in 30 min of observation, we obtain the first direct detection of a previously known companion, at an angular distance of about 40 mas and with a flux ratio of 2.05e-2 \pm 0.16e-2. Due to the limited u,v plane coverage, its position can, however, not be unambiguously determined. Conclusions. After only a few months of operation, PIONIER has already achieved one of the best dynamic ranges world-wide for multi-aperture interferometers. A dynamic range up to about 1:500 is demonstrated, but significant improvements are still required to reach the ultimate goal of directly detecting hot giant extrasolar planets.

As the brightest O-type X-ray source in the Carina nebula, HD 93250 (O4 III(fc)) is X-ray overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales of years suggesting either a single star, a binary system seen pole-one or a very long period and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 0.8+/-0.1 and a separation of 1.5+/-0.2 x 10E-03 arcsec. At the distance of Carina, this corresponds to a projected physical distance of 3.5 A.U. While a quantitative investigation would require a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively explain both its X-ray flux and hardness and its non-thermal radio emission in the framework of a colliding wind scenario. We also discuss various observational biases. We show that, due to line-blending of two similar spectral components, HD 93250 could have a period as short as one to several years despite the lack of measurable radial velocity variations.

We present the earliest results of our NACO/VLT sparse aperture masking (SAM) campaign to search for binarity in a sample of 60 O-type stars. We detect Delta Ks < 5 mag companions for 20-25% of our targets with separations in the range 30-100 mas (typically, 40 - 200 A.U.). Most of these companions were unknown, shedding thus new light on the multiplicity properties of massive stars in a separation and brightness regime that has been difficult to explore so far. Adding detections from other techniques (spectroscopy, interferometry, speckle, lucky imaging, AO), the fraction of O stars with at least one companion is 85% (51/60 targets). This is the largest multiplicity fraction ever found.

We report on near-IR interferometric observations of the outburst of the recurrent nova T Pyx. We obtained near-IR observations of T Pyx at dates ranging from t=2.37d to t=48.2d after the outburst, with the CLASSIC recombiner, located at the CHARA array, and with the PIONIER and AMBER recombiners, located at the VLTI array. These data are supplemented with near-IR photometry and spectra obtained at Mount Abu, India. Slow expansion velocities were measured (<300km/s) before t=20d (assuming D=3.5kpc). From t=28d on, the AMBER and PIONIER continuum visibilities (K and H band, respectively) are best simulated with a two component model consisting of an unresolved source plus an extended source whose expansion velocity onto the sky plane is lower than 700km/s. The expansion of the Brgamma line forming region, as inferred at t=28d and t=35d is slightly larger, implying velocities in the range 500-800km/s, still strikingly lower than the velocities of 1300-1600km/s inferred from the Doppler width of the line. Moreover, a remarkable pattern was observed in the Brgamma differential phases. A semi-quantitative model using a bipolar flow with a contrast of 2 between the pole and equator velocities, an inclination of i=15^{\circ} and a position angle P.A.=110^{\circ} provides a good match to the AMBER observables (spectra, differential visibilities and phases). At t=48d, a PIONIER dataset confirms the two component nature of the H band emission, consisting of an unresolved stellar source and an extended region whose appearance is circular and symmetric within error bars.These observations are most simply interpreted within the frame of a bipolar model, oriented nearly face-on. This finding has profound implications for the interpretation of past, current and future observations of the expanding nebula.

PIONIER stands for Precision Integrated-Optics Near-infrared Imaging ExpeRiment. It combines four 1.8m Auxilliary Telescopes or four 8m Unit Telescopes of the Very Large Telescope Interferometer (ESO, Chile) using an integrated optics combiner. The instrument has been integrated at IPAG starting in December 2009 and commissioned at the Paranal Observatory in October 2010. It provides scientific observations since November 2010. In this paper, we detail the instrumental concept, we describe the standard operational modes and the data reduction strategy. We present the typical performance and discuss how to improve them. This paper is based on laboratory data obtained during the integrations at IPAG, as well as on-sky data gathered during the commissioning at VLTI. We illustrate the imaging capability of PIONIER on the binaries deltaSco and HIP11231. PIONIER provides 6 visibilities and 3 independent closure phases in the H band, either in a broadband mode or with a low spectral dispersion (R=40), using natural light (i.e. unpolarized). The limiting magnitude is Hmag=7 in dispersed mode under median atmospheric conditions (seeing<1", tau0>3ms) with the 1.8m Auxiliary Telescopes. We demonstrate a precision of 0.5deg on the closure phases. The precision on the calibrated visibilities ranges from 3 to 15% depending on the atmospheric conditions. PIONIER has been installed and successfully tested as a visitor instrument for the VLTI. It permits high angular resolution imaging studies at an unprecedented level of sensitivity. The successful combination of the four 8m Unit Telescopes in March 2011 demonstrates that VLTI is ready for 4-telescope operation.

Although debris discs have been detected around a significant number of main-sequence stars, only a few of them are known to harbour hot dust in their inner part where terrestrial planets may have formed. Thanks to infrared interferometric observations, it is possible to obtain a direct measurement of these regions, which are of prime importance for preparing future exo-Earth characterisation missions. In this context, we have resolved the exozodiacal dust disc around Vega with the help of infrared stellar interferometry and estimated the integrated H-band flux originating from the first few AUs of the debris disc. Using precise H-band interferometric measurements obtained with the 3-telescope IOTA/IONIC interferometer (Mount Hopkins, Arizona), thorough modelling of both interferometric data (squared visibility and closure phase) and spectral energy distribution was performed to constrain the nature of the near-infrared excess emission. The most straightforward scenario consists in a compact dust disc producing a thermal emission that is largely dominated by small grains located between 0.1 and 0.3 AU from Vega and accounting for 1.23 +/- 0.45% of the near-infrared stellar flux for our best-fit model. This flux ratio is shown to vary slightly with the geometry of the model used to fit our interferometric data (variations within +/-0.19%). Initially revealed by K-band CHARA/FLUOR observations, the presence of hot exozodiacal dust in the vicinity of Vega is confirmed by our H-band IOTA/IONIC measurements at the 3-sigma level. Whereas the origin of the dust is still uncertain, its presence and the possible connection with the outer disc suggest that the Vega system is currently undergoing major dynamical perturbations.

The pre-main sequence (PMS) star ABDorA is the main component of the quadruple system ABDoradus. The precise determination of the mass and photometry of the close companion to ABDorA, ABDorC, has provided an important benchmark for calibration of theoretical evolutionary models of low-mass stars. The limiting factor to the precision of this calibration is the age of the system, as both the mass and luminosity of ABDorA and C are well monitored by other ongoing programs. In this paper we present VLTI/AMBER observations of ABDorA which provide a direct measurement of the size of this star, 0.96+/-0.06 Rsun. The latter estimate, combined with other fundamental parameters also measured for this star, allows a precise test of PMS evolutionary models using both H-R diagrams and mass-radius relationships. We have found that our radius measurement is larger than that predicted by the models, which we interpret as an evidence of the oversizing produced by the strong magnetic activity of ABDorA. Considering, at least partially, this magnetic effect, theoretical isochrones have been used to derive constraints to the age of ABDorA, favouring an age about 40-50 Myr for this system. Older ages are not completely excluded by our data.

We present multi-instruments and multi-wavelengths observations of the famous LBV star Pistol Star. These observations are part of a larger program about early O stars at different metallicities. The Pistol star has been claimed as the most massive star known, with 250 solar masses. We present the preliminary results based on X-Shooter spectra, as well as the observations performed with the VLTI-AMBER and the VLT-NACO adaptive optics. The X-shooter spectrograph allows to obtain simultaneously a spectrum from the UV to the K-band with a resolving power of $\sim$15000. The preliminary results obtained indicate that Pistol Star has similar properties of Eta Car, including shells of matter, but also the binarity. Other objects of the program, here briefly presented, were selected for their particular nature: early O stars with mass discrepancies between stellar evolution models and observations, discrepancies with the wind momentum luminosity relation.

First results of near-IR adaptive optics (AO)-assisted imaging, interferometry, and spectroscopy of this Luminous Blue Variable (LBV) are presented. They suggest that the Pistol Star is at least double. If the association is physical, it would reinforce questions concerning the importance of multiplicity for the formation and evolution of extremely massive stars.

Aims. We search for low-mass companions in the innermost region (<300 mas, i.e., 6 AU) of the beta Pic planetary system. Methods. We obtained interferometric closure phase measurements in the K-band with the VLTI/AMBER instrument used in its medium spectral resolution mode. Fringe stabilization was provided by the FINITO fringe tracker. Results. In a search region of between 2 and 60 mas in radius, our observations exclude at 3 sigma significance the presence of companions with K-band contrasts greater than 5e-3 for 90% of the possible positions in the search zone (i.e., 90% completeness). The median 1-sigma error bar in the contrast of potential companions within our search region is 1.2e-3. The best fit to our data set using a binary model is found for a faint companion located at about 14.4 mas from beta Pic, which has a contrast of 1.8e-3 \pm 1.1e-3 (a result consistent with the absence of companions). For angular separations larger than 60 mas, both time smearing and field-of-view limitations reduce the sensitivity. Conclusions. We can exclude the presence of brown dwarfs with masses higher than 29 Mjup (resp. 47 Mjup) at a 50% (resp. 90%) completeness level within the first few AUs around beta Pic. Interferometric closure phases offer a promising way to directly image low-mass companions in the close environment of nearby young stars.

Z CMa is a young binary system consisting of an Herbig primary and a FU Ori companion. Both components seem to be surrounded by active accretion disks and a jet was associated to the Herbig B0. In Nov. 2008, K. Grankin discovered that Z CMa was exhibiting an outburst with an amplitude larger than any photometric variations recorded in the last 25 years. To study the innermost regions in which the outburst occurs and understand its origin, we have observed both binary components with AMBER/VLTI across the Br{\gamma} emission line in Dec. 2009 in medium and high spectral resolution modes. Our observations show that the Herbig Be, responsible for the increase of luminosity, also produces a strong Br{\gamma} emission, and they allow us to disentangle from various origins by locating the emission at each velocities through the line. Considering a model of a Keplerian disk alone fails at reproducing the asymmetric spectro-astrometric measurements, suggesting a major contribution from an outflow.

We present the third release of the AMBER data reduction software by the JMMC. This software is based on core algorithms optimized after several years of operation. An optional graphic interface in a high level language allows the user to control the process step by step or in a completely automatic manner. Ongoing improvement is the implementation of a robust calibration scheme, making use of the full calibration sets available during the night. The output products are standard OI-FITS files, which can be used directly in high level software like model fitting or image reconstruction tools. The software performances are illustrated on a full data set of calibrators observed with AMBER during 5 years taken in various instrumental setup.

Accretion is a fundamental process in star formation. Although the time evolution of accretion remains a matter of debate, observations and modelling studies suggest that episodic outbursts of strong accretion may dominate the formation of the protostar. Observing young stellar objects during these elevated accretion states is crucial to understanding the origin of unsteady accretion. ZCMa is a pre-main-sequence binary system composed of an embedded Herbig Be star, undergoing photometric outbursts, and a FU Orionis star. The Herbig Be component recently underwent its largest optical photometric outburst detected so far. We aim to constrain the origin of this outburst by studying the emission region of the HI Brackett gamma line, a powerful tracer of accretion/ejection processes on the AU-scale in young stars. Using the AMBER/VLTI instrument at spectral resolutions of 1500 and 12 000, we performed spatially and spectrally resolved interferometric observations of the hot gas emitting across the Brackett gamma emission line, during and after the outburst. From the visibilities and differential phases, we derive characteristic sizes for the Brackett gamma emission and spectro-astrometric measurements across the line, with respect to the continuum. We find that the line profile, the astrometric signal, and the visibilities are inconsistent with the signature of either a Keplerian disk or infall of matter. They are, instead, evidence of a bipolar wind, maybe partly seen through a disk hole inside the dust sublimation radius. The disappearance of the Brackett gamma emission line after the outburst suggests that the outburst is related to a period of strong mass loss rather than a change of the extinction along the line of sight. Based on these conclusions, we speculate that the origin of the outburst is an event of enhanced mass accretion, similar to those occuring in EX Ors and FU Ors.

We obtained near-infrared long-baseline interferometry of IRC+10420 with the AMBER instrument of ESO's Very Large Telescope Interferometer (VLTI) in low and high spectral resolution (HR) mode to probe the photosphere and the innermost circumstellar environment of this rapidly evolving yellow hypergiant. In the HR observations, the visibilities show a noticeable drop across the Brackett gamma (BrG) line on all three baselines, and we found differential phases up to -25 degrees in the redshifted part of the BrG line and a non-zero closure phase close to the line center. The calibrated visibilities were corrected for AMBER's limited field-of-view to appropriately account for the flux contribution of IRC+10420's extended dust shell. We derived FWHM Gaussian sizes of 1.05 +/- 0.07 and 0.98 +/- 0.10 mas for IRC+10420's continuum-emitting region in the H and K bands, respectively, and the BrG-emitting region can be fitted with a geometric ring model with a diameter of 4.18 +0.19/-0.09 mas, which is approximately 4 times the stellar size. The geometric model also provides some evidence that the BrG line-emitting region is elongated towards a position angle of 36 degrees, well aligned with the symmetry axis of the outer reflection nebula. The HR observations were further analyzed by means of radiative transfer modeling using CMFGEN and the 2-D Busche & Hillier codes. Our spherical CMFGEN model poorly reproduces the observed line shape, blueshift, and extension, definitively showing that the IRC+10420 outflow is asymmetric. Our 2-D radiative transfer modeling shows that the blueshifted BrG emission and the shape of the visibility across the emission line can be explained with an asymmetric bipolar outflow with a high density contrast from pole to equator (8-16), where the redshifted light is substantially diminished.

The innermost parts of dusty debris disks around main sequence stars are currently poorly known due to the high contrast and small angular separation with their parent stars. Using near-infrared interferometry, we aim to detect the signature of hot dust around the nearby A4 V star Fomalhaut, which has already been suggested to harbor a warm dust population in addition to a cold dust ring located at about 140 AU. Archival data obtained with the VINCI instrument at the VLTI are used to study the fringe visibility of the Fomalhaut system at projected baseline lengths ranging from 4 m to 140 m in the K band. A significant visibility deficit is observed at short baselines with respect to the expected visibility of the sole stellar photosphere. This is interpreted as the signature of resolved circumstellar emission, producing a relative flux of 0.88% +/- 0.12% with respect to the stellar photosphere. While our interferometric data cannot directly constrain the morphology of the excess emission source, complementary data from the literature allow us to discard an off-axis point-like object as the source of circumstellar emission. We argue that the thermal emission from hot dusty grains located within 6 AU from Fomalhaut is the most plausible explanation for the detected excess. Our study also provides a revised limb-darkened diameter for Fomalhaut (2.223 +/- 0.022 mas), taking into account the effect of the resolved circumstellar emission.

While it is well known that most massive stars are found to be part of binary or multiple systems, an accurate characterization of the statistical properties of these multiple objects is still lacking. In the present talk, we will review the current status of the field, emphasizing the need of using complementarity techniques to cover the large parameter space. We will also describe what we think is the place of interferometry in this context.

We measured the radii of 7 low and very low-mass stars using long baseline interferometry with the VLTI interferometer and its VINCI and AMBER near-infrared recombiners. We use these new data, together with literature measurements, to examine the luminosity- radius and mass-radius relations for K and M dwarfs. The precision of the new interferometric radii now competes with what can be obtained for double-lined eclipsing binaries. Interferometry provides access to much less active stars, as well as to stars with much better measured distances and luminosities, and therefore complements the information obtained from eclipsing systems. The radii of magnetically quiet late-K to M dwarfs match the predictions of stellar evolution models very well, providing direct confirmation that magnetic activity explains the discrepancy that was recently found for magnetically active eclipsing systems. The radii of the early K dwarfs are well reproduced for a mixing length parameter that approaches the solar value, as qualitatively expected.

We discuss the spin-orbit orientation of the Fomalhaut planetary system composed of a central A4V star, a debris disk, and a recently discovered planetary companion. We use spectrally resolved, near-IR long baseline interferometry to obtain precise spectro-astrometric measurements across the brG absorption line. The achieved astrometric accuracy of 3 nu-as and the spectral resolution R=1500 from the AMBER/VLTI instrument allow us to spatially and spectrally resolve the rotating photosphere. We find a position angle PAstar=65deg pm 3deg for the stellar rotation axis, perfectly perpendicular with the literature measurement for the disk position angle (PAdisk=156deg pm 0.3deg). This is the first time such test can be performed for a debris disk, and in a non-eclipsing system. Additionally, our measurements suggest unexpected backward-scattering properties for the circumstellar dust grains. Our observations validate the standard scenario for star and planet formation, in which the angular momentum of the planetary systems are expected to be collinear with the stellar spins.

Diffuse envelopes around Mira variables are among the most important sources influencing the chemical evolution of galaxies. However they represent an observational challenge because of their complex spectral features and their rapid temporal variability. We constrained the exact brightness distribution of the Mira star TLep with a model-independent analysis. We obtained single-epoch interferometric observations with a dataset continuous in the spectral domain (1.5-2.4mum) and in the spatial domain (baselines ranging from 11 to 96m). We performed a model independent image reconstruction for each spectral bin using the MIRA software. We completed the analysis by modeling the data with a simple star+layer model inspired from the images. Reconstructed images confirm the general picture of a central star partially obscured by the surrounding molecular shell of changing opacity. At 1.7mum, the shell becomes optically thin, with corresponding emission appearing as a ring circling the star. This is the first direct evidence of the spherical morphology of the molecular shell. Model fitting confirmed a spherical layer of constant size and changing opacity over the wavelengths. Rough modeling points to a continuum opacity within the shell, in addition to the CO and H2O features. Accordingly, it appeared impossible to model the data by a photosphere alone in any of the spectral bins.

Aims. In this paper we model, in a self-consistent way, polarimetric, photometric, spectrophotometric and interferometric observations of the classical Be star $\zeta$ Tauri. Our primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. We have carried out detailed three-dimensional, NLTE radiative transfer calculations using the radiative transfer code HDUST. For the input for the code we have used the most up-to-date research on Be stars to include a physically realistic description for the central star and the circumstellar disc. We adopt a rotationally deformed, gravity darkened central star, surrounded by a disc whose unperturbed state is given by a steady-state viscous decretion disc model. We further assume that disc is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disc model for $\zeta$ Tauri and a rigorous solution of the radiative transfer, we have obtained a very good fit of the time-average properties of the disc. This provides strong theoretical evidence that the viscous decretion disc model is the mechanism responsible for disc formation. With the global oscillation model we have successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations of the H$\alpha$ and Br$\gamma$ lines. This result convincingly demonstrates that the oscillation pattern in the disc is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

At the Very Large Telescope Interferometer, the purpose of the fringe-tracker FINITO is to stabilize the optical path differences between the beams, allowing longer integration times on the scientific instruments AMBER and MIDI. Our goal is to demonstrate the potential of FINITO for providing H-band interferometric visibilities, simultaneously and in addition to its normal fringe-tracking role. We use data obtained during the commissioning of the Reflective Memory Network Recorder at the Paranal observatory. This device has permitted the first recording of all relevant real-time data needed for a proper data-reduction. We show that post-processing the FINITO data allows valuable scientific visibilities to be measured. Over the several hours of our engineering experiment, the intrinsic transfer function is stable at the level of 2%. Such stability would lead to robust measurements of science stars even without the observation of a calibration star within a short period of time. We briefly discuss the current limitations and the potential improvements.

The VLTI Spectro Imager project aims to perform imaging with a temporal resolution of 1 night and with a maximum angular resolution of 1 milliarcsecond, making best use of the Very Large Telescope Interferometer capabilities. To fulfill the scientific goals (see Garcia et. al.), the system requirements are: a) combining 4 to 6 beams; b) working in spectral bands J, H and K; c) spectral resolution from R= 100 to 12000; and d) internal fringe tracking on-axis, or off-axis when associated to the PRIMA dual-beam facility. The concept of VSI consists on 6 sub-systems: a common path distributing the light between the fringe tracker and the scientific instrument, the fringe tracker ensuring the co-phasing of the array, the scientific instrument delivering the interferometric observables and a calibration tool providing sources for internal alignment and interferometric calibrations. The two remaining sub-systems are the control system and the observation support software dedicated to the reduction of the interferometric data. This paper presents the global concept of VSI science path including the common path, the scientific instrument and the calibration tool. The scientific combination using a set of integrated optics multi-way beam combiners to provide high-stability visibility and closure phase measurements are also described. Finally we will address the performance budget of the global VSI instrument. The fringe tracker and scientific spectrograph will be shortly described.

The VLTI Spectro Imager (VSI) was proposed as a second-generation instrument of the Very Large Telescope Interferometer providing the ESO community with spectrally-resolved, near-infrared images at angular resolutions down to 1.1 milliarcsecond and spectral resolutions up to R=12000. Targets as faint as K=13 will be imaged without requiring a brighter nearby reference object. The unique combination of high-dynamic-range imaging at high angular resolution and high spectral resolution enables a scientific program which serves a broad user community and at the same time provides the opportunity for breakthroughs in many areas of astrophysic including: probing the initial conditions for planet formation in the AU-scale environments of young stars; imaging convective cells and other phenomena on the surfaces of stars; mapping the chemical and physical environments of evolved stars, stellar remnants, and stellar winds; and disentangling the central regions of active galactic nuclei and supermassive black holes. VSI will provide these new capabilities using technologies which have been extensively tested in the past and VSI requires little in terms of new infrastructure on the VLTI. At the same time, VSI will be able to make maximum use of new infrastructure as it becomes available; for example, by combining 4, 6 and eventually 8 telescopes, enabling rapid imaging through the measurement of up to 28 visibilities in every wavelength channel within a few minutes. The current studies are focused on a 4-telescope version with an upgrade to a 6-telescope one. The instrument contains its own fringe tracker and tip-tilt control in order to reduce the constraints on the VLTI infrastructure and maximize the scientific return.

Nowadays, compact sources like surfaces of nearby stars, circumstellar environments of stars from early stages to the most evolved ones and surroundings of active galactic nuclei can be investigated at milli-arcsecond scales only with the VLT in its interferometric mode. We propose a spectro-imager, named VSI (VLTI spectro-imager), which is capable to probe these sources both over spatial and spectral scales in the near-infrared domain. This instrument will provide information complementary to what is obtained at the same time with ALMA at different wavelengths and the extreme large telescopes.

Aims: We present J, H, K interferometry with a spectral resolution of 35 for the Mira variable S Orionis. We aim at measuring the diameter variation as a function of wavelength that is expected due to molecular layers lying above the continuum-forming photosphere. Methods: Visibility data of S Ori were obtained at phase 0.78 with the VLTI/AMBER instrument using the fringe tracker FINITO at 29 spectral channels between 1.29 and 2.32 mu. Apparent uniform disk (UD) diameters were computed for each spectral channel. In addition, the visibility data were directly compared to predictions by recent self-excited dynamic model atmospheres. Results: S Ori shows significant variations in the visibility values as a function of spectral channel that can only be described by a clear variation in the apparent angular size with wavelength. The closure phase values are close to zero at all spectral channels, indicating the absence of asymmetric intensity features. The apparent UD angular diameter is smallest at about 1.3 and 1.7 mu and increases by a factor of ~1.4 around 2.0 mu. The minimum UD angular diameter is 8.1 pm 0.5 mas, corresponding to ~420 R_sun. The S Ori visibility data and the apparent UD variations can be explained reasonably well by a dynamic atmosphere model that includes molecular layers. Conclusions: The measured visibility and UD diameter variations with wavelength resemble and generally confirm the predictions by recent dynamic model atmospheres. [abridged]

Our goal is to demonstrate the potential of the interferometric AMBER instrument linked with the Very Large Telescope Interferometer (VLTI) fringe-tracking facility FINITO to derive high-precision stellar diameters. We use commissioning data obtained on the bright single star V3879 Sgr. Locking the interferometric fringes with FINITO allows us to record very low contrast fringes on the AMBER camera. By fitting the amplitude of these fringes, we measure the diameter of the target in three directions simultaneously with an accuracy of 25 micro-arcseconds. We showed that V3879 Sgr has a round photosphere down to a sub-percent level. We quickly reached this level of accuracy because the technique used is independent from absolute calibration (at least for baselines that fully span the visibility null). We briefly discuss the potential biases found at this level of precision. The proposed AMBER+FINITO instrumental setup opens several perspectives for the VLTI in the field of stellar astrophysics, like measuring with high accuracy the oblateness of fast rotating stars or detecting atmospheric starspots.

Incoming optical interferometers will allow spectro-imaging at high angular resolution. Non-homothetic Fizeau concept combines good sensitivity and high spectral resolution capabilities. However, one critical issue is the design of the beam recombination scheme, at the heart of the instrument. We tackle the possibility of reducing the number of pixels that are coding the fringes by compressing the pupil plane. Shrinking the number of pixels -- which drastically increases with the number of recombined telescopes -- is indeed a key issue that enables to reach higher limiting magnitude, but also allows to lower the required spectral resolution and fasten the fringes reading process. By means of numerical simulations, we study the performances of existing estimators of the visibility with respect to the compression process. We show that, not only the model based estimator lead to better signal to noise ratio (SNR) performances than the Fourier ones, but above all it is the only one which prevent from introducing baseline mixing biases in the visibilities as the pupil plane compression rate increases. Furthermore, we show that moderate compression allows to keep the visibilities SNR unaffected. In the light of these conclusions, we propose an optimized pupil arrangements for 6 and 8 beam recombiners.

VLTi Spectro-Imager (VSI) is a proposition for a second generation VLTI instrument which is aimed at providing the ESO community with the capability of performing image synthesis at milli-arcsecond angular resolution. VSI provides the VLTI with an instrument able to combine 4 telescopes in a baseline version and optionally up to 6 telescopes in the near-infrared spectral domain with moderate to high spectral resolution. The instrument contains its own fringe tracker in order to relax the constraints onto the VLTI infrastructure. VSI will do imaging at the milli-arcsecond scale with spectral resolution of: a) the close environments of young stars probing the initial conditions for planet formation; b) the surfaces of stars; c) the environment of evolved stars, stellar remnants and stellar winds, and d) the central region of active galactic nuclei and supermassive black holes. The science cases allowed us to specify the astrophysical requirements of the instrument and to define the necessary studies of the science group for phase A.

Our objective is to prove that integrated optics (IO) is not only a good concept for astronomical interferometry but also a working technique with high performance. We used the commissioning data obtained with the dedicated K-band integrated optics two-telescope beam combiner which now replaces the fiber coupler MONA in the VLTI/VINCI instrument. We characterize the behaviour of this IO device and compare its properties to other single mode beam combiner like the previously used MONA fiber coupler. The IO combiner provides a high optical throughput, a contrast of 89% with a night-to-night stability of a few percent. Even if a dispersive phase is present, we show that it does not bias the measured Fourier visibility estimate. An upper limit of 0.005 for the cross-talk between linear polarization states has been measured. We take advantage of the intrinsic contrast stability to test a new astronomical prodecure for calibrating diameters of simple stars by simultaneously fitting the instrumental contrast and the apparent stellar diameters. This method reaches an accuracy with diameter errors of the order of previous ones but without the need of an already known calibrator. These results are an important step of integrated optics and paves the road to incoming imaging interferometer projects.