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.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

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.

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.

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.

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.

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.

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

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

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

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

GRAVITY is the four-beam, near- infrared, AO-assisted, fringe tracking, astrometric and imaging instrument for the Very Large Telescope Interferometer (VLTI). It is requiring the development of one of the most complex instrument software systems ever built for an ESO instrument. Apart from its many interfaces and interdependencies, one of the most challenging aspects is the overall performance and stability of this complex system. The three infrared detectors and the fast reflective memory network (RMN) recorder contribute a total data rate of up to 20 MiB/s accumulating to a maximum of 250 GiB of data per night. The detectors, the two instrument Local Control Units (LCUs) as well as the five LCUs running applications under TAC (Tools for Advanced Control) architecture, are interconnected with fast Ethernet, RMN fibers and dedicated fiber connections as well as signals for the time synchronization. Here we give a simplified overview of all subsystems of GRAVITY and their interfaces and discuss two examples of high-level applications during observations: the acquisition procedure and the gathering and merging of data to the final FITS file.