We present the Mid-infrared stellar Diameters and Fluxes compilation Catalogue (MDFC) dedicated to long-baseline interferometry at mid-infrared wavelengths (3-13 mum). It gathers data for half a million stars, i.e. nearly all the stars of the Hipparcos-Tycho catalogue whose spectral type is reported in the SIMBAD database. We cross-match 26 databases to provide basic information, binarity elements, angular diameter, magnitude and flux in the near and mid-infrared, as well as flags that allow us to identify the potential calibrators. The catalogue covers the entire sky with 465 857 stars, mainly dwarfs and giants from B to M spectral types closer than 18 kpc. The smallest reported values reach 0.16 muJy in L and 0.1 muJy in N for the flux, and 2 microarcsec for the angular diameter. We build 4 lists of calibrator candidates for the L- and N-bands suitable with the Very Large Telescope Interferometer (VLTI) sub- and main arrays using the MATISSE instrument. We identify 1 621 candidates for L and 44 candidates for N with the Auxiliary Telescopes (ATs), 375 candidates for both bands with the ATs, and 259 candidates for both bands with the Unit Telescopes (UTs). Predominantly cool giants, these sources are small and bright enough to belong to the primary lists of calibrator candidates. In the near future, we plan to measure their angular diameter with 1% accuracy.

We present in this paper the general formalism and data processing steps used in the MATISSE data reduction software, as it has been developed by the MATISSE consortium. The MATISSE instrument is the mid-infrared new generation interferometric instrument of the Very Large Telescope Interferometer (VLTI). It is a 2-in-1 instrument with 2 cryostats and 2 detectors: one 2k x 2k Rockwell Hawaii 2RG detector for L\&M-bands, and one 1k x 1k Raytheon Aquarius detector for N-band, both read at high framerates, up to 30 frames per second. MATISSE is undergoing its first tests in laboratory today.

We analyse interferometric data obtained for Regulus with AMBER (Astronomical Multi- BEam combineR) at high spectral resolution ($\lambda/\delta\lambda \approx 12000$) across the Br$\gamma$ spectral line. The study of the photocentre displacement allows us to constrain a large number of stellar parameters -- equatorial radius $R_{\rm eq}$, equatorial velocity $V_{\rm eq}$, inclination $i$, rotation-axis position angle $PA_{\rm rot}$, and flattening -- with an estimation of gravity-darkening coefficient $\beta$ using previously published theoretical results. We use the Simulation Code of Interferometric-observations for ROtators and CirCumstellar Objects (SCIROCCO), a semi-analytical algorithm dedicated to fast rotators. We chose Regulus because it is a very well-known edge-on star, for which an alternative approach is needed to check the previously published results. Our analysis showed that a significant degeneracy of solution is present. By confronting the results obtained by differential interferometry with those obtained by conventional long-base interferometry, we obtain similar results (within the uncertainties), thereby validating our approach, where $V_{eq}$ and $i$ are found separately. From the photocentre displacement, we can independently deduce $PA_{rot}$. We use two minimization methods to restrict observed stellar parameters via a fast rotator model: a non-stochastic method ($\chi^2$ fit) and a stochastic one (Markov Chain Monte Carlo method), in order to check whether the correct global minimum is achieved particularly with respect to the degeneracies of the gravity darkening parameter $\beta$, where we demonstrate, using a quantitative analysis of parameters, that the estimate of $\beta$ is easier for stars with an inclination angle of around $45^\circ$.

The rapidly rotating primary component of Regulus A system has been observed, for the first time, using the technique of differential interferometry at high spectral resolution. The observations have been performed across the Br$_\gamma$ spectral line with the VLTI/AMBER focal instrument in high spectral resolution mode (R $\approx$ 12000) at $\approx$ 80-130m (projected on the sky) Auxiliary Telescopes triplet baseline configurations. We confirm, within the uncertainties, the results previously obtained using the techniques of classical long-baseline interferometry, although the question of anomalous gravity darkening remains open for the future study.

The full scientific potential of the VLTI with its second generation instruments MATISSE and GRAVITY require fringe tracking up to magnitudes K>14 with the UTs and K>10 with the ATs. The GRAVITY fringe tracker (FT) will be limited to K~10.5 with UTs and K~7.5 with ATs, for fundamental conceptual reasons: the flux of each telescope is distributed among 3 cophasing pairs and then among 5 spectral channels for coherencing. To overcome this limit we propose a new FT concept, called Hierarchical Fringe Tracker (HFT) that cophase pairs of apertures with all the flux from two apertures and only one spectral channel. When the pair is cophased, most of the flux is transmitted as if it was produced by an unique single mode beam and then used to cophase pairs of pairs and then pairs of groups. At the deeper level, the flux is used in an optimized dispersed fringe device for coherencing. On the VLTI such a system allows a gain of about 3 magnitudes over the GRAVITY FT. On interferometers with more apertures such as CHARA (6 telescopes) or a future Planet Formation Imager (12 to 20 telescopes), the HFT would be even more decisive, as its performance does not decrease with the number of apertures. It would allow building a PFI reaching a coherent magnitude H~10 with 16 apertures with diameters smaller than 2 m. We present the HFT concept, the first steps of its feasibility demonstration from computer simulations and the optical design of a 4 telescopes HFT prototype.

MATISSE (Multi AperTure mid-Infrared SpectroScopic Experiment) is the next generation spectro-interferometer at the European Southern Observatory VLTI operating in the spectral bands L, M and N, and combining four beams from the unit and auxiliary telescopes. MATISSE is now fully integrated at the Observatoire de la C\^ote d'Azur in Nice (France), and has entered very recently its testing phase in laboratory. This paper summarizes the equations describing the MATISSE signal and the associated sources of noise. The specifications and the expected performances of the instrument are then evaluated taking into account the current characteristics of the instrument and the VLTI infrastructure, including transmission and contrast degradation budgets. In addition, we present the different MATISSE simulation tools that will be made available to the future users.

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.

We analyse the main physical parameters and the circumstellar environment of the young Herbig Be star HD 98922. We present AMBER/VLTI high spectral resolution (R =12000) interferometric observations across the Br$\gamma$ line, accompanied by UVES high-resolution spectroscopy and SINFONI-AO assisted near-infrared integral field spectroscopic data. To interpret our observations, we develop a magneto-centrifugally driven disc-wind model. Our analysis of the UVES spectrum shows that HD 98922 is a young (~5x10^5 yr) Herbig Be star (SpT=B9V), located at a distance of 440(+60-50) pc, with a mass accretion rate of ~9+/-3x10^(-7) M_sun yr^(-1). SINFONI K-band AO-assisted imaging shows a spatially resolved circumstellar disc-like region (~140 AU in diameter) with asymmetric brightness distribution. Our AMBER/VLTI UT observations indicate that the Br$\gamma$ emitting region (radius ~0.31+/-0.04 AU) is smaller than the continuum emitting region (inner dust radius ~0.7+/-0.2 AU), showing significant non-zero V-shaped differential phases (i.e. non S-shaped, as expected for a rotating disc). The value of the continuum-corrected pure Br$\gamma$ line visibility at the longest baseline (89 m) is ~0.8+/-0.1, i.e. the Br$\gamma$ emitting region is partially resolved. Our modelling suggests that the observed Br$\gamma$ line-emitting region mainly originates from a disc wind with a half opening angle of 30deg, and with a mass-loss rate of ~2x10(-7) M_sun yr^(-1). The observed V-shaped differential phases are reliably reproduced by combining a simple asymmetric continuum disc model with our Br$\gamma$ disc-wind model. The Br$\gamma$ emission of HD 98922 can be modelled with a disc wind that is able to approximately reproduce all interferometric observations if we assume that the intensity distribution of the dust continuum disc is asymmetric.

The limiting magnitude is a key issue for optical interferometry. Pairwise fringe trackers based on the integrated optics concepts used for example in GRAVITY seem limited to about K=10.5 with the 8m Unit Telescopes of the VLTI, and there is a general "common sense" statement that the efficiency of fringe tracking, and hence the sensitivity of optical interferometry, must decrease as the number of apertures increases, at least in the near infrared where we are still limited by detector readout noise. Here we present a Hierarchical Fringe Tracking (HFT) concept with sensitivity at least equal to this of a two apertures fringe trackers. HFT is based of the combination of the apertures in pairs, then in pairs of pairs then in pairs of groups. The key HFT module is a device that behaves like a spatial filter for two telescopes (2TSF) and transmits all or most of the flux of a cophased pair in a single mode beam. We give an example of such an achromatic 2TSF, based on very broadband dispersed fringes analyzed by grids, and show that it allows piston measures from very broadband fringes with only 3 to 5 pixels per fringe tracker. We show the results of numerical simulation indicating that our device is a good achromatic spatial filter and allowing a first evaluation of its coupling efficiency, which is similar to this of a single mode fiber on a single aperture. Our very preliminary results indicate that HFT has a good chance to be a serious candidate for the most sensitive fringe tracking with the VLTI and also interferometers with much larger number of apertures. On the VLTI the first rough estimate of the magnitude gain with regard to the GRAVITY internal FT is between 2.5 and 3.5 magnitudes in K, with a decisive impact on the VLTI science program for AGNs, Young stars and planet forming disks.

Unveiling the structure of the Broad Line Region (BLR) of AGNs is critical to understand the quasar phenomenon. Resolving a few BLRs by optical interferometry will bring decisive information to confront, complement and calibrate the reverberation mapping technique, basis of the mass-luminosity relation in quasars. BLRs are much smaller than the angular resolution of the VLT and Keck interferometers and they can be resolved only by differential interferometry very accurate measurements of differential visibility and phase as a function of wavelength. The latter yields the photocenter variation with wavelength, and constrains the size, position and velocity law of various regions of the BLR. AGNs are below the magnitude limit for spectrally resolved interferometry set by currently available fringe trackers. A new "blind" observation method and a data processing based on the accumulation of 2D Fourier power and cross spectra permitted us to obtain the first spectrally resolved interferometric observation of a BLR, on the K=10 quasar 3C273. A careful bias analysis is still in progress, but we report strong evidence that, as the baseline increases, the differential visibility decreases in the Pa-alpha line. Combined with a differential phase smaller than 3 degree, this yields an angular equivalent radius of the BLR larger than 0.4 milliarcseconds, or 1000 light days at the distance of 3C273, much larger than the reverberation mapping radius of 300 light days. Explaining the coexistence of these two different sizes, and possibly structures and mechanisms, implies very new insights into the BLR of 3C273

The UX Ori type variables (named after the prototype of their class) are intermediate-mass pre-main sequence objects. One of the most likely causes of their variability is the obscuration of the central star by orbiting dust clouds. We investigate the structure of the circumstellar environment of the UX~Ori star V1026 Sco (HD 142666) and test whether the disk inclination is large enough to explain the UX Ori variability. We observed the object in the low-resolution mode of the near-infrared interferometric VLTI/AMBER instrument and derived H- and K-band visibilities and closure phases. We modeled our AMBER observations, published Keck Interferometer observations, archival MIDI/VLTI visibilities, and the spectral energy distribution using geometric and temperature-gradient models. Employing a geometric inclined-ring disk model, we find a ring radius of 0.15 +- 0.06 AU in the H band and 0.18 +- 0.06 AU in the K band. The best-fit temperature-gradient model consists of a star and two concentric, ring-shaped disks. The inner disk has a temperature of 1257^{+133}_{-53} K at the inner rim and extends from 0.19 +- 0.01 AU to 0.23 +- 0.02 AU. The outer disk begins at 1.35^{+0.19}_{-0.20} AU and has an inner temperature of 334^{+35}_{-17} K. The derived inclination of 48.6^{+2.9}_{-3.6}deg approximately agrees with the inclination derived with the geometric model (49 +- 5deg in the K band and 50 +- 11deg in the H band). The position angle of the fitted geometric and temperature-gradient models are 163 +- 9deg (K band; 179 +- 17deg in the H band) and 169.3^{+4.2}_{-6.7}deg, respectively. The narrow width of the inner ring-shaped model disk and the disk gap might be an indication for a puffed-up inner rim shadowing outer parts of the disk. The intermediate inclination of ~50deg is consistent with models of UX Ori objects where dust clouds in the inclined disk obscure the central star.

Pupil-remapping is a new high-dynamic range imaging technique that has recently demonstrated feasibility on sky. The current prototypes present however deceiving limiting magnitude, restricting the current use to the brightest stars in the sky. We propose to combine pupil-remapping with spatio-spectral encoding, a technique first applied to the VEGA/CHARA interferometer. The result is an instrument proposal, called "Dividing Interferometer for Stars Characterizations and Observations" (DISCO). The idea is to take profit of wavelength multiplexing when using a spectrograph in order to pack as much as possible the available information, yet providing a potential boost of 1.5 magnitude if used in existing prototypes. We detail in this paper the potential of such a concept.

We study the sub-AU-scale circumstellar environment of the Herbig Ae star HD144432 with near-infrared (NIR) VLTI/AMBER observations to investigate the structure of its inner dust disk. The interferometric observations were carried out with the AMBER instrument in the H and K band. We interpret the measured H- and K-band visibilities, the near- and mid-infrared visibilities from the literature, and the SED of HD144432 by using geometric ring models and ring-shaped temperature-gradient disk models with power-law temperature distributions. We derived a K-band ring-fit radius of 0.17 \pm 0.01 AU and an H-band radius of 0.18 \pm 0.01 AU (for a distance of 145 pc). This measured K-band radius of \sim0.17 AU lies in the range between the dust sublimation radius of \sim0.13 AU (predicted for a dust sublimation temperature of 1500 K and gray dust) and the prediction of models including backwarming (\sim0.27 AU). We found that an additional extended halo component is required in both the geometric and temperature-gradient modeling. In the best temperature- gradient model, the disk consists of two components. The inner part of the disk is a thin ring with an inner radius of \sim0.21 AU, a temperature of \sim1600 K, and a ring thickness \sim0.02 AU. The outer part extends from \sim1 AU to \sim10 AU with an inner temperature of \sim400 K. We find that the disk is nearly face-on with an inclination angle of < 28 degree. Our temperature-gradient modeling suggests that the NIR excess is dominated by emission from a narrow, bright rim located at the dust sublimation radius, while an extended halo component contributes \sim6% to the total flux at 2 {\mu}m. The MIR model emission has a two-component structure with \sim20% flux from the inner ring and the rest from the outer part. This two-component structure suggests a disk gap, which is possibly caused by the shadow of a puffed-up inner rim.

Circumstellar disks and outflows play a fundamental role in star formation. Infrared spectro-interferometry allows the inner accretion-ejection region to be resolved. We measured interferometric visibilities, wavelength-differential phases, and closure phases of MWC 297 with a spectral resolution of 12000. To interpret our MWC 297 observations, we employed disk-wind models. The measured continuum visibilities confirm previous results that the continuum-emitting region of MWC 297 is remarkably compact. We derive a continuum ring-fit radius of ~2.2 mas (~0.56 AU at a distance of 250 pc), which is ~5.4 times smaller than the 3 AU dust sublimation radius expected for silicate grains (in the absence of radiation-shielding material). The strongly wavelength-dependent and asymmetric Br gamma-emitting region is more extended (~2.7 times) than the continuum-emitting region. At the center of the Br gamma line, we derive a Gaussian fit radius of ~6.3 mas HWHM (~1.6 AU). To interpret the observations, we employ a magneto-centrifugally driven disk-wind model consisting of an accretion disk, which emits the observed continuum radiation, and a disk wind, which emits the Br gamma line. The calculated wavelength-dependent model intensity distributions and Br gamma line profiles are compared with the observations (i.e., K-band spectrum, visibilities, differential phases, and closure phases). The closest fitting model predicts a continuum-emitting disk with an inner radius of ~0.3 AU and a disk wind ejection region with an inner radius of ~0.5 AU (~17.5 stellar radii). We obtain a disk-wind half-opening angle (the angle between the rotation axis and the innermost streamline of the disk wind) of ~80 degrees, which is larger than in T Tau models, and a disk inclination angle of ~20 degrees (i.e., almost pole-on).

The observable quantities in optical interferometry, which are the modulus and the phase of the complex visibility, may be corrupted by parasitic fringes superimposed on the genuine fringe pattern. These fringes are due to an interference phenomenon occurring from straylight effects inside an interferometric instrument. We developed an analytical approach to better understand this phenomenon when straylight causes crosstalk between beams. We deduced that the parasitic interference significantly affects the interferometric phase and thus the associated observables including the differential phase and the closure phase. The amount of parasitic flux coupled to the piston between beams appears to be very influential in this degradation. For instance, considering a point-like source and a piston ranging from $\lambda/500$ to $\lambda/5$ in L band ($\lambda=3.5\:\mu$m), a parasitic flux of about 1\% of the total flux produces a parasitic phase reaching at most one third of the intrinsic phase. The piston, which can have different origins (instrumental stability, atmospheric perturbations, ...), thus amplifies the effect of parasitic interference. According to specifications of piston correction in space or at ground level (respectively $\lambda/500\approx 2$nm and $\lambda/30\approx 100$nm), the detection of hot Jupiter-like planets, one of the most challenging aims for current ground-based interferometers, limits parasitic radiation to about 5\% of the incident intensity. This was evaluated by considering different types of hot Jupiter synthetic spectra. Otherwise, if no fringe tracking is used, the detection of a typical hot Jupiter-like system with a solar-like star would admit a maximum level of parasitic intensity of 0.01\% for piston errors equal to $\lambda$/15. If the fringe tracking specifications are not precisely observed, it thus appears that the allowed level of parasitic intensity dramatically decreases and may prevent the detection. In parallel, the calibration of the parasitic phase by a reference star, at this accuracy level, seems very difficult. Moreover, since parasitic phase is an object-dependent quantity, the use of a hypothetical phase abacus, directly giving the parasitic phase from a given parasitic flux level, is also impossible. Some instrumental solutions, implemented at the instrument design stage for limiting or preventing this parasitic interference, appears to be crucial and are presented in this paper.

While doing optical study in an instrument similar to the interferometers dedicated to the Very Large Telescope (VLT), we have to take care of the pupil and focus conjugations. Modules with artificial sources are designed to simulate the stellar beams, in terms of collimation and pupil location. They constitute alignment and calibration tools. In this paper, we present such a module in which the pupil mask is not located in a collimated beam thus introducing Fresnel diffraction. We study the instrumental contrast taking into account the spatial coherence of the source, and the pupil diffraction. The considered example is MATISSE, but this study can apply to any other instrument concerned with Fresnel diffraction.

This paper presents a study of the atmospheric refraction and its effect on the light coupling efficiency in an instrument using single-mode optical fibers. We show the analytical approach which allowed us to assess the need to correct the refraction in J- and H-bands while observing with an 8-m Unit Telescope. We then developed numerical simulations to go further in calculations. The hypotheses on the instrumental characteristics are those of AMBER (Astronomical Multi BEam combineR), the near infrared focal beam combiner of the Very Large Telescope Interferometric mode (VLTI), but most of the conclusions can be generalized to other single-mode instruments. We used the software package caos (Code for Adaptive Optics Systems) to take into account the atmospheric turbulence effect after correction by the ESO system MACAO (Multi-Application Curvature Adaptive Optics). The opto-mechanical study and design of the system correcting the atmospheric refraction on AMBER is then detailed. We showed that the atmospheric refraction becomes predominant over the atmospheric turbulence for some zenith angles z and spectral conditions: for z larger than 30{\circ} in J-band for example. The study of the optical system showed that it allows to achieve the required instrumental performance in terms of throughput in J- and H-bands. First observations in J-band of a bright star, alpha Cir star, at more than 30{\circ} from zenith clearly showed the gain to control the atmospheric refraction in a single mode instrument, and validated the operating law.

This paper describes the design goals and engineering efforts that led to the realization of AMBER (Astronomical Multi BEam combineR) and to the achievement of its present performance. On the basis of the general instrumental concept, AMBER was decomposed into modules whose functions and detailed characteristics are given. Emphasis is put on the spatial filtering system, a key element of the instrument. We established a budget for transmission and contrast degradation through the different modules, and made the detailed optical design. The latter confirmed the overall performance of the instrument and defined the exact implementation of the AMBER optics. The performance was assessed with laboratory measurements and commissionings at the VLTI, in terms of spectral coverage and resolution, instrumental contrast higher than 0.80, minimum magnitude of 11 in K, absolute visibility accuracy of 1%, and differential phase stability of 1E-3 rad over one minute.

Long baseline interferometry is now a mature technique in the optical domain. Current interferometers are however highly limited in number of sub apertures and concepts are being developed for future generations of very large optical arrays and especially with the goal of direct imaging. In this paper, we study the effects of introducing single-mode fibers in direct imaging optical interferometers. We show how the flexibility of optical fibers is well adapted to the pupil densification scheme. We study the effects of the truncation of the gaussian beams in the imaging process, either in the Fizeau mode or in the densified pupil mode or in the densified image mode. Finally, in the pupil densification configuration, we identify an optimum of the diaphragm width. This optimum maximizes the on-axis irradiance and corresponds to a trade-off between the loss of transmission and the efficiency of the densification.

In this work, we present the first AMBER observations, of the Wolf-Rayet and O (WR+O) star binary system y2 Velorum. The AMBER instrument was used with the telescopes UT2, UT3, and UT4 on baselines ranging from 46m to 85m. It delivered spectrally dispersed visibilities, as well as differential and closure phases, with a resolution R = 1500 in the spectral band 1.95-2.17 micron. We interpret these data in the context of a binary system with unresolved components, neglecting in a first approximation the wind-wind collision zone flux contribution. We show that the AMBER observables result primarily from the contribution of the individual components of the WR+O binary system. We discuss several interpretations of the residuals, and speculate on the detection of an additional continuum component, originating from the free-free emission associated with the wind-wind collision zone (WWCZ), and contributing at most to the observed K-band flux at the 5% level. The expected absolute separation and position angle at the time of observations were 5.1&plusmn;0.9mas and 66&plusmn;15&deg; respectively. However, we infer a separation of 3.62+0.11-0.30 mas and a position angle of 73+9-11&deg;. Our analysis thus implies that the binary system lies at a distance of 368+38-13 pc, in agreement with recent spectrophotometric estimates, but significantly larger than the Hipparcos value of 258+41-31 pc.

AMBER is the first near infrared focal instrument of the VLTI. It combines three telescopes and produces spectrally resolved interferometric measures. This paper discusses some preliminary results of the first scientific observations of AMBER with three Unit Telescopes at medium (1500) and high (12000) spectral resolution. We derive a first set of constraints on the structure of the circumstellar material around the Wolf Rayet Gamma2 Velorum and the LBV Eta Carinae.