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.

One of the major limitations of adaptive optics (AO) corrected image post-processing is the lack of knowledge on the system point spread function (PSF). The PSF is not always available as a direct imaging on isolated point like objects such as stars. Its prediction using AO telemetry also suffers from serious limitations and requires complex and yet not fully operational algorithms. A very attractive solution consists in a direct PSF estimation using the scientific images themselves thanks to blind or myopic post-processing approaches. We demonstrate that such approaches suffer from severe limitations when a joint restitution of object and PSF parameters is performed. As an alternative we propose here a marginalized PSF identification that overcomes this limitation. Then the PSF is used for image post-processing. Here we focus on deconvolution, a post-processing technique to restore the object, given the image and the PSF. We show that the PSF estimated by marginalisation provides good quality deconvolution. The full process of marginalized PSF estimation and deconvolution constitutes a successful blind deconvolution technique. It is tested on simulated data to measure its performance. It is also tested on experimental adaptive optics images of the asteroid 4-Vesta by VLT/SPHERE/Zimpol to demonstrate application to on-sky data.

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

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

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.

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.

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

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.

Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of approximately 20 m/pixel scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta's albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of approximately 20-30 km. Contour extraction shows an accuracy of approximately 1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30-40m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north-south crater dichotomy, provided that they operate at the diffraction limit.

Ground-based astronomy is severely limited by the atmospheric turbulence, resulting in a large Point-Spread-Function (PSF) and poor imaging resolution. Even imaging with Adaptive Optics (AO) cannot completely correct the aberrated wavefront, and a residual turbulence still corrupts the observation. Thus the consequences of the turbulence on the PSF is of first interest when building any ground-based telescope. The Power Spectral Density (PSD) of a spatially stationary turbulent phase carries all the information needed for describing the long-exposure PSF. We then develop an analytical description of the long-exposure PSF as a series expansion of the aberrated phase PSD. Our description of the PSF given the PSD of the phase is a simple theoretical way to describe the impact of turbulence on the PSF. We also show accordance with previous papers when restricting our model to its first expansion order. Finally we derive applications of our formula to some particular cases, such as Kolmogorov or von-K\'arm\'an models, or the AO correction impact on the PSF.

The MEDUSAE method (Multispectral Exoplanet Detection Using Simultaneous Aberration Estimation) is dedicated to the detection of exoplanets and disks features in multispectral high-contrast images. The concept of MEDUSAE is to retrieve both the speckle field and the object map via a stochastic approach to inverse problem (taking into account the statistics of the noise) in the Bayesian framework (using parametric regularization). One fundamental aspect of MEDUSAE is that the model of the coronagraphic PSF is analytic and parametrized by the optical path difference which, contrary to the phase, is achromatic. The speckle field is thus estimated by a phase retrieval, using the spectral diversity to disentangle the planetary signal from the residual starlight. The object map is restored via a non-myopic deconvolution under adequate regularization. The basis of this MEDUSAE method have been previously published and validated on an inverse crime. In this communication, we present its application to realistic simulated data, in preparation for real data application. The solution we proposed to attempt bypassing the main differences between the model used for the inversion and the real data is not sufficient: it is now necessary to make the model of the coronagraphic PSF more realistic.

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 G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Whether is ground-based or space-based, any optical instrument suffers from some amount of optical geometric distortion. Recently, the diffraction-limited image quality afforded by space-based telescopes and by Adaptive Optics (AO) corrected instruments on ground based-telescope, have increased the relative importance of the error terms induced by optical distortions. In particular, variable distortions present in Multi-Conjugated Adaptive Optics (MCAO) data are limiting the astrometry and photometry accuracy of such high resolution instruments. From there, the ability to deal with those phenomenon had become a critical issue for high-precision studies. We present in this paper an optimal method of distortion correction for high angular resolution images. Based on a prior-knowledge of the static distortion the method aims to correct the dynamical distortions specifically for each observation set and each frame. The method follows an inverse problem approach based on the work done by Gratadour et al. 2005 on image re-centering, and we aim to generalized it to any kind of distortion mode. The complete formalism of a Weighted Least Square minimization, as well as a detailed characterization of the error budget are presented. In particular we study the influence of different parameters such as the number of frames and the density of the field (sparse or crowed images), of the noise level, and of the aliasing effect. Finally, we show the first application of the method on real observations collected with the Gemini MCAO instrument, GeMS/GSAOI. The performance as well as the gain brought by this method are presented.

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\%

A novel method nicknamed ELASTIC is proposed for the alignment of multiple-aperture telescopes, in particular segmented telescopes. It only needs the acquisition of two diversity images of an unresolved source, and is based on the computation of a modified, frequency-shifted, cross-spectrum. It provides a polychromatic large range tip/tilt estimation with the existing hardware and an inexpensive noniterative unsupervised algorithm. Its performance is studied and optimized by means of simulations. They show that with 5000 photo-electrons/sub-aperture/frame and 1024x1024 pixel images, residues are within the capture range of interferometric phasing algorithms such as phase diversity. The closed-loop alignment of a 6 sub-aperture mirror provides an experimental demonstration of the effectiveness of the method. Author accepted version. Final version is Copyright 2017 Optical Society of America. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modifications of the content of this paper are prohibited.

51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, recently discovered by direct imaging. Being only 0.5" away from its host star it is well suited for spectroscopic analysis using integral field spectrographs. We aim to refine the atmospheric properties of this and to further constrain the architecture of the system by searching for additional companions. Using the SPHERE instrument at the VLT we extend the spectral coverage of the planet to the complete Y- to H-band range and provide photometry in the K12-bands (2.11, 2.25 micron). The object is compared to other cool and peculiar dwarfs. Furthermore, the posterior probability distributions of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. For probing the innermost region for additional companions, archival VLT-NACO (L') SAM data is used. We present the first spectrophotometric measurements in the Y- and K-bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity and cloud sedimentation parameter f_sed. We find that the atmosphere is highly super-solar (Fe/H~1.0) with an extended, thick cloud cover of small particles. The model radius and surface gravity suggest planetary masses of about 9 M_jup. The evolutionary model only provides a lower mass limit of >2 M_jup (for pure hot-start). The cold-start model cannot explain the planet's luminosity. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at Solar System scales and exclude brown-dwarf companions more massive than 20 M_jup beyond separations of ~2.5 au and giant planets more massive than 2 M_jup beyond 9 au.

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.

The advent of Wide Field Adaptive Optics (WFAO) systems marks the beginning of a new era in high spatial resolution imaging. The newly commissioned Gemini South Multi-Conjugate Adaptive Optics System (GeMS) combined with the infrared camera Gemini South Adaptive Optics Imager (GSAOI), delivers quasi diffraction-limited images over a field of 2 arc-minutes across. However, despite this excellent performance, some variable residues still limit the quality of the analyses. In particular, distortions severely affect GSAOI and become a critical issue for high-precision astrometry and photometry. In this paper, we investigate an optimal way to correct for the distortion following an inverse problem approach. Formalism as well as applications on GeMS data are presented.

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

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.

COupled SLope and scIntillation Detection And Ranging (CO-SLIDAR) is a recent profiling method of the vertical distribution of atmospheric turbulence strength ($C_n^2$ profile). It takes advantage of correlations of slopes and of scintillation, both measured with a Shack-Hartmann wavefront sensor on a binary star. In this paper, we present the improved CO-SLIDAR reconstruction method of the $C_n^2$ profile and the first on-sky results of the CO-SLIDAR profiler. We examine CO-SLIDAR latest performance in simulation, taking into account the detection noise bias and estimating error bars along with the turbulence profile. The estimated $C_n^2$ profiles demonstrate the accuracy of the CO-SLIDAR method, showing sensitivity to both low and high altitude turbulent layers. CO-SLIDAR is tested on-sky for the first time, on the $1.5$~m MeO (M\'etrologie Optique) telescope at Observatoire de la Cote d'Azur (France). The reconstructed profiles are compared to turbulence profiles estimated from meteorological data and a good agreement is found. We discuss CO-SLIDAR's contribution in the $C_n^2$ profilers' landscape and we propose some improvements of the instrument.

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)

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.

High-contrast imaging for the detection and characterization of exoplanets relies on the instrument's capability to block out the light of the host star. Some current post-processing methods for calibrating out the residual speckles use information redundancy offered by multispectral imaging but do not use any prior information on the origin of these speckles. We investigate whether additional information on the system and image formation process can be used to more finely exploit the multispectral information. We developed an inversion method in a Bayesian framework that is based on an analytical imaging model to estimate both the speckles and the object map. The model links the instrumental aberrations to the speckle pattern in the image focal plane, distinguishing between aberrations upstream and downstream of the coronagraph. We propose and validate several numerical techniques to handle the difficult minimization problems of phase retrieval and achieve a contrast of 10^6 at 0.2 arcsec from simulated images, in the presence of photon noise. This opens up the the possibility of tests on real data where the ultimate performance may override the current techniques if the instrument has good and stable coronagraphic imaging quality. This paves the way for new astrophysical exploitations or even new designs for future instruments.

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

We present the last results of our PSF reconstruction (PSF-R) project for the Keck-II and Gemini-North AO systems in natural guide star mode. Our initial tests have shown that the most critical aspects of PSF-R are the determination of the system static aberrations and the optical turbulence parameters, and we have set up a specific observation campaign on the two systems to explore this. We demonstrate that deformable mirror based seeing monitor works well, and 10% accuracy is easily obtained. Phase diversity has been demonstrated to work on sky sources. Besides, residual phase stationarity is an important assumption in PSF-R, and we demonstrate here that it is basically true. As a result of these tests and verifications, we have been able for the first time to obtain a very good PSF reconstruction for the Keck-II system, in bright natural guide star mode.

Cn2 profile monitoring usually makes use of wavefront slope correlations or of scintillation pattern correlations. Wavefront slope correlations provide sensitivity to layers close to the receiving plane. In addition, scintillation correlations allow a better sensitivity to high turbulence layers. Wavefront slope and scintillation correlations are therefore complementary. Slopes and scintillation being recorded simultaneously with a Shack-Hartmann wavefront sensor (SHWFS), we propose here to exploit their correlation to retrieve the Cn2 profile. The measurement method named COupled SLodar scIDAR (CO-SLIDAR) uses correlations of SHWFS data from two separated stars. A maximum-likelihood method is developed to estimate precisely the positions and intensities corresponding to each SHWFS spot, which are used as inputs for CO-SLIDAR. First results are presented using SHWFS real data from a binary star.

The galactic central black hole Sgr A* exhibits outbursts of radiation in the near infrared (so-called IR flares). One model of these events consists in a hotspot orbiting on the innermost stable circular orbit (ISCO) of the hole. These outbursts can be used as a probe of the central gravitational potential. One main scientific goal of the second generation VLTI instrument GRAVITY is to observe these flares astrometrically. Here, the astrometric precision of GRAVITY is investigated in imaging mode, which consists in analysing the image computed from the interferometric data. The capability of the instrument to put in light the motion of a hotspot orbiting on the ISCO of our central black hole is then discussed. We find that GRAVITY's astrometric precision for a single star in imaging mode is smaller than the Schwarzschild radius of Sgr A*. The instrument can also demonstrate that a body orbiting on the last stable orbit of the black hole is indeed moving. It yields a typical size of the orbit, if the source is as bright as m_K=14. These results show that GRAVITY allows one to study the close environment of Sgr A*. Having access to the ISCO of the central massive black hole probably allows constraining general relativity in its strong regime. Moreover, if the hotspot model is appropriate, the black hole spin can be constrained.

We present the results of the fourth Optical/IR Interferometry Imaging Beauty Contest. The contest consists of blind imaging of test data sets derived from model sources and distributed in the OI-FITS format. The test data consists of spectral data sets on an object "observed" in the infrared with spectral resolution. There were 4 different algorithms competing this time: BSMEM the Bispectrum Maximum Entropy Method by Young, Baron & Buscher; RPR the Recursive Phase Reconstruction by Rengaswamy; SQUEEZE a Markov Chain Monte Carlo algorithm by Baron, Monnier & Kloppenborg; and, WISARD the Weak-phase Interferometric Sample Alternating Reconstruction Device by Vannier & Mugnier. The contest model image, the data delivered to the contestants and the rules are described as well as the results of the image reconstruction obtained by each method. These results are discussed as well as the strengths and limitations of each algorithm.

This paper reports on H-band interferometric observations of Betelgeuse made at the three-telescope interferometer IOTA. We image Betelgeuse and its asymmetries to understand the spatial variation of the photosphere, including its diameter, limb darkening, effective temperature, surrounding brightness, and bright (or dark) star spots. We used different theoretical simulations of the photosphere and dusty environment to model the visibility data. We made images with parametric modeling and two image reconstruction algorithms: MIRA and WISARD. We measure an average limb-darkened diameter of 44.28 +/- 0.15 mas with linear and quadratic models and a Rosseland diameter of 45.03 +/- 0.12 mas with a MARCS model. These measurements lead us to derive an updated effective temperature of 3600 +/- 66 K. We detect a fully-resolved environment to which the silicate dust shell is likely to contribute. By using two imaging reconstruction algorithms, we unveiled two bright spots on the surface of Betelgeuse. One spot has a diameter of about 11 mas and accounts for about 8.5% of the total flux. The second one is unresolved (diameter < 9 mas) with 4.5% of the total flux. Resolved images of Betelgeuse in the H band are asymmetric at the level of a few percent. The MOLsphere is not detected in this wavelength range. The amount of measured limb-darkening is in good agreement with model predictions. The two spots imaged at the surface of the star are potential signatures of convective cells.

Since the discovery of its dusty disk in 1984, Beta Pictoris has become the prototype of young early-type planetary systems, and there are now various indications that a massive Jovian planet is orbiting the star at ~ 10 AU. However, no planets have been detected around this star so far. Our goal was to investigate the close environment of Beta Pic, searching for planetary companion(s). Deep adaptive-optics L'-band images of Beta Pic were recorded using the NaCo instrument at the Very Large Telescope. A faint point-like signal is detected at a projected distance of ~ 8 AU from the star, within the North-East side of the dust disk. Various tests were made to rule out with a good confidence level possible instrumental or atmospheric artifacts. The probability of a foreground or background contaminant is extremely low, based in addition on the analysis of previous deep Hubble Space Telescope images. The object L'=11.2 apparent magnitude would indicate a typical temperature of ~1500 K and a mass of ~ 8 Jovian masses. If confirmed, it could explain the main morphological and dynamical peculiarities of the Beta Pic system. The present detection is unique among A-stars by the proximity of the resolved planet to its parent star. Its closeness and location inside the Beta Pic disk suggest a formation process by core accretion or disk instabilities rather than a binary-like formation process.

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.

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

This paper is an H band interferometric examination of Arcturus, a star frequently used as a spatial and spectral calibrator. Using the IOTA 3 telescope interferometer, we performed spectro-interferometric observation (R~35) of Arcturus. Atmospheric models and prescriptions were fitted to the data to derive the brightness distribution of the photosphere. Image reconstruction was also obtained using two software algorithms: Wisard and Mira. An achromatic power law proved to be a good model of the brightness distribution, with a limb darkening compatible with the one derived from atmospheric model simulations using our Marcs model. A Rosseland diameter of 21.05+/-0.21 was derived, corresponding to an effective temperature of T_eff=4295+/-26 K. No companion was detected from the closure phases, with an upper limit on the brightness ratio of 8e-4 at 1AU. Dynamic range at such distance from the photosphere was established at 1.5e-4 (1sigma rms). An upper limit of 1.7e-3 was also derived for the level of brightness asymmetries present on the photosphere.

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.

We present K, L and M diffraction-limited images of NGC 1068 obtained with NAOS+CONICA at VLT/YEPUN over a 3.5" field around the central engine. Hot dust (Tcol = 550-650 K) is found in three different regions : (a) in the true nucleus, seen as a slightly NS elongated, core of extremely hot dust, "resolved" in K and L with respective diameters of ~5 pc and 8.5 pc ; (b) along the NS direction, as a "spiral arm" and a southern tongue ; (c) as a set of parallel elongated nodules ("wave-like") bracketting the jet. Several structures observed on radio maps, mid-IR or HST UV-visible maps are seen, so that a precise registration can be done from UV to 6 cm. These results do support the current interpretion that source (a) corresponds to emission from dust near sublimation temperature delimiting the walls of the cavity in the central obscuring torus. Structure (b) is thought to be a mixture of hot dust and active star forming regions along a micro spiral structure that could trace the tidal mechanism bringing matter to the central engine. Structure c)which was not known, exhibits too high a temperature for "classical'' grains ; it is most probably the signature of transiently heated very small dust grains (VSG) : "nano-diamonds", which are resistant and can form in strong UV field or in shocks, are very attractive candidates. The "waves'' can be condensations triggered by jet induced shocks, as predicted by recent models. First estimates, based on a simple VSG model and on a detailed radiative transfer model, do agree with those interpretations, both qualitatively and quantitatively.