We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

Abridged: Debris disks are valuable systems to study dust properties. Because they are optically thin at all wavelengths, we have direct access to the properties of dust grains. One very promising technique to study them is to measure their phase function. Disks that are highly inclined are promising targets as a wider range of scattering angles can be probed. The phase function is usually either inferred by comparing the observations to synthetic disk models assuming a parametrized phase function, or estimating it from the surface brightness of the disk. We argue here that the latter approach can be biased due to projection effects leading to an increase in column density along the major axis of a non flat disk. We present a novel approach to account for those column density effects. The method remains model dependent, as one still requires a disk model to estimate the density variations as a function of the scattering angle. This method allows us however to estimate the shape of the phase function without having to invoke any parametrized form. We apply our method to SPHERE/ZIMPOL observations of HR\,4796 and highlight the differences with previous measurements. Our modelling results suggest that the disk is not vertically flat at optical wavelengths. We discuss some of the caveats of the approach, mostly that our method remains blind to real local increase of the dust density, and that it cannot yet be readily applied to angular differential imaging observations. Similarly to previous studies on HR\,4796, we still cannot reconcile the full picture using a given scattering theory to explain the shape of the phase function, a long lasting problem for debris disks. Nonetheless, we argue that similar effects as the ones highlighted in this study can also bias the determination of the phase function in total intensity.

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc.

Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions.

Context. Since 1995 and the first discovery of an exoplanet orbiting a main-sequence star, 4000 exoplanets have been discovered using several techniques. However, only a few of these exoplanets were detected through direct imaging. Indeed, the imaging of circumstellar environments requires high-contrast imaging facilities and accurate control of wavefront aberrations. Ground-based planet imagers such as VLT/SPHERE or Gemini/GPI have already demonstrated great performance. However, their limit of detection is hampered by suboptimal correction of aberrations unseen by adaptive optics (AO). Aims. Instead of focusing on the phase minimization of the pupil plane as in standard AO, we aim to directly minimize the stellar residual light in the SPHERE science camera behind the coronagraph to improve the contrast as close as possible to the inner working angle. Methods. We propose a dark hole (DH) strategy optimized for SPHERE. We used a numerical simulation to predict the global improvement of such a strategy on the overall performance of the instrument for different AO capabilities and particularly in the context of a SPHERE upgrade. Then, we tested our algorithm on the internal source with the AO in closed loop. Results. We demonstrate that our DH strategy can correct for aberrations of phase and amplitude. Moreover, this approach has the ability to strongly reduce the diffraction pattern induced by the telescope pupil and the coronagraph, unlike methods operating at the pupil plane. Our strategy enables us to reach a contrast of 5e-7 at 150 mas from the optical axis in a few minutes using the SPHERE internal source. This experiment establishes the grounds for implementing the algorithm on sky in the near future.

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible.

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.

Context. The wind driven halo is a feature observed within the images delivered by the latest generation of ground-based instruments equipped with an extreme adaptive optics system and a coronagraphic device, such as SPHERE at the VLT. This signature appears when the atmospheric turbulence conditions are varying faster than the adaptive optics loop can correct. The wind driven halo shows as a radial extension of the point spread function along a distinct direction (sometimes referred to as the butterfly pattern). When present, it significantly limits the contrast capabilities of the instrument and prevents the extraction of signals at close separation or extended signals such as circumstellar disks. This limitation is consequential because it contaminates the data a substantial fraction of the time: about 30% of the data produced by the VLT/SPHERE instrument are affected by the wind driven halo.Aims. This paper reviews the causes of the wind driven halo and presents a method to analyze its contribution directly from the scientific images. Its effect on the raw contrast and on the final contrast after post-processing is demonstrated.Methods. We used simulations and on-sky SPHERE data to verify that the parameters extracted with our method are capable of describing the wind driven halo present in the images. We studied the temporal, spatial and spectral variation of these parameters to point out its deleterious effect on the final contrast.Results. The data driven analysis we propose does provide information to accurately describe the wind driven halo contribution in the images. This analysis justifies why this is a fundamental limitation to the final contrast performance reached.Conclusions. With the established procedure, we will analyze a large sample of data delivered by SPHERE in order to propose, in the future, post-processing techniques tailored to remove the wind driven halo.

Tiny meteoroids entering the Earth's atmosphere and inducing meteor showers have long been thought to originate partly from cometary dust. Together with other dust particles, they form a huge cloud around the Sun, the zodiacal cloud. From our previous studies of the zodiacal light, as well as other independent methods (dynamical studies, infrared observations, data related to Earth's environment), it is now established that a significant fraction of dust particles entering the Earth's atmosphere comes from Jupiter-family comets (JFCs). This paper relies on our understanding of key properties of the zodiacal cloud and of comet 67P/Churyumov-Gerasimenko, extensively studied by the Rosetta mission to a JFC. The interpretation, through numerical and experimental simulations of zodiacal light local polarimetric phase curves, has recently allowed us to establish that interplanetary dust is rich in absorbing organics and consists of fluffy particles. The ground-truth provided by Rosetta presently establishes that the cometary dust particles are rich in organic compounds and consist of quite fluffy and irregular aggregates. Our aims are as follows: (1) to make links, back in time, between peculiar micrometeorites, tiny meteoroids, interplanetary dust particles, cometary dust particles, and the early evolution of the Solar System, and (2) to show how detailed studies of such meteoroids and of cometary dust particles can improve the interpretation of observations of dust in protoplanetary and debris disks. Future modeling of dust in such disks should favor irregular porous particles instead of more conventional compact spherical particles.

Point spread function (PSF) reconstruction (PSF-R) is a well established technique to determine reliably and accurately the PSF from Adaptive Optics (AO) control loop data. We have successfully applied this technique to improve the precision on photometry and astrometry to observation of NGC6121 obtained with SPHERE/ZIMPOL as it will be presented in a forthcoming letter. Firstly, we present the methodology we followed to reconstruct the PSF combining pupil-plane and focal-plane measurements using using our PSF-R method PRIME (Beltramo-Martin et al. 2019), with upgrade of both the model and best-fitting steps compared to previous papers. Secondly, we highlight that PRIME allows to maintain the PSF fitting residual below 0.2% over 2 hours of observation and using only 30 s of AO telemetry, which may have important consequences for telemetry storage for PSF-R purpose on future 30-40 m class telescopes. Finally, we deploy PRIME in a more realistic regime using faint stars so as to identify the precision needed on the initial guess parameters to ensure the convergence towards the optimal solution.

Precise photometric and astrometric measurements on astronomical images require an accurate knowledge of the Point Spread Function (PSF). When the PSF cannot be modelled directly from the image, PSF-reconstruction techniques become the only viable solution. So far, however, their performance on real observations has rarely been quantified. Aims. In this Letter, we test the performance of a novel hybrid technique, called PRIME, on Adaptive Optics-assisted SPHERE/ZIMPOL observations of the Galactic globular cluster NGC6121. Methods. PRIME couples PSF-reconstruction techniques, based on control-loop data and direct image fitting performed on the only bright point-like source available in the field of view of the ZIMPOL exposures, with the aim of building the PSF model. Results. By exploiting this model, the magnitudes and positions of the stars in the field can be measured with an unprecedented precision, which surpasses that obtained by more standard methods by at least a factor of four for on-axis stars and by up to a factor of two on fainter, off-axis stars. Conclusions. Our results demonstrate the power of PRIME in recovering precise magnitudes and positions when the information directly coming from astronomical images is limited to only a few point-like sources and, thus, paving the way for a proper analysis of future Extremely Large Telescope observations of sparse stellar fields or individual extragalactic objects.

We performed observations of the Sco-Cen F star HD 117214 aiming at a search for planetary companions and the characterization of the debris disk structure. HD 117214 was observed with the SPHERE subsystems IRDIS, IFS and ZIMPOL at optical and near-IR wavelengths using angular and polarimetric differential imaging techniques. This provided the first images of scattered light from the debris disk with a spatial resolution reaching 25 mas and an inner working angle $< 0.1''$. With the observations with IRDIS and IFS we derive detection limits for substellar companions. The geometrical parameters of the detected disk are constrained by fitting 3D models for the scattering of an optically thin dust disk. Investigating the possible origin of the disk gap, we introduced putative planets therein and modeled the planet-disk and planet-planet dynamical interactions. The obtained planetary architectures are compared with the detection limit curves. The debris disk has an axisymmetric ring structure with a radius of $0.42(\pm 0.01)''$ or $\sim45$ au and an inclination of $71(\pm 2.5)^\circ$ and exhibits a $0.4''$ ($\sim40$ au) wide inner cavity. From the polarimetric data, we derive a polarized flux contrast for the disk of $(F_{\rm pol})_{\rm disk}/F_{\rm \ast}> (3.1 \pm 1.2)\cdot 10^{-4}$ in the RI band. The fractional scattered polarized flux of the disk is eight times smaller than the fractional infrared flux excess. This ratio is similar to the one obtained for the debris disk HIP 79977 indicating that dust radiation properties are not very different between these two disks. Inside the disk cavity we achieve the high sensitivity limits on planetary companions with a mass down to $\sim 4 M_{\rm J}$ at projected radial separations between $0.2''$ and $0.4''$. We can exclude the stellar companions at a radial separation larger than 75 mas from the star.

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

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.

At Paranal Observatory, the least predictable parameter affecting the short-term scheduling of astronomical observations is the optical turbulence, especially the seeing, coherence time and ground layer fraction. These are critical variables driving the performance of the instruments of the Very Large Telescope (VLT), especially those fed with adaptive optics systems. Currently, the night astronomer does not have a predictive tool to support him/her in decision-making at night. As most service-mode observations at the VLT last less than two hours, it is critical to be able to predict what will happen in this time frame, to avoid time losses due to sudden changes in the turbulence conditions, and also to enable more aggressive scheduling. We therefore investigate here the possibility to forecast the turbulence conditions over the next two hours. We call this "turbulence nowcasting", analogously with weather nowcasting, a term already used in meteorology coming from the contraction of "now" and "forecasting". We present here the results of a study based on historical data of the Paranal Astronomical Site Monitoring combined with ancillary data, in a machine learning framework. We show the strengths and shortcomings of such an approach, and present some perspectives in the context of the Extremely Large Telescope.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

Observations of circumstellar environments to look for the direct signal of exoplanets and the scattered light from disks has significant instrumental implications. In the past 15 years, major developments in adaptive optics, coronagraphy, optical manufacturing, wavefront sensing and data processing, together with a consistent global system analysis have enabled a new generation of high-contrast imagers and spectrographs on large ground-based telescopes with much better performance. One of the most productive is the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) designed and built for the ESO Very Large Telescope (VLT) in Chile. SPHERE includes an extreme adaptive optics system, a highly stable common path interface, several types of coronagraphs and three science instruments. Two of them, the Integral Field Spectrograph (IFS) and the Infra-Red Dual-band Imager and Spectrograph (IRDIS), are designed to efficiently cover the near-infrared (NIR) range in a single observation for efficient young planet search. The third one, ZIMPOL, is designed for visible (VIR) polarimetric observation to look for the reflected light of exoplanets and the light scattered by debris disks. This suite of three science instruments enables to study circumstellar environments at unprecedented angular resolution both in the visible and the near-infrared. In this work, we present the complete instrument and its on-sky performance after 4 years of operations at the VLT.

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AU Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morphology at the time of our observations. The loop has a projected width of 1.5 au and rises 2.3 au above the disk midplane. We also explored TESS photometric observations of AU Mic that are consistent with evidence of two starspot complexes in the system. The likely co-alignment of the stellar and disk rotational axes breaks degeneracies in detailed spot modeling, indicating that AU Mic's projected magnetic field axis is offset from its rotational axis. We speculate that small grains in AU Mic's disk could be sculpted by a time-dependent wind that is influenced by this offset magnetic field axis, analogous to co-rotating Solar interaction regions that sculpt and influence the inner and outer regions of our own Heliosphere. Alternatively, if the observed spot modulation is indicative of a significant mis-alignment of the stellar and disk rotational axes, we suggest the disk could still be sculpted by the differential equatorial versus polar wind that it sees with every stellar rotation.

Debris disks are the natural by-products of the planet formation process. Scattered or polarized light observations are mostly sensitive to small dust grains that are released from the grinding down of bigger planetesimals. High angular resolution observations at optical wavelengths can provide key constraints on the radial and azimuthal distribution of the small dust grains. These constraints can help us better understand where most of the dust grains are released upon collisions. We present SPHERE/ZIMPOL observations of the debris disk around HR 4796 A, and model the radial profiles along several azimuthal angles of the disk with a code that accounts for the effect of stellar radiation pressure. This enables us to derive an appropriate description for the radial and azimuthal distribution of the small dust grains. Even though we only model the radial profiles along (or close to) the semi-major axis of the disk, our best-fit model is not only in good agreement with our observations but also with previously published datasets (from near-IR to sub-mm wavelengths). We find that the reference radius is located at $76.4\pm0.4$ au, and the disk has an eccentricity of $0.076_{-0.010}^{+0.016}$, with the pericenter located on the front side of the disk (north of the star). We find that small dust grains must be preferentially released near the pericenter to explain the observed brightness asymmetry. Even though parent bodies spend more time near the apocenter, the brightness asymmetry implies that collisions happen more frequently near the pericenter of the disk. Our model can successfully reproduce the shape of the outer edge of the disk, without having to invoke an outer planet shepherding the debris disk. With a simple treatment of the effect of the radiation pressure, we conclude that the parent planetesimals are located in a narrow ring of about $3.6$ au in width.

Spectro-photometry of debris disks in total intensity and polarimetry can provide new insight into the properties of the dust grains therein (size distribution and optical properties). We aim to constrain the morphology of the highly inclined debris disk HD 32297. We also intend to obtain spectroscopic and polarimetric measurements to retrieve information on the particle size distribution within the disk for certain grain compositions. We observed HD 32297 with SPHERE in Y, J, and H bands in total intensity and in J band in polarimetry. The observations are compared to synthetic models of debris disks and we developed methods to extract the photometry in total intensity overcoming the data-reduction artifacts, namely the self-subtraction. The spectro-photometric measurements averaged along the disk mid-plane are then compared to model spectra of various grain compositions. These new images reveal the very inner part of the system as close as 0.15". The disk image is mostly dominated by the forward scattering making one side (half-ellipse) of the disk more visible, but observations in total intensity are deep enough to also detect the back side for the very first time. The images as well as the surface brightness profiles of the disk rule out the presence of a gap as previously proposed. We do not detect any significant asymmetry between the northeast and southwest sides of the disk. The spectral reflectance features a "gray to blue" color which is interpreted as the presence of grains far below the blowout size. The presence of sub-micron grains in the disk is suspected to be the result of gas drag and/or "avalanche mechanisms". The blue color of the disk could be further investigated with additional total intensity and polarimetric observations in K and H bands respectively to confirm the spectral slope and the fraction of polarization.

Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

Second-generation exoplanet imagers using extreme adaptive optics and coronagraphy have demonstrated their great potential for studying close circumstellar environments and for detecting new companions and helping to understand their physical properties. However, at very small angular separation, their performance in contrast is limited by several factors: diffraction by the complex telescope pupil not perfectly canceled by the coronagraph, residual dynamic wavefront errors, chromatic wavefront errors, and wavefront errors resulting from noncommon path aberrations (NCPAs). In a previous work, we demonstrated the use of a Zernike wavefront sensor called ZELDA for sensing NCPAs in VLT/SPHERE and their compensation. In the present work, we move to the next step with the on-sky validation of NCPA compensation with ZELDA. We start by reproducing previous results on the internal source and show that the amount of aberration integrated between 1 and 15 cycles/pupil is decreased by a factor of five, which translates into a gain in raw contrast of between 2 and 3 below 300 mas. On sky, we demonstrate that NCPA compensation works in closed loop, leading to an attenuation of the amount of aberration by a factor of approximately two. However, we identify a loss of sensitivity for the sensor that is only partly explained by the difference in Strehl ratio between the internal and on-sky measurements. Coronagraphic imaging on sky is improved in raw contrast by a factor of 2.5 at most in the ExAO-corrected region. We use coronagraphic image reconstruction based on a detailed model of the instrument to demonstrate that both internal and on-sky raw contrasts can be precisely explained, and we establish that the observed performance after NCPA compensation is no longer limited by an improper compensation for aberration but by the current apodized-pupil Lyot coronagraph design. [abridged]

High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwarf was discovered with the VLT/SPHERE instrument in 2017, which orbits at $\sim$ 11 au around HD 206893. Its mass was estimated between 12 and 50 $M_{Jup}$ from evolutionary models and its photometry. However, given the significant uncertainty on the age of the system and the peculiar spectrophotometric properties of the companion, this mass is not well constrained. We aim at constraining the orbit and dynamical mass of HD 206893 B. We combined radial velocity data obtained with HARPS spectra and astrometric data obtained with the high contrast imaging VLT/SPHERE and VLT/NaCo instruments, with a time baseline less than three years. We then combined those data with astrometry data obtained by Hipparcos and Gaia with a time baseline of 24 years. We used a MCMC approach to estimate the orbital parameters and dynamical mass of the brown dwarf from those data. We infer a period between 21 and 33{\deg} and an inclination in the range 20-41{\deg} from pole-on from HD 206893 B relative astrometry. The RV data show a significant RV drift over 1.6 yrs. We show that HD 206893 B cannot be the source of this observed RV drift as it would lead to a dynamical mass inconsistent with its photometry and spectra and with Hipparcos and Gaia data. An additional inner (semimajor axis in the range 1.4-2.6 au) and massive ($\sim$ 15 $M_{Jup}$) companion is needed to explain the RV drift, which is compatible with the available astrometric data of the star, as well as with the VLT/SPHERE and VLT/NaCo nondetection.

In a multi-wavelength study of thermal emission and scattered light images we analyse the dust properties and structure of the debris disc around the A1-type main sequence star 49~Cet. As a basis for this study, we present new scattered light images of the debris disc known to possess both a high amount of dust and gas. The outer region of the disc is revealed in former coronagraphic H-band and our new Y-band images from the Very Large Telescope SPHERE instrument. We use the knowledge of the disc's radial extent inferred from ALMA observations and the grain size distribution found by SED fitting to generate semi-dynamical dust models of the disc. We compare the models to scattered light and thermal emission data and find that a disc with a maximum of the surface density at 110~au and shallow edges can describe both thermal emission and scattered light observations. This suggests that grains close to the blow-out limit and large grains stem from the same planetesimal population and are mainly influenced by radiation pressure. The influence of inwards transport processes could not be analysed in this study.

Various structures are visible within Spectro-Polarimetric High-contrast Exoplanet REsearch instrument (SPHERE) images that are not always straightforward to interpret. In this article we present a review of these features and demonstrate their origin using simulations. We also identify which expected or unexpected features are limiting the contrast reached by the instrument and how they may be tackled. This vision paves the way to designing a future upgrade of the SPHERE instrument and the next generation of high-contrast instruments such as those planned for the Extremely Large Telescope (ELT).

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts.

The scattering properties of the dust originating from debris discs are still poorly known. The analysis of scattered light is however a powerful remote-sensing tool to understand the physical properties of dust particles orbiting other stars. Scattered light is indeed widely used to characterise the properties of cometary dust in the solar system. We aim to measure the morphology and scattering properties of the dust from the debris ring around HR4796A in polarised optical light. We obtained high-contrast polarimetric images of HR4796A in the wavelength range 600-900nm with the SPHERE / ZIMPOL instrument on the Very Large Telescope. We measured for the first time the polarised phase function of the dust in a debris system over a wide range of scattering angles in the optical. We confirm that it is incompatible with dust particles being compact spheres under the assumption of the Mie theory, and propose alternative scenarios compatible with the observations, such as particles with irregular surface roughness or aggregate particles.

Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of $\sim$85 and $\sim$115\,au, and with a mutual inclination of about 5$\degb$. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.

Aims: We want to detect and quantify observables related to accretion processes occurring locally in circumstellar disks, which could be attributed to young forming planets. We focus on objects known to host protoplanet candidates and/or disk structures thought to be the result of interactions with planets. Methods: We analyzed observations of 6 young stars (age $3.5-10$ Myr) and their surrounding environments with the SPHERE/ZIMPOL instrument on the VLT in the H$\alpha$ filter (656 nm) and a nearby continuum filter (644.9 nm). Results: We re-detect the known accreting M-star companion HD142527 B with the highest published signal to noise to date in both H$\alpha$ and the continuum. We derive new astrometry ($r = 62.8^{+2.1}_{-2.7}$ mas and $\text{PA} = (98.7\,\pm1.8)^\circ$) and photometry ($\Delta$N_Ha=$6.3^{+0.2}_{-0.3}$ mag, $\Delta$B_Ha=$6.7\pm0.2$ mag and $\Delta$Cnt_Ha=$7.3^{+0.3}_{-0.2}$ mag) for the companion in agreement with previous studies, and estimate its mass accretion rate ($\dot{M}\approx1-2\,\times10^{-10}\,M_\odot\text{ yr}^{-1}$). A faint point-like source around HD135344 B (SAO206462) is also investigated, but a second deeper observation is required to reveal its nature. No other companions are detected. In the framework of our assumptions we estimate detection limits at the locations of companion candidates around HD100546, HD169142 and MWC758 and calculate that processes involving H$\alpha$ fluxes larger than $\sim8\times10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}>10^{-10}-10^{-12}\,M_\odot\text{ yr}^{-1}$) can be excluded. Furthermore, flux upper limits of $\sim10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}<10^{-11}-10^{-12}\,M_\odot \text{ yr}^{-1}$) are estimated within the gaps identified in the disks surrounding HD135344B and TW Hya.

The inclination distribution of the Kuiper belt provides unique constraints on its origin and dynamical evolution, motivating vertically resolved observations of extrasolar planetesimal belts. We present ALMA observations of millimeter emission in the near edge-on planetesimal belt around $\beta$ Pictoris, finding that the vertical distribution is significantly better described by the sum of two Gaussians compared to a single Gaussian. This indicates that, as for the Kuiper belt, the inclination distribution of $\beta$ Pic's belt is better described by the sum of dynamically hot and cold populations rather than a single component. The hot and cold populations have RMS inclinations of 8.9$^{+0.7}_{-0.5}$ and 1.1$^{+0.5}_{-0.5}$ degrees. We also report that an axisymmetric belt model provides a good fit to new and archival ALMA visibilities, and confirm that the midplane is misaligned with respect to $\beta$ Pic b's orbital plane. However, we find no significant evidence for either the inner disk tilt observed in scattered light and CO emission or the South-West/North-East (SW/NE) asymmetry previously reported for millimeter emission. Finally, we consider the origin of the belt's inclination distribution. Secular perturbations from $\beta$ Pic b are unlikely to provide sufficient dynamical heating to explain the hot population throughout the belt's radial extent, and viscous stirring from large bodies within the belt alone cannot reproduce the two populations observed. This argues for an alternative or additional scenario, such as planetesimals being born with high inclinations, or the presence of a `$\beta$ Pic c' planet, potentially migrating outwards near the belt's inner edge.

Circumbinary disks are common around post-asymptotic giant branch (post-AGB) stars with a stellar companion on orbital time scales of a few 100 to few 1000 days. The presence of a disk is usually inferred from the system's spectral energy distribution (SED), and confirmed, for a sub-sample, by interferometric observations. We used the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument on the Very Large Telescope to obtain extreme adaptive optics assisted scattered light images of the post-AGB binary system AR Puppis. Data have been obtained in the V, I, and H bands. Our observations have produced the first resolved images of AR Puppis' circumbinary disk and confirm its edge-on orientation. In our high angular-resolution and high dynamic-range images we identify several structural components such as a dark mid-plane, the disk surface, and arc-like features. We discuss the nature of these components and use complementary photometric monitoring to relate them to the orbital phase of the binary system. Because the star is completely obscured by the disk at visible wavelengths, we conclude that the long-term photometric variability of the system must be caused by variable scattering, not extinction, of star light by the disk over the binary orbit. Finally, we discuss how the short disk life times and fast evolution of the host stars compared to the ages at which protoplanetary disks are typically observed make systems like AR Puppis valuable extreme laboratories to study circumstellar disk evolution and constrain the time scale of dust grain growth during the planet formation process.

(Abr.) Laser guide stars employed at astronomical observatories provide artificial wavefront reference sources to help correct (in part) the impact of atmospheric turbulence on astrophysical observations. Following the recent commissioning of the 4 Laser Guide Star Facility (4LGSF) on UT4 at the VLT, we characterize the spectral signature of the uplink beams from the 22W lasers to assess the impact of laser scattering from the 4LGSF on science observations. We use the MUSE optical integral field spectrograph to acquire spectra at a resolution of R~3000 of the uplink laser beams over the wavelength range of 4750\AA\ to 9350\AA. We report the first detection of laser-induced Raman scattering by N2, O2, CO2, H2O and (tentatively) CH4 molecules in the atmosphere above the astronomical observatory of Cerro Paranal. In particular, our observations reveal the characteristic spectral signature of laser photons -- but 480\AA\ to 2210\AA\ redder than the original laser wavelength of 5889.959\AA\ -- landing on the 8.2m primary mirror of UT4 after being Raman-scattered on their way up to the sodium layer. Laser-induced Raman scattering is not unique to the observatory of Cerro Paranal, but common to any astronomical telescope employing a laser-guide-star (LGS) system. It is thus essential for any optical spectrograph coupled to a LGS system to handle thoroughly the possibility of a Raman spectral contamination via a proper baffling of the instrument and suitable calibrations procedures. These considerations are particularly applicable for the HARMONI optical spectrograph on the upcoming Extremely Large Telescope. At sites hosting multiple telescopes, laser collision prediction tools also ought to account for the presence of Raman emission from the uplink laser beam(s) to avoid the unintentional contamination of observations acquired with telescopes in the vicinity of a LGS system.

Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs.

The Coronagraph Instrument (CGI) for NASA's Wide Field Infrared Survey Telescope (WFIRST) will constitute a dramatic step forward for high-contrast imaging, integral field spectroscopy, and polarimetry of exoplanets and circumstellar disks, aiming to improve upon the sensitivity of current ground-based direct imaging facilities by 2-3 orders of magnitude. Furthermore, CGI will serve as a pathfinder for future exo-Earth imaging and characterization missions by demonstrating wavefront control, coronagraphy, and spectral retrieval in a new contrast regime, and by validating instrument and telescope models at unprecedented levels of precision. To achieve this jump in performance, it is critical to draw on the experience of ground-based high-contrast facilities. We discuss several areas of relevant commonalities, including: wavefront control, post-processing of integral field unit data, and calibration and observing strategies.

We present new observations of the edge-on debris disk around HD 15115 (F star at 48.2 pc) obtained in the near-IR. We search for observational evidence for a second inner planetesimal ring in the system. We obtained total intensity and polarimetric data in the broad bands J and H and processed the data with differential imaging techniques achieving an angular resolution of about 40 mas. We observe an axisymmetric planetesimal belt with a radius of $\sim$2$''$, an inclination of $85.8^{\circ} \pm 0.7^{\circ}$ and position angle of 278.9$^{\circ} \pm 0.1^{\circ}$. A grid of models describing the spatial distribution of the grains in the disk is generated to constrain the geometric parameters of the disk and to explore the presence of a second belt. We perform a photometric analysis of the data and compare disk brightness in two bands in scattered and in polarized light. The analysis shows that the west side is $\sim$2.5 times brighter in total intensity than the east side in both bands, while for polarized light in the J band this ratio is only 1.25. The maximum polarization fraction is 15--20% at $r\sim$2.5$''$. We also find that the J - H color of the disk appears to be red for the radial separations $r\lesssim2''$ and is getting bluer for the larger separations. This apparent change of disk color from red to blue with an increasing radial separation could be explained by the decreasing average grain size with distance. The presence of an additional inner belt slightly inclined with respect to the main planetesimal belt is suspected from the polarized intensity image but the analysis and modeling presented here cannot establish a firm conclusion due to the faintness of the disk and its high inclination.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

HD~105 is a nearby, pre-main sequence G0 star hosting a moderately bright debris disc ($L_{\rm dust}/L_{\star} \sim 2.6\times10^{-4}$). HD~105 and its surroundings might therefore be considered an analogue of the young Solar System. We refine the stellar parameters based on an improved Gaia parallax distance, identify it as a pre-main sequence star {with an age of 50~$\pm$~16~Myr}. The circumstellar disc was marginally resolved by \textit{Herschel}/PACS imaging at far-infrared wavelengths. Here we present an archival ALMA observation at 1.3~mm, revealing the extent and orientation of the disc. We also present \textit{HST}/NICMOS and VLT/SPHERE near-infrared images, where we recover the disc in scattered light at the $\geq$~5-$\sigma$ level. This was achieved by employing a novel annular averaging technique, and is the first time this has been achieved for a disc in scattered light. Simultaneous modelling of the available photometry, disc architecture, and detection in scattered light allow better determination of the disc's architecture, and dust grain minimum size, composition, and albedo. We measure the dust albedo to lie between 0.19 and 0.06, the lower value being consistent with Edgeworth-Kuiper belt objects.

The latest generation of high contrast instruments dedicated to exoplanets and circumstellar disks imaging are equipped with extreme adaptive optics and coronagraphs to reach contrasts of up to 10^--4 at a few tenths of arc-seconds in the near infrared. The resulting image shows faint features, only revealed with this combination, such as the wind driven halo. The wind driven halo is due to the lag between the adaptive optics correction and the turbulence speed over the telescope pupil. However we observe an asymmetry of this wind driven halo that was not expected when the instrument was designed. In this letter, we describe and demonstrate the physical origin of this asymmetry and support our explanation by simulating the asymmetry with an end-to-end approach. From this work, we found out that the observed asymmetry is explained by the interference between the AO-lag error and scintillation effects, mainly originating from the fast jet stream layer located at about 12 km in altitude. Now identified and interpreted, this effect can be taken into account for further design of high contrast imaging simulators, next generation or upgrade of high contrast instruments, predictive control algorithms for adaptive optics or image post-processing techniques.

The first observations of laser guide star photons Raman-scattered by air molecules above the Very Large Telescope (VLT) were reported in June 2017. The initial detection came from the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph, following the installation of the 4 Laser Guide Star Facility (4LGSF) on the Unit Telescope 4 (UT4) of the VLT. In this Letter, we delve further into the symbiotic relationship between the 4LGSF laser guide star system, the UT4 telescope, and MUSE by monitoring the spectral contamination of MUSE observations by Raman photons over a 27 month period. This dataset reveals that dust particles deposited on the primary and tertiary mirrors of UT4 -- responsible for a reflectivity loss of ~8% at 6000{\AA} -- contribute (60$\pm5)% to the laser line fluxes detected by MUSE. The flux of Raman lines, contaminating scientific observations acquired with optical spectrographs, thus provides a new, non-invasive means to monitor the evolving scatter properties of the mirrors of astronomical telescopes equipped with laser guide star systems.

We describe the Zurich Imaging Polarimeter (ZIMPOL), the visual focal plane subsystem of the SPHERE "VLT planet finder", which pushes the limits of current AO systems to shorter wavelengths, higher spatial resolution, and much improved polarimetric performance. We provide new benchmarks for the performance of high contrast instruments, in particular for polarimetric differential imaging. We have analyzed SPHERE/ZIMPOL point spread functions and measure the peak surface brightness, the encircled energy, and the full width half maximum (FWHM) for different wavelengths, atmospheric conditions, star brightness, and instrument modes. Coronagraphic images are described and analized and the performance for different coronagraphs is compared with tests for the binary alpha Hyi with a separation of 92 mas and a contrast of 6 mag. For the polarimetric mode we made the instrument calibrations using zero polarization and high polarization standard stars and here we give a recipe for the absolute calibration of polarimetric data. The data show a small <1 mas but disturbing differential polarimetric beam shifts, which can be explained as Goos-H\"ahnchen shifts from the inclined mirrors, and we discuss how to correct this effect. The polarimetric sensitivity is investigated with non-coronagraphic and deep, coronagraphic observations of the dust scattering around the symbiotic Mira variable R Aqr. SPHERE/ZIMPOL achieves imaging performances in the visual range with unprecedented characteristics, in particular very high spatial resolution and very high polarimetric contrast. This instrument opens up many new research opportunities for the detailed investigation of circumstellar dust, in scattered and therefore polarized light, for the investigation of faint companions, and for the mapping of circumstellar Halpha emission.

The low wind effect is a phenomenon disturbing the phase of the wavefront in the pupil of a large telescope obstructed by spiders, in the absence of wind. It can be explained by the radiative cooling of the spiders, creating air temperature inhomogeneities across the pupil. Because it is unseen by traditional adaptive optics (AO) systems, thus uncorrected, it significantly degrades the quality of AO-corrected images. We provide a statistical analysis of the strength of this effect as seen by VLT/SPHERE after 4 years of operations. We analyse its dependence upon the wind and temperature conditions. We describe the mitigation strategy implemented in 2017: a specific coating with low thermal emissivity in the mid-infrared was applied on the spiders of Unit Telescope 3. We quantify the improvement in terms of image quality, contrast and wave front error using both focal plane images and measured phase maps.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

Proxima Centauri, our closest stellar neighbour, is a low-mass M5 dwarf orbiting in a triple system. An Earth-mass planet with an 11 day period has been discovered around this star. The star's mass has been estimated only indirectly using a mass-luminosity relation, meaning that large uncertainties affect our knowledge of its properties. To refine the mass estimate, an independent method has been proposed: gravitational microlensing. By taking advantage of the close passage of Proxima Cen in front of two background stars, it is possible to measure the astrometric shift caused by the microlensing effect due to these close encounters and estimate the gravitational mass of the lens (Proxima Cen). Microlensing events occurred in 2014 and 2016 with impact parameters, the closest approach of Proxima Cen to the background star, of 1\farcs6 $\pm$ 0\farcs1 and 0\farcs5 $\pm$ 0\farcs1, respectively. Accurate measurements of the positions of the background stars during the last two years have been obtained with HST/WFC3, and with VLT/SPHERE from the ground. The SPHERE campaign started on March 2015, and continued for more than two years, covering 9 epochs. The parameters of Proxima Centauri's motion on the sky, along with the pixel scale, true North, and centering of the instrument detector were readjusted for each epoch using the background stars visible in the IRDIS field of view. The experiment has been successful and the astrometric shift caused by the microlensing effect has been measured for the second event in 2016. We used this measurement to derive a mass of 0.150$^{\textrm{+}0.062}_{-0.051}$ (an error of $\sim$ 40\%) \MSun for Proxima Centauri acting as a lens. This is the first and the only currently possible measurement of the gravitational mass of Proxima Centauri.

SPHERE is the high-contrast exoplanet imager and spectrograph installed at the Unit Telescope 3 of the Very Large Telescope. After more than two years of regular operations, we analyse statistically the performance of the adaptive optics system and its dependence on the atmospheric conditions above the Paranal observatory, as measured by the suite of dedicated instruments which are part of the Astronomical Site Monitor and as estimated by the SPHERE real-time calculator. We also explain how this information can be used to schedule the observations in order to yield the best data quality and to guide the astronomer when processing his/her data.

The major source of noise in high-contrast imaging is the presence of slowly evolving speckles that do not average with time. The temporal stability of the point-spread-function (PSF) is therefore critical to reach a high contrast with extreme adaptive optics (xAO) instruments. Understanding on which timescales the PSF evolves and what are the critical parameters driving the speckle variability allow to design an optimal observing strategy and data reduction technique to calibrate instrumental aberrations and reveal faint astrophysical sources. We have obtained a series of 52 min, AO-corrected, coronagraphically occulted, high-cadence (1.6Hz), H-band images of the star HR 3484 with the SPHERE (Spectro-Polarimeter High-contrast Exoplanet REsearch instrument on the VLT. This is a unique data set from an xAO instrument to study its stability on timescales as short as one second and as long as several tens of minutes. We find different temporal regimes of decorrelation. We show that residuals from the atmospheric turbulence induce a fast, partial decorrelation of the PSF over a few seconds, before a transition to a regime with a linear decorrelation with time, at a rate of several tens parts per million per second (ppm/s). We analyze the spatial dependence of this decorrelation, within the well-corrected radius of the adaptive optics system and show that the linear decorrelation is faster at short separations. Last, we investigate the influence of the distance to the meridian on the decorrelation.

Circumstellar environments are now routinely observed by dedicated high-contrast imagers on large, ground-based observatories. These facilities combine extreme adaptive optics and coronagraphy to achieve unprecedented sensitivities for exoplanet detection and spectral characterization. However, non-common path aberrations (NCPA) in these coronagraphic systems represent a critical limitation for the detection of giant planets with a contrast lower than a few $10^{-6}$ at very small separations ($<$0.3$^{\prime\prime}$) from their host star. In 2013 we proposed ZELDA, a Zernike wavefront sensor to measure these residual quasi-static phase aberrations and a prototype was installed in SPHERE, the exoplanet imager for the VLT. In 2016, we demonstrated the ability of our sensor to provide a nanometric calibration and compensation for these aberrations on an internal source in the instrument, resulting in a contrast gain of 10 at 0.2$^{\prime\prime}$ in coronagraphic images. However, initial on-sky tests in 2017 did not show a tangible gain in contrast when calibrating the NCPA internally and then applying the correction on sky. In this communication, we present recent on-sky measurements to demonstrate the potential of our sensor for the NCPA compensation during observations and quantify the contrast gain in coronagraphic data.

Context. WR104 is an emblematic dusty Wolf-Rayet star and the prototypical member of a subgroup hosting spirals that are mainly observable with high-angular resolution techniques. Previous aperture masking observations showed that WR104 is likely an interacting binary star at the end of its life. However, several aspects of the system are still unknown. This includes the opening angle of the spiral, the dust formation locus, and the link between the central binary star and a candidate companion star detected with the Hubble Space Telescope (HST) at 1". Aims. Our aim was to directly image the dusty spiral or "pinwheel" structure around WR104 for the first time and determine its physical properties at large spatial scales. We also wanted to address the characteristics of the candidate companion detected by the HST. Methods. For this purpose, we used SPHERE and VISIR at the Very Large Telescope to respectively image the system in the near-and mid-infrared. Both instruments furnished an excellent view of the system at the highest angular resolution a single, ground-based telescope can provide. Based on these direct images, we then used analytical and radiative transfer models to determine several physical properties of the system. Results. Employing a different technique than previously used, our new images have allowed us to confirm the presence of the dust pinwheel around the central star. We have also detected up to 5 revolutions of the spiral pattern of WR104 in the K-band for the first time. The circumstellar dust extends up to 2 arcsec from the central binary star in the N-band, corresponding to the past 20 years of mass loss. Moreover, we found no clear evidence of a shadow of the first spiral coil onto the subsequent ones, which likely points to a dusty environment less massive than inferred in previous studies. We have also confirmed that the stellar candidate companion previously detected by the HST is gravitationally bound to WR104 and herein provide information about its nature and orbital elements.

Knowledge of the Earth's atmospheric optical turbulence is critical for astronomical instrumentation. Not only does it enable performance verification and optimisation of existing systems but it is required for the design of future instruments. As a minimum this includes integrated astro-atmospheric parameters such as seeing, coherence time and isoplanatic angle, but for more sophisticated systems such as wide field adaptive optics enabled instrumentation the vertical structure of the turbulence is also required. Stereo-SCIDAR is a technique specifically designed to characterise the Earth's atmospheric turbulence with high altitude resolution and high sensitivity. Together with ESO, Durham University has commissioned a Stereo-SCIDAR instrument at Cerro Paranal, Chile, the site of the Very Large Telescope (VLT), and only 20~km from the site of the future Extremely Large Telescope (ELT). Here we provide results from the first 18 months of operation at ESO Paranal including instrument comparisons and atmospheric statistics. Based on a sample of 83 nights spread over 22 months covering all seasons, we find the median seeing to be 0.64" with 50% of the turbulence confined to an altitude below 2 km and 40% below 600 m. The median coherence time and isoplanatic angle are found as 4.18 ms and 1.75" respectively. A substantial campaign of inter-instrument comparison was also undertaken to assure the validity of the data. The Stereo-SCIDAR profiles (optical turbulence strength and velocity as a function of altitude) have been compared with the Surface-Layer SLODAR, MASS-DIMM and the ECMWF weather forecast model. The correlation coefficients are between 0.61 (isoplanatic angle) and 0.84 (seeing).

In the present study we aim to investigate the circumstellar environment of the spectroscopic binary T Tauri star CS Cha. From unresolved mid- to far-infrared photometry it is predicted that CS Cha hosts a disk with a large cavity. In addition, SED modeling suggests significant dust settling, pointing towards an evolved disk that may show signs of ongoing or completed planet formation. We observed CS Cha with the high contrast imager VLT/SPHERE in polarimetric differential imaging mode to resolve the circumbinary disk in near infrared scattered light. These observations were followed-up by VLT/NACO L-band observations and complemented by archival VLT/NACO K-band and HST/WFPC2 I-band data. We resolve the compact circumbinary disk around CS Cha for the first time in scattered light. We find a smooth, low inclination disk with an outer radius of $\sim$55 au (at 165 pc). We do not detect the inner cavity but find an upper limit for the cavity size of $\sim$15 au. Furthermore, we find a faint co-moving companion with a projected separation of 210 au from the central binary outside of the circumbinary disk. The companion is detected in polarized light and shows an extreme degree of polarization (13.7$\pm$0.4 \% in J-band). The companion's J- and H-band magnitudes are compatible with masses of a few M$_\mathrm{Jup}$. However, K-, L- and I-band data draw this conclusion into question. We explore with radiative transfer modeling whether an unresolved circum-companion disk can be responsible for the high polarization and complex photometry. We find that the set of observations is best explained by a heavily extincted low mass ($\sim 20 \mathrm{M}_\mathrm{Jup}$) brown dwarf or high mass planet with an unresolved disk and dust envelope.

Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris disks around low-mass stars, especially when it comes to spatially resolved observations. We present new VLT/SPHERE IRDIS Dual-Polarization Imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA\,7. Combined with additional Angular Differential Imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. We model the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and perform simple N-body simulations. We find that the dust density distribution peaks at 25 au, with a very shallow outer power-law slope, and that the disk has an inclination of 13 degrees with a position angle of 90 degrees East of North. We also report low signal-to-noise detections of an outer belt at a distance of ~52 au from the star, of a spiral arm in the Southern side of the star, and of a possible dusty clump at 3.9 au. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at 7 au and another belt at 25 au. We report the detections of several unexpected features in the disk around TWA\,7. A yet undetected 100 M$_\oplus$ planet with a semi-major axis at 20-30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.

The nearby and young M star AU Mic is surrounded by a debris disk in which we previously identified a series of large-scale arch-like structures that have never been seen before in any other debris disk and that move outward at high velocities. We initiated a monitoring program with the following objectives: 1) track the location of the structures and better constrain their projected speeds, 2) search for new features emerging closer in, and ultimately 3) understand the mechanism responsible for the motion and production of the disk features. AU Mic was observed at 11 different epochs between August 2014 and October 2017 with the IR camera and spectrograph of SPHERE. These high-contrast imaging data were processed with a variety of angular, spectral, and polarimetric differential imaging techniques to reveal the faintest structures in the disk. We measured the projected separations of the features in a systematic way for all epochs. We also applied the very same measurements to older observations from the Hubble Space Telescope (HST) with the visible cameras STIS and ACS. The main outcomes of this work are 1) the recovery of the five southeastern broad arch-like structures we identified in our first study, and confirmation of their fast motion (projected speed in the range 4-12 km/s); 2) the confirmation that the very first structures observed in 2004 with ACS are indeed connected to those observed later with STIS and now SPHERE; 3) the discovery of two new very compact structures at the northwest side of the disk (at 0.40" and 0.55" in May 2015) that move to the southeast at low speed; and 4) the identification of a new arch-like structure that might be emerging at the southeast side at about 0.4" from the star (as of May 2016). Abridged.

HD95086 (A8V, 17Myr) hosts a rare planetary system for which a multi-belt debris disk and a giant planet of 4-5MJup have been directly imaged. Our study aims to characterize the physical and orbital properties of HD95086b, search for additional planets at short and wide orbits and image the cold outer debris belt in scattered light. We used HARPS at the ESO 3.6m telescope to monitor the radial velocity of HD95086 over 2 years and investigate the existence of giant planets at less than 3au orbital distance. With the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE at VLT, we imaged the faint circumstellar environment beyond 10au at six epochs between 2015 and 2017. We do not detect additional giant planets around HD95086. We identified the nature (bound companion or background contaminant) of all point-like sources detected in the IRDIS field of view. None of them correspond to the ones recently discovered near the edge of the cold outer belt by ALMA. HD95086b is resolved for the first time in J-band with IFS. Its near-infrared spectral energy distribution is well fitted by a few dusty and/or young L7-L9 dwarf spectral templates. The extremely red 1-4um spectral distribution is typical of low-gravity objects at the L/T spectral type transition. The planet's orbital motion is resolved between January 2015 and May 2017. Together with past NaCo measurements properly re-calibrated, our orbital fitting solutions favor a retrograde low to moderate-eccentricity orbit e=0.2 (0.0 to 0.5), with a semi-major axis 52au corresponding to orbital periods of 288$ yrs and an inclination that peaks at i = 141deg, which is compatible with a planet-disk coplanar configuration. Finally, we report the detection in polarimetric differential imaging of the cold outer debris belt between 100 and 300au, consistent in radial extent with recent ALMA 1.3mm resolved observations.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

A large number of systems harboring a debris disk show evidence for a double belt architecture. One hypothesis for explaining the gap between the belts is the presence of one or more planets dynamically carving it. This work aims to investigate this scenario in systems harboring two components debris disks. All the targets in the sample were observed with the SPHERE instrument which performs high-contrast direct imaging. Positions of the inner and outer belts were estimated by SED fitting of the infrared excesses or, when available, from resolved images of the disk. Very few planets have been observed so far in debris disks gaps and we intended to test if such non-detections depend on the observational limits of the present instruments. This aim is achieved by deriving theoretical predictions of masses, eccentricities and semi-major axes of planets able to open the observed gaps and comparing such parameters with detection limits obtained with SPHERE. The relation between the gap and the planet is due to the chaotic zone around the orbit of the planet. The radial extent of this zone depends on the mass ratio between the planet and the star, on the semi-major axis and on the eccentricity of the planet and it can be estimated analytically. We apply the formalism to the case of one planet on a circular or eccentric orbit. We then consider multi-planetary systems: 2 and 3 equal-mass planets on circular orbits and 2 equal-mass planets on eccentric orbits in a packed configuration. We then compare each couple of values (M,a), derived from the dynamical analysis of single and multiple planetary models, with the detection limits obtained with SPHERE. Our results show that the apparent lack of planets in gaps between double belts could be explained by the presence of a system of two or more planets possibly of low mass and on an eccentric orbits whose sizes are below the present detection limits.

We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the $JHK$ bands, possibly correlated with the presence of the shadows. Analysis of archival VLTI/PIONIER $H$-band visibilities constrain the orientation of the inner disk to $i = 18.2\deg^{+3.4}_{-4.1}$ and ${\rm PA} = 57.3\deg \pm 5.7\deg$, consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere ($H/r \lesssim 0.2$). Possible mechanisms that may cause asymmetric variations in the optical depth ($\Delta\tau\lesssim1$) through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon.

We present observations of the known edge-on debris disk around HIP 79977 (HD 146897, F star in Upper Sco, 123 pc), taken with the ZIMPOL differential polarimeter of the SPHERE instrument in the Very Broad Band filter ($\lambda_c=735$ nm, $\Delta\lambda=290$ nm) with a spatial resolution of about 25 mas. We measure the polarization flux along and perpendicular to the disk spine of the highly inclined disk for projected separations between 0.2" (25 AU) and 1.6" (200 AU) and investigate the diagnostic potential of such data with model simulations. The polarized flux contrast ratio for the disk is $F_{pol}/F_\ast= (5.5 \pm 0.9) 10^{-4}$. The surface brightness reaches a maximum of 16.2 mag arcsec$^{-2}$ at a separation of $0.2"-0.5"$ along the disk spine with a maximum surface brightness contrast of 7.64 mag arcsec$^{-2}$. The polarized flux has a minimum near the star $<0.2"$ because no or only little polarization is produced by forward or backward scattering in the disk section lying in front of or behind the star. The data are modeled as a circular dust belt with an inclination $i=85(\pm 1.5)^\circ$ and a radius between $r_0$ = 60 AU and 90 AU. The radial density dependence is described by $(r/r_0)^{\alpha}$ with a steep power law index $\alpha=5$ inside $r_0$ and a more shallow index $\alpha=-2.5$ outside $r_0$. The scattering asymmetry factor lies between $g$ = 0.2 and 0.6 adopting a scattering angle-dependence for the fractional polarization as for Rayleigh scattering. Our data are qualitatively very similar to the case of AU Mic and they confirm that edge-on debris disks have a polarization minimum at a position near the star and a maximum near the projected separation of the main debris belt. The comparison of the polarized flux contrast ratio $F_{pol}/F_{\ast}$ with the fractional infrared excess provides strong constraints on the scattering albedo of the dust.

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.

R Aqr is a symbiotic binary system consisting of a mira variable, a hot companion with a spectacular jet outflow, and an extended emission line nebula. We have used R Aqr as test target for the visual camera subsystem ZIMPOL, which is part of the new extreme adaptive optics (AO) instrument SPHERE at the Very Large Telescope (VLT). We compare our observations with data from the Hubble Space Telescope (HST) and illustrate the complementarity of the two instruments. We determine from the Halpha emission the position, size, geometric structure, and line fluxes of the jet source and the clouds in the innermost region (<2") of R Aqr and determine Halpha emissivities mean density, mass, recombination time scale, and other cloud parameters. Our data resolve for the first time the R Aqr binary and we measure for the jet source a relative position 46+/-1 mas West of the mira. The central jet source is the strongest Halpha component. North east and south west from the central source there are many clouds with very diverse structures. We see in the SW a string of bright clouds arranged in a zig-zag pattern and, further out, more extended bubbles. In the N and NE we see a bright, very elongated filamentary structure and faint perpendicular "wisps" further out. Some jet clouds are also detected in the ZIMPOL [OI] and He I filters, as well as in the HST line filters for Halpha, [OIII], [NII], and [OI]. We determine jet cloud parameters and find a very well defined anti-correlation between cloud density and distance to the central binary. Future Halpha observations will provide the orientation of the orbital plane of the binary and allow detailed hydrodynamical investigations of this jet outflow and its interaction with the wind of the red giant companion.

HR4796A is surrounded by a debris disc, observed in scattered light as an inclined ring. Past observations raised several questions. First, a strong brightness asymmetry detected in polarized reflected light recently challenged our understanding of scattering by the dust particles in this system. Secondly, the morphology of the ring strongly suggests the presence of planets, although no planets have been detected to date. We obtained high-angular resolution coronagraphic images of the circumstellar environment around HR4796A with VLT/SPHERE during the commissioning of the instrument in May 2014 and during guaranteed-time observations in February 2015. The observations reveal for the first time the entire ring of dust, including the semi-minor axis that was previously hidden either behind the coronagraphic spot or in the speckle noise. We determine empirically the scattering phase function of the dust in the H band from 13.6deg to 166.6deg. It shows a prominent peak of forward scattering, never detected before, for scattering angles below 30deg. We analyse the reflectance spectra of the disc from the 0.95 to 1.6 microns, confirming the red colour of the dust, and derive detection limits on the presence of planetary mass objects. We confirm which side of the disc is inclined towards the Earth. The analysis of the phase function suggests that the dust population is dominated by particles much larger than the observation wavelength, of about 20 microns. Compact Mie grains of this size are incompatible with the spectral energy distribution of the disc, however the observed rise in scattering efficiency beyond 50deg points towards aggregates which could reconcile both observables. We do not detect companions orbiting the star but our high-contrast observations provide the most stringent constraints yet on the presence of planets responsible for the morphology of the dust.

The development of adaptive optics (AO) played a major role in modern astronomy over the last three decades. By compensating for the atmospheric turbulence, these systems enable to reach the diffraction limit on large telescopes. In this review, we will focus on high contrast applications of adaptive optics, namely, imaging the close vicinity of bright stellar objects and revealing regions otherwise hidden within the turbulent halo of the atmosphere to look for objects with a contrast ratio lower than 10^-4 with respect to the central star. Such high-contrast AO-corrected observations have led to fundamental results in our current understanding of planetary formation and evolution as well as stellar evolution. AO systems equipped three generations of instruments, from the first pioneering experiments in the nineties, to the first wave of instruments on 8m-class telescopes in the years 2000, and finally to the extreme AO systems that have recently started operations. Along with high-contrast techniques, AO enables to reveal the circumstellar environment: massive protoplanetary disks featuring spiral arms, gaps or other asymmetries hinting at on-going planet formation, young giant planets shining in thermal emission, or tenuous debris disks and micron-sized dust leftover from collisions in massive asteroid-belt analogs. After introducing the science case and technical requirements, we will review the architecture of standard and extreme AO systems, before presenting a few selected science highlights obtained with recent AO instruments.

The debris disk of HIP73145 has been detected in scattered light in the near-IR, and at far-IR wavelengths before, but no substructure has been seen so far. Detection of such substructures in combination with detailed modeling can hint at the presence of perturbing planetary bodies, or reveal other mechanisms acting to replenish gas and dust reservoirs and forming structures such as spirals or rings. We obtained multiwavelength images with SPHERE in the near-IR in the H2 and H3 bands with the IRDIS camera and a 0.95-1.35 micron spectral cube with the IFS. Data were acquired in pupil-tracking mode, thus allowing for angular differential imaging. The SPHERE standard suite of angular differential imaging algorithms was applied. ALMA Band 6 observations complement the SPHERE data. We detect a bright ring of scattered light plus more structures inside, at least one of them forming a secondary, concentric ring with the first. This is the first detection of this disk in total-intensity scattered light. A second object is detected in the field at high contrast but concluded to be a background star. Forward modeling yields information on the primary parameters of the disk and confirms that the detected substructures are not due to the data analysis approach, which sometimes leads to spurious structures. We detect a series of concentric rings in the disk around HIP73145. This is one of the rare cases where multiple components are necessary to fit the SED and are also detected in scattered light. The presence of such ring structures somewhat questions the nature of the object as a pure debris disk, but the gas and dust content would presumably offer sufficient explanations for such structures to form.

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.

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.

The recent discovery of an earth-like planet around Proxima Centauri has drawn much attention to this star and its environment. We performed a series of observations of Proxima Centauri using SPHERE, the planet finder instrument installed at the ESO Very Large Telescope UT3, using its near infrared modules, IRDIS and IFS. No planet was directly detected but we set upper limits on the mass up to 7 au exploiting the AMES-COND models. Our IFS observations reveal that no planet more massive than ~6-7 M Jup can be present within 1 au. The dual band imaging camera IRDIS also enables us to probe larger separations than the other techniques like the radial velocity or astrometry. We obtained mass limits of the order of 4 M Jup at separations of 2 au or larger representing the most stringent mass limits at separations larger than 5 au available at the moment. We also did an attempt to estimate the radius of possible planets around Proxima using the reflected light. Since the residual noise for this observations are dominated by photon noise and thermal background, longer exposures in good observing conditions could further improve the achievable contrast limit.

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.

Whether BP Piscium (BP Psc) is either a pre-main sequence T Tauri star at d ~ 80 pc, or a post-main sequence G giant at d ~ 300 pc is still not clear. As a first-ascent giant, it is the first to be observed with a molecular and dust disk. Alternatively, BP Psc would be among the nearest T Tauri stars with a protoplanetary disk (PPD). We investigate whether the disk geometry resembles typical PPDs, by comparing polarimetric images with radiative transfer models. Our VLT/SPHERE/ZIMPOL observations allow us to perform Polarimetric Differential Imaging; Reference Star Differential Imaging; and Richardson-Lucy deconvolution. We present the first visible light polarization and intensity images of the disk of BP Psc. Our deconvolution confirms the disk shape as detected before, mainly showing the southern side of the disk. In polarized intensity the disk is imaged at larger detail and also shows the northern side, giving it the typical shape of high inclination flared disks. We explain the observed disk features by retrieving the large-scale geometry with MCMax radiative transfer modeling, which yields a strongly flared model, atypical for disks of T Tauri stars.

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 the current results of the astrometric characterization of the VLT planet finder SPHERE over 2 years of on-sky operations. We first describe the criteria for the selection of the astrometric fields used for calibrating the science data: binaries, multiple systems, and stellar clusters. The analysis includes measurements of the pixel scale and the position angle with respect to the North for both near-infrared subsystems, the camera IRDIS and the integral field spectrometer IFS, as well as the distortion for the IRDIS camera. The IRDIS distortion is shown to be dominated by an anamorphism of 0.60+/-0.02% between the horizontal and vertical directions of the detector, i.e. 6 mas at 1". The anamorphism is produced by the cylindrical mirrors in the common path structure hence common to all three SPHERE science subsystems (IRDIS, IFS, and ZIMPOL), except for the relative orientation of their field of view. The current estimates of the pixel scale and North angle for IRDIS are 12.255+/-0.009 milliarcseconds/pixel for H2 coronagraphic images and -1.75+/-0.08 deg. Analyses of the IFS data indicate a pixel scale of 7.46+/-0.02 milliarcseconds/pixel and a North angle of -102.18+/-0.13 deg. We finally discuss plans for providing astrometric calibration to the SPHERE users outside the instrument consortium.

In 2015, we initiated a survey of Scorpius-Centaurus A-F stars that are predicted to host warm-inner and cold-outer belts of debris similar to the case of the system HR~8799. The survey aims to resolve the disks and detect planets responsible for the disk morphology. In this paper, we study the F-type star HIP~67497 and present a first-order modelisation of the disk in order to derive its main properties. We used the near-infrared integral field spectrograph (IFS) and dual-band imager IRDIS of VLT/SPHERE to obtain angular-differential imaging observations of the circumstellar environnement of HIP~67497. We removed the stellar halo with PCA and TLOCI algorithms. We modeled the disk emission with the GRaTeR code. We resolve a ring-like structure that extends up to $\sim$450 mas ($\sim$50 au) from the star in the IRDIS and IFS data. It is best reproduced by models of a non-eccentric ring with an inclination of $80\pm1^{\circ}$, a position angle of $-93\pm1^{\circ}$, and a semi-major axis of $59\pm3$ au. We also detect an additional, but fainter, arc-like structure with a larger extension (0.65 arcsec) South of the ring that we model as a second belt of debris at $\sim$130 au. We detect 10 candidate companions at separations $\geq$1''. We estimate the mass of putative perturbers responsible for the disk morphology and compare it to our detection limits. Additional data are needed to find those perturbers, and to relate our images to large-scale structures seen with HST/STIS.

We studied the well known circumstellar disk around the Herbig Ae/Be star HD 97048 with high angular resolution to reveal undetected structures in the disk, which may be indicative of disk evolutionary processes such as planet formation. We used the IRDIS near-IR subsystem of the extreme adaptive optics imager SPHERE at the ESO/VLT to study the scattered light from the circumstellar disk via high resolution polarimetry and angular differential imaging. We imaged the disk in unprecedented detail and revealed four ring-like brightness enhancements and corresponding gaps in the scattered light from the disk surface with radii between 39 au and 341 au. We derived the inclination and position angle as well as the height of the scattering surface of the disk from our observational data. We found that the surface height profile can be described by a single power law up to a separation ~270 au. Using the surface height profile we measured the scattering phase function of the disk and found that it is well consistent with theoretical models of compact dust aggregates. We discuss the origin of the detected features and find that low mass (< 1 M_Jup) nascent planets are a possible explanation.

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.

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.

Debris disks offer valuable insights into the latest stages of circumstellar disk evolution, and can possibly help us to trace the outcomes of planetary formation processes. In the age range 10 to 100\,Myr, most of the gas is expected to have been removed from the system, giant planets (if any) must have already been formed, and the formation of terrestrial planets may be on-going. Pluto-sized planetesimals, and their debris released in a collisional cascade, are under their mutual gravitational influence, which may result into non-axisymmetric structures in the debris disk. High angular resolution observations are required to investigate these effects and constrain the dynamical evolution of debris disks. Furthermore, multi-wavelength observations can provide information about the dust dynamics by probing different grain sizes. Here we present new VLT/SPHERE and ALMA observations of the debris disk around the 40\,Myr-old solar-type star HD\,61005. We resolve the disk at unprecedented resolution both in the near-infrared (in scattered and polarized light) and at millimeter wavelengths. Thanks to the new observations, we propose a solution for both the radial and azimuthal distribution of the dust grains in the debris disk. We find that the disk has a moderate eccentricity ($e \sim 0.1$) and that the dust density is two times larger at the pericenter compared to the apocenter. With no giant planets detected in our observations, we investigate alternative explanations besides planet-disk interactions to interpret the inferred disk morphology. We postulate that the morphology of the disk could be the consequence of a massive collision between $\sim$\,1000\,km-sized bodies at $\sim$\,61\,au. If this interpretation holds, it would put stringent constraints on the formation of massive planetesimals at large distances from the star.

Transition disks correspond to a short stage between the young protoplanetary phase and older debris phase. Along this evolutionary sequence, the gas component disappears leaving room for a dust-dominated environment where already-formed planets signpost their gravitational perturbations. We endeavor to study the very inner region of the well-known and complex debris, but still gas-rich disk, around HD 141569A using the exquisite high-contrast capability of SPHERE at the VLT. Recent near-infrared (IR) images suggest a relatively depleted cavity within ~200 au, while former mid-IR data indicate the presence of dust at separations shorter than ~100 au. We obtained multi-wavelength images in the near-IR in J, H2, H3 and Ks bands with the IRDIS camera and a 0.95-1.35 micrometers spectral data cube with the IFS. Data were acquired in pupil-tracking mode, thus allowing for angular differential imaging. We discovered several new structures inside 1", of which the most prominent is a bright ring with sharp edges (semi-major axis: 0.4") featuring a strong north-south brightness asymmetry. Other faint structures are also detected from 0.4" to 1" in the form of concentric ringlets and at least one spiral arm. Finally, the VISIR data at 8.6 micrometers suggests the presence of an additional dust population closer in. Besides, we do not detect companions more massive than 1-3 mass of Jupiter. The performance of SPHERE allows us to resolve the extended dust component, which was previously detected at thermal and visible wavelengths, into very complex patterns with strong asymmetries ; the nature of these asymmetries remains to be understood. Scenarios involving shepherding by planets or dust-gas interactions will have to be tested against these observations.

Context. It has recently been proposed that the surface composition of icy main-belt asteroids (B-,C-,Cb-,Cg-,P-,and D-types) may be consistent with that of chondritic porous interplanetary dust particles (CPIDPs). Aims. In the light of this new association, we re-examine the surface composition of a sample of asteroids belonging to the Themis family in order to place new constraints on the formation and evolution of its parent body. Methods. We acquired NIR spectral data for 15 members of the Themis family and complemented this dataset with existing spectra in the visible and mid-infrared ranges to perform a thorough analysis of the composition of the family. Assuming end-member minerals and particle sizes (<2\mum) similar to those found in CPIDPs, we used a radiative transfer code adapted for light scattering by small particles to model the spectral properties of these asteroids. Results. Our best-matching models indicate that most objects in our sample possess a surface composition that is consistent with the composition of CP IDPs.We find ultra-fine grained Fe-bearing olivine glasses to be among the dominant constituents. We further detect the presence of minor fractions of Mg-rich crystalline silicates. The few unsuccessfully matched asteroids may indicate the presence of interlopers in the family or objects sampling a distinct compositional layer of the parent body. Conclusions. The composition inferred for the Themis family members suggests that the parent body accreted from a mixture of ice and anhydrous silicates (mainly amorphous) and subsequently underwent limited heating. By comparison with existing thermal models that assume a 400km diameter progenitor, the accretion process of the Themis parent body must have occurred relatively late (>4Myr after CAIs) so that only moderate internal heating occurred in its interior, preventing aqueous alteration of the outer shell.

[Abridged] Context. The young systems PZ Tel and HD 1160, hosting known low-mass companions, were observed during the commissioning of the new planet finder SPHERE with several imaging and spectroscopic modes. Aims. We aim to refine the physical properties and architecture of both systems. Methods. We use SPHERE commissioning data and REM observations, as well as literature and unpublished data from VLT/SINFONI, VLT/NaCo, Gemini/NICI, and Keck/NIRC2. Results. We derive new photometry and confirm the nearly daily photometric variability of PZ Tel A. Using literature data spanning 38 yr, we show that the star also exhibits a long-term variability trend. The 0.63-3.8 mic SED of PZ Tel B allows us to revise its properties: spectral type M7+/-1, Teff=2700+/-100 K, log(g)<4.5 dex, log(L/L_Sun)=-2.51+/-0.10 dex, and mass 38-72 MJ. The 1-3.8 mic SED of HD 1160 B suggests a massive brown dwarf or a low-mass star with spectral type M5.5-7.0, Teff=3000+/-100 K, [M/H]=-0.5-0.0 dex, log(L/L_Sun)=-2.81+/-0.10 dex, and mass 39-168 MJ. We confirm the deceleration and high eccentricity (e>0.66) of PZ Tel B. For e<0.9, the inclination, longitude of the ascending node, and time of periastron passage are well constrained. The system is seen close to an edge-on geometry. We reject other brown dwarf candidates outside 0.25" for both systems, and massive giant planets (>4 MJ) outside 0.5" for the PZ Tel system. We also show that K1-K2 color can be used with YJH low-resolution spectra to identify young L-type companions, provided high photometric accuracy (<0.05 mag) is achieved. Conclusions. SPHERE opens new horizons in the study of young brown dwarfs and giant exoplanets thanks to high-contrast imaging capabilities at optical and near-infrared wavelengths, as well as high signal-to-noise spectroscopy in the near-infrared from low (R~30-50) to medium resolutions (R~350).

The system of four planets around HR8799 offers a unique opportunity to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R~30) spectra of HR8799 d and e with the SPHERE instrument (paper III). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these two planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2{\sigma}) the whole set of spectrophotometric datapoints available for HR8799 d and e for Teff = 1200 K, log g in the range 3.0-4.5, and M/H=+0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. We conclude that an enhanced content in dust and decreased CIA of H2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $\mu$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $\mu$m, 1.667 $\mu$m), $K1K2$ (2.110 $\mu$m, 2.251 $\mu$m), and broadband $J$ (1.245 $\mu$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $\mu$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Lo\`eve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

GJ758 B is a brown dwarf companion to a nearby (15.76 pc) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (~600K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the VLT. The data was obtained in Y-, J-, H-, and Ks-bands with the dual-band imaging (DBI) mode of IRDIS, providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new 8 photometric points for an extended comparison of GJ758 B with empirical objects and 4 families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison object can accurately represent the observed near-IR fluxes of GJ758 B. From comparison to atmospheric models, we attribute a Teff = 600K $\pm$ 100K, but we find that no atmospheric model can adequately fit all the fluxes of GJ758 B. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches, Least-Square Monte Carlo and Markov Chain Monte Carlo. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU. [abridged]

Context. Some of transition disks show asymmetric structures in thermal sub-millimetre emission and optical scattered light. These structures can be the result of planet(s) or companions embedded in the disk. Aims. We aim to detect and analyse the scattered light of the transition disk J160421.7-213028, identify disk structures, and compare the results with previous observations of this disk at other wavelengths. Methods. We obtained and analysed new polarised intensity observations of the transition disk J160421.7-213028 with VLT/SPHERE using the visible light instrument ZIMPOL at $R'$ band (0.626$\mu$m). We probe the disk gap down to a radius of confidence of 0.1'' (${\sim}15$ AU at 145 pc). We interpret the results in the context of dust evolution when planets interact with the parental disk. Results. We observe a gap from 0.1 to 0.3'' (${\sim}15$ to 40AU) and a bright annulus as previously detected by HiCIAO $H$ band observations at $1.65\mu$m. The radial width of the annulus is around $40$AU, and its centre is at ${\sim}61$AU from the central star. The peak of the reflected light at 0.626$\mu$m is located 20 AU inward of the cavity detected in the submillimetre. In addition, we detect a dip at a position angle of ${\sim}46.2 \pm 5.4^\circ$. A dip was also detected with HiCIAO, but located at ${\sim}85^\circ$. If the dip observed with HiCIAO is the same, this suggests an average dip rotation of ${\sim}12^\circ/$year, which is inconsistent with the local Keplerian angular velocity of $\sim$0.8$^\circ$/yr at $\sim$61AU. Conclusions. The spatial discrepancy in the radial emission in J160421.7-213028 at different wavelengths is consistent with dust filtration at the outer edge of a gap carved by a massive planet. The dip rotation can be interpreted by fast variability of the inner disk and/or the presence of a warp or circumplanetary material of a planet at ${\sim}9.6$ AU.

HD~106906AB is so far the only young binary system around which a planet has been imaged and a debris disk evidenced thanks to a strong IR excess. As such, it represents a unique opportunity to study the dynamics of young planetary systems. We aim at further investigating the close (tens of au scales) environment of the HD~106906AB system. We used the extreme AO fed, high contrast imager SPHERE recently installed on the VLT to observe HD~106906. Both the IRDIS imager and the Integral Field Spectrometer were used. We discovered a very inclined, ring-like disk at a distance of 65~au from the star. The disk shows a strong brightness asymmetry with respect to its semi-major axis. It shows a smooth outer edge, compatible with ejection of small grains by the stellar radiation pressure. We show furthermore that the planet's projected position is significantly above the disk's PA. Given the determined disk inclination, it is not excluded though that the planet could still orbit within the disk plane if at a large separation (2000--3000 au). We identified several additional point sources in the SPHERE/IRDIS field-of-view, that appear to be background objects. We compare this system with other debris disks sharing similarities, and we briefly discuss the present results in the framework of dynamical evolution.

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.

CONTEXT. Little is known about the planetary systems around single white dwarfs although there is strong evidence that they do exist. AIMS. We performed a pilot study with the extreme-AO system on the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) on the Very Large Telescopes (VLT) to look for giant planets around a young white dwarf, GD 50. METHODS. We were awarded science verification time on the new ESO instrument SPHERE. Observations were made with the InfraRed Dual-band Imager and Spectrograph in classical imaging mode in H band. RESULTS. Despite the faintness of the target (14.2 mag in R band), the AO loop was closed and a strehl of 37\% was reached in H band. No objects were detected around GD 50. We achieved a 5-sigma contrast of 6.2, 8.0 and 8.25 mags at 0{\farcs}2, 0{\farcs}4 and 0{\farcs}6 and beyond, respectively. We exclude any substellar objects more massive than 4.0 M$_\textrm{J}$ at 6.2 AU, 2.9 M$_\textrm{J}$ at 12.4 AU and 2.8 M$_\textrm{J}$ at 18.6 AU and beyond. This rivals the previous upper limit set by Spitzer. We further show that SPHERE is the most promising instrument available to search for close-in substellar objects around nearby white dwarfs.

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

HR 4796A is surrounded by a well-structured and very bright circumstellar disc shaped like an annulus with many interesting features: very sharp inner and outer edges, brightness asymmetries, centre offset, and suspected distortions in the ring. We aim to constrain the properties of the dust surrounding the star HR 4796A, in particular the grain size and composition. We also want to confirm and refine the morphological parameters derived from previous scattered light observations, and reveal the dust spatial extent in regions unexplored so far due to their proximity to the star. We have obtained new images in polarised light of the binary system HR 4796A and B in the Ks and Lp band with the NaCo instrument at the Very Large Telescope (VLT). In addition, we revisit two archival data sets obtained in the Lp band with that same instrument and at 2.2 microns with the NICMOS instrument on the Hubble Space Telescope. We analyse these observations with simulations using the radiative transfer code MCFOST to investigate the dust properties. We explore a grid of models with various dust compositions and sizes in a Bayesian approach. We detect the disc in polarised light in the Ks band and reveal for the first time the innermost regions down to $0.3"$ along the semi-minor axis. We measure a polarised fraction of $29\%\pm8\%$ in the two disc ansae, with a maximum occurring more than $13^\circ$ westwards from the ansae. A very pronounced brightness asymmetry between the north-west and south-east side is detected. This contradicts the asymmetry previously reported in all images of the disc in unpolarised light at wavelengths smaller than or equal to 2.2 microns and is inconsistent with the predicted scattered light from spherical grains using the Mie theory. Our modelling suggests the north-west side is most likely inclined towards the Earth, contrary to previous conclusions.

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

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.

Imaging of debris disks has found evidence for both eccentric and offset disks. One hypothesis is that these provide evidence for massive perturbers that sculpt the observed structures. One such disk was recently observed in the far-IR by the Herschel Space Observatory around $\zeta^2$ Ret. In contrast with previously reported systems, the disk is significantly eccentric, and the system is Gyr-old. We aim to investigate the long-term evolution of eccentric structures in debris disks caused by a perturber on an eccentric orbit. Both analytical predictions and numerical N-body simulations are used to investigate the observable structures that could be produced by eccentric perturbers. The long-term evolution of the disk geometry is examined, with particular application to the $\zeta^2$ Ret system. In addition, synthetic images of the disk are produced for comparison with Herschel observations. We show that an eccentric companion can produce both the observed offsets and eccentric disks. Such effects are not immediate and we characterise the timescale required for the disk to develop to an eccentric state. For the case of $\zeta^2$ Ret, we place limits on the mass and orbit of the companion required to produce the observations. Synthetic images show that the pattern observed around $\zeta^2$ Ret can be produced by an eccentric disk seen close to edge-on, and allow us to bring additional constraints on the disk parameters of our model (disk flux, extent). We determine that eccentric planets or stellar companions can induce long-lived eccentric structures in debris disks. Observations of such eccentric structures provide potential evidence of the presence of such a companion in a planetary system. We consider the example of $\zeta^2$ Ret, whose observed eccentric disk can be explained by a distant companion at tens of AU, on an eccentric orbit ($e_p\gtrsim 0.3$).

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.

Over the past five years, radial-velocity and transit techniques have revealed a new population of Earth-like planets with masses of a few Earth masses. Their very close orbit around their host star requires an exquisite inner working angle to be detected in direct imaging and sets a challenge for direct imagers that work in the visible range, like SPHERE / ZIMPOL. Among all known exoplanets with less than twenty-five Earth-masses we first predict the best candidate for direct imaging. Our primary objective is then to provide the best instrument setup and observing strategy for detecting such a peculiar object with ZIMPOL. As a second step, we aim at predicting its detectivity. Using exoplanet properties constrained by radial velocity measurements, polarimetric models and the diffraction propagation code CAOS, we estimate the detection sensitivity of ZIMPOL for such a planet in different observing modes of the instrument. We show how observing strategies can be optimized to yield the best detection performance on a specific target. In our current knowledge of exoplanetary systems, alpha Centauri B b is the most promising target with less than twenty-five Earth-masses for ZIMPOL. With a gaseous Rayleigh-scattering atmosphere and favorable inclinations, the planet could be detected in about four hours of observing time, using the four-quadrant phase-mask coronograph in the I band. However, if alpha Centauri B b should display unfavorable polarimetric and reflective properties similar to that of our Moon, it is around 50 times fainter than the best sensitivity of ZIMPOL. Alpha Centauri B is a primary target for SPHERE. Dedicated deep observations specifically targeting the radial velocity-detected planet can lead to a detection if the polarimetric properties of the planet are favorable.

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.

Context. It has long been suggested that circumstellar disks surrounding young stars may be the signposts of planets, and still more since the recent discoveries of embedded substellar companions. The planet-disk interaction may create, according to models, large structures, gaps, rings or spirals, in the disk. In that sense, the Herbig star HD 142527 is particularly compelling as, its massive disk displays intriguing asymmetries that suggest the existence of a dynamical peturber of unknown nature. Aims. Our goal was to obtain deep thermal images of the close circumstellar environment of HD 142527 to re-image the reported close-in structures (cavity, spiral arms) of the disk and to search for stellar and substellar companions that could be connected to their presence. Results. The circumstellar environment of HD 142527 is revealed at an unprecedented spatial resolution down to the sub arcsecond level for the first time at 3.8 microns. Our images reveal important radial and azimuthal asymmetries which invalidate an elliptical shape for the disk as previously proposed. It rather suggests a bright inhomogeneous spiral arm plus various fainter spiral arms. We also confirm an inner cavity down to 30 AU and two important dips at position angles of 0 and 135 deg. The detection performance in angular differential imaging enables the exploration of the planetary mass regime for projected physical separations as close as 40 AU. The use of our detection map together with Monte Carlo simulations sets stringent constraints on the presence of planetary mass, brown dwarf or stellar companions as a function of the semi-major axis. They severely constrain the presence of massive giant planets with semi-major axis beyond 50AU, i.e. probably within the large disk's cavity that radially extends up to 145 AU or even further outside.

Direct imaging of circumstellar disks requires high-contrast and high-resolution techniques. The angular differential imaging (ADI) technique is one of them, initially developed for point-like sources but now increasingly applied to extended objects. This new field of application raises many questions because the disk images reduced with ADI depend strongly on the amplitude of field rotation and the ADI data reduction strategy. Both of them directly affect the disk observable properties. Our aim is to characterize the applicability and biases of some ADI data reduction strategies for different disk morphologies. A particular emphasis is placed on parameters mostly used for disks: their surface brightness, their width for a ring, and local features such as gaps or asymmetries. We first present a general method for predicting and quantifying those biases. In a second step we illustrate them for some widely used ADI algorithms applied to typical debris disk morphologies: inclined rings with various inner/outer slopes and width. Last, our aim is also to propose improvements of classical ADI to limit the biases on extended objects. Simulated disks seen under various observing conditions were used to reduce ADI data and quantify the resulting biases. These conclusions complements previous results from NaCo L' real-disk images of HR4796A. ADI induces flux losses on disks. This makes this technique appropriate only for low- to medium-inclination disks. A theoretical criterion is derived to predict the amount of flux loss for a given disk morphology, and quantitative estimates of the biases are given in some specific configurations. These biases alter the disk observable properties, such as the slopes of their surface brightness or their radial/azimuthal extent. Additionally, this work demonstrates that ADI can very easily create artificial features without involving astrophysical processes.

Small-angle coronagraphy is technically and scientifically appealing because it enables the use of smaller telescopes, allows covering wider wavelength ranges, and potentially increases the yield and completeness of circumstellar environment - exoplanets and disks - detection and characterization campaigns. However, opening up this new parameter space is challenging. Here we will review the four posts of high contrast imaging and their intricate interactions at very small angles (within the first 4 resolution elements from the star). The four posts are: choice of coronagraph, optimized wavefront control, observing strategy, and post-processing methods. After detailing each of the four foundations, we will present the lessons learned from the 10+ years of operations of zeroth and first-generation adaptive optics systems. We will then tentatively show how informative the current integration of second-generation adaptive optics system is, and which lessons can already be drawn from this fresh experience. Then, we will review the current state of the art, by presenting world record contrasts obtained in the framework of technological demonstrations for space-based exoplanet imaging and characterization mission concepts. Finally, we will conclude by emphasizing the importance of the cross-breeding between techniques developed for both ground-based and space-based projects, which is relevant for future high contrast imaging instruments and facilities in space or on the ground.

HR4796 is a young, early A-type star harbouring a well structured debris disk, shaped as a ring with sharp inner edges. It forms with the M-type star HR4796B a binary system, with a proj. sep. ~560 AU. Our aim is to explore the surroundings of HR4796A and B, both in terms of extended or point-like structures. Adaptive optics images at L'-band were obtained with NaCo in Angular Differential Mode and with Sparse Aperture Masking (SAM). We analyse the data as well as the artefacts that can be produced by ADI reduction on an extended structure with a shape similar to that of HR4796A dust ring. We determine constraints on the presence of companions using SAM and ADI on HR4796A, and ADI on HR4796B. We also performed dynamical simulations of a disk of planetesimals and dust produced by collisions, perturbed by a planet located close to the disk outer edge. The disk ring around HR4796A is well resolved. We highlight the potential effects of ADI reduction of the observed disk shape and surface brightness distribution, and side-to-side asymmetries. No planet is detected around the star, with masses as low as 3.5 M_Jup at 0.5" (58 AU) and less than 3 M_Jup in the 0.8-1" range along the semi-major axis. We exclude massive brown dwarfs at separations as close as 60 mas (4.5 AU) from the star thanks to SAM data. The detection limits obtained allow us to exclude a possible close companion to HR4796A as the origin of the offset of the ring center with respect to the star; they also allow to put interesting constraints on the (mass, separation) of any planet possibly responsible for the inner disk steep edge. Using detailed dynamical simulations, we show that a giant planet orbiting outside the ring could sharpen the disk outer edge and reproduce the STIS images published by Schneider et al. (2009).

Direct imaging of circumstellar disks at high angular resolution is mandatory to provide morphological information that bring constraints on their properties, in particular the spatial distribution of dust. New techniques combining observing strategy and data processing now allow very high contrast imaging with 8-m class ground-based telescopes (10^-4 to 10^-5 at ~1") and complement space telescopes while improving angular resolution at near infrared wavelengths. We carried out a program at the VLT with NACO to image known debris disks with higher angular resolution in the near IR than ever before in order to study morphological properties and ultimately to detect signpost of planets. The observing method makes use of advanced techniques: Adaptive Optics, Coronagraphy and Differential Imaging, a combination designed to directly image exoplanets with the upcoming generation of "planet finders" like GPI (Gemini Planet Imager) and SPHERE (Spectro-Polarimetric High contrast Exoplanet REsearch). Applied to extended objects like circumstellar disks, the method is still successful but produces significant biases in terms of photometry and morphology. We developed a new model-matching procedure to correct for these biases and hence to bring constraints on the morphology of debris disks. From our program, we present new images of the disk around the star HD 32297 obtained in the H (1.6mic) and Ks (2.2mic) bands with an unprecedented angular resolution (~65 mas). The images show an inclined thin disk detected at separations larger than 0.5-0.6". The modeling stage confirms a very high inclination (i=88{\deg}) and the presence of an inner cavity inside r_0~110AU. We also found that the spine (line of maximum intensity along the midplane) of the disk is curved and we attributed this feature to a large anisotropic scattering factor (g~0.5, valid for an non-edge on disk). Abridged ...

Context. We detected in 2009 a giant, close-by planet orbiting {\beta} Pic, a young star surrounded with a disk, extensively studied for more than 20 years. We showed that if located on an inclined orbit, the planet could explain several peculiarities of {\beta} Pictoris system. However, the available data did not permit to measure the inclination of {\beta} Pic b with respect to the disk, and in particular to establish in which component of the disk - the main, extended disk or the inner inclined component/disk-, the planet was located. Comparison between the observed planet position and the disk orientation measured on previous imaging data was not an option because of potential biases in the measurements. Aims. Our aim is to measure precisely the planet location with respect to the dust disk using a single high resolution image, and correcting for systematics or errors that degrades the precision of the disk and planet relative position measurements. Methods. We gathered new NaCo data at Ks band, with a set-up optimized to derive simultaneously the orientation(s) of the disk(s) and that of the planet. Results. We show that the projected position of {\beta} Pic b is above the midplane of the main disk. With the current data and knowledge on the system, this implies that {\beta} Pic b cannot be located in the main disk. The data rather suggest the planet being located in the inclined component.