Direct imaging of Earth-like planets from space requires dedicated observatories, combining large segmented apertures with instruments and techniques such as coronagraphs, wavefront sensors, and wavefront control in order to reach the high contrast of 10^10 that is required. The complexity of these systems would be increased by the segmentation of the primary mirror, which allows for the larger diameters necessary to image Earth-like planets but also introduces specific patterns in the image due to the pupil shape and segmentation and making high-contrast imaging more challenging. Among these defects, the phasing errors of the primary mirror are a strong limitation to the performance. In this paper, we focus on the wavefront sensing of segment phasing errors for a high-contrast system, using the COronagraphic Focal plane wave-Front Estimation for Exoplanet detection (COFFEE) technique. We implemented and tested COFFEE on the High-contrast imaging for Complex Aperture Telescopes (HiCAT) testbed, in a configuration without any coronagraph and with a classical Lyot coronagraph, to reconstruct errors applied on a 37 segment mirror. We analysed the quality and limitations of the reconstructions. We demonstrate that COFFEE is able to estimate correctly the phasing errors of a segmented telescope for piston, tip, and tilt aberrations of typically 100nm RMS. We also identified the limitations of COFFEE for the reconstruction of low-order wavefront modes, which are highly filtered by the coronagraph. This is illustrated using two focal plane mask sizes on HiCAT. We discuss possible solutions, both in the hardware system and in the COFFEE optimizer, to mitigate these issues.

Extremely Large Telescopes have overwhelmingly opted for the Pyramid wavefront sensor (PyWFS) over the more widely used Shack-Hartmann WaveFront Sensor (SHWFS) to perform their Single Conjugate Adaptive Optics (SCAO) mode. The PyWFS, a sensor based on Fourier filtering, has proven to be highly successful in many astronomy applications. However, it exhibits non-linearity behaviors that lead to a reduction of its sensitivity when working with non-zero residual wavefronts. This so-called Optical Gains (OG) effect, degrades the close loop performance of SCAO systems and prevents accurate correction of Non-Common Path Aberrations (NCPA). In this paper, we aim at computing the OG using a fast and agile strategy in order to control the PyWFS measurements in adaptive optics closed loop systems. Using a novel theoretical description of the PyFWS, which is based on a convolutional model, we are able to analytically predict the behavior of the PyWFS in closed-loop operation. This model enables us to explore the impact of residual wavefront error on particular aspects such as sensitivity and associated OG. The proposed method relies on the knowledge of the residual wavefront statistics and enables automatic estimation of the current OG. End-to-End numerical simulations are used to validate our predictions and test the relevance of our approach. We demonstrate, using on non-invasive strategy, that our method provides an accurate estimation of the OG. The model itself only requires AO telemetry data to derive statistical information on atmospheric turbulence. Furthermore, we show that by only using an estimation of the current Fried parameter r_0 and the basic system-level characteristics, OGs can be estimated with an accuracy of less than 10%. Finally, we highlight the importance of OG estimation in the case of NCPA compensation. The proposed method is applied to the PyWFS.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

We present observations of the young multiple system UX Tauri to look for circumstellar disks and for signs of dynamical interactions. We obtained SPHERE/IRDIS deep differential polarization images in the J and H bands. We also used ALMA archival CO data. Large extended spirals are well detected in scattered light coming out of the disk of UX Tau A. The southern spiral forms a bridge between UX Tau A and C. These spirals, including the bridge connecting the two stars, all have a CO (3-2) counterpart seen by ALMA. The disk of UX Tau C is detected in scattered light. It is much smaller than the disk of UX Tau A and has a major axis along a different position angle, suggesting a misalignment. We performed PHANTOM SPH hydrodynamical models to interpret the data. The scattered light spirals, CO emission spirals and velocity patterns of the rotating disks, and the compactness of the disk of UX Tau C all point to a scenario in which UX Tau A has been perturbed very recently (about 1000 years) by the close passage of UX Tau C.

High-contrast imaging (HCI) observations of exoplanets can be limited by the island effect (IE). On the current generation of telescopes, the IE becomes a severe problem when the ground wind speed is below a few meters per second. This is referred to as the low-wind effect (LWE). The LWE severely distorts the point spread function (PSF), significantly lowering the Strehl ratio and degrading the contrast. In this article, we aim to show that the focal-plane wavefront sensing (FPWFS) algorithm, Fast and Furious (F&F), can be used to measure and correct the IE/LWE. We deployed the algorithm on the SCExAO HCI instrument at the Subaru Telescope using the internal near-infrared camera in H-band. We tested F&F with the internal source, and it was deployed on-sky to test its performance with the full end-to-end system and atmospheric turbulence. The performance of the algorithm was evaluated by two metrics based on the PSF quality: 1) the Strehl ratio approximation ($SRA$), and 2) variance of the normalized first Airy ring ($VAR$). Random LWE phase screens with a peak-to-valley wavefront error between 0.4 $\mu$m and 2 $\mu$m were all corrected to a $SRA$ $>$90\% and an $VAR\lessapprox0.05$. Furthermore, the on-sky results show that F&F is able to improve the PSF quality during very challenging atmospheric conditions (1.3-1.4'' seeing at 500 nm). Closed-loop tests show that F&F is able to improve the $VAR$ from 0.27 to 0.03 and therefore significantly improve the symmetry of the PSF. Simultaneous observations of the PSF in the {optical} ($\lambda = $ 750 nm, $\Delta \lambda =$ 50 nm) show that during these tests we were correcting aberrations common to the optical and NIR paths within SCExAO. Going forward, the algorithm is suitable for incorporation into observing modes, which will enable PSFs of higher quality and stability during science observations.

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.

Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

Current designs for all three extremely large telescopes show the overwhelming adoption of the pyramid wavefront sensor (P-WFS) as the WFS of choice for adaptive optics (AO) systems sensing on natural guide stars (NGS) or extended objects. The key advantages of the P-WFS over the Shack-Hartmann are known and are mainly provided by the improved sensitivity (fainter NGS) and reduced sensitivity to spatial aliasing. However, robustness and tolerances of the P-WFS for the ELTs are not currently well understood. In this paper, we present simulation results for the single-conjugate AO mode of HARMONI, a visible and near-infrared integral field spectrograph for the European Extremely Large Telescope. We first explore the wavefront sensing issues related to the telescope itself; namely the island effect (i.e. differential piston) and M1 segments phasing errors. We present mitigation strategies to the island effect and their performance. We then focus on some performance optimisation aspects of the AO design to explore the impact of the RTC latency and the optical gain issues, which will in particular affect the high-contrast mode of HARMONI. Finally, we investigate the influence of the quality of glass pyramid prism itself, and of optical aberrations on the final AO performance. By relaxing the tolerances on the fabrication of the prism, we are able to reduce hardware costs and simplify integration. We show the importance of calibration (i.e. updating the control matrix) to capture any displacement of the telescope pupil and rotation of the support structure for M4. We also show the importance of the number of pixels used for wavefront sensing to relax tolerances of the pyramid prism. Finally, we present a detailed optical design of the pyramid prism, central element of the P-WFS.

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.

Advanced AO systems will likely utilise Pyramid wave-front sensors (PWFS) over the traditional Shack-Hartmann sensor in the quest for increased sensitivity, peak performance and ultimate contrast. Here, we wish to bring knowledge and quantify the PWFS theoretical limits as a means to highlight its properties and use cases. We explore forward models for the PWFS in the spatial-frequency domain for they prove quite useful since a) they emanate directly from physical-optics (Fourier) diffraction theory; b) provide a straightforward path to meaningful error breakdowns, c) allow for reconstruction algorithms with $O (n\,log(n))$ complexity for large-scale systems and d) tie in seamlessly with decoupled (distributed) optimal predictive dynamic control for performance and contrast optimisation. All these aspects are dealt with here. We focus on recent analytical PWFS developments and demonstrate the performance using both analytic and end-to-end simulations. We anchor our estimates with observed on-sky contrast on existing systems and then show very good agreement between analytical and Monte-Carlo estimates for the PWFS. For a potential upgrade of existing high-contrast imagers on 10\,m-class telescopes with visible or near-infrared PWFS, we show under median conditions at Paranal a contrast improvement (limited by chromatic and scintillation effects) of 2x-5x by replacing the wave-front sensor alone at large separations close to the AO control radius where aliasing dominates, and factors in excess of 10x by coupling distributed control with the PWFS over most of the AO control region, from small separations starting with the Inner Working Angle of typically 1-2 $\lambda/D$ to the AO correction edge (here 20 $\lambda/D$).

The Low Wind Effect (LWE) refers to a phenomenon that occurs when the wind speed inside a telescope dome drops below $3$m/s creating a temperature gradient near the telescope spider. This produces phase discontinuities in the pupil plane that are not detected by traditional Adaptive Optics (AO) systems such as the pyramid wavefront sensor or the Shack-Hartmann. Considering the pupil as divided in 4 quadrants by regular spiders, the phase discontinuities correspond to piston, tip and tilt aberrations in each quadrant of the pupil. Uncorrected, it strongly decreases the ability of high contrast imaging instruments utilizing coronagraphy to detect exoplanets at small angular separations. Multiple focal plane wavefront sensors are currently being developed and tested on the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument at Subaru Telescope: Among them, the Zernike Asymmetric Pupil (ZAP) wavefront sensor already showed on-sky that it could measure the LWE induced aberrations in focal plane images. The Fast and Furious algorithm, using previous deformable mirror commands as temporal phase diversity, showed in simulations its efficiency to improve the wavefront quality in the presence of LWE. A Neural Network algorithm trained with SCExAO telemetry showed promising PSF prediction on-sky. The Linearized Analytic Phase Diversity (LAPD) algorithm is a solution for multi-aperture cophasing and is studied to correct for the LWE aberrations by considering the Subaru Telescope as a 4 sub-aperture instrument. We present the different algorithms, show the latest results and compare their implementation on SCExAO/SUBARU as real-time wavefront sensors for the LWE compensation.

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.

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.

Direct imaging of exo-Earths and search for life is one of the most exciting and challenging objectives for future space observatories. Segmented apertures in space will be required to reach the needed large diameters beyond the capabilities of current or planned launch vehicles. These apertures present additional challenges for high-contrast coronagraphy, not only in terms of static phasing but also in terms of their stability. The Pair-based Analytical model for Segmented Telescope Imaging from Space (PASTIS) was developed to model the effects of segment-level optical aberrations on the final image contrast. In this paper, we extend the original PASTIS propagation model from a purely analytical to a semi-analytical method, in which we substitute the use of analytical images with numerically simulated images. The inversion of this model yields a set of orthonormal modes that can be used to determine segment-level wavefront tolerances. We present results in the case of segment-level piston error applied to the baseline coronagraph design of LUVOIR A, with minimum and maximum wavefront error constraint between 56 pm and 290 pm per segment. The analysis is readily generalizable to other segment-level aberrations modes, and can also be expanded to establish stability tolerances for these missions.

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]

Quasi-static aberrations in coronagraphic systems are the ultimate limitation to the capabilities of exoplanet imagers both ground-based and space-based. These aberrations - which can be due to various causes such as optics alignment or moving optical parts during the observing sequence - create light residuals called speckles in the focal plane that might be mistaken for a planets. For ground-based instruments, the presence of residual turbulent wavefront errors due to partial adaptive optics correction causes an additional difficulty to the challenge of measuring aberrations in the presence of a coronagraph. In this paper, we present an extension of COFFEE, the coronagraphic phase diversity, to the estimation of quasi-static aberrations in the presence of adaptive optics-corrected residual turbulence. We perform realistic numerical simulations to assess the performance that can be expected on an instrument of the current generation. We perform the first experimental validation in the laboratory which demonstrates that quasistatic aberrations can be corrected during the observations by means of coronagraphic phase diversity.

In this paper, we describe Fourier-based Wave Front Sensors (WFS) as linear integral operators, characterized by their Kernel. In a first part, we derive the dependency of this quantity with respect to the WFS's optical parameters: pupil geometry, filtering mask, tip/tilt modulation. In a second part we focus the study on the special case of convolutional Kernels. The assumptions required to be in such a regime are described. We then show that these convolutional kernels allow to drastically simplify the WFS's model by summarizing its behavior in a concise and comprehensive quantity called the WFS's Impulse Response. We explain in particular how it allows to compute the sensor's sensitivity with respect to the spatial frequencies. Such an approach therefore provides a fast diagnostic tool to compare and optimize Fourier-based WFSs. In a third part, we develop the impact of the residual phases on the sensor's impulse response, and show that the convolutional model remains valid. Finally, a section dedicated to the Pyramid WFS concludes this work, and illustrates how the slopes maps are easily handled by the convolutional model.

Jets are rarely associated with pre-main-sequence intermediate-mass stars. Optical and near-IR observations of jet-driving sources are often hindered by the presence of a natal envelope. Jets around partly embedded sources are a useful diagnostic to constrain the geometry of the concealed protoplanetary disk. In fact, the jet-driving mechanisms are affected by both spatial anisotropies and episodic variations at the (sub-)au scale from the star. We obtained a rich set of high-contrast VLT/SPHERE observations from 0.6 micron to 2.2 micron of the young intermediate-mass star RY Tau. Given the proximity to the Sun of this source, our images have the highest spatial resolution ever obtained for an atomic jet. Optical observations in polarized light show no sign of the protoplanetary disk detected by ALMA. Instead, we observed a diffuse signal resembling a remnant envelope with an outflow cavity. The jet is detected in four spectral lines. The jet appears to be wiggling and its radial width increasing with the distance is complementary to the shape of the outflow cavity suggesting a strong jet/envelope interaction. Through the estimated tangential velocity, we revealed a possible connection between the launching time of the jet sub-structures and the stellar activity of RY Tau. RY Tau is at an intermediate stage toward the dispersal of the natal envelope. This source shows episodic increases of mass accretion/ejection similarly to other known intermediate-mass stars. The amount of observed jet wiggle is consistent with the presence of a precessing disk warp or misaligned inner disk that would be induced by an unseen planetary/sub-stellar companion at sub-/few-au scales. The high disk mass of RY Tau and of two other jet-driving intermediate-mass stars, HD163296 and MWC480, suggests that massive, full disks are more efficient at launching prominent jets.

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.

Wavefront sensors encode phase information of an incoming wavefront into an intensity pattern that can be measured on a camera. Several kinds of wavefront sensors (WFS) are used in astronomical adaptive optics. Amongst them, Fourier-based wavefront sensors perform a filtering operation on the wavefront in the focal plane. The most well known example of a WFS of this kind is the Zernike wavefront sensor, and the pyramid wavefront sensor (PWFS) also belongs to this class. Based on this same principle, new WFSs can be proposed such as the n-faced pyramid (which ultimately becomes an axicone) or the flattened pyramid, depending on whether the image formation is incoherent or coherent. In order to test such novel concepts, the LOOPS adaptive optics testbed hosted at the Laboratoire d'Astrophysique de Marseille has been upgraded by adding a Spatial Light Modulator (SLM). This device, placed in a focal plane produces high-definition phase masks that mimic otherwise bulk optic devices. In this paper, we first present the optical design and upgrades made to the experimental setup of the LOOPS bench. Then, we focus on the generation of the phase masks with the SLM and the implications of having such a device in a focal plane. Finally, we present the first closed-loop results in either static or dynamic mode with different WFS applied on the SLM.

Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered light observations with VLT/SPHERE of the transition disks J1608 and J1852. We complement our datasets with ALMA observations and with unresolved photometric observations covering a wide range of wavelengths. We performed radiative transfer modeling to analyze the morphology of the disks and compare the results with a sample of 20 other transition disks observed with both SPHERE and ALMA. Results. The scattered light image of J1608 reveals a very inclined disk, with two bright lobes and a large cavity. J1852 shows an inner ring extending beyond the coronagraphic radius up to 15au, a gap and a second ring at 42au. Our radiative transfer model of J1608 indicates that the millimeter-sized grains are less extended vertically and radially than the micron-sized grains, indicating advanced settling and radial drift. We find good agreement with the observations of J1852 with a similar model, but due to the low inclination of the system, the model remains partly degenerate. The analysis of 22 transition disks shows that, in general, the cavities observed in scattered light are smaller than the ones detected at millimeter wavelengths. Conclusions. The analysis of a sample of transition disks indicates that the small grains can flow inward of the region where millimeter grains are trapped. While 15 out of the 22 cavities in our sample could be explained by a planet of less than 13 Jupiter masses, the others either require the presence of a more massive companion or of several low-mass planets.

The European Extremely Large Telescope will see first lights by the end of 2024. With a diameter of almost 40 meters, it will be the biggest optical telescope ever built from the ground. The ELT will open a brand new window in a sensitivity / spatial angular resolution parameter space. To take the full benefit of the scientific potential of this giant, all the instruments will be equipped with Adaptive Optics (AO), providing the sharpest images. This paper provides a quick overview of the AO capabilities of the future instruments to be deployed at the ELT, highlight some of the expected performance and describe a couple of technical challenges that are still to tackle for an optimal scientific use. This paper has been presented at the "Societe Francaise d'Astronomie et Astrophysique" symposium in Bordeaux 2018, it is then naturally biased toward the French contribution for the ELT.

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 segmentation of the telescope pupil (by spiders & the segmented M4) create areas of phase isolated by the width of the spiders on the wavefront sensor (WFS), breaking the spatial continuity of the wavefront. The poor sensitivity of the Pyramid WFS (PWFS) to differential piston leads to badly seen and therefore uncontrollable differential pistons. In close loop operation, differential pistons between segments will settle around integer values of the average sensing wavelength. The differential pistons typically range from one to ten times the sensing wavelength and vary rapidly over time, leading to extremely poor performance. In addition, aberrations created by atmospheric turbulence will contain large amounts of differential piston between the segments. Removing piston contribution over each of the DM segments leads to poor performance. In an attempt to reduce the impact of unwanted differential pistons that are injected by the AO correction, we compare three different approaches. We first limit ourselves to only use the information measured by the PWFS, in particular by reducing the modulation. We show that using this information sensibly is important but will not be sufficient. We discuss possible ways of improvement by using prior information. A second approach is based on phase closure of the DM commands and assumes the continuity of the correction wavefront over the entire unsegmented pupil. The last approach is based on the pair-wise slaving of edge actuators and shows the best results. We compare the performance of these methods using realistic end-to-end simulations. We find that pair-wise slaving leads to a small increase of the total wavefront error, only adding between 20-45 nm RMS in quadrature for seeing conditions between 0.45-0.85 arcsec. Finally, we discuss the possibility of combining the different proposed solutions to increase robustness.

In this paper we present HighRes: a laboratory demonstration of a 3U CubeSat with a deployable primary mirror that has the potential of achieving high-resolution imaging for Earth Observation. The system is based on a Cassegrain telescope with a segmented primary mirror composed of 4 petals that form an effective aperture of 300 mm. The design provides diffraction limited performance over the entire field-of-view and allows for a panchromatic ground-sampling distance of less than 1 m at an altitude of 350 km. The alignment and co-phasing of the mirror segments is performed by focal plane sharpening and is validated through rigorous numerical simulations. The opto-mechanical design of the prototype and its laboratory demonstration are described and measurements from the on-board metrology sensors are presented. This data verifies that the performance of the mirror deployment and manipulation systems is sufficient for co-phasing. In addition, it is shown that the mirrors can be driven to any target position with an accuracy of 25 nm using closed-loop feedback between the mirror motors and the on-board metrology.

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.

Context. The next generation of space-borne instruments dedicated to the direct detection of exoplanets requires unprecedented levels of wavefront control precision. Coronagraphic wavefront sensing techniques for these instruments must measure both the phase and amplitude of the optical aberrations using the scientific camera as a wavefront sensor. Aims. In this paper, we develop an extension of coronagraphic phase diversity to the estimation of the complex electric field, that is, the joint estimation of phase and amplitude. Methods. We introduced the formalism for complex coronagraphic phase diversity. We have demonstrated experimentally on the Tr\`es Haute Dynamique testbed at the Observatoire de Paris that it is possible to reconstruct phase and amplitude aberrations with a subnanometric precision using coronagraphic phase diversity. Finally, we have performed the first comparison between the complex wavefront estimated using coronagraphic phase diversity (which relies on time-modulation of the speckle pattern) and the one reconstructed by the self-coherent camera (which relies on the spatial modulation of the speckle pattern). Results. We demonstrate that coronagraphic phase diversity retrieves complex wavefront with subnanometric precision with a good agreement with the reconstruction performed using the self-coherent camera. Conclusions. This result paves the way to coronagraphic phase diversity as a coronagraphic wave-front sensor candidate for very high contrast space missions.

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.

The SPHERE (spectro-photometric exoplanet research) extreme-AO planet hunter saw first light at the VLT observatory on Mount Paranal in May 2014 after ten years of development. Great efforts were put into modelling its performance, particularly in terms of achievable contrast, and to budgeting instrumental features such as wave front errors and optical transmission to each of the instrument's three focal planes, the near infrared dual imaging camera IRDIS, the near infrared integral field spectrograph IFS and the visible polarimetric camera ZIMPOL. In this paper we aim at comparing predicted performance with measured performance. In addition to comparing on-sky contrast curves and calibrated transmission measurements, we also compare the PSD-based wave front error budget with in-situ wave front maps obtained thanks to a Zernike phase mask, ZELDA, implemented in the infrared coronagraph wheel. One of the most critical elements of the SPHERE system is its high-order deformable mirror, a prototype 40x40 actuator piezo stack design developed in parallel with the instrument itself. The development was a success, as witnessed by the instrument performance, in spite of some bad surprises discovered on the way. The devastating effects of operating without taking properly into account the loss of several actuators and the thermally and temporally induced variations in the DM shape will be analysed, and the actions taken to mitigate these defects through the introduction of specially designed Lyot stops and activation of one of the mirrors in the optical train will be described.

The imaging and spectroscopy of habitable worlds will require large-aperture space-based telescopes, to increase the collecting area and the angular resolution. These large telescopes will necessarily use segmented primaries to fit in a rocket. However, these massively segmented mirrors make high-contrast performance very difficult to achieve and stabilize, compared to more common monolithic primaries. Despite space telescopes operating in a friendlier environment than ground-based telescopes, remaining vibrations and resonant modes on the segments can still deteriorate the performance. In this context, we present the Pair-based Analytical model for Segmented Telescopes Imaging from Space (PASTIS) that enables the establishment of a comprehensive error budget, both in term of segment alignment and stability. Using this model, one may evaluate the influence of the segment cophasing and surface quality evolution on the final images and contrasts, and set up requirements for any given mission. One can also identify the dominant modes of a given geometry for a given coronagraphic instrument and design the feedback control systems accordingly. In this paper, we first develop and validate this analytical model by comparing its outputs to the images and contrasts predicted by an end-to-end simulation. We show that the contrasts predicted using PASTIS are accurate enough compared to the end-to-end propagation results, at the exo-Earth detection level. Second, we develop a method for a fast and efficient error budget in term of segment manufacturing and alignment that takes into account the disparities of the segment effects on the final performance. This technique is then applied on a specific aperture to provide static and quasi-static requirements on each segment for local aberrations. Finally we discuss potential application of this new technique to future missions.

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.

HARMONI is a visible and near-infrared integral field spectrograph equipped with two complementary adaptive optics systems, fully integrated within the instrument. A Single Conjugate AO (SCAO) system offers high performance for a limited sky coverage and a Laser Tomographic AO (LTAO) system provides AO correction with a very high sky-coverage. While the deformable mirror performing real-time correction of the atmospheric disturbances is located within the telescope itself, the instrument contains a suite of state-of-the-art and innovative wavefront sensor systems. Laser guide star sensors (LGSS) are located at the entrance of the instrument and fed by a dichroic beam splitter, while the various natural guide star sensors for LTAO and SCAO are located close to the science focal plane. We present opto-mechanical architecture and design at PDR level for these wavefront sensor systems.

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.

HARMONI is a visible and NIR integral field spectrograph, providing the E-ELT's core spectroscopic capability at first light. HARMONI will work at the diffraction limit of the E-ELT, thanks to a Classical and a Laser Tomographic AO system. In this paper, we present the system choices that have been made for these SCAO and LTAO modules. In particular, we describe the strategy developed for the different Wave-Front Sensors: pyramid for SCAO, the LGSWFS concept, the NGSWFS path, and the truth sensor capabilities. We present first potential implementations. And we asses the first system performance.

We present a proof of concept for a new algorithm which can be used to detect exoplanets in high contrast images. The algorithm properly combines mutliple observations acquired during different nights, taking into account the orbital motion of the planet. Methods. We simulate SPHERE/IRDIS time series of observations in which we blindly inject planets on random orbits, at random level of S/N, below the detection limit (down to S/N 1.5). We then use an optimization algorithm to guess the orbital parameters, and take into account the orbital motion to properly recombine the different images, and eventually detect the planets. We show that an optimization algorithm can indeed be used to find undetected planets in temporal sequences of images, even if they are spread over orbital time scales. As expected, the typical gain in S/N ratio is sqrt(n), n being the number of observations combined. We find that the K-Stacker algorithm is able de-orbit and combine the images to reach a level of performance similar to what could be expected if the planet was not moving. We find recovery rates of 50% at S/N=5. We also find that the algorithm is able to determine the position of the planet in individual frames at one pixel precision, even despite the fact that the planet itself is below the detection limit in each frame. Our simulations show that K-Stacker can be used to increase the contrast limit of current exoplanet imaging instruments and to discover fainter bodies. We also suggest that the ability of K-Stacker to determine the position of the planet in every image of the time serie could be used as part of a new observing strategy in which long exposures would be broken into shorter ones spread over months. This could make possible to determine the orbital parameters of a planet without requiring multiple high S/N > 5 detections.

The low wind effect (LWE) refers to a characteristic set of quasi-static wavefront aberrations seen consistently by the SPHERE instrument when dome-level wind speeds drop below 3 m/s. This effect produces bright low-order speckles in the stellar PSF, which severely limit the contrast performance of SPHERE under otherwise optimal observing conditions. In this paper we propose the Fast & Furious (F&F) phase diversity algorithm as a viable software-only solution for real-time LWE compensation, which would utilise image sequences from the SPHERE differential tip-tilt sensor (DTTS). We evaluated the closed-loop performance of F&F on the MITHIC high-contrast test-bench under a variety of conditions emulating LWE-affected DTTS images, in order to assess the expected performance of an on-sky implementation of F&F in SPHERE. The algorithm was found to be capable of returning such LWE-affected images to Strehl ratios of greater than 90% within five iterations, for all appropriate laboratory test cases. These results are highly representative of predictive simulations, and demonstrate the stability of the algorithm against a wide range of factors including low image signal-to-noise ratio (S/N), small image field of view, and amplitude errors. It was also found in simulation that closed-loop stability can be preserved down to image S/N as low as five while still improving overall wavefront quality, allowing for reliable operation even on faint targets. The Fast & Furious algorithm is an extremely promising solution for real-time compensation of the LWE, which can operate simultaneously with science observations and may be implemented in SPHERE without requiring additional hardware. The robustness and relatively large effective dynamic range of F&F also make it suitable for general wavefront optimisation applications, including the co-phasing of segmented ELT-class telescopes.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

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.

HD$\,$169142 is an excellent target to investigate signs of planet-disk interaction due to the previous evidence of gap structures. We performed J-band (~1.2{\mu}m) polarized intensity imaging of HD169142 with VLT/SPHERE. We observe polarized scattered light down to 0.16" (~19 au) and find an inner gap with a significantly reduced scattered light flux. We confirm the previously detected double ring structure peaking at 0.18" (~21 au) and 0.56" (~66 au), and marginally detect a faint third gap at 0.70"-0.73" (~82-85 au). We explore dust evolution models in a disk perturbed by two giant planets, as well as models with a parameterized dust size distribution. The dust evolution model is able to reproduce the ring locations and gap widths in polarized intensity, but fails to reproduce their depths. It, however, gives a good match with the ALMA dust continuum image at 1.3 mm. Models with a parameterized dust size distribution better reproduce the gap depth in scattered light, suggesting that dust filtration at the outer edges of the gaps is less effective. The pile-up of millimeter grains in a dust trap and the continuous distribution of small grains throughout the gap likely require a more efficient dust fragmentation and dust diffusion in the dust trap. Alternatively, turbulence or charging effects might lead to a reservoir of small grains at the surface layer that is not affected by the dust growth and fragmentation cycle dominating the dense disk midplane. The exploration of models shows that extracting planet properties such as mass from observed gap profiles is highly degenerate.

We build on a long-standing tradition in astronomical adaptive optics (AO) of specifying performance metrics and error budgets using linear systems modeling in the spatial-frequency domain. Our goal is to provide a comprehensive tool for the calculation of error budgets in terms of residual temporally filtered phase power spectral densities and variances. In addition, the fast simulation of AO-corrected point spread functions (PSFs) provided by this method can be used as inputs for simulations of science observations with next-generation instruments and telescopes, in particular to predict post-coronagraphic contrast improvements for planet finder systems. We extend the previous results and propose the synthesis of a distributed Kalman filter to mitigate both aniso-servo-lag and aliasing errors whilst minimizing the overall residual variance. We discuss applications to (i) analytic AO-corrected PSF modeling in the spatial-frequency domain, (ii) post-coronagraphic contrast enhancement, (iii) filter optimization for real-time wavefront reconstruction, and (iv) PSF reconstruction from system telemetry. Under perfect knowledge of wind velocities, we show that $\sim$60 nm rms error reduction can be achieved with the distributed Kalman filter embodying anti- aliasing reconstructors on 10 m class high-order AO systems, leading to contrast improvement factors of up to three orders of magnitude at few ${\lambda}/D$ separations ($\sim1-5{\lambda}/D$) for a 0 magnitude star and reaching close to one order of magnitude for a 12 magnitude star.

We propose and apply two methods to estimate pupil plane phase discontinuities for two realistic scenarios on VLT and Keck. The methods use both Phase Diversity and a form of image sharpening. For the case of VLT, we simulate the `low wind effect' (LWE) which is responsible for focal plane errors in the SPHERE system in low wind and good seeing conditions. We successfully estimate the simulated LWE using both methods, and show that they are complimentary to one another. We also demonstrate that single image Phase Diversity (also known as Phase Retrieval with diversity) is also capable of estimating the simulated LWE when using the natural de-focus on the SPHERE/DTTS imager. We demonstrate that Phase Diversity can estimate the LWE to within 30 nm RMS WFE, which is within the allowable tolerances to achieve a target SPHERE contrast of 10$^{-6}$. Finally, we simulate 153 nm RMS of piston errors on the mirror segments of Keck and produce NIRC2 images subject to these effects. We show that a single, diverse image with 1.5 waves (PV) of focus can be used to estimate this error to within 29 nm RMS WFE, and a perfect correction of our estimation would increase the Strehl ratio of a NIRC2 image by 12\%

We have observed the protoplanetary disk of the well-known young Herbig star HD 142527 using ZIMPOL Polarimetric Differential Imaging with the VBB (Very Broad Band, ~600-900nm) filter. We obtained two datasets in May 2015 and March 2016. Our data allow us to explore dust scattering around the star down to a radius of ~0.025" (~4au). The well-known outer disk is clearly detected, at higher resolution than before, and shows previously unknown sub-structures, including spirals going inwards into the cavity. Close to the star, dust scattering is detected at high signal-to-noise ratio, but it is unclear whether the signal represents the inner disk, which has been linked to the two prominent local minima in the scattering of the outer disk, interpreted as shadows. An interpretation of an inclined inner disk combined with a dust halo is compatible with both our and previous observations, but other arrangements of the dust cannot be ruled out. Dust scattering is also present within the large gap between ~30 and ~140au. The comparison of the two datasets suggests rapid evolution of the inner regions of the disk, potentially driven by the interaction with the close-in M-dwarf companion, around which no polarimetric signal is detected.

The use of Fourier methods in wave-front reconstruction can significantly reduce the computation time for large telescopes with a high number of degrees of freedom. However, Fourier algorithms for discrete data require a rectangular data set which conform to specific boundary requirements, whereas wave-front sensor data is typically defined over a circular domain (the telescope pupil). Here we present an iterative Gerchberg routine modified for the purposes of discrete wave-front reconstruction which adapts the measurement data (wave-front sensor slopes) for Fourier analysis, fulfilling the requirements of the Fast Fourier Transform (FFT) and providing accurate reconstruction. The routine is used in the adaptation step only and can be coupled to any other Wiener-like or least-squares method. We compare simulations using this method with previous Fourier methods and show an increase in performance in terms of Strehl ratio and a reduction in noise propagation for a 40x40 SPHERE-like adaptive optics system. For closed loop operation with minimal iterations the Gerchberg method provides an improvement in Strehl, from 95.4% to 96.9% in K-band. This corresponds to ~40 nm improvement in rms, and avoids the high spatial frequency errors present in other methods, providing an increase in contrast towards the edge of the correctable band.

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.

In this article, we compare a set of Wave Front Sensors (WFS) based on Fourier filtering technique. In particular, this study explores the "class of pyramidal WFS" defined as the 4 faces pyramid WFS, all its recent variations (6, 8 faces, the flattened PWFS, etc.) and also some new WFSs as the flattened cone WFS or the 3 faces pyramid WFS. In the first part, we describe such a sensors class thanks to the optical parameters of the Fourier filtering mask and the modulation parameters. In the second part, we use the unified formalism to create a set of performance criteria: size of the signal on the detector, efficiency of incoming flux, sensitivity, linear range and chromaticity. In the third part, we show the influence of the previous optical and modulation parameters on these performance criteria. This exhaustive study allows to know how to optimize the sensor regarding to performance specifications. We show in particular that the number of faces has influence on the number of pixels required to do the wave front sensing but no influence on the sensitivity and linearity range. To modify these criteria, we show that the modulation radius and the apex angle are much more relevant. Moreover we observe that the time spent on edges or faces during a modulation cycle allows to adjust the trade-off between sensitivity and linearity range.

The ultimate performance of coronagraphic high contrast exoplanet imaging systems such as SPHERE or GPI is limited by quasi-static aberrations. These aberrations produce speckles that can be mistaken for planets in the image. In order to design instruments, correct quasi-static aberrations or analyse data, the expression of the point spread function of a coronagraphic instrument in the presence of residual turbulence is most useful. Here we derive an analytic expression for this point spread function that is an extension to coronagraphic imaging of Roddier's expression for imaging through turbulence. We give a physical interpretation of its structure, we validate it by numerical simulations and we show that it is computationally efficient. Finally, we incorporate this imaging model into a coronagraphic phase diversity method (COFFEE) and validate by simulations that it allows wave-front reconstruction in the presence of residual turbulence. The preliminary results, which give a sub-nanometric precision in the case of a SPHERE-like system, strongly suggest that quasi-static aberrations could be calibrated during observations by this method.

We introduce in this article a general formalism for Fourier based wave front sensing. To do so, we consider the filtering mask as a free parameter. Such an approach allows to unify sensors like the Pyramid Wave Front Sensor (PWFS) and the Zernike Wave Front Sensor (ZWFS). In particular, we take the opportunity to generalize this two sensors in terms of sensors' class where optical quantities as, for instance, the apex angle for the PWFS or the depth of the Zernike mask for the ZWFS become free parameters. In order to compare all the generated sensors of this two classes thanks to common performance criteria, we firstly define a general phase-linear quantity that we call meta-intensity. Analytical developments allow then to split the perfectly phase-linear behavior of a WFS from the non-linear contributions making robust and analytic definitions of the sensitivity and the linearity range possible. Moreover, we define a new quantity called the SD factor which characterizes the trade-off between these two antagonist quantities. These developments are generalized for modulation device and polychromatic light. A non-exhaustive study is finally led on the two classes allowing to retrieve the usual results and also make explicit the influence of the optical parameters introduced above.

Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.

We present scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 0.63, 0.79, 1.24 and 1.62 micron. We also present H2/H3-band ADI observations. Three distinct radial depressions in the polarized intensity distribution are seen, around 85, 21, and 6~au. The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average towards the South and darker towards the North-West. The ADI observations yielded no signifiant detection of point sources in the disk. Our observations have a linear spatial resolution of 1 to 2au, similar to that of recent ALMA dust continuum observations. The sub-micron sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas. We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling. This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations. We find that the gas surface density in the "gaps" is reduced by 50% to 80% relative to an unperturbed model. Should embedded planets be responsible for carving the gaps then their masses are at most a few 10 Mearth. The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets. If forming planetary bodies have undergone collapse and are in the "detachted phase" then they may be directly observable with future facilities such as METIS at the E-ELT.

Sky-coverage in laser-assisted AO observations largely depends on the system's capability to guide on the faintest natural guide-stars possible. Here we give an up-to-date status of our natural guide-star processing tailored to the European-ELT's visible and near-infrared (0.47 to 2.45 {\mu}m) integral field spectrograph - Harmoni. We tour the processing of both the isoplanatic and anisoplanatic tilt modes using the spatio-angular approach whereby the wave-front is estimated directly in the pupil plane avoiding a cumbersome explicit layered estimation on the 35-layer profiles we're currently using. Taking the case of Harmoni, we cover the choice of wave-front sensors, the number and field location of guide-stars, the optimised algorithms to beat down angular anisoplanatism and the performance obtained with different temporal controllers under split high-order/low-order tomography or joint tomography. We consider both atmospheric and far greater telescope wind buffeting disturbances. In addition we provide the sky-coverage estimates thus obtained.

Warm or massive gas giant planets, brown dwarfs, and debris disks around nearby stars are now routinely observed by dedicated high-contrast imaging instruments on large, ground-based observatories. These facilities include extreme adaptive optics (ExAO) and state-of-the-art coronagraphy to achieve unprecedented sensitivities for exoplanet detection and spectral characterization. However, differential aberrations between the ExAO sensing path and the science path 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\as) from their host star. In our previous work, we proposed a wavefront sensor based on Zernike phase contrast methods to circumvent this issue and measure these quasi-static aberrations at a nanometric level. We present the design, manufacturing and testing of ZELDA, a prototype that was installed on VLT/SPHERE during its reintegration in Chile. Using the internal light source of the instrument, we performed measurements in the presence of Zernike or Fourier modes introduced with the deformable mirror. Our experimental and simulation results are consistent, confirming the ability of our sensor to measure small aberrations (<50 nm rms) with nanometric accuracy. We then corrected the long-lived non-common path aberrations in SPHERE based on ZELDA measurements. We estimated a contrast gain of 10 in the coronagraphic image at 0.2\as, reaching the raw contrast limit set by the coronagraph in the instrument. The simplicity of the design and its phase reconstruction algorithm makes ZELDA an excellent candidate for the on-line measurements of quasi-static aberrations during the observations. The implementation of a ZELDA-based sensing path on the current and future facilities (ELTs, future space missions) could ease the observation of the cold gaseous or massive rocky planets around nearby stars.

We present new near-infrared photometric measurements of the core of the young massive cluster NGC 3603 obtained with extreme adaptive optics. The data were obtained with the SPHERE instrument mounted on ESO Very Large Telescope, and cover three fields in the core of this cluster. We applied a correction for the effect of extinction to our data obtained in the J and K broadband filters and estimated the mass of detected sources inside the field of view of SPHERE/IRDIS, which is 13.5"x13.5". We derived the mass function (MF) slope for each spectral band and field. The MF slope in the core is unusual compared to previous results based on Hubble space telescope (HST) and very large telescope (VLT) observations. The average slope in the core is estimated as -1.06^{+0.26}_{-0.26} for the main sequence stars with 3.5 Msun < M < 120 Msun.Thanks to the SPHERE extreme adaptive optics, 814 low-mass stars were detected to estimate the MF slope for the pre-main sequence stars with 0.6 Msun< M < 3.5 Msun , Gamma = -0.54^{+0.11}_{-0.11} in the K-band images in two fields in the core of the cluster. For the first time, we derive the mass function of the very core of the NGC 3603 young cluster for masses in the range 0.6 - 120 Msun. Previous studies were either limited by crowding, lack of dynamic range, or a combination of both.

Spectral characterization of young, giant exoplanets detected by direct imaging is one of the tasks of the new generation of high-contrast imagers. For this purpose, the VLT/SPHERE instrument includes a unique long-slit spectroscopy (LSS) mode coupled with Lyot coronagraphy in its infrared dual-band imager and spectrograph (IRDIS). The performance of this mode is intrinsically limited by the use of a non-optimal coronagraph, but in a previous work we demonstrated that it could be significantly improved at small inner-working angles using the stop-less Lyot coronagraph (SLLC). We now present the development, testing, and validation of the first SLLC prototype for VLT/SPHERE. Based on the transmission profile previously proposed, the prototype was manufactured using microdots technology and was installed inside the instrument in 2014. The transmission measurements agree well with the specifications, except in the very low transmissions (<5% in amplitude). The performance of the SLLC is tested in both imaging and spectroscopy using data acquired on the internal source. In imaging, we obtain a raw contrast gain of a factor 10 at 0.3" and 5 at 0.5" with the SLLC. Using data acquired with a focal-plane mask, we also demonstrate that no Lyot stop is required to reach the full performance, which validates the SLLC concept. Comparison with a realistic simulation model shows that we are currently limited by the internal phase aberrations of SPHERE. In spectroscopy, we obtain a gain of ~1 mag in a limited range of angular separations. Simulations show that although the main limitation comes from phase errors, the performance in the non-SLLC case is very close to the ultimate limit of the LSS mode. Finally, we obtain the very first on-sky data with the SLLC, which appear extremely promising for the future scientific exploitation of an apodized LSS mode in SPHERE.

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

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]

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.

This year, a second generation of coronagraphs dedicated to high-contrast direct imaging of exoplanets is starting operations. Among them, SPHERE, installed at the focus of the UT3 Very Large Telescope, reaches unprecedented contrast ratios up to $10^{-6}$ -$ 10^{-7}$, using eXtreme Adaptive Optics and the Angular Differential Imaging (ADI) techniques. In this paper, we present a new method called Keplerian-Stacker that improves the detection limit of high contrast instruments like SPHERE, by up to a factor of 10. It consists of observing a star on a long enough period to let a hypothetical planet around that star move along its orbit. Even if in each individual observation taken during one night, we do not detect anything, we show that it is possible, using an optimization algorithm, to re-center the images according to keplerian motions (ex: 10-100 images taken over a long period of typically 1-10 years) and detect planets otherwise unreachable. This method can be used in combination with the ADI technics (or possibly any other high contrast data reduction method) to improve the Signal to Noise Ratio in each individual image, and to further improve the global detection limit. It also directly provides orbital parameters of the detected planets, as a by-product of the optimization algorithm.

We propose a new type of Wave Front Sensor (WFS) derived from the Pyramid WFS (PWFS). This new WFS, called the Flattened Pyramid-WFS (FPWFS), has a reduced Pyramid angle in order to optically overlap the four pupil images into an unique intensity. This map is then used to derive the phase information. In this letter this new WFS is compared to three existing WFSs, namely the PWFS, the Modulated PWFS (MPWFS) and the Zernike WFS (ZWFS) following tests about sensitivity, linearity range and low photon flux behavior. The FPWFS turns out to be more linear than a modulated pyramid for the high-spatial order aberrations but it provides an improved sensitivity compared to the non-modulated pyramid. The noise propagation may even be as low as the ZWFS for some given radial orders. Furthermore, the pixel arrangement being more efficient than for the PWFS, the FPWFS seems particularly well suited for high-contrast applications.

We present the performance of the Integral Field Spectrograph (IFS) of SPHERE, the high-contrast imager for the ESO VLT telescope designed to perform imaging and spectroscopy of extrasolar planets, obtained from tests performed at the Institute de Plan\'etologie et d'Astrophysique de Grenoble facility during the integration phase of the instrument.} {The tests were performed using the instrument software purposely prepared for SPHERE. The output data were reduced applying the SPHERE data reduction and handling software, adding an improved spectral deconvolution procedure. To this aim, we prepared an alternative procedure for the spectral subtraction exploiting the principal components analysis algorithm. Moreover, a simulated angular differential imaging procedure was also implemented to estimate how the instrument performed once this procedure was applied at telescope. The capability of the IFS to faithfully retrieve the spectra of the detected faint companions was also considered.} {We found that the application of the updated version of the spectral deconvolution procedure alone, when the algorithm throughput is properly taken into account, gives us a $5\sigma$ limiting contrast of the order of 5$\times$$10^{-6}$ or slightly better. The further application of the angular differential imaging procedure on these data should allow us to improve the contrast by one order of magnitude down to around 7$\times$$10^{-7}$ at a separation of 0.3 arcsec. The application of a principal components analysis procedure that simultaneously uses spectral and angular data gives comparable results. Finally, we found that the reproducibility of the spectra of the detected faint companions is greatly improved when angular differential imaging is applied in addition to the spectral deconvolution.

The new planet finder for the Very Large Telescope (VLT), the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE), just had its first light in Paranal. A dedicated instrument for the direct detection of planets, SPHERE, is composed of a polametric camera in visible light, the Zurich IMager POLarimeter (ZIMPOL), and two near-infrared sub-systems: the Infra-Red Dual-beam Imager and Spectrograph (IRDIS), a multi-purpose camera for imaging, polarimetry, and long-slit spectroscopy, and the integral field spectrograph (IFS), an integral field spectrograph. We present the results obtained from the analysis of data taken during the laboratory integration and validation phase, after the injection of synthetic planets. Since no continuous field rotation could be performed in the laboratory, this analysis presents results obtained using reduction techniques that do not use the angular differential imaging (ADI) technique. To perform the simulations, we used the instrumental point spread function (PSF) and model spectra of L and T-type objects scaled in contrast with respect to the host star. We evaluated the expected error in astrometry and photometry as a function of the signal to noise of companions, after spectral differential imaging (SDI) reduction for IRDIS and spectral deconvolution (SD) or principal component analysis (PCA) data reductions for IFS. We deduced from our analysis, for example, that $\beta$Picb, a 12~Myr old planet of $\sim$10~\MJ and semi-major axis of 9--10 AU, would be detected with IRDIS with a photometric error of 0.16~mag and with a relative astrometric position error of 1.1~mas. With IFS, we could retrieve a spectrum with error bars of about 0.15~mag on each channel and astrometric relative position error of 0.6~mas. For a fainter object such as HR8799d, a 13~\MJ planet at a distance of 27~AU, IRDIS could obtain a relative astrometric error of 3~mas.

The second-generation instrument SPHERE, dedicated to high-contrast imaging, will soon be in operation on the European Very Large Telescope. Such an instrument relies on an extreme adaptive optics system coupled with a coronagraph that suppresses most of the diffracted stellar light. However, the coronagraph performance is strongly limited by quasi-static aberrations that create long-lived speckles in the scientific image plane, which can easily be mistaken for planets. The ultimate performance is thus limited by the unavoidable differential aberrations between the wave-front sensor and the scientific camera, which have to be estimated andcompensated for. In this paper, we use the COFFEE approach to measure and compensate for SPHERE's quasi-static aberrations. COFFEE (for COronagraphic Focal-plane wave-Front Estimation for Exoplanet detection), which consists in an extension of phase diversity to coronagraphic imaging, estimates the quasi-static aberrations, including the differential ones, using only two focal plane images recorded by the scientific camera. In this paper, we use coronagraphic images recorded from SPHERE's infrared detector IRDIS to estimate the aberrations upstream of the coronagraph, which are then compensated for using SPHERE's extreme adaptive optics loop SAXO. We first validate the ability of COFFEE to estimate high-order aberrations by estimating a calibrated influence function pattern introduced upstream of the coronagraph. We then use COFFEE in an original iterative compensation process to compensate for the estimated aberrations, leading to a contrast improvement by a factor that varies from 1.4 to 4.7 between 2l/D and 15l/D on IRDIS. The performance of the compensation process is also evaluated through simulations. An excellent match between experimental results and these simulations is found.

The estimation and compensation of quasi-static aberrations is mandatory to reach the ultimate performance of high-contrast imaging systems. COFFEE is a focal plane wave-front sensing method that consists in the extension of phase diversity to high-contrast imaging systems. Based on a Bayesian approach, it estimates the quasi-static aberrations from two focal plane images recorded from the scientific camera itself. In this paper, we present COFFEE's extension which allows an estimation of low and high order aberrations with nanometric precision for any coronagraphic device. The performance is evaluated by realistic simulations, performed in the SPHERE instrument framework. We develop a myopic estimation that allows us to take into account an imperfect knowledge on the used diversity phase. Lastly, we evaluate COFFEE's performance in a compensation process, to optimize the contrast on the detector, and show it allows one to reach the 10^-6 contrast required by SPHERE at a few resolution elements from the star. Notably, we present a non-linear energy minimization method which can be used to reach very high contrast levels (better than 10^-7 in a SPHERE-like context)

Until now, just a few extrasolar planets (~30 out of 860) have been found through the direct imaging method. This number should greatly improve when the next generation of High Contrast Instruments like Gemini Planet Imager (GPI) at Gemini South Telescope or SPHERE at VLT will became operative at the end of this year. In particular, the Integral Field Spectrograph (IFS), one of the SPHERE subsystems, should allow a first characterization of the spectral type of the found extrasolar planets. Here we present the results of the last performance tests that we have done on the IFS instrument at the Institut de Planetologie et d'Astrophysique de Grenoble (IPAG) in condition as similar as possible to the ones that we will find at the telescope. We have found that we should be able to reach contrast down to 5x10$^{-7}$ and make astrometry at sub-mas level with the instrument in the actual conditions. A number of critical issues have been identified. The resolution of these problems could allow to further improve the performance of the instrument.

Next year the second generation instrument SPHERE will begin science operations at the Very Large Telecope (ESO). This instrument will be dedicated to the search for exoplanets through the direct imaging techniques, with the new generation extreme adaptive optics. In this poster, we present the performances of one of the focal instruments, the Infra-Red Dual-beam Imaging and Spectroscopy (IRDIS). All the results have been obtained with tests in laboratory, simulating the observing conditions in Paranal. We tested several configurations using the sub-system Integral Field Spectrograph (IFS) in parallel and simulating long coronographic exposures on a star, calibrating instrumental ghosts, checking the performance of the adaptive optics system and reducing data with the consortium pipeline. The contrast one can reach with IRDIS is of the order of 2\times 10^{-6}$ at 0.5 arcsec separation from the central star.

The final performance of current and future instruments dedicated to exoplanet detection and characterization (such as SPHERE on the European Very Large Telescope, GPI on Gemini North, or future instruments on Extremely Large Telescopes) is limited by uncorrected quasi-static aberrations. These aberrations create long-lived speckles in the scientific image plane, which can easily be mistaken for planets. Common adaptive optics systems require dedicated components to perform wave-front analysis. The ultimate wave-front measurement performance is thus limited by the unavoidable differential aberrations between the wavefront sensor and the scientific camera. To reach the level of detectivity required by high-contrast imaging, these differential aberrations must be estimated and compensated for. In this paper, we characterize and experimentally validate a wave-front sensing method that relies on focal-plane data. Our method, called COFFEE (for COronagraphic Focal-plane wave-Front Estimation for Exoplanet detection), is based on a Bayesian approach, and it consists in an extension of phase diversity to high-contrast imaging. It estimates the differential aberrations using only two focal-plane coronagraphic images recorded from the scientific camera itself. In this paper, we first present a thorough characterization of COFFEE's performance by means of numerical simulations. This characterization is then compared with an experimental validation of COFFEE using an in-house adaptive optics bench and an apodized Roddier & Roddier phase mask coronagraph. An excellent match between experimental results and the theoretical study is found. Lastly, we present a preliminary validation of COFFEE's ability to compensate for the aberrations upstream of a coronagraph.

Observing sequences have shown that the major noise source limitation in high-contrast imaging is due to the presence of quasi-static speckles. The timescale on which quasi-static speckles evolve, is determined by various factors, among others mechanical or thermal deformations. Understanding of these time-variable instrumental speckles, and especially their interaction with other aberrations, referred to as the pinning effect, is paramount for the search of faint stellar companions. The temporal evolution of quasi-static speckles is for instance required for a quantification of the gain expected when using angular differential imaging (ADI), and to determine the interval on which speckle nulling techniques must be carried out. Following an early analysis of a time series of adaptively corrected, coronagraphic images obtained in a laboratory condition with the High-Order Test bench (HOT) at ESO Headquarters, we confirm our results with new measurements carried out with the SPHERE instrument during its final test phase in Europe. The analysis of the residual speckle pattern in both direct and differential coronagraphic images enables the characterization of the temporal stability of quasi-static speckles. Data were obtained in a thermally actively controlled environment reproducing realistic conditions encountered at the telescope. The temporal evolution of the quasi-static wavefront error exhibits linear power law, which can be used to model quasi-static speckle evolution in the context of forthcoming high-contrast imaging instruments, with implications for instrumentation (design, observing strategies, data reduction). Such a model can be used for instance to derive the timescale on which non-common path aberrations must be sensed and corrected. We found in our data that quasi-static wavefront error increases with ~0.7 angstrom per minute.

The phase diversity technique is a useful tool to measure and pre-compensate for quasi-static aberrations, in particular non-common path aberrations, in an adaptive optics corrected imaging system. In this paper, we propose and validate by simulations an extension of the phase diversity technique that uses long exposure adaptive optics corrected images for sensing quasi-static aberrations during the scientific observation, in particular for high-contrast imaging. The principle of the method is that, for a sufficiently long exposure time, the residual turbulence is averaged into a convolutive component of the image and that phase diversity estimates the sole static aberrations of interest. The advantages of such a procedure, compared to the processing of short-exposure image pairs, are that the separation between static aberrations and turbulence-induced ones is performed by the long-exposure itself and not numerically, that only one image pair must be processed, that the estimation benefits from the high SNR of long-exposure images, and that only the static aberrations of interest are to be estimated. Long-exposure phase diversity can also be used as a phasing sensor for a segmented aperture telescope. Thus, it may be particularly useful for future planet finder projects such as EPICS on the European ELT.

Non-Common Path Aberrations (NCPA) are one of the main limitations for extreme Adaptive Optics (AO) system. NCPA prevent extreme AO systems to achieve their ultimate performance. These static aberrations are unseen by the wave front sensor and therefore not corrected in closed loop. We present experimental results validating new procedures of measurement and pre-compensation of the NCPA on the AO bench at ONERA. The measurement procedure is based on new refined algorithms of phase diversity. The pre-compensation procedure makes use of a pseudo-closed loop scheme to overcome the AO wavefront sensor model uncertainties. Strehl Ratio obtained in the images reaches 98.7% @ 632.8 nm. This result allows us to be confident of achieving the challenging performance required for extrasolar planet direct observation.

The direct imaging from the ground of extrasolar planets has become today a major astronomical and biological focus. This kind of imaging requires simultaneously the use of a dedicated high performance Adaptive Optics [AO] system and a differential imaging camera in order to cancel out the flux coming from the star. In addition, the use of sophisticated post-processing techniques is mandatory to achieve the ultimate detection performance required. In the framework of the SPHERE project, we present here the development of a new technique, based on Maximum A Posteriori [MAP] approach, able to estimate parameters of a faint companion in the vicinity of a bright star, using the multi-wavelength images, the AO closed-loop data as well as some knowledge on non-common path and differential aberrations. Simulation results show a 10^-5 detectivity at 5sigma for angular separation around 15lambda/D with only two images.