Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

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

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

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.

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

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.

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.

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

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.

The estimation of atmospheric turbulence parameters is of relevance for: a) site evaluation & characterisation; b) prediction of the point spread function; c) live assessment of error budgets and optimisation of adaptive optics performance; d) optimisation of fringe trackers for long baseline optical interferometry. The ubiquitous deployment of Shack-Hartmann wavefront sensors in large telescopes makes them central for atmospheric turbulence parameter estimation via adaptive optics telemetry. Several methods for the estimation of the Fried parameter and outer scale have been developed, most of which are based on the fitting of Zernike polynomial coefficients variances reconstructed from the telemetry. The non-orthogonality of Zernike polynomial derivatives introduces modal cross coupling, which affects the variances. Furthermore, the finite resolution of the sensor introduces aliasing. In this article the impact of these effects on atmospheric turbulence parameter estimation is addressed with simulations. It is found that cross coupling is the dominant bias. An iterative algorithm to overcome it is presented. Simulations are conducted for typical ranges of the outer scale (4 to 32m), Fried parameter (10 cm) and noise in the variances (signal-to-noise ratio of 10 and above). It is found that, using the algorithm, both parameters are recovered with sub-percent accuracy.

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.

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.

V923 Sco is a bright ($V$ = 5.91), nearby ($\pi$ = 15.46$\pm$0.40 mas) southern eclipsing binary. Because both components are slow rotators, the minimum masses of the components are known with 0.2% precision from spectroscopy. The system seems ideal for very precise mass, radius, and luminosity determinations and, owing to its proximity and long orbital period ($\sim$ 34.8 days), promises to be resolved with long-baseline interferometry. The principal aim is very accurate determinations of absolute stellar parameters for both components of the eclipsing binary and a model-independent determination of the distance.} New high-precision photometry of both eclipses of V923 Sco with the MOST satellite was obtained. The system was spatially resolved with the VLTI AMBER, PIONIER, and GRAVITY instruments at nine epochs. Combining the projected size of the spectroscopic orbit (in km) and visual orbit (in mas) the distance to the system is derived. Simultaneous analysis of photometric, spectroscopic, and interferometric data was performed to obtain a robust determination of the absolute parameters. Very precise absolute parameters of the components were derived in spite of the parameter correlations. The primary component is found to be overluminous for its mass. Combining spectroscopic and interferometric observations enabled us to determine the distance to V923 Sco with better than 0.2% precision, which provides a stringent test of Gaia parallaxes. It is shown that combining spectroscopic and interferometric observations of nearby eclipsing binaries can lead to extremely accurate parallaxes and stellar parameters.

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.

Shack-Hartmann wavefront sensing relies on accurate spot centre measurement. Several algorithms were developed with this aim, mostly focused on precision, i.e. minimizing random errors. In the solar and extended scene community, the importance of the accuracy (bias error due to peak-locking, quantisation or sampling) of the centroid determination was identified and solutions proposed. But these solutions only allow partial bias corrections. To date, no systematic study of the bias error was conducted. This article bridges the gap by quantifying the bias error for different correlation peak-finding algorithms and types of sub-aperture images and by proposing a practical solution to minimize its effects. Four classes of sub-aperture images (point source, elongated laser guide star, crowded field and solar extended scene) together with five types of peak-finding algorithms (1D parabola, the centre of gravity, Gaussian, 2D quadratic polynomial and pyramid) are considered, in a variety of signal-to-noise conditions. The best performing peak-finding algorithm depends on the sub-aperture image type, but none is satisfactory to both bias and random errors. A practical solution is proposed that relies on the anti-symmetric response of the bias to the sub-pixel position of the true centre. The solution decreases the bias by a factor of ~7 to values of < 0.02 pix. The computational cost is typically twice of current cross-correlation algorithms.

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.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

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.

Assessing the quality of aperture synthesis maps is relevant for benchmarking image reconstruction algorithms, for the scientific exploitation of data from optical long-baseline interferometers, and for the design/upgrade of new/existing interferometric imaging facilities. Although metrics have been proposed in these contexts, no systematic study has been conducted on the selection of a robust metric for quality assessment. This article addresses the question: what is the best metric to assess the quality of a reconstructed image? It starts by considering several metrics, and selecting a few based on general properties. Then, a variety of image reconstruction cases is considered. The observational scenarios are phase closure and phase referencing at the Very Large Telescope Interferometer (VLTI), for a combination of two, three, four and six telescopes. End-to-end image reconstruction is accomplished with the MiRA software, and several merit functions are put to test. It is found that convolution by an effective point spread function is required for proper image quality assessment. The effective angular resolution of the images is superior to naive expectation based on the maximum frequency sampled by the array. This is due to the prior information used in the aperture synthesis algorithm and to the nature of the objects considered. The l1 norm is the most robust of all considered metrics, because being linear it is less sensitive to image smoothing by high regularisation levels. For the cases considered, this metric allows the implementation of automatic quality assessment of reconstructed images, with a performance similar to human selection.

Z CMa is a complex pre-main sequence binary with a current separation of 100 mas, known to consist of an FU Orionis star (SE component) and an embedded Herbig Be star (NW component). Immediately when the late-2008 outburst of Z CMa was announced to the community, we initiated a high angular resolution imaging campaign with VLT/NaCo, Keck/NIRC2, VLT/SINFONI, and Keck/OSIRIS which aimed at characterizing the outburst of both components of the system in the near-infrared. We confirm that the NW star dominates the system flux in the 1.1-3.8 microns range and is responsible for the photometric outburst. We extract the first medium-resolution (R=2000-4000) near-infrared (1.1-2.4 microns) spectra of the individual components during and after the outburst. The SE component has a spectrum typical of FU Orionis objects. The NW component spectrum is characteristic of embedded outbursting protostars and EX Or objects. It displays numerous emission lines during the outburst whose intensity correlates with the system activity. In particular, we find a correlation between the Brackett gamma equivalent width and the system brightness. The bluing of the continuum of the NW component along with the absolute flux and color-variation of the system during the outburst suggests that the outburst was caused by a complex interplay between a variation of the extinction in the line of sight of the NW component on one hand, and the emission of shocked regions close to the NW component on the other. We confirm the recently reported wiggling of the SE component jet from [Fe II] line emission. We find a point-like structure associated with a peak emission at 2.098 microns coincidental with the clump or arm seen in broadband polarization differential imaging as well as additional diffuse emission along a PA=214 degrees. The origin of these two structures is unclear and deserves further investigation.

FU Orionis objects are a class of young stars with important bursts in luminosity and which show evidence of accretion and ejection activity. It is generally accepted that they are surrounded by a Keplerian circumstellar disk and an infalling envelope. The outburst would occurs because of a sudden increase in the accretion rate. We aim at studying the regions closer to the central star in order to observe the signs of the accretion/ejection activity. We present optical observations of the Halpha line using the Integral Field Spectrograph OASIS, at the William Herschel Telescope, combined with Adaptive Optics. Since this technique gives the spectral information for both spatial directions, we carried out a two-dimensional spectro-astrometric study of the signal. We measured a clear spectro-astrometric signal in the North-South direction. The cross-correlation between the spectra showed a spatial distribution in velocity suggestive of scattering by a disk surrounding the star. This would be one of the few spatial inferences of a disk observed in a FU Orionis object. However, in order to fully understand the observed structure, higher angular and spectral resolution observations are required. V1515 Cyg appears now as an important object to be observed with a new generation of instruments to increase our knowledge about the disk and outflows structure in FU Orionis objects.

The Columbia - U. de Chile CO Survey of the Southern Milky Way is used for separating the CO(1-0) emission of the fourth Galactic quadrant within the solar circle into its dominant components, giant molecular clouds (GMCs). After the subtraction of an axisymmetric model of the CO background emission in the inner Southern Galaxy, 92 GMCs are identified, and for 87 of them the two-fold distance ambiguity is solved. Their total molecular mass is M(H2) = 1.14 +/- 0.05 x 10^8 Msolar accounting for around 40% of the molecular mass estimated from an axisymmetric analysis of the H2 volume density in the Galactic disk (Bronfman et al. 1988b) M(H2) disk = 3.03 x 10^8 Msolar. The large scale spiral structure in the Southern Galaxy, within the solar circle, is traced by the GMCs in our catalog; 3 spiral arm segments: the Centaurus, Norma, and 3-kpc expanding arm are analyzed. After fitting a logarithmic spiral arm model to the arms, tangent directions at 310{\deg}, 330{\deg}, and 338{\deg}, respectively, are found, consistent with previous values from the literature. A complete CS(2-1) survey toward IRAS point-like sources with FIR colors characteristic of UC HII regions is used to estimate the massive star formation rate per unit H2 mass (MSFR), and the massive star formation efficiency for GMCs. The average MSFR for GMCs is 0.41 +/- 0.06 Lsolar/Msolar, and for the most massive clouds in the Norma arm it is 0.58 +/- 0.09 Lsolar/Msolar. Massive star formation efficiencies of GMCs are on average 3% of their available molecular mass.

Active pre-main-sequence binaries with separations of around ten stellar radii present a wealth of phenomena unobserved in common systems. The study of these objects is extended from Classical T Tauri stars to the Herbig Ae star HD 104237. Spectro-interferometry with the VLTI/AMBER is presented. It is found that the K-band continuum squared visibilities are compatible with a circumbinary disc with a radius of ~0.5 AU. However, a significant fraction (~50 per cent) of the flux is unresolved and not fully accounted by the stellar photospheres. The stars probably don't hold circumstellar discs, in addition to the circumbinary disk, due to the combined effects of inner magnetospheric truncation and outer tidal truncation. This unresolved flux likely arises in compact structures inside the tidally disrupted circumbinary disc. Most ($\gtrsim 90$ per cent) of the Br gamma line emission is unresolved. The line-to-continuum spectro-astrometry shifts in time, along the direction of the Ly alpha jet known to be driven by the system. The shift is anti-correlated with the Br gamma equivalent width. It is shown that the unresolved Br gamma emission cannot originate in the jet but instead is compatible with stellar emission from the orbiting binary components. The increase in the absolute value of the equivalent width of the line takes place at periastron passage; it could arise in an accretion burst, a flare or in the increase in effective size of the emission region by the interaction of the magnetospheres. The binary longitude of the ascending node is found to be $\Omega=(235\pm3)\degr$ and the orbit retrograde. The origin of the jet is revisited. The tidal disruption of the circumstellar disks creates difficulties to ejection models that rely on stellar magnetosphere and disc coupling. A scenario of a stellar wind collimated by a circumbinary disc wind is suggested.

Molecular counterparts to atomic jets have been detected within 1000 AU of young stars. Reproducing them is a challenge for proposed ejection models. We explore whether molecules may survive in an MHD disk wind invoked to reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. The coupled ionization, chemical and thermal evolution along dusty flow streamlines is computed for a prescribed MHD disk wind solution, using a method developed for magnetized shocks in the interstellar medium. Irradiation by wind-attenuated coronal X-rays and FUV photons from accretion hot spots is included, with self-shielding of H2 and CO. Disk accretion rates of 5e-6, 1e-6 and 1e-7 solar masses per year are considered, representative of low-mass young protostars (Class 0), evolved protostars (Class I) and very active T Tauri stars (Class II). The disk wind has an onion-like thermo-chemical structure, with streamlines launched from larger radii having lower temperature and ionisation, and higher H2 abundance. The coupling between charged and neutral fluids is sufficient to eject molecules from the disk out to 9 AU. The launch radius beyond which most H2 survives moves outward with evolutionary stage. CO survives in the Class 0 but is significantly photodissociated in the Class I/II. Balance between ambipolar heating and molecular cooling establishes an asymptotic temperature 700-3000 K, with cooler jets at earlier protostellar stages. Endothermic formation of H2O is efficient with abundances up to 1e-4, while CH+ and SH+ can exceed 1e-6 in the Class I/II winds. A centrifugal MHD disk wind launched from beyond 0.2-1 AU can produce molecular jets/winds up to speeds 100 km/s in young low-mass stars. The model predicts a high ratio H2/CO and an increase of molecular launch radius, temperature, and flow width as the source evolves, in agreement with current observed trends.

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys.

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.

The Z CMa binary is understood to undergo both FU Orionis (FUOR) and EX Orionis (EXOR) type outbursts. While the SE component has been spectro- scopically identified as an FUOR, the NW component, a Herbig Be star, is the source of the EXOR outbursts. The system has been identified as the source of a large outflow, however, previous studies have failed to identify the driver. Here we present adaptive optics (AO) assisted [FeII] spectro-images which reveal for the first time the presence of two jets. Observations made using OSIRIS at the Keck Observatory show the Herbig Be star to be the source of the parsec-scale outflow, which within 2'' of the source shows signs of wiggling and the FUOR to '' be driving a ~ 0.4 jet. The wiggling of the Herbig Be star's jet is evidence for an additional companion which could in fact be generating the EXOR outbursts, the last of which began in 2008 (Grankin & Artemenko 2009). Indeed the dy- namical scale of the wiggling corresponds to a time-scale of 4-8 years which is in agreement with the time-scale of these outbursts. The spectro-images also show a bow-shock shaped feature and possible associated knots. The origin of this structure is as of yet unclear. Finally interesting low velocity structure is also observed. One possibility is that it originates in a wide-angle outflow launched from a circumbinary disk.

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.

Classically, optical and near-infrared interferometry have relied on closure phase techniques to produce images. Such techniques allow us to achieve modest dynamic ranges. In order to test the feasibility of next generation optical interferometers in the context of the VLTI-spectro-imager (VSI), we have embarked on a study of image reconstruction and analysis. Our main aim was to test the influence of the number of telescopes, observing nights and distribution of the visibility points on the quality of the reconstructed images. Our results show that observations using six Auxiliary Telescopes (ATs) during one complete night yield the best results in general and is critical in most science cases; the number of telescopes is the determining factor in the image reconstruction outcome. In terms of imaging capabilities, an optical, six telescope VLTI-type configuration and ~200 meter baseline will achieve 4 mas spatial resolution, which is comparable to ALMA and almost 50 times better than JWST will achieve at 2.2 microns. Our results show that such an instrument will be capable of imaging, with unprecedented detail, a plethora of sources, ranging from complex stellar surfaces to microlensing events.

One of the aims of next generation optical interferometric instrumentation is to be able to make use of information contained in the visibility phase to construct high dynamic range images. Radio and optical interferometry are at the two extremes of phase corruption by the atmosphere. While in radio it is possible to obtain calibrated phases for the science objects, in the optical this is currently not possible. Instead, optical interferometry has relied on closure phase techniques to produce images. Such techniques allow only to achieve modest dynamic ranges. However, with high contrast objects, for faint targets or when structure detail is needed, phase referencing techniques as used in radio interferometry, should theoretically achieve higher dynamic ranges for the same number of telescopes. Our approach is not to provide evidence either for or against the hypothesis that phase referenced imaging gives better dynamic range than closure phase imaging. Instead we wish to explore the potential of this technique for future optical interferometry and also because image reconstruction in the optical using phase referencing techniques has only been performed with limited success. We have generated simulated, noisy, complex visibility data, analogous to the signal produced in radio interferometers, using the VLTI as a template. We proceeded with image reconstruction using the radio image reconstruction algorithms contained in AIPS IMAGR (CLEAN algorithm). Our results show that image reconstruction is successful in most of our science cases, yielding images with a 4 milliarcsecond resolution in K band. (abridged)

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.

Asteroseismology provides us with a unique opportunity to improve our understanding of stellar structure and evolution. Recent developments, including the first systematic studies of solar-like pulsators, have boosted the impact of this field of research within Astrophysics and have led to a significant increase in the size of the research community. In the present paper we start by reviewing the basic observational and theoretical properties of classical and solar-like pulsators and present results from some of the most recent and outstanding studies of these stars. We centre our review on those classes of pulsators for which interferometric studies are expected to provide a significant input. We discuss current limitations to asteroseismic studies, including difficulties in mode identification and in the accurate determination of global parameters of pulsating stars, and, after a brief review of those aspects of interferometry that are most relevant in this context, anticipate how interferometric observations may contribute to overcome these limitations. Moreover, we present results of recent pilot studies of pulsating stars involving both asteroseismic and interferometric constraints and look into the future, summarizing ongoing efforts concerning the development of future instruments and satellite missions which are expected to have an impact in this field of research.

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.

Most of the observed emission lines and continuum excess from young accreting low mass stars (Classical T Tauri stars -- CTTSs) take place in the star-disk or inner disk region. These regions have a complex emission topology still largely unknown. In this paper the magnetospheric accretion and inner wind contributions to the observed permitted He and H near infrared (NIR) lines of the bright southern CTTS RU Lupi are investigated for the first time. Previous optical observations of RU Lupi showed a large H-alpha profile, due to the emission from a wind in the line wings, and a micro-jet detected in forbidden lines. We extend this analysis to NIR lines through seeing-limited high spectral resolution spectra taken with VLT/ISAAC, and adaptive optics (AO) aided narrow-band imaging and low spectral resolution spectroscopy with VLT/NACO. Using spectro-astrometric analysis we investigate the presence of extended emission down to very low spatial scales (a few AU). The HeI 10830 line presents a P Cygni profile whose absorption feature indicates the presence of an inner stellar wind. Moreover the spectro-astrometric analysis evidences the presence of an extended emission superimposed to the absorption feature and likely coming from the micro-jet detected in the optical. On the contrary, the origin of the Hydrogen Paschen and Brackett lines is difficult to address. We tried tentatively to explain the observed line profiles and flux ratios with both accretion and wind models showing the limits of both approaches. The lack of spectro-astrometric signal indicates that the HI emission is either compact or symmetric. Our analysis confirms the sensitivity of the HeI line to the presence of faint extended emission regions in the close proximity of the star.

Preparing observing proposals is a classical exercise in an observational astronomer's life. In this paper, we summarize the general requirements to write a good proposal for the ESO telescopes, and eventually obtain telescope time. We detail in particular the different proposal types, how they are evaluated by ESO, what they should contain and what to avoid. The advice given in the present paper are applicable to observing proposals for the AMBER and MIDI instruments of the VLTI, but also to other ESO instruments, and to some extent to other observatories.

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.

VITRUV is a second generation spectro-imager for the PRIMA enabled Very Large Telescope Interferometer. By combining simultaneously up to 8 telescopes VITRUV makes the VLTI up to 6 times more efficient. This operational gain allows two novel scientific methodologies: 1) massive surveys of sizes; 2) routine interferometric imaging. The science cases presented concentrate on the qualitatively new routine interferometric imaging methodology. The science cases are not exhaustive but complementary to the PRIMA reference mission. The focus is on: a) the close environment of young stars probing for the initial conditions of planet formation and disk evolution; b) the surfaces of stars tackling dynamos, activity, pulsation, mass-loss and evolution; c) revealing the origin of the extraordinary morphologies of Planetary Nebulae and related stars; d) studying the accretion-ejection structures of stellar black-holes (microquasars) in our galaxy; e) unveiling the different interacting components (torus, jets, BLRs) of Active Galactic Nuclei; and f) probing the environment of nearby supermassive black-holes and relativistic effects in the Galactic Center black-hole.

The VITRUV project has the objective to deliver milli-arcsecond spectro-images of the environment of compact sources like young stars, active galaxies and evolved stars to the community. This instrument of the VLTI second generation based on the integrated optics technology is able to combine from 4 to 8 beams from the VLT telescopes. Working primarily in the near infrared, it will provide intermediate to high spectral resolutions and eventually polarization analysis. This paper summarizes the result from the concept study led within the Joint Research Activity advanced instruments of the OPTICON program.

Motivated by the first detections of rotation signatures in the DG Tau jet (Bacciotti et al. 2002), we examine possible biases affecting the relation between detected rotation signatures and true azimuthal velocity for self-similar MHD disc winds, taking into account projection, convolution as well as excitation gradients effects. We find that computed velocity shifts are systematically smaller than the true underlying rotation curve. When outer slower streamlines dominate the emission, we predict observed shifts increasing with transverse distance to the jet axis, opposite to the true rotation profile. Determination of the full transverse rotation profile thus requires high angular resolution observations (< 5 AU) on an object with dominant inner faster streamlines. Comparison of our predictions with HST/STIS observations of DG Tau clearly shows that self-similar, warm MHD disc wind models with lambda = 13 and an outer radius of the disc ~ 3 AU are able to reproduce detected velocity shifts, while cold disc wind models (lambda > 50) are ruled out for the medium-velocity component in the DG Tau jet.

We investigate the emissivity properties of the main near-IR transitions of the Fe+ ion in the conditions prevailing in the inner regions of jets from young stars, based on a simplified 16-level atom model. We present new diagnostic diagrams involving prominent near-IR line ratios that allow us to constrain the electronic density, temperature, and Fe gas phase abundance ratio, independently of the heating process. Comparison with recent near-IR observations of a sample of HH objects indicates gas phase Fe abundances ranging from 15-50 % up to 100 % of the solar value in agreement with the moderate depletions previously derived from optical line ratios or shock models. Hence, it appears that Fe-bearing dust is efficiently destroyed in stellar jets. We then use our Fe+ emissivity model to predict near-IR [Fe II] emission maps for self-similar, cold MHD disk wind models. We show that observations in [Fe II] with AMBER on the VLTI could severely constrain the MHD solution and the inner launch radius of the jet. We also compare theoretical predictions with recent observations in the near-IR [Fe II] lines of the L1551-IRS5 and DG Tau jets. The cold disk wind model reproduces quite well the two velocity components observed at -100 and -300 km/s, although the high velocity component appears overestimated by a factor of 1.5 in the DG Tau jet. However, the model predicts too little emission at intermediate-velocity and insufficient densities. Denser disk winds with stronger heating at the jet base, which have been invoked for optical jets, also appear needed in younger, embedded Class I jet sources.

Motivated by recent subarcsecond resolution observations of jets from T Tauri stars, we extend the work of \citep{Safier93a, Safier93b} by computing the thermal and ionization structure of self-similar, magnetically-driven, atomic disk winds heated by ambipolar diffusion. Improvements over his work include: (1) new magnetized cold jet solutions consistent with the underlying accretion disk \citep{Ferreira97}; (2) a more accurate treatment of ionization and ion-neutral momentum exchange rates; and (3) predictions for spatially resolved forbidden line emission (maps, long-slit spectra, and line ratios), presented in a companion paper, \cite{Garcia2001b}. As in \citep{Safier93a}, we obtain jets with a temperature plateau around $10^4$ K, but ionization fractions are revised downward by a factor of 10-100. This is due to previous omission of thermal speeds in ion-neutral momentum-exchange rates and to different jet solutions. The physical origin of the hot temperature plateau is outlined. In particular we present three analytical criteria for the presence of a hot plateau, applicable to any given MHD wind solution where ambipolar diffusion and adiabatic expansion are the dominant heating and cooling terms. We finally show that, for solutions favored by observations, the jet thermal structure remains consistent with the usual approximations used for MHD jet calculations (thermalized, perfectly conducting, single hydromagnetic cold fluid calculations).

Thermal and ionization structures of self-similar, magnetically-driven disk winds obtained in a companion paper \cite{Garcia2001a} are used to compute a series of jet synthetic observations. These include spatially resolved forbidden line emission maps, long-slit spectra, as well as line ratios. Line profiles and jet widths appear to be good tracers of the wind dynamics and collimation, whereas line ratios essentially trace gas excitation conditions. All the above diagnostics are confronted to observations of T Tauri star microjets. Convolution by the observing beam is shown to be essential for a meaningful test of the models. We find that jet widths and qualitative variations in line profiles with both distance and line tracers are well reproduced. A low-velocity [\ion{O}{i}] component is also obtained, originating from the disk wind base. However, this component is too weak, predicted maximum velocities are too high and electronic and total densities are too low. Denser and slower magnetized winds, launched from disks with warm chromospheres, might resolve these discrepancies.

We summarize results on the thermal and ionization structure of self-similar, magnetically-driven, atomic disk winds heated by ambipolar diffusion. We improve upon earlier work by Safier by considering (1) new MHD solutions consistent with underlying cold keplerian disk equilibrium, (2) a more accurate treatment of the micro-physics, and (3) predictions for spatially resolved forbidden line emission (maps, long-slit spectra). The temperature plateau $\simeq 10^4$ K found earlier is recovered, but ionization fractions are revised downward by a factor of 10, due to previous omission of thermal speeds in ion-neutral momentum-exchange rates. The physical origin of the temperature plateau is outlined. Predictions are then compared with T~Tauri star observations, with emphasis on the necessity of suitable beam convolution. Jet widths and variations in line profiles with distance and line tracer are well reproduced. However, predicted maximum velocities are too high, total densities too low, and the low-velocity [OI] component is too weak. Denser, slower MHD winds from warm disks might resolve these discrepancies.

Using time resolved 2-dimensional aperture photometry we have established that the optical candidate for PSR 1509-58 does not pulse. Our pulsed upper limits (m_V = 24.3 and m_B = 25.7) put severe constraints on this being the optical counterpart. Furthermore the colours of the candidate star are consistent with a main sequence star at a distance of 2-4 kpc. The probability of a chance coincidence with a normal star and the difficulty of explaining the lack of pulsed emission leads us to conclude that this object is an intermediate field star.