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

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

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

While exoplanets are now routinely detected, the detection of small bodies in extrasolar systems remains challenging. Since the discovery of sporadic events interpreted as exocomets (Falling Evaporating Bodies) around $\beta$ Pic in the early 80s, only $\sim$20 stars have been reported to host exocomet-like events. We aim to expand the sample of known exocomet-host stars, as well as to monitor the hot-gas environment around stars with previously known exocometary activity. We have obtained high-resolution optical spectra of a heterogeneous sample of 117 main-sequence stars in the spectral type range from B8 to G8. The data have been collected in 14 observing campaigns expanding over 2 years from both hemispheres. We have analysed the Ca ii K&H and Na i D lines in order to search for non-photospheric absorptions originated in the circumstellar environment, and for variable events that could be caused by outgassing of exocomet-like bodies. We have detected non-photospheric absorptions towards 50% of the sample, attributing a circumstellar origin to half of the detections (i.e. 26% of the sample). Hot circumstellar gas is detected in the metallic lines inspected via narrow stable absorptions, and/or variable blue-/red-shifted absorption events. Such variable events were found in 18 stars in the Ca ii and/or Na i lines; 6 of them are reported in the context of this work for the first time. In some cases the variations we report in the Ca ii K line are similar to those observed in $\beta$ Pic. While we do not find a significant trend with the age or location of the stars, we do find that the probability of finding CS gas in stars with larger vsin i is higher. We also find a weak trend with the presence of near-infrared excess, and with anomalous ($\lambda$ Boo-like) abundances, but this would require confirmation by expanding the sample.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates.

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

Beyond the choice of wavefront control systems or coronographs, advanced data processing methods play a crucial role in disentangling potential planetary signals from bright quasi-static speckles. Among these methods, angular differential imaging (ADI) for data sets obtained in pupil tracking mode (ADI sequences) is one of the foremost research avenues, considering the many observing programs performed with ADI-based techniques and the associated discoveries. Inspired by the field of econometrics, here we propose a new detection algorithm for ADI sequences, deriving from the regime-switching model first proposed in the 1980s. The proposed model is very versatile as it allows the use of PSF-subtracted data sets (residual cubes) provided by various ADI-based techniques, separately or together, to provide a single detection map. The temporal structure of the residual cubes is used for the detection as the model is fed with a concatenated series of pixel-wise time sequences. The algorithm provides a detection probability map by considering two possible regimes for concentric annuli, the first one accounting for the residual noise and the second one for the planetary signal in addition to the residual noise. The algorithm performance is tested on data sets from two instruments, VLT/NACO and VLT/SPHERE. The results show an overall better performance in the receiver operating characteristic space when compared with standard signal-to-noise-ratio maps for several state-of-the-art ADI-based post-processing algorithms.

Context. HD72946 is a bright and nearby solar-type star hosting a low-mass companion at long period (P~16 yr) detected with the radial velocities (RV) method. The companion has a minimum mass of 60.4+/-2.2 MJ and might be a brown dwarf. Its expected semi-major axis of ~243 mas makes it a suitable target for further characterization with high-contrast imaging, in particular to measure its inclination, mass, and spectrum and thus definitely establish its substellar nature. Aims. We aim to further characterize the orbit, atmosphere, and physical nature of HD72946B. Methods. We present high-contrast imaging data in the near-infrared with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. We also use proper motion measurements of the star from Hipparcos and Gaia. Results. The SPHERE data reveal a point source with a contrast of ~9 mag at a projected separation of ~235 mas. No other point sources are detected in the field of view. By jointly fitting the RV, imaging, and proper motion data, we constrain all the orbital parameters of HD72946B and assess a dynamical mass of 72.4+/-1.6 MJ and a semi-major axis of 6.45$^{+0.08}_{-0.07}$ au. Empirical comparison of its IFS spectrum to template dwarfs indicates a spectral type of L5.0+/-1.5. The J-H3 color is close to the expectations of the DUSTY models and suggests a cloudy atmosphere. Comparison with atmospheric models of the spectrophotometry suggests an effective temperature of ~1700 K. The bolometric luminosity (log(L/LS)=-4.11+/-0.10 dex) and dynamical mass of HD72946B are more compatible with evolutionary models for an age range of ~0.9-3 Gyr. The formation mechanism of the companion is currently unclear as it appears slightly away from the bulk of model predictions. HD72946B is currently the closest benchmark brown dwarf companion to a solar-type star with imaging, RV, and proper motion measurements.

METIS is one of the three first-light instruments planned for the ELT, mainly dedicated to high contrast imaging in the mid-infrared. On the SPHERE high-contrast instrument currently installed at the VLT, we observe that one of the main contrast limitations is the wind driven halo, due to the limited AO running speed with respect to the atmospheric turbulence temporal evolution. From this observation, we extrapolate this signature to the ELT/METIS instrument, which is equipped with a single conjugated adaptive optics system and with several coronagraphic devices. By making use of an analytic AO simulator, we compare the amount of wind driven halo observed with SPHERE and with METIS, under the same turbulence conditions.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

This paper is framed within a large project devoted to studying the presence of circumstellar material around main sequence stars, and looking for exocometary events. The work concentrates on HR 10 (A2 IV/V), known for its conspicuous variability in the circumstellar narrow absorption features of Ca II K and other lines, so far interpreted as $\beta$ Pic-like phenomena, within the falling evaporating body scenario. The main goal of this paper is to carry out a thorough study of HR 10 to find the origin of the observed variability, determine the nature of the star, its absolute parameters, and evolutionary status. Interferometric near-infrared (NIR) observations, multi-epoch high-resolution optical spectra spanning a time baseline of more than 32 years, and optical and NIR photometry, together with theoretical modelling, were used to tackle the above objectives. Our results reveal that HR 10 is a binary. The narrow circumstellar absorption features superimposed on the photospheric Ca II K lines -- and lines of other species -- can be decomposed into two or more components, the two deep ones tracing the radial velocity of the individual stars, which implies that their origin cannot be ascribed to transient exocometary events, their variability being fully explained by the binarity of the object. There does not appear to be transient events associated with potential exocomets. Each individual star holds its own circumstellar shell and there are no traces of a circumbinary envelope. The combined use of the interferometric and radial velocity data leads to a complete spectrometric and orbital solution for the binary, the main parameters being: an orbital period of 747.6 days, eccentricities of the orbits around the centre of mass 0.25 (HR 10-A), 0.21 (HR 10-B) and a mass ratio of q=M$_{\rm B}$/M$_{\rm A}$=0.72-0.84. The stars are slightly off the main sequence, the binary being $\sim530$ Myr old.

We present the first observational evidence for a circumplanetary disk around the protoplanet PDS~70~b, based on a new spectrum in the $K$ band acquired with VLT/SINFONI. We tested three hypotheses to explain the spectrum: Atmospheric emission from the planet with either (1) a single value of extinction or (2) variable extinction, and (3) a combined atmospheric and circumplanetary disk model. Goodness-of-fit indicators favour the third option, suggesting circumplanetary material contributing excess thermal emission --- most prominent at $\lambda \gtrsim 2.3 \mu$m. Inferred accretion rates ($\sim 10^{-7.8}$--$10^{-7.3} M_J$ yr$^{-1}$) are compatible with observational constraints based on the H$\alpha$ and Br$\gamma$ lines. For the planet, we derive an effective temperature of 1500--1600 K, surface gravity $\log(g)\sim 4.0$, radius $\sim 1.6 R_J$, mass $\sim 10 M_J$, and possible thick clouds. Models with variable extinction lead to slightly worse fits. However, the amplitude ($\Delta A_V \gtrsim 3$mag) and timescale of variation ($\lesssim$~years) required for the extinction would also suggest circumplanetary material.

Direct imaging of exoplanets is a challenging task as it requires to reach a high contrast at very close separation to the star. Today, the main limitation in the high-contrast images is the quasi-static speckles that are created by residual instrumental aberrations. They have the same angular size as planetary companions and are often brighter, hence hindering our capability to detect exoplanets. Dedicated observation strategies and signal processing techniques are necessary to disentangle these speckles from planetary signals. The output of these methods is a detection map in which the value of each pixel is related to a probability of presence of a planetary signal. The detection map found in the literature relies on the assumption that the residual noise is Gaussian. However, this is known to lead to higher false positive rates, especially close to the star. In this paper, we re-visit the notion of detection map by analyzing the speckle noise distribution, namely the Modified Rician distribution. We use non-asymptotic analysis of the sum of random variables to show that the tail of the distribution of the residual noise decays as an exponential distribution, hence explaining the high false detection rate obtained with the Gaussian assumption. From this analysis, we introduce a novel time domain detection map and we demonstrate its capabilities and the relevance of our approach through experiments on real data. We also provide an empirical rule to determine detection threshold providing a good trade off between true positive and false positive rates for exoplanet detection.

Transition discs are prime targets to look for protoplanets and study planet-disc interactions. We present VLT/SINFONI observations of PDS~70, a transition disc with a recently claimed embedded protoplanet. We take advantage of the angular and spectral diversity present in our data for an optimal PSF modeling and subtraction using principal component analysis (PCA). We report the redetection of PDS 70 b, both the front and far side of the outer disc edge, and the detection of several extended features in the annular gap. We compare spectral differential imaging applied before (PCA-SADI), and after (PCA-ASDI) angular differential imaging. Our tests suggest that PCA-SADI better recovers extended features, while PCA-ASDI is more sensitive to point sources. We adapted the negative fake companion (NEGFC) technique to infer the astrometry of the companion, and derived $r = 193.5 \pm 4.9 \mathrm{mas}$ and PA = 158.7deg $\pm$ 3.0deg. We used both NEGFC and ANDROMEDA to infer the $K$-band spectro-photometry of the protoplanet, and found results consistent with recent VLT/SPHERE observations, except for their 2018/02 epoch measurement in the $K2$ filter. Finally, we derived an upper limit of $\dot{M_b} < 1.26 \times 10^{-7} \big[ \frac{5 M_{\rm Jup}}{M_b} \big] \big[ \frac{R_b}{R_{\rm Jup}}\big] M_{\rm Jup} $ yr$^{-1}$ for the accretion rate of the companion based on an adaptation of PCA-SADI/PCA-ASDI around the Br$\gamma$ line.

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

The ELT is a 39m large, ground-based optical and near- to mid-infrared telescope under construction in the Chilean Atacama desert. Operation is planned to start around the middle of the next decade. All first light instruments will come with wavefront sensing devices that allow control of the ELT's intrinsic M4 and M5 wavefront correction units, thus building an adaptive optics (AO) system. To take advantage of the ELT's optical performance, full diffraction-limited operation is required and only a high performance AO system can deliver this. Further technically challenging requirements for the AO come from the exoplanet research field, where the task to resolve the very small angular separations between host star and planet, has also to take into account the high-contrast ratio between the two objects. We present in detail the results of our simulations and their impact on high-contrast imaging in order to find the optimal wavefront sensing device for the METIS instrument. METIS is the mid-infrared imager and spectrograph for the ELT with specialised high-contrast, coronagraphic imaging capabilities, whose performance strongly depends on the AO residual wavefront errors. We examined the sky and target sample coverage of a generic wavefront sensor in two spectral regimes, visible and near-infrared, to pre-select the spectral range for the more detailed wavefront sensor type analysis. We find that the near-infrared regime is the most suitable for METIS. We then analysed the performance of Shack-Hartmann and pyramid wavefront sensors under realistic conditions at the ELT, did a balancing with our scientific requirements, and concluded that a pyramid wavefront sensor is the best choice for METIS. For this choice we additionally examined the impact of non-common path aberrations, of vibrations, and the long-term stability of the SCAO system including high-contrast imaging performance.

Reference star differential imaging (RDI) is a powerful strategy for high contrast imaging. Using example observations taken with the vortex coronagraph mode of Keck/NIRC2 in $L^\prime$ band, we demonstrate that RDI provides improved sensitivity to point sources at small angular separations compared to angular differential imaging (ADI). Applying RDI to images of the low-mass stellar companions HIP 79124 C (192 mas separation, $\Delta L^\prime$=4.01) and HIP 78233 B (141 mas separation, $\Delta L^\prime$=4.78), the latter a first imaging detection, increases the significance of their detections by up to a factor of 5 with respect to ADI. We compare methods for reference frames selection and find that pre-selection of frames improves detection significance of point sources by up to a factor of 3. In addition, we use observations of the circumstellar disks around MWC 758 and 2MASS J16042165-2130284 to show that RDI allows for accurate mapping of scattered light distributions without self-subtraction artifacts.

The quest for other habitable worlds and the search for life among them are major goals of modern astronomy. One way to make progress towards these goals is to obtain high-quality spectra of a large number of exoplanets over a broad range of wavelengths. While concepts currently investigated in the United States are focused on visible/NIR wavelengths, where the planets are probed in reflected light, a compelling alternative to characterize planetary atmospheres is the mid-infrared waveband (5-20~$\mu$m). Indeed, mid-infrared observations provide key information on the presence of an atmosphere, the surface conditions (e.g., temperature, pressure, habitability), and the atmospheric composition in important species such as H$_2$O, CO$_2$, O$_3$, CH$_4$, and N$_2$O. This information is essential to investigate the potential habitability of exoplanets and to make progress towards the search for life in the universe. Obtaining high-quality mid-infrared spectra of exoplanets from the ground is however extremely challenging due to the overwhelming brightness and turbulence of Earth's atmosphere. In this paper, we present a concept of space-based mid-infrared interferometer that can tackle this observing challenge and discuss the main technological developments required to launch such a sophisticated instrument.

We present the most sensitive direct imaging and radial velocity (RV) exploration of $\epsilon$ Eridani to date. $\epsilon$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 years, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7$\mu$m) with the vortex coronagraph recently installed in W.M. Keck Observatory's infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. $\epsilon$ Eridani b has a mass of $0.78^{+0.38}_{-0.12}$ $M_{Jup}$ and is orbiting $\epsilon$ Eridani at about $3.48\pm 0.02$ AU with a period of $7.37 \pm 0.07$ years. The eccentricity of $\epsilon$ Eridani b's orbit is $0.07^{+0.06}_{-0.05}$, an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet-disk interactions and prospects for future detection and characterization with the James Webb Space Telescope.

Distinct gap features in the nearest protoplanetary disk, TW Hya (distance of 59.5$\pm$0.9 pc), may be signposts of ongoing planet formation. We performed long-exposure thermal infrared coronagraphic imaging observations to search for accreting planets especially within dust gaps previously detected in scattered light and submm-wave thermal emission. Three nights of observations with the Keck/NIRC2 vortex coronagraph in $L^\prime$ (3.4-4.1$\mu$m) did not reveal any statistically significant point sources. We thereby set strict upper limits on the masses of non-accreting planets. In the four most prominent disk gaps at 24, 41, 47, and 88 au, we obtain upper mass limits of 1.6-2.3, 1.1-1.6, 1.1-1.5, and 1.0-1.2 Jupiter masses ($M_J$) assuming an age range of 7-10 Myr for TW Hya. These limits correspond to the contrast at 95\% completeness (true positive fraction of 0.95) with a 1\% chance of a false positive within $1^{\prime\prime}$ of the star. We also approximate an upper limit on the product of planet mass and planetary accretion rate of $M_p\dot{M}\lesssim10^{-8} M_J^2/yr$ implying that any putative $\sim0.1 M_J$ planet, which could be responsible for opening the 24 au gap, is presently accreting at rates insufficient to build up a Jupiter mass within TW Hya's pre-main sequence lifetime.

The circumstellar disk of the Herbig Fe star HD 142527 is host to several remarkable features including a warped inner disk, a 120 au-wide annular gap, a prominent dust trap and several spiral arms. A low-mass companion, HD 142527 B, was also found orbiting the primary star at $\sim$14 au. This study aims to better characterize this companion, which could help explain its impact on the peculiar geometry of the disk. We observed the source with VLT/SINFONI in $H$+$K$ band in pupil-tracking mode. Data were post-processed with several algorithms based on angular differential imaging (ADI). HD 142527 B is conspicuously re-detected in most spectral channels, which enables us to extract the first medium-resolution spectrum of a low-mass companion within 0.1'' from its central star. Fitting our spectrum with both template and synthetic spectra suggests that the companion is a young M2.5$\pm$1.0 star with an effective temperature of $3500\pm100$ K, possibly surrounded with a hot (1700 K) circum-secondary environment. Pre-main sequence evolutionary tracks provide a mass estimate of $0.34\pm0.06 M_{\odot}$, independent of the presence of a hot environment. However, the estimated stellar radius and age do depend on that assumption; we find a radius of $1.37 \pm 0.05 R_{\odot}$ (resp. $1.96 \pm 0.10 R_{\odot}$) and an age of $1.8^{+1.2}_{-0.5}$ Myr (resp. $0.75 \pm 0.25$ Myr) in the case of the presence (resp. absence) of a hot environment contributing in $H$+$K$. Our new values for the mass and radius yield a mass accretion rate of $\sim$5 $\times 10^{-9} M_{\odot}$ yr$^{-1}$ (2-3% that of the primary). Our results illustrate thus the potential for SINFONI+ADI to characterize faint close-in companions. The new spectral type makes HD 142527 B a twin of the well known TW Hya T-Tauri star, and the revision of its mass to higher values further supports its role in shaping the disk.

[abridged] The presence of large amounts of dust in the habitable zones of nearby stars is a significant obstacle for future exo-Earth imaging missions. We executed an N band nulling interferometric survey to determine the typical amount of such exozodiacal dust around a sample of nearby main sequence stars. The majority of our data have been analyzed and we present here an update of our ongoing work. We find seven new N band excesses in addition to the high confidence confirmation of three that were previously known. We find the first detections around Sun-like stars and around stars without previously known circumstellar dust. Our overall detection rate is 23%. The inferred occurrence rate is comparable for early type and Sun-like stars, but decreases from 71% [+11%/-20%] for stars with previously detected mid- to far-infrared excess to 11% [+9%/-4%] for stars without such excess, confirming earlier results at high confidence. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal luminosity function of the dust, we find upper limits on the median dust level around all stars without previously known mid to far infrared excess of 11.5 zodis at 95% confidence level. The corresponding upper limit for Sun-like stars is 16 zodis. An LBTI vetted target list of Sun-like stars for exo-Earth imaging would have a corresponding limit of 7.5 zodis. We provide important new insights into the occurrence rate and typical levels of habitable zone dust around main sequence stars. Exploiting the full range of capabilities of the LBTI provides a critical opportunity for the detailed characterization of a sample of exozodiacal dust disks to understand the origin, distribution, and properties of the dust.

The NIRC2 vortex coronagraph is an instrument on Keck II designed to directly image exoplanets and circumstellar disks at mid-infrared bands $L^\prime$ (3.4-4.1 $\mu$m) and $M_s$ (4.55-4.8 $\mu$m). We analyze imaging data and corresponding adaptive optics telemetry, observing conditions, and other metadata over a three year time period to characterize the performance of the instrument and predict the detection limits of future observations. We systematically process images from 359 observations of 304 unique stars to subtract residual starlight (i.e., the coronagraphic point spread function) of the target star using two methods: angular differential imaging (ADI) and reference star differential imaging (RDI). We find that for the typical parallactic angle (PA) rotation of our dataset ($\sim$10$^{\circ}$), RDI provides gains over ADI for angular separations smaller than 0.25$^{\prime\prime}$. Furthermore, we find a power-law relation between the angular separation from the host star and the minimum PA rotation required for ADI to outperform RDI, with a power-law index of -1.18$\pm$0.08. Finally, we use random forest models to estimate ADI and RDI post-processed detection limits a priori. These models, which we provide publicly on a website, explain 70%-80% of the variance in ADI detection limits and 30%-50% of the variance in RDI detection limits. Averaged over a range of angular separations, our models predict both ADI and RDI contrast to within a factor of 2. These results illuminate important factors in high-contrast imaging observations with the NIRC2 vortex coronagraph, help improve observing strategies, and inform future upgrades to the hardware.

One of the long-term goals of exoplanet science is the (atmospheric) characterization of a large sample (>100) of terrestrial planets to assess their potential habitability and overall diversity. Hence, it is crucial to quantitatively evaluate and compare the scientific return of various mission concepts. Here we discuss the exoplanet yield of a space-based mid-infrared (MIR) nulling interferometer. We use Monte-Carlo simulations, based on the observed planet population statistics from the Kepler mission, to quantify the number and properties of detectable exoplanets (incl. potentially habitable planets) and we compare the results to those for a large aperture optical/NIR space telescope. We investigate how changes in the underlying technical assumptions (sensitivity and spatial resolution) impact the results and discuss scientific aspects that influence the choice for the wavelength coverage and spectral resolution. Finally, we discuss the advantages of detecting exoplanets at MIR wavelengths, summarize the current status of some key technologies, and describe what is needed in terms of further technology development to pave the road for a space-based MIR nulling interferometer for exoplanet science.

Proxima b is our nearest potentially rocky exoplanet and represents a formidable opportunity for exoplanet science and possibly astrobiology. With an angular separation of only 35~mas (or 0.05~AU) from its host star, Proxima b is however hardly observable with current imaging telescopes and future space-based coronagraphs. One way to separate the photons of the planet from those of its host star is to use an interferometer that can easily resolve such spatial scales. In addition, its proximity to Earth and its favorable contrast ratio compared with its host M dwarf (approximately 10$^{-5}$ at 10 microns) makes it an ideal target for a space-based nulling interferometer with relatively small apertures. In this paper, we present the motivation for observing this planet in the mid-infrared (5-20 microns) and the corresponding technological challenges. Then, we describe the concept of a space-based infrared interferometer with relatively small ($<$1m in diameter) apertures that can measure key details of Proxima b, such as its size, temperature, climate structure, as well as the presence of important atmospheric molecules such as H$_2$O, CO$_2$, O$_3$, and CH$_4$. Finally, we illustrate the concept by showing realistic observations using synthetic spectra of Proxima b computed with coupled climate chemistry models.

Exo-zodiacal dust, exozodi for short, is warm (~300K) or hot (up to ~2000K) dust found in the inner regions of planetary systems around main sequence stars. In analogy to our own zodiacal dust, it may be located in or near the habitable zone or closer in, down to the dust sublimation distance. The study of the properties, distribution, and evolution of exozodis can inform about the architecture and dynamics of the innermost regions of planetary systems, close to their habitable zones. On the other hand, the presence of large amounts of exo-zodiacal dust may be an obstacle for future space missions aiming to image Earth-like exoplanets. The dust can be the most luminous component of extrasolar planetary systems, but predominantly emits in the near- to mid-infrared where it is outshone by the host star. Interferometry provides a unique method of separating the dusty from the stellar emission. We discuss the prospects of exozodi observations with the next generation VLTI instruments and summarize critical instrument specifications.

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

The Optimal Optical CoronagraphWorkshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. This contribution is the final part of a series of three papers summarizing the outcomes of the workshop, and presents an overview of novel optical technologies and systems that are implemented or considered for high-contrast imaging instruments on both ground-based and space telescopes. The overall objective of high contrast instruments is to provide direct observations and characterizations of exoplanets at contrast levels as extreme as 10^-10. We list shortcomings of current technologies, and identify opportunities and development paths for new technologies that enable quantum leaps in performance. Specifically, we discuss the design and manufacturing of key components like advanced deformable mirrors and coronagraphic optics, and their amalgamation in "adaptive coronagraph" systems. Moreover, we discuss highly integrated system designs that combine contrast-enhancing techniques and characterization techniques (like high-resolution spectroscopy) while minimizing the overall complexity. Finally, we explore extreme implementations using all-photonics solutions for ground-based telescopes and dedicated huge apertures for space telescopes.

The Optimal Optical Coronagraph (OOC) Workshop held at the Lorentz Center in September 2017 in Leiden, the Netherlands, gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this second installment of a series of three papers summarizing the outcomes of the OOC workshop (see also~\citenum{ruane2018,snik2018}), we present an overview of common path wavefront sensing/control and Coherent Differential Imaging techniques, highlight the latest results, and expose their relative strengths and weaknesses. We layout critical milestones for the field with the aim of enhancing future ground/space based high contrast imaging platforms. Techniques like these will help to bridge the daunting contrast gap required to image a terrestrial planet in the zone where it can retain liquid water, in reflected light around a G type star from space.

The Optimal Optical Coronagraph (OOC) Workshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this first installment of a series of three papers summarizing the outcomes of the OOC workshop, we present an overview of design methods and optical performance metrics developed for coronagraph instruments. The design and optimization of coronagraphs for future telescopes has progressed rapidly over the past several years in the context of space mission studies for Exo-C, WFIRST, HabEx, and LUVOIR as well as ground-based telescopes. Design tools have been developed at several institutions to optimize a variety of coronagraph mask types. We aim to give a broad overview of the approaches used, examples of their utility, and provide the optimization tools to the community. Though it is clear that the basic function of coronagraphs is to suppress starlight while maintaining light from off-axis sources, our community lacks a general set of standard performance metrics that apply to both detecting and characterizing exoplanets. The attendees of the OOC workshop agreed that it would benefit our community to clearly define quantities for comparing the performance of coronagraph designs and systems. Therefore, we also present a set of metrics that may be applied to theoretical designs, testbeds, and deployed instruments. We show how these quantities may be used to easily relate the basic properties of the optical instrument to the detection significance of the given point source in the presence of realistic noise.

We present observations of the nearby (D$\sim$100\,pc) Herbig star HD~163296 taken with the vortex coronograph at Keck/NIRC2 in the L' band (3.7~$\mu$m), to search for planetary mass companions in the ringed disc surrounding this pre-main sequence star. The images reveal an arc-like region of scattered light from the disc surface layers that is likely associated with the first bright ring detected with ALMA in the $\lambda$=1.3mm dust continuum at $\sim$65~au. We also detect a point-like source at $\sim$0\farcs5 projected separation in the North-East direction, close to the inner edge of the second gap in the millimetre images. Comparing the point source photometry with the atmospheric emission models of non-accreting giant planets, we obtain a mass of 6--7~M$_J$ for a putative protoplanet, assuming a system age of 5~Myr. Based on the contrast at a 95\% level of completeness calculated on the emission-free regions of our images, we set upper limits for the masses of giant planets of 8--15~M$_J$, 4.5--6.5~M$_J$ and 2.5-4.0~M$_J$ at the locations of the first, second and third gap in the millimetre dust continuum, respectively. Further deep, high resolution thermal IR imaging of the HD~163296 system are warranted, to confirm the presence and nature of the point source and to better understand the structure of the dust disc.

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

The HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey searches for dust near the habitable zones (HZs) around nearby, bright main sequence stars. We use nulling interferometry in N band to suppress the bright stellar light and to probe for low levels of HZ dust around the 30 stars observed so far. Our overall detection rate is 18%, including four new detections, among which are the first three around Sun-like stars and the first two around stars without any previously known circumstellar dust. The inferred occurrence rates are comparable for early type and Sun-like stars, but decrease from 60 (+16/-21)% for stars with previously detected cold dust to 8 (+10/-3)% for stars without such excess, confirming earlier results at higher sensitivity. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal excess luminosity function, we put upper limits on the median HZ dust level of 13 zodis (95% confidence) for a sample of stars without cold dust and of 26 zodis when focussing on Sun-like stars without cold dust. However, our data suggest that a more complex luminosity function may be more appropriate. For stars without detectable LBTI excess, our upper limits are almost reduced by a factor of two, demonstrating the strength of LBTI target vetting for future exo-Earth imaging missions. Our statistics are so far limited and extending the survey is critical to inform the design of future exo-Earth imaging surveys.

A warm/hot dust component (at temperature $>$ 300K) has been detected around $\sim$ 20% of stars. This component is called "exozodiacal dust" as it presents similarities with the zodiacal dust detected in our Solar System, even though its physical properties and spatial distribution can be significantly different. Understanding the origin and evolution of this dust is of crucial importance, not only because its presence could hamper future detections of Earth-like planets in their habitable zones, but also because it can provide invaluable information about the inner regions of planetary systems. In this review, we present a detailed overview of the observational techniques used in the detection and characterisation of exozodiacal dust clouds ("exozodis") and the results they have yielded so far, in particular regarding the incidence rate of exozodis as a function of crucial parameters such as stellar type and age, or the presence of an outer cold debris disc. We also present the important constraints that have been obtained, on dust size distribution and spatial location, by using state-of-the-art radiation transfer models on some of these systems. Finally, we investigate the crucial issue of how to explain the presence of exozodiacal dust around so many stars (regardless of their ages) despite the fact that such dust so close to its host star should disappear rapidly due to the coupled effect of collisions and stellar radiation pressure. Several potential mechanisms have been proposed to solve this paradox and are reviewed in detail in this paper. The review finishes by presenting the future of this growing field.

Debris discs have often been described as gas-poor discs as the gas-to-dust ratio is expected to be considerably lower than in primordial,protoplanetary discs. However, recent observations have confirmed the presence of a non-negligible amount of cold gas in the circumstellar (CS) debris discs around young main-sequence stars.This cold gas has been suggested to be related to the outgassing of planetesimals and cometary-like objects. The aim of the paper is to investigate the presence of hot gas in the surroundings of stars bearing cold-gas debris discs. High-resolution optical spectra of all currently known cold-gas-bearing debris-disc systems, with the exception of $\beta$ Pic and Fomalhaut, have been obtained from different observatories.We have analysed the Ca II H & K and the Na I D lines searching for non-photospheric absorptions of CS origin, usually attributed to cometary-like activity. Narrow, stable Ca II and/or Na I absorption features have been detected superimposed to the photospheric lines in 10 out of the 15 observed cold-gas-bearing debris disc.Features are found at the radial velocity of the stars, or slightly blue- or red-shifted, and/or at the velocity of the local interstellar medium (ISM). Some stars also present transient variable events or absorptions extended towards red wavelengths. These are the first detections of such Ca II features in 7 out of the 15 observed stars. In some of these stars, results suggest that the stable and variable absorptions arise from relatively hot gas located in the CS close-in environment. This hot gas is detected in at least ~80%, of edge-on cold-gas-bearing debris discs, while in only ~10% of the discs seen close to face-on. We interpret this as a geometrical effect, and suggest that the non-detection of hot gas absorptions is due to the disc inclination rather than to the absence of the hot-gas component.

Post-processing algorithms play a key role in pushing the detection limits of high-contrast imaging (HCI) instruments. State-of-the-art image processing approaches for HCI enable the production of science-ready images relying on unsupervised learning techniques, such as low-rank approximations, for generating a model PSF and subtracting the residual starlight and speckle noise. In order to maximize the detection rate of HCI instruments and survey campaigns, advanced algorithms with higher sensitivities to faint companions are needed, especially for the speckle-dominated innermost region of the images. We propose a reformulation of the exoplanet detection task (for ADI sequences) that builds on well-established machine learning techniques to take HCI post-processing from an unsupervised to a supervised learning context. In this new framework, we present algorithmic solutions using two different discriminative models: SODIRF (random forests) and SODINN (neural networks). We test these algorithms on real ADI datasets from VLT/NACO and VLT/SPHERE HCI instruments. We then assess their performances by injecting fake companions and using receiver operating characteristic analysis. This is done in comparison with state-of-the-art ADI algorithms, such as ADI-PCA. This study shows the improved sensitivity vs specificity trade-off of the proposed approach. At the diffraction limit, SODINN improves the true positive rate by a factor ranging from ~2 to ~10 (depending on the dataset and angular separation) with respect to ADI-PCA when working at the same false positive level. The proposed supervised detection framework outperforms state-of-the-art techniques in the task of discriminating planet signal from speckles. In addition, it offers the possibility of re-processing existing HCI databases to maximize their scientific return and potentially improve the demographics of directly imaged exoplanets.

Transition disks offer the extraordinary opportunity to look for newly born planets and investigate the early stages of planet formation. In this context we observed the Herbig A5 star MWC 758 with the L band vector vortex coronagraph installed in the near-infrared camera and spectrograph NIRC2 at the Keck II telescope, with the aim of unveiling the nature of the spiral structure by constraining the presence of planetary companions in the system. Our high-contrast imaging observations show a bright (delta L=7.0+/-0.3 mag) point-like emission, south of MWC 758 at a deprojected separation of about 20 au (r=0.111+/- 0. 004 arcsec) from the central star. We also recover the two spiral arms (south-east and north-west), already imaged by previous studies in polarized light, and discover a third one to the south-west of the star. No additional companions were detected in the system down to 5 Jupiter masses beyond 0.6 arcsec from the star. We propose that the bright L band emission could be caused by the presence of an embedded and accreting protoplanet, although the possibility of it being an asymmetric disk feature cannot be excluded. The spiral structure is probably not related to the protoplanet candidate, unless on an inclined and eccentric orbit, and it could be due to one (or more) yet undetected planetary companions at the edge of or outside the spiral pattern. Future observations and additional simulations will be needed to shed light on the true nature of the point-like source and its link with the spiral arms.

As planning for the next generation of high contrast imaging instruments (e.g. WFIRST, HabEx, and LUVOIR, TMT-PFI, EELT-EPICS) matures, and second-generation ground-based extreme adaptive optics facilities (e.g. VLT-SPHERE, Gemini-GPI) are halfway through their principal surveys, it is imperative that the performance of different designs, post-processing algorithms, observing strategies, and survey results be compared in a consistent, statistically robust framework. In this paper, we argue that the current industry standard for such comparisons -- the contrast curve -- falls short of this mandate. We propose a new figure of merit, the "performance map," that incorporates three fundamental concepts in signal detection theory: the true positive fraction (TPF), false positive fraction (FPF), and detection threshold. By supplying a theoretical basis and recipe for generating the performance map, we hope to encourage the widespread adoption of this new metric across subfields in exoplanet imaging.

We report the results of high-angular-resolution observations that search for exozodiacal light in a sample of main sequence stars and sub-giants. Using the "jouvence" of the fiber linked unit for optical recombination (JouFLU) at the center for high angular resolution astronomy (CHARA) telescope array, we have observed a total of 44 stars. Out of the 44 stars, 33 are new stars added to the initial, previously published survey of 42 stars performed at CHARA with the fiber linked unit for optical recombiation (FLUOR). Since the start of the survey extension, we have detected a K-band circumstellar excess for six new stars at the ~ 1\% level or higher, four of which are known or candidate binaries, and two for which the excess could be attributed to exozodiacal dust. We have also performed follow-up observations of 11 of the stars observed in the previously published survey and found generally consistent results. We do however detect a significantly larger excess on three of these follow-up targets: Altair, $\upsilon$ And and $\kappa$ CrB. Interestingly, the last two are known exoplanet host stars. We perform a statistical analysis of the JouFLU and FLUOR samples combined, which yields an overall exozodi detection rate of $21.7^{+5.7}_{-4.1}\%$. We also find that the K-band excess in FGK-type stars correlates with the existence of an outer reservoir of cold ($\lesssim 100\,$K) dust at the $99\%$ confidence level, while the same cannot be said for A-type stars.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical high-contrast imaging. Our package relies on the extensive python stack of scientific libraries and aims to provide a flexible framework for high-contrast data and image processing. In this paper, we describe the capabilities of VIP related to processing image sequences acquired using the angular differential imaging (ADI) observing technique. VIP implements functionalities for building high-contrast data processing pipelines, encompass- ing pre- and post-processing algorithms, potential sources position and flux estimation, and sensitivity curves generation. Among the reference point-spread function subtraction techniques for ADI post-processing, VIP includes several flavors of principal component analysis (PCA) based algorithms, such as annular PCA and incremental PCA algorithm capable of processing big datacubes (of several gigabytes) on a computer with limited memory. Also, we present a novel ADI algorithm based on non-negative matrix factorization (NMF), which comes from the same family of low-rank matrix approximations as PCA and provides fairly similar results. We showcase the ADI capabilities of the VIP library using a deep sequence on HR8799 taken with the LBTI/LMIRCam and its recently commissioned L-band vortex coronagraph. Using VIP we investigated the presence of additional companions around HR8799 and did not find any significant additional point source beyond the four known planets. VIP is available at http://github.com/vortex-exoplanet/VIP and is accompanied with Jupyter notebook tutorials illustrating the main functionalities of the library.

A vortex coronagraph is now available for high contrast observations with the Keck/NIRC2 instrument at L band. Reaching the optimal performance of the coronagraph requires fine control of the wavefront incident on the phase mask. In particular, centering errors can lead to significant stellar light leakage that degrades the contrast performance and prevents the observation of faint planetary companions around the observed stars. It is thus critical to correct for the possible slow drift of the star image from the phase mask center, generally due to mechanical flexures induced by temperature and/or gravity field variation, or to misalignment between the optics that rotate in pupil tracking mode. A control loop based on the QACITS algorithm for the vortex coronagraph has thus been developed and deployed for the Keck/NIRC2 instrument. This algorithm executes the entire observing sequence, including the calibration steps, initial centering of the star on the vortex center and stabilisation during the acquisition of science frames. On-sky data show that the QACITS control loop stabilizes the position of the star image down to 2.4 mas rms at a frequency of about 0.02 Hz. However, the accuracy of the estimator is probably limited by a systematic error due to a misalignment of the Lyot stop with respect to the entrance pupil, estimated to be on the order of 4.5 mas. A method to reduce the amplitude of this bias down to 1 mas is proposed. The QACITS control loop has been successfully implemented and provides a robust method to center and stabilize the star image on the vortex mask. In addition, QACITS ensures a repeatable pointing quality and significantly improves the observing efficiency compared to manual operations. It is now routinely used for vortex coronagraph observations at Keck/NIRC2, providing contrast and angular resolution capabilities suited for exoplanet and disk imaging.

We present the first scattered-light images of the debris disk around 49 ceti, a ~40 Myr A1 main sequence star at 59 pc, famous for hosting two massive dust belts as well as large quantities of atomic and molecular gas. The outer disk is revealed in reprocessed archival Hubble Space Telescope NICMOS F110W images, as well as new coronagraphic H band images from the Very Large Telescope SPHERE instrument. The disk extends from 1.1" (65 AU) to 4.6" (250 AU), and is seen at an inclination of 73degr, which refines previous measurements at lower angular resolution. We also report no companion detection larger than 3 M_Jup at projected separations beyond 20 AU from the star (0.34"). Comparison between the F110W and H-band images is consistent with a grey color of 49 ceti's dust, indicating grains larger than >2microns. Our photometric measurements indicate a scattering efficiency / infrared excess ratio of 0.2-0.4, relatively low compared to other characterized debris disks. We find that 49 ceti presents morphological and scattering properties very similar to the gas-rich HD 131835 system. From our constraint on the disk inclination we find that the atomic gas previously detected in absorption must extend to the inner disk, and that the latter must be depleted of CO gas. Building on previous studies, we propose a schematic view of the system describing the dust and gas structure around 49 ceti and hypothetic scenarios for the gas nature and origin.

HD 141569 A is a pre-main sequence B9.5 Ve star surrounded by a prominent and complex circumstellar disk, likely still in a transition stage from protoplanetary to debris disk phase. Here, we present a new image of the third inner disk component of HD 141569 A made in the L' band (3.8 micron) during the commissioning of the vector vortex coronagraph recently installed in the near-infrared imager and spectrograph NIRC2 behind the W.M. Keck Observatory Keck II adaptive optics system. We used reference point spread function subtraction, which reveals the innermost disk component from the inner working distance of $\simeq 23$ AU and up to $\simeq 70$ AU. The spatial scale of our detection roughly corresponds to the optical and near-infrared scattered light, thermal Q, N and 8.6 micron PAH emission reported earlier. We also see an outward progression in dust location from the L'-band to the H-band (VLT/SPHERE image) to the visible (HST/STIS image), likely indicative of dust blowout. The warm disk component is nested deep inside the two outer belts imaged by HST NICMOS in 1999 (respectively at 406 and 245 AU). We fit our new L'-band image and spectral energy distribution of HD 141569 A with the radiative transfer code MCFOST. Our best-fit models favor pure olivine grains, and are consistent with the composition of the outer belts. While our image shows a putative very-faint point-like clump or source embedded in the inner disk, we did not detect any true companion within the gap between the inner disk and the first outer ring, at a sensitivity of a few Jupiter masses.

An optical vortex coronagraph has been implemented within the NIRC2 camera on the Keck II telescope and used to carry out on-sky tests and observations. The development of this new L'-band observational mode is described, and an initial demonstration of the new capability is presented: a resolved image of the low-mass companion to HIP79124, which had previously been detected by means of interferometry. With HIP79124 B at a projected separation of 186.5 mas, both the small inner working angle of the vortex coronagraph and the related imaging improvements were crucial in imaging this close companion directly. Due to higher Strehl ratios and more relaxed contrasts in L' band versus H band, this new coronagraphic capability will enable high-contrast small-angle observations of nearby young exoplanets and disks on a par with those of shorter-wavelength extreme adaptive optics coronagraphs.

Uncovering the ingredients and the architecture of planetary systems is a very active field of research that has fuelled many new theories on giant planet formation, migration, composition, and interaction with the circumstellar environment. We aim at discovering and studying new such systems, to further expand our knowledge of how low-mass companions form and evolve. We obtained high-contrast H-band images of the circumstellar environment of the F5V star HD206893, known to host a debris disc never detected in scattered light. These observations are part of the SPHERE High Angular Resolution Debris Disc Survey (SHARDDS) using the InfraRed Dual-band Imager and Spectrograph (IRDIS) installed on VLT/SPHERE. We report the detection of a source with a contrast of 3.6x10^{-5} in the H-band, orbiting at a projected separation of 270 milliarcsecond or 10 au, corresponding to a mass in the range 24 to 73 Mjup for an age of the system in the range 0.2 to 2 Gyr. The detection was confirmed ten months later with VLT/NaCo, ruling out a background object with no proper motion. A faint extended emission compatible with the disc scattered light signal is also observed. The detection of a low-mass companion inside a massive debris disc makes this system an analog of other young planetary systems such as beta Pictoris, HR8799 or HD95086 and requires now further characterisation of both components to understand their interactions.

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

We present the first resolved image of the debris disk around the 16+/-8 Myr old star, HD 114082. The observation was made in the H-band using the SPHERE instrument. The star is at a distance of 92+/-6 pc in the Lower Centaurus Crux association. Using a Markov Chain Monte Carlo analysis, we determined that the debris is likely in the form of a dust ring with an inner edge of 27.7+2.8/-3.5 au, position angle -74+0.5/-1.5 deg, and an inclination with respect to the line of sight of 6.7+3.8/-0.4 deg. The disk imaged in scattered light has a surface density declining with radius like ~r^(-4), steeper than expected for grain blowout by radiation pressure. We find only marginal evidence (2 sigma) of eccentricity, and rule out planets more massive than 1.0 Mjup orbiting within 1 au of the ring's inner edge, since such a planet would have disrupted the disk. The disk has roughly the same fractional disk luminosity (Ldisk/L*=3.3x10^(-3)) as HR4796A and Beta Pictoris, however it was not detected by previous instrument facilities most likely because of its small angular size (radius~0.4"), low albedo (~0.2) and low scattering efficiency far from the star due to high scattering anisotropy. With the arrival of extreme adaptive optics systems like SPHERE and GPI, the morphology of smaller, fainter and more distant debris disks are being revealed, providing clues to planet-disk interactions in young protoplanetary systems.

Context. The Annular Groove Phase Mask (AGPM) is one possible implementation of the vector vortex coronagraph, where the helical phase ramp is produced by a concentric subwavelength grating. For several years, we have been manufacturing AGPMs by etching gratings into synthetic diamond substrates using inductively coupled plasma etching. Aims. We aim to design, fabricate, optimize, and evaluate new L-band AGPMs that reach the highest possible coronagraphic performance, for applications in current and forthcoming infrared high-contrast imagers. Methods. Rigorous coupled wave analysis (RCWA) is used for designing the subwavelength grating of the phase mask. Coronagraphic performance evaluation is performed on a dedicated optical test bench. The experimental results of the performance evaluation are then used to accurately determine the actual profile of the fabricated gratings, based on RCWA modeling. Results. The AGPM coronagraphic performance is very sensitive to small errors in etch depth and grating profile. Most of the fabricated components therefore show moderate performance in terms of starlight rejection (a few 100:1 in the best cases). Here we present new processes for re-etching the fabricated components in order to optimize the parameters of the grating and hence significantly increase their coronagraphic performance. Starlight rejection up to 1000:1 is demonstrated in a broadband L filter on the coronagraphic test bench, which corresponds to a raw contrast of about 1e-5 at two resolution elements from the star for a perfect input wave front on a circular, unobstructed aperture. Conclusions. Thanks to their exquisite performance, our latest L-band AGPMs are good candidates for installation in state-of-the-art and future high-contrast thermal infrared imagers, such as METIS for the E-ELT.

HR8799 is orbited by at least four giant planets, making it a prime target for the recently commissioned Spectro-Polarimetric High-contrast Exoplanet REsearch (VLT/SPHERE). As such, it was observed on five consecutive nights during the SPHERE science verification in December 2014. We aim to take full advantage of the SPHERE capabilities to derive accurate astrometric measurements based on H-band images acquired with the Infra-Red Dual-band Imaging and Spectroscopy (IRDIS) subsystem, and to explore the ultimate astrometric performance of SPHERE in this observing mode. We also aim to present a detailed analysis of the orbital parameters for the four planets. We report the astrometric positions for epoch 2014.93 with an accuracy down to 2.0 mas, mainly limited by the astrometric calibration of IRDIS. For each planet, we derive the posterior probability density functions for the six Keplerian elements and identify sets of highly probable orbits. For planet d, there is clear evidence for nonzero eccentricity ($e \simeq 0.35$), without completely excluding solutions with smaller eccentricities. The three other planets are consistent with circular orbits, although their probability distributions spread beyond $e = 0.2$, and show a peak at $e \simeq 0.1$ for planet e. The four planets have consistent inclinations of about $30\deg$ with respect to the sky plane, but the confidence intervals for the longitude of ascending node are disjoint for planets b and c, and we find tentative evidence for non-coplanarity between planets b and c at the $2 \sigma$ level.

We present an intensive monitoring of high-resolution spectra of the Ca {\sc ii} K line in the A7IV shell star $\Phi$ Leo at very short (minutes, hours), short (night to night), and medium (weeks, months) timescales. The spectra show remarkable variable absorptions on timescales of hours, days, and months. The characteristics of these sporadic events are very similar to most that are observed toward the debris disk host star $\beta$ Pic, which are commonly interpreted as signs of the evaporation of solid, comet-like bodies grazing or falling onto the star. Therefore, our results suggest the presence of solid bodies around $\Phi$ Leo. To our knowledge, with the exception of $\beta$ Pic, our monitoring has the best time resolution at the mentioned timescales for a star with events attributed to exocomets. Assuming the cometary scenario and considering the timescales of our monitoring, our results indicate that $\Phi$ Leo presents the richest environment with comet-like events known to date, second only to $\beta$ Pic.

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

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

Small inner working angle coronagraphs are essential to benefit from the full potential of large and future extremely large ground-based telescopes, especially in the context of the detection and characterization of exoplanets. Among existing solutions, the vortex coronagraph stands as one of the most effective and promising solutions. However, for focal-plane coronagraph, a small inner working angle comes necessarily at the cost of a high sensitivity to pointing errors. This is the reason why a pointing control system is imperative to stabilize the star on the vortex center against pointing drifts due to mechanical flexures, that generally occur during observation due for instance to temperature and/or gravity variations. We have therefore developed a technique called QACITS (Quadrant Analysis of Coronagraphic Images for Tip-tilt Sensing), which is based on the analysis of the coronagraphic image shape to infer the amount of pointing error. It has been shown that the flux gradient in the image is directly related to the amount of tip-tilt affecting the beam. The main advantage of this technique is that it does not require any additional setup and can thus be easily implemented on all current facilities equipped with a vortex phase mask. In this paper, we focus on the implementation of the QACITS sensor at Keck/NIRC2, where an L-band AGPM has been recently commissioned (June and October 2015), successfully validating the QACITS estimator in the case of a centrally obstructed pupil. The algorithm has been designed to be easily handled by any user observing in vortex mode, which is available for science in shared risk mode since 2016B.

For several years, we have been developing vortex phase masks based on sub-wavelength gratings, known as Annular Groove Phase Masks. Etched onto diamond substrates, these AGPMs are currently designed to be used in the thermal infrared (ranging from 3 to 13 {\mu}m). Our AGPMs were first installed on VLT/NACO and VLT/VISIR in 2012, followed by LBT/LMIRCam in 2013 and Keck/NIRC2 in 2015. In this paper, we review the development, commissioning, on-sky performance, and early scientific results of these new coronagraphic modes and report on the lessons learned. We conclude with perspectives for future developments and applications.

The Annular Groove Phase Mask (AGPM) is a vectorial vortex phase mask. It acts as a half-wave plate with a radial fast axis orientation operating in the mid infrared domain. When placed at the focus of a telescope element provides a continuous helical phase ramp for an on axis sources, which creates the orbital angular momentum. Thanks to that phase, the intensity of the central source is canceled by a down-stream pupil stop, while the off axis sources are not affected. However due to experimental conditions the nulling is hardly perfect. To improve the null, a Mach-Zehnder interferometer containing Dove prisms differently oriented can be proposed to sort out light based on its orbital angular momentum (OAM). Thanks to the differential rotation of the beam, a {\pi} phase shift is achieved for the on axis light affected by a non zero OAM. Therefore the contrast between the star and its faint companion is enhanced. Nevertheless, due the Dove prisms birefringence, the performance of the interferometer is relatively poor. To solve this problem, we propose to add a birefringent wave-plate in each arm to compensate this birefringence. In this paper, we will develop the mathematical model of the wave front using the Jones formalism. The performance of the interferometer is at first computed for the simple version without the birefringent plate. Then the effect of the birefringent plate is be mathematically described and the performance is re-computed.

In this paper, we present the infrared coronagraphic test bench of the University of Li\`ege named VODCA (Vortex Optical Demonstrator for Coronagraphic Applications). The goal of the bench is to assess the performances of the Annular Groove Phase Masks (AGPMs) at near- to mid-infrared wavelengths. The AGPM is a subwavelength grating vortex coronagraph of charge two (SGVC2) made out of diamond. The bench is designed to be completely achromatic and will be composed of a super continuum laser source emitting in the near to mid-infrared, several parabolas, diaphragms and an infrared camera. This way, we will be able to test the different AGPMs in the M, L, K and H bands. Eventually, the bench will also allow the computation of the incident wavefront aberrations on the coronagraph. A reflective Lyot stop will send most of the stellar light to a second camera to perform low-order wavefront sensing. This second system coupled with a deformable mirror will allow the correction of the wavefront aberrations. We also aim to test other pre- and/or post-coronagraphic concepts such as optimal apodization.

The mid-infrared region is well suited for exoplanet detection thanks to the reduced contrast between the planet and its host star with respect to the visible and near-infrared wavelength regimes. This contrast may be further improved with Vector Vortex Coronagraphs (VVCs), which allow us to cancel the starlight. One flavour of the VVC is the AGPM (Annular Groove Phase Mask), which adds the interesting properties of subwavelength gratings (achromaticity, robustness) to the already known properties of the VVC. In this paper, we present the optimized designs, as well as the expected performances of mid-IR AGPMs etched onto synthetic diamond substrates, which are considered for the E-ELT/METIS instrument.

Debris discs are typically revealed through excess emission at infrared wavelengths. Most discs exhibit excess at mid- and far-infrared wavelengths, analogous to the solar system's Asteroid and Edgeworth-Kuiper belts. Recently, stars with strong (1 per cent) excess at near-infrared wavelengths were identified through interferometric measurements. Using the HIgh Precision Polarimetric Instrument (HIPPI), we examined a sub-sample of these hot dust stars (and appropriate controls) at parts-per-million sensitivity in SDSS g' (green) and r' (red) filters for evidence of scattered light. No detection of strongly polarized emission from the hot dust stars is seen. We therefore rule out scattered light from a normal debris disk as the origin of this emission. A wavelength dependent contribution from multiple dust components for hot dust stars is inferred from the dispersion (difference in polarization angle in red and green) of southern stars. Contributions of 17 ppm (green) and 30 ppm (red) are calculated, with strict 3 sigma upper limits of 76 and 68 ppm, respectively. This suggests weak hot dust excesses consistent with thermal emission, although we cannot rule out contrived scenarios, e.g. dust in a spherical shell or face on discs. We also report on the nature of the local interstellar medium, obtained as a byproduct of the control measurements. Highlights include the first measurements of the polarimetric colour of the local interstellar medium and discovery of a southern sky region with a polarization per distance thrice the previous maximum. The data suggest the wavelength of maximum polarization is bluer than typical.

Data processing constitutes a critical component of high-contrast exoplanet imaging. Its role is almost as important as the choice of a coronagraph or a wavefront control system, and it is intertwined with the chosen observing strategy. Among the data processing techniques for angular differential imaging (ADI), the most recent is the family of principal component analysis (PCA) based algorithms. PCA serves, in this case, as a subspace projection technique for constructing a reference point spread function (PSF) that can be subtracted from the science data for boosting the detectability of potential companions present in the data. Unfortunately, when building this reference PSF from the science data itself, PCA comes with certain limitations such as the sensitivity of the lower dimensional orthogonal subspace to non-Gaussian noise. Inspired by recent advances in machine learning algorithms such as robust PCA, we aim to propose a localized subspace projection technique that surpasses current PCA-based post-processing algorithms in terms of the detectability of companions at near real-time speed, a quality that will be useful for future direct imaging surveys. We used randomized low-rank approximation methods recently proposed in the machine learning literature, coupled with entry-wise thresholding to decompose an ADI image sequence locally into low-rank, sparse, and Gaussian noise components (LLSG). This local three-term decomposition separates the starlight and the associated speckle noise from the planetary signal, which mostly remains in the sparse term. We tested the performance of our new algorithm on a long ADI sequence obtained on beta Pictoris with VLT/NACO. Compared to a standard PCA approach, LLSG decomposition reaches a higher signal-to-noise ratio and has an overall better performance in the receiver operating characteristic space. (abridged).

The Large Binocular Telescope Interferometer (LBTI) is a versatile instrument designed for high-angular resolution and high-contrast infrared imaging (1.5-13 microns). In this paper, we focus on the mid-infrared (8-13 microns) nulling mode and present its theory of operation, data reduction, and on-sky performance as of the end of the commissioning phase in March 2015. With an interferometric baseline of 14.4 meters, the LBTI nuller is specifically tuned to resolve the habitable zone of nearby main-sequence stars, where warm exozodiacal dust emission peaks. Measuring the exozodi luminosity function of nearby main-sequence stars is a key milestone to prepare for future exoEarth direct imaging instruments. Thanks to recent progress in wavefront control and phase stabilization, as well as in data reduction techniques, the LBTI demonstrated in February 2015 a calibrated null accuracy of 0.05% over a three-hour long observing sequence on the bright nearby A3V star beta Leo. This is equivalent to an exozodiacal disk density of 15 to 30 zodi for a Sun-like star located at 10pc, depending on the adopted disk model. This result sets a new record for high-contrast mid-infrared interferometric imaging and opens a new window on the study of planetary systems.

Small inner working angle coronagraphs, like the vortex phase mask, are essential to exploit the full potential of ground-based telescopes in the context of exoplanet detection and characterization. However, the drawback of this attractive feature is a high sensitivity to pointing errors, which degrades the performance of the coronagraph. We propose a tip-tilt retrieval technique based on the analysis of the final coronagraphic image, hereafter called Quadrant Analysis of Coronagraphic Images for Tip-tilt Sensing (QACITS). Under the assumption of small phase aberrations, we show that the behaviour of the vortex phase mask can be simply described from the entrance pupil to the Lyot stop plane by Zernike polynomials. This convenient formalism is used to establish the theoretical basis of the QACITS technique. Simulations have been performed to demonstrate the validity and limits of the technique, including the case of a centrally obstructed pupil. The QACITS technique principle is further validated by experimental results in the case of an unobstructed circular aperture. The typical configuration of the Keck telescope (24% central obstruction) has been simulated with additional high order aberrations. In these conditions, our simulations show that the QACITS technique is still adapted to centrally obstructed pupils and performs tip-tilt retrieval with a precision of $5 \times 10^{-2}$ {\lambda}/D when wavefront errors amount to {\lambda}/14 rms and $10^{-2}$ {\lambda}/D for {\lambda}/70 rms errors (with {\lambda} the wavelength and D the pupil diameter). The implementation of the QACITS technique is based on the analysis of the scientific image and does not require any modification of the original setup. Current facilities equipped with a vortex phase mask can thus directly benefit from this technique to improve the contrast performance close to the axis.

We present methods for optimizing pupil and focal plane optical elements that improve the performance of vortex coronagraphs on telescopes with obstructed or segmented apertures. Phase-only and complex masks are designed for the entrance pupil, focal plane, and the plane of the Lyot stop. Optimal masks are obtained using both analytical and numerical methods. The latter makes use of an iterative error reduction algorithm to calculate "correcting" optics that mitigate unwanted diffraction from aperture obstructions. We analyze the achieved performance in terms of starlight suppression, contrast, off-axis image quality, and chromatic dependence. Manufacturing considerations and sensitivity to aberrations are also discussed. This work provides a path to joint optimization of multiple coronagraph planes to maximize sensitivity to exoplanets and other faint companions.

The vortex coronagraph is an optical instrument that precisely removes on-axis starlight allowing for high contrast imaging at small angular separation from the star, thereby providing a crucial capability for direct detection and characterization of exoplanets and circumstellar disks. Telescopes with aperture obstructions, such as secondary mirrors and spider support structures, require advanced coronagraph designs to provide adequate starlight suppression. We introduce a phase-only Lyot-plane optic to the vortex coronagraph that offers improved contrast performance on telescopes with complicated apertures. Potential solutions for the European Extremely Large Telescope (E-ELT) are described and compared. Adding a Lyot-plane phase mask relocates residual starlight away from a region of the image plane thereby reducing stellar noise and improving sensitivity to off-axis companions. The phase mask is calculated using an iterative phase retrieval algorithm. Numerically, we achieve a contrast on the order of $10^{-6}$ for a companion with angular displacement as small as $4~\lambda/D$ with an E-ELT type aperture. Even in the presence of aberrations, improved performance is expected compared to either a conventional vortex coronagraph or optimized pupil plane phase element alone.

The Large Binocular Telescope Interferometer (LBTI) is a strategic instrument of the LBT designed for high-sensitivity, high-contrast, and high-resolution infrared (1.5-13 $\mu$m) imaging of nearby planetary systems. To carry out a wide range of high-spatial resolution observations, it can combine the two AO-corrected 8.4-m apertures of the LBT in various ways including direct (non-interferometric) imaging, coronagraphy (APP and AGPM), Fizeau imaging, non-redundant aperture masking, and nulling interferometry. It also has broadband, narrowband, and spectrally dispersed capabilities. In this paper, we review the performance of these modes in terms of exoplanet science capabilities and describe recent instrumental milestones such as first-light Fizeau images (with the angular resolution of an equivalent 22.8-m telescope) and deep interferometric nulling observations.

Context. The direct detection of exoplanets with high-contrast imaging requires advanced data processing methods to disentangle potential planetary signals from bright quasi-static speckles. Among them, angular differential imaging (ADI) permits potential planetary signals with a known rotation rate to be separated from instrumental speckles that are either statics or slowly variable. The method presented in this paper, called ANDROMEDA for ANgular Differential OptiMal Exoplanet Detection Algorithm is based on a maximum likelihood approach to ADI and is used to estimate the position and the flux of any point source present in the field of view. Aims. In order to optimize and experimentally validate this previously proposed method, we applied ANDROMEDA to real VLT/NaCo data. In addition to its pure detection capability, we investigated the possibility of defining simple and efficient criteria for automatic point source extraction able to support the processing of large surveys. Methods. To assess the performance of the method, we applied ANDROMEDA on VLT/NaCo data of TYC-8979-1683-1 which is surrounded by numerous bright stars and on which we added synthetic planets of known position and flux in the field. In order to accommodate the real data properties, it was necessary to develop additional pre-processing and post-processing steps to the initially proposed algorithm. We then investigated its skill in the challenging case of a well-known target, $\beta$ Pictoris, whose companion is close to the detection limit and we compared our results to those obtained by another method based on principal component analysis (PCA). Results. Application on VLT/NaCo data demonstrates the ability of ANDROMEDA to automatically detect and characterize point sources present in the image field. We end up with a robust method bringing consistent results with a sensitivity similar to the recently published algorithms, with only two parameters to be fine tuned. Moreover, the companion flux estimates are not biased by the algorithm parameters and do not require a posteriori corrections. Conclusions. ANDROMEDA is an attractive alternative to current standard image processing methods that can be readily applied to on-sky data.

We present an analytic model to estimate the capabilities of space missions dedicated to the search for biosignatures in the atmosphere of rocky planets located in the habitable zone of nearby stars. Relations between performance and mission parameters such as mirror diameter, distance to targets, and radius of planets, are obtained. Two types of instruments are considered: coronagraphs observing in the visible, and nulling interferometers in the thermal infrared. Missions considered are: single-pupil coronagraphs with a 2.4 m primary mirror, and formation flying interferometers with 4 x 0.75 m collecting mirrors. The numbers of accessible planets are calculated as a function of {\eta}earth. When Kepler gives its final estimation for {\eta}earth, the model will permit a precise assessment of the potential of each instrument. Based on current estimations, {\eta}earth = 10% around FGK stars and 50% around M stars, the coronagraph could study in spectroscopy only ~1.5 relevant planets, and the interferometer ~14.0. These numbers are obtained under the major hypothesis that the exozodiacal light around the target stars is low enough for each instrument. In both cases, a prior detection of planets is assumed and a target list established. For the long-term future, building both types of spectroscopic instruments, and using them on the same targets, will be the optimal solution because they provide complementary information. But as a first affordable space mission, the interferometer looks the more promising in term of biosignature harvest.

Exozodiacal dust is warm or hot dust found in the inner regions of planetary systems orbiting main sequence stars, in or around their habitable zones. The dust can be the most luminous component of extrasolar planetary systems, but predominantly emits in the near- to mid-infrared where it is outshone by the host star. Interferometry provides a unique method of separating this dusty emission from the stellar emission. The visitor instrument PIONIER at the Very Large Telescope Interferometer (VLTI) has been used to search for hot exozodiacal dust around a large sample of nearby main sequence stars. The results of this survey are summarised: 9 out of 85 stars show excess exozodiacal emission over the stellar photospheric emission.

We report on the first nulling interferometric observations with the Large Binocular Telescope Interferometer (LBTI), resolving the N' band (9.81 - 12.41 um) emission around the nearby main-sequence star eta Crv (F2V, 1-2 Gyr). The measured source null depth amounts to 4.40% +/- 0.35% over a field-of-view of 140 mas in radius (~2.6\,AU at the distance of eta Corvi) and shows no significant variation over 35{\deg} of sky rotation. This relatively low null is unexpected given the total disk to star flux ratio measured by Spitzer/IRS (~23% across the N' band), suggesting that a significant fraction of the dust lies within the central nulled response of the LBTI (79 mas or 1.4 AU). Modeling of the warm disk shows that it cannot resemble a scaled version of the Solar zodiacal cloud, unless it is almost perpendicular to the outer disk imaged by Herschel. It is more likely that the inner and outer disks are coplanar and the warm dust is located at a distance of 0.5-1.0 AU, significantly closer than previously predicted by models of the IRS spectrum (~3 AU). The predicted disk sizes can be reconciled if the warm disk is not centrosymmetric, or if the dust particles are dominated by very small grains. Both possibilities hint that a recent collision has produced much of the dust. Finally, we discuss the implications for the presence of dust at the distance where the insolation is the same as Earth's (2.3 AU).

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

One possible solution to achieve high contrast direct imaging at a small inner working angle (IWA) is to use a vector vortex coronagraph (VVC), which provides a continuous helical phase ramp in the focal plane of the telescope with a phase singularity in its center. Such an optical vortex is characterized by its topological charge, i.e., the number of times the phase accumulates 2{\pi} radians along a closed path surrounding the singularity. Over the past few years, we have been developing a charge-2 VVC induced by rotationally symmetric subwavelength gratings (SGVC2), also known as the Annular Groove Phase Mask (AGPM). Since 2013, several SGVC2s (or AGPMs) were manufactured using synthetic diamond substrate, then validated on dedicated optical benches, and installed on 10-m class telescopes. Increasing the topological charge seems however mandatory for cancelling the light of bright stars which will be partially resolved by future Extremely Large Telescopes in the near-infrared. In this paper, we first detail our motivations for developing an SGVC4 (charge 4) dedicated to the near-infrared domain. The challenge lies in the design of the pattern which is unrealistic in the theoretically perfect case, due to state-of-the-art manufacturing limitations. Hence, we propose a new realistic design of SGVC4 with minimized discontinuities and optimized phase ramp, showing conclusive improvements over previous works in this field. A preliminary validation of our concept is given based on RCWA simulations, while full 3D finite-difference time-domain simulations (and eventually laboratory tests) will be required for a final validation.

In recent years, phase mask coronagraphy has become increasingly efficient in imaging the close environment of stars, enabling the search for exoplanets and circumstellar disks. Coronagraphs are ideally suited instruments, characterized by high dynamic range imaging capabilities, while preserving a small inner working angle. The AGPM (Annular Groove Phase Mask, Mawet et al. 2005) consists of a vector vortex induced by a rotationally symmetric subwavelength grating. This technique constitutes an almost unique solution to the achromatization at longer wavelengths (mid-infrared). For this reason, we have specially conceived a mid-infrared AGPM coronagraph for the forthcoming upgrade of VISIR, the mid-IR imager and spectrograph on the VLT at ESO (Paranal), in collaboration with members of the VISIR consortium. The implementation phase of the VISIR Upgrade Project is foreseen for May-August 2012, and the AGPM installed will cover the 11-13.2 {\mu}m spectral range. In this paper, we present the entire fabrication process of our AGPM imprinted on a diamond substrate. Diamond is an ideal material for mid-infrared wavelengths owing to its high transparency, small dispersion, extremely low thermal expansion and outstanding mechanical and chemical properties. The design process has been performed with an algorithm based on the rigorous coupled wave analysis (RCWA), and the micro-fabrication has been carried out using nano-imprint lithography and reactive ion etching. A precise grating profile metrology has also been conducted using cleaving techniques. Finally, we show the deposit of fiducials (i.e. centering marks) with Aerosol Jet Printing (AJP). We conclude with the ultimate coronagraph expected performances.

We present the first observations obtained with the L'-band AGPM vortex coronagraph recently installed on LBTI/LMIRCam. The AGPM (Annular Groove Phase Mask) is a vector vortex coronagraph made from diamond subwavelength gratings. It is designed to improve the sensitivity and dynamic range of high-resolution imaging at very small inner working angles, down to 0.09 arcseconds in the case of LBTI/LMIRCam in the L' band. During the first hours on sky, we observed the young A5V star HR\,8799 with the goal to demonstrate the AGPM performance and assess its relevance for the ongoing LBTI planet survey (LEECH). Preliminary analyses of the data reveal the four known planets clearly at high SNR and provide unprecedented sensitivity limits in the inner planetary system (down to the diffraction limit of 0.09 arcseconds).

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

(abridged) Vortex coronagraphs are among the most promising solutions to perform high contrast imaging at small angular separations. They feature a very small inner working angle, a clear 360 degree discovery space, have demonstrated very high contrast capabilities, are easy to implement on high-contrast imaging instruments, and have already been extensively tested on the sky. Since 2005, we have been designing, developing and testing an implementation of the charge-2 vector vortex phase mask based on concentric subwavelength gratings, referred to as the Annular Groove Phase Mask (AGPM). Science-grade mid-infrared AGPMs were produced in 2012 for the first time, using plasma etching on synthetic diamond substrates. They have been validated on a coronagraphic test bench, showing broadband peak rejection up to 500:1 in the L band, which translates into a raw contrast of about $6\times 10^{-5}$ at $2 \lambda/D$. Three of them have now been installed on world-leading diffraction-limited infrared cameras (VLT/NACO, VLT/VISIR and LBT/LMIRCam). During the science verification observations with our L-band AGPM on NACO, we observed the beta Pictoris system and obtained unprecedented sensitivity limits to planetary companions down to the diffraction limit ($0.1''$). More recently, we obtained new images of the HR 8799 system at L band during the AGPM first light on LMIRCam. After reviewing these first results obtained with mid-infrared AGPMs, we will discuss the short- and mid-term goals of the on-going VORTEX project, which aims to improve the performance of our vortex phase masks for future applications on second-generation high-contrast imagers and on future extremely large telescopes (ELTs).

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

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

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

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

In this paper, we review the impact of small sample statistics on detection thresholds and corresponding confidence levels (CLs) in high contrast imaging at small angles. When looking close to the star, the number of resolution elements decreases rapidly towards small angles. This reduction of the number of degrees of freedom dramatically affects CLs and false alarm probabilities. Naively using the same ideal hypothesis and methods as for larger separations, which are well understood and commonly assume Gaussian noise, can yield up to one order of magnitude error in contrast estimations at fixed CL. The statistical penalty exponentially increases towards very small inner working angles. Even at 5-10 resolution elements from the star, false alarm probabilities can be significantly higher than expected. Here we present a rigorous statistical analysis which ensures robustness of the CL, but also imposes a substantial limitation on corresponding achievable detection limits (thus contrast) at small angles. This unavoidable fundamental statistical effect has a significant impact on current coronagraphic and future high contrast imagers. Finally, the paper concludes with practical recommendations to account for small number statistics when computing the sensitivity to companions at small angles and when exploiting the results of direct imaging planet surveys.

Very few debris discs have been imaged in scattered light at wavelengths beyond 3 microns because the thermal emission from both the sky and the telescope is generally too strong with respect to the faint emission of a debris disc. We present here the first analysis of a high angular resolution image of the disc of beta Pictoris at 3.8 microns. Our primary objective is to probe the innermost parts of the beta Pictoris debris disc and describe its morphology. We performed extensive forward modelling to correct for the biases induced by angular differential imaging on extended objects and derive the physical parameters of the disc. This work relies on a new analysis of seven archival datasets of beta Pictoris observed with VLT/NaCo in the Lp band, associated with disc forward modelling to correct for the biases induced by that technique. The disc is detected above a 5 sigma level between 0.4" and 3.8". The two extensions have a similar brightness within error bars. We confirm an asymmetry previously observed at larger distances from the star and at shorter wavelengths: the isophotes are more widely spaced on the north-west side (above the disc apparent midplane) than on the south-east side. This is interpreted as a small inclination of the disc combined with anisotropic scattering. Our best-fit model has an inclination of 86 degrees with an anisotropic Henyey- Greenstein coefficient of 0.36. This interpretation is supported by a new asymmetry detected in the disc: the disc is significantly bowed towards the north-west within 3" (above the apparent midplane). We also detect a possible new asymmetry within 1", but at this stage we cannot discern between a real feature and an underlying speckle.

Context. The orbit of the giant planet discovered around beta Pic is slightly inclined with respect to the outer parts of the debris disc, which creates a warp in the inner debris disc. This inclination might be explained by gravitational interactions with other planets. Aims. We aim to search for additional giant planets located at smaller angular separations from the star. Methods. We used the new L'-band AGPM coronagraph on VLT/NACO, which provides an exquisite inner working angle. A long observing sequence was obtained on beta Pic in pupil-tracking mode. To derive sensitivity limits, the collected images were processed using a principal-component analysis technique specifically tailored to angular differential imaging. Results. No additional planet is detected down to an angular separation of 0.2" with a sensitivity better than 5 MJup. Meaningful upper limits (< 10 MJup) are derived down to an angular separation of 0.1", which corresponds to 2 AU at the distance of beta Pic.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

(Abridged) Dust is expected to be ubiquitous in extrasolar planetary systems owing to the dynamical activity of minor bodies. Inner dust populations are, however, still poorly known because of the high contrast and small angular separation with respect to their host star. We aim to determine the level of near-infrared exozodiacal dust emission around a sample of 42 nearby main sequence stars with spectral types ranging from A to K and to investigate its correlation with various stellar parameters and with the presence of cold dust belts. We use high-precision K-band visibilities obtained with the FLUOR interferometer on the shortest baseline of the CHARA array. The calibrated visibilities are compared with the expected visibility of the stellar photosphere to assess whether there is an additional, fully resolved circumstellar emission. Near-infrared circumstellar emission amounting to about 1% of the stellar flux is detected around 13 of our 42 target stars. Follow-up observations showed that one of them (eps Cep) is associated with a stellar companion, while another one was detected around what turned out to be a giant star (kap CrB). The remaining 11 excesses found around single main sequence stars are most probably associated with hot circumstellar dust, yielding an overall occurrence rate of 28+8-6% for our (biased) sample. We show that the occurrence rate of bright exozodiacal discs correlates with spectral type, K-band excesses being more frequent around A-type stars. It also correlates with the presence of detectable far-infrared excess emission in the case of solar-type stars. This study provides new insight into the phenomenon of bright exozodiacal discs, showing that hot dust populations are probably linked to outer dust reservoirs in the case of solar-type stars. For A-type stars, no clear conclusion can be made regarding the origin of the detected near-infrared excesses.

(abridged) Infrared excesses associated with debris disk host stars detected so far peak at wavelengths around ~100{\mu}m or shorter. However, six out of 31 excess sources in the Herschel OTKP DUNES have been seen to show significant - and in some cases extended - excess emission at 160{\mu}m, which is larger than the 100{\mu}m excess. This excess emission has been suggested to stem from debris disks colder than those known previously. Using several methods, we re-consider whether some or even all of the candidates may be associated with unrelated galactic or extragalactic emission and conclude that it is highly unlikely that none of the candidates represents a true circumstellar disk. For true disks, both the dust temperatures inferred from the SEDs and the disk radii estimated from the images suggest that the dust is nearly as cold as a blackbody. This requires the grains to be larger than ~100{\mu}m, regardless of their material composition. To explain the dearth of small grains, we explore several conceivable scenarios: transport-dominated disks, disks of low dynamical excitation, and disks of unstirred primordial macroscopic grains. Our qualitative analysis and collisional simulations rule out the first two of these scenarios, but show the feasibility of the third one. We show that such disks can survive for gigayears, largely preserving the primordial size distribution. They should be composed of macroscopic solids larger than millimeters, but smaller than kilometers in size. Thus planetesimal formation, at least in the outer regions of the systems, has stopped before "cometary" or "asteroidal" sizes were reached.

[Abridged] Debris disks are extrasolar analogs to the solar system planetesimal belts. The star Fomalhaut harbors a cold debris belt at 140 AU as well as evidence of a warm dust component, which is suspected of being a bright analog to the solar system's zodiacal dust. Interferometric observations obtained with the VLTI and the KIN have identified near- and mid-infrared excesses attributed to hot and warm exozodiacal dust in the inner few AU of the star. We performed parametric modeling of the exozodiacal disk using the GRaTeR radiative transfer code to reproduce the interferometric data, complemented by mid- to far-infrared measurements. A detailed treatment of sublimation temperatures was introduced to explore the hot population at the sublimation rim. We then used an analytical approach to successively testing several source mechanisms. A good fit to the data is found by two distinct dust populations: (1) very small, hence unbound, hot dust grains confined in a narrow region at the sublimation rim of carbonaceous material; (2) bound grains at 2 AU that are protected from sublimation and have a higher mass despite their fainter flux level. We propose that the hot dust is produced by the release of small carbon grains following the disruption of aggregates that originate from the warm component. A mechanism, such as gas braking, is required to further confine the small grains for a long enough time. In situ dust production could hardly be ensured for the age of the star, so the observed amount of dust must be triggered by intense dynamical activity. Fomalhaut may be representative of exozodis that are currently being surveyed worldwide. We propose a framework for reconciling the hot exozodi phenomenon with theoretical constraints: the hot component of Fomalhaut is likely the tip of the iceberg since it could originate from a warm counterpart residing near the ice line.

Nulling interferometry aims to detect faint objects close to bright stars. Its principle is to produce a destructive interference along the line-of-sight so that the stellar flux is rejected, while the flux of the off-axis source can be transmitted. In practice, various instrumental perturbations can degrade the nulling performance. Any imperfection in phase, amplitude, or polarization produces a spurious flux that leaks to the interferometer output and corrupts the transmitted off-axis flux. One of these instrumental pertubations is the crosstalk phenomenon, which occurs because of multiple parasitic reflections inside transmitting optics, and/or diffraction effects related to beam propagation along finite size optics. It can include a crosstalk of a beam with itself, and a mutual crosstalk between different beams. This can create a parasitic interference pattern, which degrades the intrinsic transmission map - or intensity response - of the interferometer. In this context, we describe how this instrumental effect impairs the performance of a Bracewell interferometer. A simple formalism is developed to derive the corresponding modified intensity response of the interferometer, as a function of the two parameters of interest: the crosstalk level (or contamination rate) and the phase shift between the primary and secondary - parasitic - beams. We then apply our mathematical approach to a few scientific cases, both analytically and using the GENIEsim simulation software, adapted to handle coherent crosstalk. Our results show that a coherent crosstalk level of about 1% implies a 20% drop of the SNR at most. Careful attention should thus be paid to reduce the crosstalk level inside an interferometric instrument and ensure an instrumental stability that provides the necessary sensitivity through calibration procedures.

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $\mu$m were obtained, complemented in some cases with observations at 70 $\mu$m, and at 250, 350 and 500 $\mu$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $\mu$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$\mu$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

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

Context. High contrast imaging has thoroughly combed through the limited search space accessible with first-generation ground-based adaptive optics instruments and the Hubble Space Telescope. Only a few objects were discovered, and many non-detections reported and statistically interpreted. The field is now in need of a technological breakthrough. Aim. Our aim is to open a new search space with first-generation systems such as NACO at the Very Large Telescope, by providing ground-breaking inner working angle (IWA) capabilities in the L' band. The L' band is a sweet spot for high contrast coronagraphy since the planet-to-star brightness ratio is favorable, while the Strehl ratio is naturally higher. Methods. An annular groove phase mask (AGPM) vector vortex coronagraph optimized for the L' band, made from diamond subwavelength gratings was manufactured and qualified in the lab. The AGPM enables high contrast imaging at very small IWA, potentially being the key to unexplored discovery space. Results. Here we present the installation and successful on-sky tests of an L'-band AGPM coronagraph on NACO. Using angular differential imaging, which is well suited to the rotational symmetry of the AGPM, we demonstrated a \Delta L' > 7.5 mag contrast from an IWA ~ 0".09 onwards, during average seeing conditions, and for total integration times of a few hundred seconds.

Coronagraphy is a powerful technique to achieve high contrast imaging and hence to image faint companions around bright targets. Various concepts have been used in the visible and near-infrared regimes, while coronagraphic applications in the mid-infrared remain nowadays largely unexplored. Vector vortex phase masks based on concentric subwavelength gratings show great promise for such applications. We aim at producing and validating the first high-performance broadband focal plane phase mask coronagraphs for applications in the mid-infrared regime, and in particular the L band with a fractional bandwidth of ~16% (3.5-4.1 \mu m). Based on rigorous coupled wave analysis, we designed an annular groove phase mask (AGPM) producing a vortex effect in the L band, and etched it onto a series of diamond substrates. The grating parameters were measured by means of scanning electron microscopy. The resulting components were then tested on a mid-infrared coronagraphic test bench. A broadband raw null depth of 2 x 10^{-3} was obtained for our best L-band AGPM after only a few iterations between design and manufacturing. This corresponds to a raw contrast of about 6 x 10^{-5} (10.5 mag) at 2\lambda/D. This result is fully in line with our projections based on rigorous coupled wave analysis modeling, using the measured grating parameters. The sensitivity to tilt and focus has also been evaluated. After years of technological developments, mid-infrared vector vortex coronagraphs finally become a reality and live up to our expectations. Based on their measured performance, our L-band AGPMs are now ready to open a new parameter space in exoplanet imaging at major ground-based observatories.

Warm debris disks are a sub-sample of the large population of debris disks, and display excess emission in the mid-IR. Around solar-type stars, very few objects show emission features in mid-IR spectroscopic observations, that are attributed to small, warm silicate dust grains. The origin of this warm dust can possibly be explained either by a collision between several bodies or by transport from an outer belt. We present and analyse new far-IR Herschel/Pacs observations, supplemented by ground-based data in the mid-IR (VLTI/Midi and VLT/Visir), for one of these rare systems: the 10-16 Myr old debris disk around HD 113766 A. We improve an existing model to account for these new observations, and better constrain the spatial distribution of the dust and its composition. We underline the limitations of SED modelling and the need for spatially resolved observations. We find that the system is best described by an inner disk located within the first AU, well constrained by the Midi data, and an outer disk located between 9-13 AU. In the inner dust belt, our previous finding of Fe-rich crystalline olivine grains still holds. We do not observe time variability of the emission features over at least a 8 years time span, in a environment subjected to strong radiation pressure. The time stability of the emission features indicates that {\mu}m-sized dust grains are constantly replenished from the same reservoir, with a possible depletion of sub-{\mu}m-sized grains. We suggest that the emission features may arise from multi-composition aggregates. We discuss possible scenarios concerning the origin of the warm dust. The compactness of the innermost regions as probed by Midi, as well as the dust composition, suggest that we are witnessing the outcomes of (at least) one collision between partially differentiated bodies, in an environment possibly rendered unstable by terrestrial planetary formation.

Cyg OB2 #9 is one of a small set of non-thermal radio emitting massive O-star binaries. The non-thermal radiation is due to synchrotron emission in the colliding-wind region. Cyg OB2 #9 was only recently discovered to be a binary system and a multi-wavelength campaign was organized to study its 2011 periastron passage. We report here on the results of the radio observations obtained in this monitoring campaign. We used the Expanded Very Large Array (EVLA) radio interferometer to obtain 6 and 20 cm continuum fluxes. The observed radio light curve shows a steep drop in flux sometime before periastron. The fluxes drop to a level that is comparable to the expected free-free emission from the stellar winds, suggesting that the non-thermal emitting region is completely hidden at that time. After periastron passage, the fluxes slowly increase. We introduce a simple model to solve the radiative transfer in the stellar winds and the colliding-wind region, and thus determine the expected behaviour of the radio light curve. From the asymmetry of the light curve, we show that the primary has the stronger wind. This is somewhat unexpected if we use the astrophysical parameters based on theoretical calibrations. But it becomes entirely feasible if we take into account that a given spectral type - luminosity class combination covers a range of astrophysical parameters. The colliding-wind region also contributes to the free-free emission, which can help to explain the high values of the spectral index seen after periastron passage. Combining our data with older Very Large Array (VLA) data allows us to derive a period P = 860.0 +- 3.7 days for this system. With this period, we update the orbital parameters that were derived in the first paper of this series.

We report on high contrast mid-infrared observations of Fomalhaut obtained with the Keck Interferometer Nuller (KIN) showing a small resolved excess over the level expected from the stellar photosphere. The measured null excess has a mean value of 0.35% +/- 0.10% between 8 and 11 microns and increases from 8 to 13 microns. Given the small field of view of the instrument, the source of this marginal excess must be contained within 2AU of Fomalhaut. This result is reminiscent of previous VLTI K-band observations, which implied the presence of a ~ 0.88% excess, and argued that thermal emission from hot dusty grains located within 6 AU from Fomalhaut was the most plausible explanation. Using a parametric 2D radiative transfer code and a Bayesian analysis, we examine different dust disk structures to reproduce both the near and mid-infrared data simultaneously. While not a definitive explanation of the hot excess of Fomalhaut, our model suggests that the most likely inner few AU disk geometry consists of a two-component structure, with two different and spatially distinct grain populations. The 2 to 11 microns data are consistent with an inner hot ring of very small (~ 10 to 300 nm) carbon-rich grains concentrating around 0.1AU. The second dust population consists of larger grains (size of a few microns to a few tens of microns) located further out in a colder region where regular astronomical silicates could survive, with an inner edge around 1AU. From a dynamical point of view, the presence of the inner concentration of sub-micron sized grains is surprising, as such grains should be expelled from the inner planetary system by radiation pressure within only a few years. This could either point to some inordinate replenishment rates (e.g. many grazing comets coming from an outer reservoir) or to the existence of some braking mechanism preventing the grains from moving out.

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

Context: Nonthermal radio emission in massive stars is expected to arise in wind-wind collisions occurring inside a binary system. One such case, the O-type star Cyg OB2 #9, was proven to be a binary only four years ago, but the orbital parameters remained uncertain. The periastron passage of 2011 was the first one to be observable under good conditions since the discovery of binarity. Aims: In this context, we have organized a large monitoring campaign to refine the orbital solution and to study the wind-wind collision. Methods: This paper presents the analysis of optical spectroscopic data, as well as of a dedicated X-ray monitoring performed with Swift and XMM. Results: In light of our refined orbital solution, Cyg OB2 #9 appears as a massive O+O binary with a long period and high eccentricity; its components (O5-5.5I for the primary and O3-4III for the secondary) have similar masses and similar luminosities. The new data also provide the first evidence that a wind-wind collision is present in the system. In the optical domain, the broad Ha line varies, displaying enhanced absorption and emission components at periastron. X-ray observations yield the unambiguous signature of an adiabatic collision because, as the stars approach periastron, the X-ray luminosity closely follows the 1/D variation expected in that case. The X-ray spectrum appears, however, slightly softer at periastron, which is probably related to winds colliding at slightly lower speeds at that time. Conclusions: It is the first time that such a variation has been detected in O+O systems, and the first case where the wind-wind collision is found to remain adiabatic even at periastron passage.

Most exoplanet imagers consist of ground-based adaptive optics coronagraphic cameras which are currently limited in contrast, sensitivity and astrometric precision, but advantageously observe in the near-IR (1- 5{\mu}m). Because of these practical limitations, our current observational aim at detecting and characterizing planets puts heavy constraints on target selection, observing strategies, data reduction, and follow-up. Most surveys so far have thus targeted young systems (1-100Myr) to catch the putative remnant thermal radiation of giant planets, which peaks in the near-IR. They also favor systems in the solar neighborhood (d<80pc), which eases angular resolution requirements but also ensures a good knowledge of the distance and proper motion, which are critical to secure the planet status, and enable subsequent characterization. Because of their youth, it is very tempting to target the nearby star forming regions, which are typically twice as far as the bulk of objects usually combed for planets by direct imaging. Probing these interesting reservoirs sets additional constraints that we review in this paper by presenting the planet search that we initiated in 2008 around the disk-bearing T Tauri star IM Lup (Lupus star forming region, 140-190pc). We show and discuss why age determination, the choice of evolutionary model for the central star and the planet, precise knowledge of the host star proper motion, relative or absolute astrometric accuracy, and patience are the key ingredients for exoplanet searches around more distant young stars. Unfortunately, most of the time, precision and perseverance are not paying off: we discovered a candidate companion around IM Lup in 2008, which we report here to be an unbound background object. We nevertheless review in details the lessons learned from our endeavor, and additionally present the best detection limits ever calculated for IM Lup.

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

Debris dust in the habitable zones of stars - otherwise known as exozodiacal dust - comes from extrasolar asteroids and comets and is thus an expected part of a planetary system. Background flux from the Solar System's zodiacal dust and the exozodiacal dust in the target system is likely to be the largest source of astrophysical noise in direct observations of terrestrial planets in the habitable zones of nearby stars. Furthermore, dust structures like clumps, thought to be produced by dynamical interactions with exoplanets, are a possible source of confusion. In this paper, we qualitatively assess the primary impact of exozodical dust on high-contrast direct imaging at optical wavelengths, such as would be performed with a coronagraph. Then we present the sensitivity of previous, current, and near-term facilities to thermal emission from debris dust at all distances from nearby solar-type stars, as well as our current knowledge of dust levels from recent surveys. Finally, we address the other method of detecting debris dust, through high-contrast imaging in scattered light. This method is currently far less sensitive than thermal emission observations, but provides high spatial resolution for studying dust structures. This paper represents the first report of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG).

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

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

Aims. We present photometric data of debris disks around HIP 103389 (HD 199260), HIP 107350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3 sigma sensitivity of a few mJy at 100 um and 160 um. In addition, we obtained Herschel/PACS photometric data at 70 um for HIP 103389. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70 um. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a very distinct range of grain sizes is implied to dominate the thermal emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented. (abridged)

We present Herschel PACS 100 and 160 micron observations of the solar-type stars alpha Men, HD 88230 and HD 210277, which form part of the FGK stars sample of the Herschel Open Time Key Programme (OTKP) DUNES (DUst around NEarby S tars). Our observations show small infrared excesses at 160 micron for all three stars. HD 210277 also shows a small excess at 100 micron, while the 100 micron fluxes of alpha Men and HD 88230 agree with the stellar photospheric predictions. We attribute these infrared excesses to a new class of cold, faint debris discs. alpha Men and HD 88230 are spatially resolved in the PACS 160 micron images, while HD 210277 is point-like at that wavelength. The projected linear sizes of the extended emission lie in the range from ~ 115 to ~ 250 AU. The estimated black body temperatures from the 100 and 160 micron fluxes are $\lesssim$ 22 K, while the fractional luminosity of the cold dust is Ldust/Lstar ~ 10E-6, close to the luminosity of the Solar-System's Kuiper belt. These debris discs are the coldest and faintest discs discovered so far around mature stars and cannot easily be explained by invoking "classical" debris disc models.

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

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

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

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

Using a vector vortex coronagraph behind the 1.5-m well-corrected subaperture (WCS) at Palomar, we detected a second object very close to \epsilon Cephei, a \delta Scuti F0 IV star. The candidate companion, ~50 times fainter than \epsilon Cephei, if physically associated, is a late-type K or early M star, and lies at an angular separation of 330 mas, or 1.1 \lambda/D for the WCS, making it the smallest angle detection ever realized with a coronagraph in terms of \lambda/D units. The projected separation of the putative companion is ~8.6 AU, most likely on a highly eccentric orbit. The recently detected near-infrared excess is thus likely not due to hot dust. Moreover, we also show that the previously reported IRAS 60 \mu m excess was due to source confusion on the galactic plane.

Dusty debris discs around main sequence stars are thought to be the result of continuous collisional grinding of planetesimals in the system. The majority of these systems are unresolved and analysis of the dust properties is limited by the lack of information regarding the dust location.vThe Herschel DUNES key program is observing 133 nearby, Sun-like stars (<20 pc, FGK spectral type) in a volume limited survey to constrain the absolute incidence of cold dust around these stars by detection of far infrared excess emission at flux levels comparable to the Edgeworth-Kuiper belt (EKB). We have observed the Sun-like star HD 207129 with Herschel PACS and SPIRE. In all three PACS bands we resolve a ring-like structure consistent with scattered light observations. Using {\alpha} Bo\"otis as a reference point spread function (PSF), we deconvolved the images, clearly resolving the inner gap in the disc at both 70 and 100 {\mu}m. We have resolved the dust-producing planetesimal belt of a debris disc at 100 {\mu}m for the first time. We measure the radial profile and fractional luminosity of the disc, and compare the values to those of discs around stars of similar age and/or spectral type, placing this disc in context of other resolved discs observed by Herschel/DUNES.

Context. With the arrival of the next generation of ground-based imaging interferometers combining from 4 to possibly 6 telescopes simultaneously, there is also a strong need for a new generation of fringe trackers able to cophase such arrays. These instruments have to be very sensitive and to provide robust operations in quickly varying observational conditions. Aims. We aim at defining the optimal characteristics of fringe sensor concepts operating with 4 or 6 telescopes. The current detector limitations impose us to consider solutions based on co-axial pairwise combination schemes. Methods. We independently study several aspects of the fringe sensing process: 1) how to measure the phase and the group delay, and 2) how to combine the telescopes in order to ensure a precise and robust fringe tracking in real conditions. Thanks to analytical developments and numerical simulations, we define the optimal fringe-sensor concepts and compute the expected performance of the 4-telescope one with our dedicated end-to-end simulation tool sim2GFT. Results. We first show that measuring the phase and the group delay by obtaining the data in several steps (i.e. by temporally modulating the optical path difference) is extremely sensitive to atmospheric turbulence and therefore conclude that it is better to obtain the fringe position with a set of data obtained simultaneously. Subsequently, we show that among all co-axial pairwise schemes, moderately redundant concepts increase the sensitivity as well as the robustness in various atmospheric or observing conditions. Merging all these results, end-to-end simulations show that our 4-telescope fringe sensor concept is able to track fringes at least 90% of the time up to limiting magnitudes of 7.5 and 9.5 for the 1.8- and 8.2-meter VLTI telescopes respectively.

When observing an extrasolar planetary system, the most luminous component after the star itself is generally the light scattered and/or thermally emitted by a population of micron-sized dust grains. These grains are expected to be continuously replenished by the collisions and evaporation of larger bodies just as in our solar zodiacal cloud. Exozodiacal clouds ("exozodis") must therefore be seriously taken into account when attempting to directly image faint Earth-like planets (exoEarths, for short). This paper summarizes the oral contributions and discussions that took place during the Satellite Meeting on exozodiacal dust disks, in an attempt to address the following two questions: Do we need to solve the exozodi question? If yes, how to best solve it?

In this paper, we review the various ways in which an infrared stellar interferometer can be used to perform direct detection of extrasolar planetary systems. We first review the techniques based on classical stellar interferometry, where (complex) visibilities are measured, and then describe how higher dynamic ranges can be achieved with nulling interferometry. The application of nulling interferometry to the study of exozodiacal discs and extrasolar planets is then discussed and illustrated with a few examples.

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

PIONIER is a 4-telescope visitor instrument for the VLTI, planned to see its first fringes in 2010. It combines four ATs or four UTs using a pairwise ABCD integrated optics combiner that can also be used in scanning mode. It provides low spectral resolution in H and K band. PIONIER is designed for imaging with a specific emphasis on fast fringe recording to allow closure-phases and visibilities to be precisely measured. In this work we provide the detailed description of the instrument and present its updated status.

In a few years, the second generation instruments of the Very Large Telescope Interferometer (VLTI) will routinely provide observations with 4 to 6 telescopes simultaneously. To reach their ultimate performance, they will need a fringe sensor capable to measure in real time the randomly varying optical paths differences. A collaboration between LAOG (PI institute), IAGL, OCA and GIPSA-Lab has proposed the Planar Optics Phase Sensor concept to ESO for the 2nd Generation Fringe Tracker. This concept is based on the integrated optics technologies, enabling the conception of extremely compact interferometric instruments naturally providing single-mode spatial filtering. It allows operations with 4 and 6 telescopes by measuring the fringes position thanks to a spectrally dispersed ABCD method. We present here the main analysis which led to the current concept as well as the expected on-sky performance and the proposed design.

We present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

About two dozen exo-solar debris systems have been spatially resolved. These debris discs commonly display a variety of structural features such as clumps, rings, belts, eccentric distributions and spiral patterns. In most cases, these features are believed to be formed, shaped and maintained by the dynamical influence of planets orbiting the host stars. In very few cases has the presence of the dynamically important planet(s) been inferred from direct observation. The solar-type star q1 Eri is known to be surrounded by debris, extended on scales of < 30''. The star is known to host at least one planet, albeit on an orbit far too small to make it responsible for structures at distances of tens to hundreds of AU. The aim of the present investigation is twofold: to determine the optical and material properties of the debris and to infer the spatial distribution of the dust, which may hint at the presence of additional planets. The photodetector array camera and spectrometer (PACS) aboard the Herschel Space Observatory allows imaging observations in the far infrared at unprecedented resolution, i.e. at better than 6'' to 12'' over the wavelength range of 60 {\mu}m to 210 {\mu}m. Together with the results from ground-based observations, these spatially resolved data can be modelled to determine the nature of the debris and its evolution more reliably than would be possible from unresolved data alone. For the first time has the q1 Eri disc been resolved at far infrared wavelengths. The PACS observations at 70, 100 and 160 {\mu}m reveal an oval image showing a disc-like structure in all bands, the size of which increases with wavelength. Assuming a circular shape yields the inclination of its equatorial plane with respect to that of the sky, i > 53deg. The results of image de-convolution indicate that i likely is larger than 63deg, where 90deg corresponds to an edge-on disc. {abridged}

The direct images of giant extrasolar planets recently obtained around several main sequence stars represent a major step in the study of planetary systems. These high-dynamic range images are among the most striking results obtained by the current generation of high angular resolution instruments, which will be superseded by a new generation of instruments in the coming years. It is therefore an appropriate time to review the contributions of high angular resolution visible/infrared techniques to the rapidly growing field of extrasolar planetary science. During the last 20 years, the advent of the Hubble Space Telescope, of adaptive optics on 4- to 10-m class ground-based telescopes, and of long-baseline infrared stellar interferometry has opened a new viewpoint on the formation and evolution of planetary systems. By spatially resolving the optically thick circumstellar discs of gas and dust where planets are forming, these instruments have considerably improved our models of early circumstellar environments and have thereby provided new constraints on planet formation theories. High angular resolution techniques are also directly tracing the mechanisms governing the early evolution of planetary embryos and the dispersal of optically thick material around young stars. Finally, mature planetary systems are being studied with an unprecedented accuracy thanks to single-pupil imaging and interferometry, precisely locating dust populations and putting into light a whole new family of long-period giant extrasolar planets.

The Fourier-Kelvin Stellar Interferometer (FKSI) is a mission concept for a nulling interferometer for the near-to-mid-infrared spectral region. FKSI is conceived as a mid-sized strategic or Probe class mission. FKSI has been endorsed by the Exoplanet Community Forum 2008 as such a mission and has been costed to be within the expected budget. The current design of FKSI is a two-element nulling interferometer. The two telescopes, separated by 12.5 m, are precisely pointed (by small steering mirrors) on the target star. The two path lengths are accurately controlled to be the same to within a few nanometers. A phase shifter/beam combiner (Mach-Zehnder interferometer) produces an output beam consisting of the nulled sum of the target planet's light and the host star's light. When properly oriented, the starlight is nulled by a factor of 10^-4, and the planet light is undiminished. Accurate modeling of the signal is used to subtract the residual starlight, permitting the detection of planets much fainter than the host star. The current version of FKSI with 0.5-m apertures and waveband 3-8 microns has the following main capabilities: (1) detect exozodiacal emission levels to that of our own solar system (1 Solar System Zodi) around nearby F, G, and K, stars; (2) characterize spectroscopically the atmospheres of a large number of known non-transiting planets; (3) survey and characterize nearby stars for planets down to 2 Earth radii from just inside the habitable zone and inward. An enhanced version of FKSI with 1-m apertures separated by 20 m and cooled to 40 K, with science waveband 5-15 microns, allows for the detection and characterization of 2 Earth-radius super-Earths and smaller planets in the habitable zone around stars within about 30 pc.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

Earth-sized planets around nearby stars are being detected for the first time by ground-based radial velocity and space-based transit surveys. This milestone is opening the path towards the definition of missions able to directly detect the light from these planets, with the identification of bio-signatures as one of the main objectives. In that respect, both ESA and NASA have identified nulling interferometry as one of the most promising techniques. The ability to study distant planets will however depend on exozodiacal dust clouds surrounding the target stars. In this paper, we assess the impact of exozodiacal dust clouds on the performance of an infrared nulling interferometer in the Emma X-array configuration. For the nominal mission architecture with 2-m aperture telescopes, we found that point-symmetric exozodiacal dust discs about 100 times denser than the solar zodiacal cloud can be tolerated in order to survey at least 150 targets during the mission lifetime. Considering modeled resonant structures created by an Earth-like planet orbiting at 1 AU around a Sun-like star, we show that the tolerable dust density for planet detection goes down to about 15 times the solar zodiacal density for face-on systems and decreases with the disc inclination. The upper limits on the tolerable exozodiacal dust density derived in this study must be considered as rather pessimistic, but still give a realistic estimation of the typical sensitivity that we will need to reach on exozodiacal discs in order to prepare the scientific programme of future Earth-like planet characterisation missions.

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

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

We present infrared interferometric observations of the inner regions of two A-star debris disks, beta Leo and zeta Lep, using the FLUOR instrument at the CHARA interferometer on both short (30 m) and long (>200 m) baselines. For the target stars, the short baseline visibilities are lower than expected for the stellar photosphere alone, while those of a check star, delta Leo, are not. We interpret this visibility offset of a few percent as a near-infrared excess arising from dust grains which, due to the instrumental field of view, must be located within several AU of the central star. For beta Leo, the near-infrared excess producing grains are spatially distinct from the dust which produces the previously known mid-infrared excess. For zeta Lep, the near-infrared excess may be spatially associated with the mid-infrared excess producing material. We present simple geometric models which are consistent with the near and mid-infrared excess and show that for both objects, the near-infrared producing material is most consistent with a thin ring of dust near the sublimation radius with typical grain sizes smaller than the nominal radiation pressure blowout radius. Finally, we discuss possible origins of the near-infrared emitting dust in the context of debris disk evolution models.

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)

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.

Characterising the circumstellar dust around nearby main sequence stars is a necessary step in understanding the planetary formation process and is crucial for future life-finding space missions such as ESA's Darwin or NASA's Terrestrial Planet Finder (TPF). Besides paving the technological way to Darwin/TPF, the space-based infrared interferometers Pegase and FKSI (Fourier-Kelvin Stellar Interferometer) will be valuable scientific precursors in that respect. In this paper, we investigate the performance of Pegase and FKSI for exozodiacal disc detection and compare the results with ground-based nulling interferometers. Besides their main scientific goal (characterising hot giant extrasolar planets), Pegase and FKSI are very efficient in assessing within a few minutes the level of circumstellar dust in the habitable zone around nearby main sequence stars. They are capable of detecting exozodiacal discs respectively 5 and 1 time as dense as the solar zodiacal cloud and they outperform any ground-based instrument. Unlike Pegase, FKSI can achieve this sensitivity for most targets of the Darwin/TPF catalogue thanks to an appropriate combination of baseline length and observing wavelength. The sensitivity of Pegase could, however, be significantly boosted by considering a shorter interferometric baseline length. These space-based interferometers would be complementary to Antarctica-based instruments in terms of sky coverage and would be ideal instruments for preparing future life-finding space missions.

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.

In this paper, we study the feasibility of obtaining near-infrared spectra of bright extrasolar planets with the 2nd generation VLTI Spectro-Imager instrument (VSI), which has the required angular resolution to resolve nearby hot Extrasolar Giant Planets (EGPs) from their host stars. Taking into account fundamental noises, we simulate closure phase measurements of several extrasolar systems using four 8-m telescopes at the VLT and a low spectral resolution (R = 100). Synthetic planetary spectra from T. Barman are used as an input. Standard chi2-fitting methods are then used to reconstruct planetary spectra from the simulated data. These simulations show that low-resolution spectra in the H and K bands can be retrieved with a good fidelity for half a dozen targets in a reasonable observing time (about 10 hours, spread over a few nights). Such observations would strongly constrain the planetary temperature and albedo, the energy redistribution mechanisms, as well as the chemical composition of their atmospheres. Systematic errors, not included in our simulations, could be a serious limitation to these performance estimations. The use of integrated optics is however expected to provide the required instrumental stability (around 10^-4 on the closure phase) to enable the first thorough characterisation of extrasolar planetary emission spectra in the near-infrared.

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.

High-precision interferometric observations of six early-type main sequence stars known to harbour cold debris discs have been obtained in the near-infrared K band with the FLUOR instrument at the CHARA Array. The measured squared visibilities are compared to the expected visibility of the stellar photospheres based on theoretical photospheric models taking into account rotational distortion, searching for potential visibility reduction at short baselines due to circumstellar emission. Our observations bring to light the presence of resolved circumstellar emission around one of the six target stars (zeta Aql) at the 5 sigma level. The morphology of the emission source cannot be directly constrained because of the sparse spatial frequency sampling of our interferometric data. Using complementary adaptive optics observations and radial velocity measurements, we find that the presence of a low-mass companion is a likely origin for the excess emission. The potential companion has a K-band contrast of four magnitudes, a most probable mass of about 0.6 Msun, and is expected to orbit between about 5.5 AU and 8 AU from its host star assuming a purely circular orbit. Nevertheless, by adjusting a physical debris disc model to the observed Spectral Energy Distribution of the zeta Aql system, we also show that the presence of hot dust within 10 AU from zeta Aql, producing a total thermal emission equal to 1.69 +- 0.31% of the photospheric flux in the K band, is another viable explanation for the observed near-infrared excess. Our re-interpretation of archival near- to far-infrared photometric measurements shows however that cold dust is not present around zeta Aql at the sensitivity limit of the IRS and MIPS instruments onboard Spitzer, and urges us to remove zeta Aql from the category of bona fide debris disc stars.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

We present a thorough study of the impact of a migrating planet on a planetesimal disk, by exploring a broad range of masses and eccentricities for the planet. We discuss the sensitivity of the structures generated in debris disks to the basic planet parameters. We perform many N-body numerical simulations, using the symplectic integrator SWIFT, taking into account the gravitational influence of the star and the planet on massless test particles. A constant migration rate is assumed for the planet. The effect of planetary migration on the trapping of particles in mean motion resonances is found to be very sensitive to the initial eccentricity of the planet and of the planetesimals. A planetary eccentricity as low as 0.05 is enough to smear out all the resonant structures, except for the most massive planets. The planetesimals also initially have to be on orbits with a mean eccentricity of less than than 0.1 in order to keep the resonant clumps visible. This numerical work extends previous analytical studies and provides a collection of disk images that may help in interpreting the observations of structures in debris disks. Overall, it shows that stringent conditions must be fulfilled to obtain observable resonant structures in debris disks. Theoretical models of the origin of planetary migration will therefore have to explain how planetary systems remain in a suitable configuration to reproduce the observed structures.

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.

We probed the first 3AU around tau Ceti and epsilon Eridani with the CHARA array (Mt Wilson, USA) in order to gauge the 2micron excess flux emanating from possible hot dust grains in the debris disks and to also resolve the stellar photospheres. High precision visibility amplitude measurements were performed with the FLUOR single mode fiber instrument and telescope pairs on baselines ranging from 22 to 241m of projected length. The short baseline observations allow us to disentangle the contribution of an extended structure from the photospheric emission, while the long baselines constrain the stellar diameter. We have detected a resolved emission around tau Cet, corresponding to a spatially integrated, fractional excess flux of 0.98 +/- 0.21 x 10^{-2} with respect to the photospheric flux in the K'-band. Around eps Eri, our measurements can exclude a fractional excess of greater than 0.6x10^{-2} (3sigma). We interpret the photometric excess around tau Cet as a possible signature of hot grains in the inner debris disk and demonstrate that a faint, physical or background, companion can be safely excluded. In addition, we measured both stellar angular diameters with an unprecedented accuracy: Theta_LD(tau Cet)= 2.015 +/- 0.011 mas and Theta_LD(eps Eri)=2.126 +/- 0.014 mas.

Detecting the presence of circumstellar dust around nearby solar-type main sequence stars is an important pre-requisite for the design of future life-finding space missions such as ESA's Darwin or NASA's Terrestrial Planet Finder (TPF). The high Antarctic plateau may provide appropriate conditions to perform such a survey from the ground. We investigate the performance of a nulling interferometer optimised for the detection of exozodiacal discs at Dome C, on the high Antarctic plateau, and compare it to the expected performance of similar instruments at temperate sites. Based on the currently available measurements of the turbulence characteristics at Dome C, we adapt the GENIEsim software (Absil et al. 2006, A&A 448) to simulate the performance of a nulling interferometer on the high Antarctic plateau. To feed a realistic instrumental configuration into the simulator, we propose a conceptual design for ALADDIN, the Antarctic L-band Astrophysics Discovery Demonstrator for Interferometric Nulling. We assume that this instrument can be placed above the 30-m high boundary layer, where most of the atmospheric turbulence originates. We show that an optimised nulling interferometer operating on a pair of 1-m class telescopes located 30 m above the ground could achieve a better sensitivity than a similar instrument working with two 8-m class telescopes at a temperate site such as Cerro Paranal. The detection of circumstellar discs about 20 times as dense as our local zodiacal cloud seems within reach for typical Darwin/TPF targets in a integration time of a few hours. Moreover, the exceptional turbulence conditions significantly relax the requirements on real-time control loops, which has favourable consequences on the feasibility of the nulling instrument.

The inner Solar System contains a cloud of small (1-100 micron) dust grains created when small bodies-asteroids, comets, and Kuiper belt objects-collide and outgas. This dust cloud, the zodiacal cloud probably has extrasolar analogs, exozodiacal clouds. Exozodiacal clouds are related to debris disks, clouds of rocks and dust orbiting main sequence stars thought to represent the debris left over from planet formation. Some debris disks appear to contain distinct inner clouds that could be considered massive exozodiacal clouds (e.g. Koerner et al. 1998, Absil et al. 2006). This white paper addresses the need for future theoretical work on the dynamics of exozodiacal clouds. This theoretical work should help us discover new planets and understand exozodiacal clouds as astrophysical noise. So far, observations of nearby stars have not provided good constraints on the structures of exozodiacal clouds. But future observations probably will demand a better theoretical understanding of these systems.

Only a handful of debris disks have been imaged up to now. Due to the need for high dynamic range and high angular resolution, very little is known about the inner planetary region, where small amounts of warm dust are expected to be found. We investigate the close neighbourhood of Vega with the help of infrared stellar interferometry and estimate the integrated K-band flux originating from the central 8 AU of the debris disk. We performed precise visibility measurements at both short (~30 m) and long (~150 m) baselines with the FLUOR beam-combiner installed at the CHARA Array (Mt Wilson, California) in order to separately resolve the emissions from the extended debris disk (short baselines) and from the stellar photosphere (long baselines). After revising Vega's K-band angular diameter (3.202+/-0.005 mas), we show that a significant deficit in squared visibility (1.88+/-0.34%) is detected at short baselines with respect to the best-fit uniform disk stellar model. This deficit can be either attributed to the presence of a low-mass stellar companion around Vega, or as the signature of the thermal and scattered emissions from the debris disk. We show that the presence of a close companion is highly unlikely, as well as other possible perturbations (stellar morphology, calibration), and deduce that we have most probably detected the presence of dust in the close neighbourhood of Vega. The resulting flux ratio between the stellar photosphere and the debris disk amounts to 1.29+/-0.19% within the FLUOR field-of-view (~7.8 AU). Finally, we complement our K-band study with archival photometric and interferometric data in order to evaluate the main physical properties of the inner dust disk. The inferred properties suggest that the Vega system could be currently undergoing major dynamical perturbations.

We have obtained high-precision interferometric measurements of Vega with the CHARA Array and FLUOR beam combiner in the K' band at projected baselines between 103m and 273m. The measured visibility amplitudes beyond the first lobe are significantly weaker than expected for a slowly rotating star characterized by a single effective temperature and surface gravity. Our measurements, when compared to synthetic visibilities and synthetic spectrophotometry from a Roche-von Zeipel gravity-darkened model atmosphere, provide strong evidence for the model of Vega as a rapidly rotating star viewed very nearly pole-on. Our best fitting model indicates that Vega is rotating at ~91% of its angular break-up rate with an equatorial velocity of 275 km/s. Together with the measured vsin(i), this velocity yields an inclination for the rotation axis of 5 degrees. For this model the pole-to-equator effective temperature difference is 2250 K, a value much larger than previously derived from spectral line analyses. The derived equatorial T_eff of 7900 K indicates Vega's equatorial atmosphere may be convective and provides a possible explanation for the discrepancy. The model has a luminosity of ~37 Lsun, a value 35% lower than Vega's apparent luminosity based on its bolometric flux and parallax, assuming a slowly rotating star. The model luminosity is consistent with the mean absolute magnitude of A0V stars. Our model predicts the spectral energy distribution of Vega as viewed from its equatorial plane; a model which may be employed in radiative models for the surrounding debris disk.

Nulling interferometry, a powerful technique for high-resolution imaging of the close neighbourhood of bright astrophysical objets, is currently considered for future space missions such as Darwin or the Terrestrial Planet Finder Interferometer (TPF-I), both aiming at Earth-like planet detection and characterization. Ground-based nulling interferometers are being studied for both technology demonstration and scientific preparation of the Darwin/TPF-I missions through a systematic survey of circumstellar dust disks around nearby stars. In this paper, we investigate the influence of atmospheric turbulence on the performance of ground-based nulling instruments, and deduce the major design guidelines for such instruments. End-to-end numerical simulations allow us to estimate the performance of the main subsystems and thereby the actual sensitivity of the nuller to faint exozodiacal disks. Particular attention is also given to the important question of stellar leakage calibration. This study is illustrated in the context of GENIE, the Ground-based European Nulling Interferometer Experiment, to be installed at the VLTI and working in the L' band. We estimate that this instrument will detect exozodiacal clouds as faint as about 50 times the Solar zodiacal cloud, thereby placing strong constraints on the acceptable targets for Darwin/TPF-I.