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.

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.

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.

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.

The performance of tomographic adaptive optics systems is intrinsically linked to the vertical profile of optical turbulence. Firstly, a sufficient number of discrete turbulent layers must be reconstructed to model the true continuous turbulence profile. Secondly over the course of an observation, the profile as seen by the telescope changes and the tomographic reconstructor must be updated. These changes can be due to the unpredictable evolution of turbulent layers on meteorological timescales as short as minutes. Here we investigate the effect of changing atmospheric conditions on the quality of tomographic reconstruction by coupling fast analytical adaptive optics simulation to a large database of 10 691 high resolution turbulence profiles measured over two years by the Stereo-SCIDAR instrument at ESO Paranal, Chile. This work represents the first investigation of these effects with a large, statistically significant sample of turbulence profiles. The statistical nature of the study allows us to assess not only the degradation and variability in tomographic error with a set of system parameters (e.g. number of layers, temporal update period) but also the required parameters to meet some error threshold. In the most challenging conditions where the profile is rapidly changing, these parameters must be far more tightly constrained in order to meet this threshold. By providing estimates of these constraints for a wide range of system geometries as well as the impact of different temporal optimisation strategies we may assist the designers of tomographic AO for the ELT to dimension their systems.

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

Here we describe a simple, efficient, and most importantly fully operational point-spread-function(PSF)-reconstruction approach for laser-assisted ground layer adaptive optics (GLAO) in the frame of the Multi Unit Spectroscopic Explorer (MUSE) Wide Field Mode. Based on clear astrophysical requirements derived by the MUSE team and using the functionality of the current ESO Adaptive Optics Facility we aim to develop an operational PSF-reconstruction (PSFR) algorithm and test it both in simulations and using on-sky data. The PSFR approach is based on a Fourier description of the GLAO correction to which the specific instrumental effects of MUSE Wide Field Mode (pixel size, internal aberrations, etc.) have been added. It was first thoroughly validated with full end-to-end simulations. Sensitivity to the main atmospheric and AO system parameters was analysed and the code was re-optimised to account for the sensitivity found. Finally, the optimised algorithm was tested and commissioned using more than one year of on-sky MUSE data. We demonstrate with an on-sky data analysis that our algorithm meets all the requirements imposed by the MUSE scientists, namely an accuracy better than a few percent on the critical PSF parameters including full width at half maximum and global PSF shape through the kurtosis parameter of a Moffat function. The PSFR algorithm is publicly available and is used routinely to assess the MUSE image quality for each observation. It can be included in any post-processing activity which requires knowledge of the PSF.

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

With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

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.

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.

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

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.

As part of the National Science Foundation funded "Gemini in the Era of MultiMessenger Astronomy" (GEMMA) program, Gemini Observatory is developing GNAO, a widefield adaptive optics (AO) facility for Gemini-North on Maunakea, the only 8m-class open-access telescope available to the US astronomers in the northern hemisphere. GNAO will provide the user community with a queue-operated Multi-Conjugate AO (MCAO) system, enabling a wide range of innovative solar system, Galactic, and extragalactic science with a particular focus on synergies with JWST in the area of time-domain astronomy. The GNAO effort builds on institutional investment and experience with the more limited block-scheduled Gemini Multi-Conjugate System (GeMS), commissioned at Gemini South in 2013. The project involves close partnerships with the community through the recently established Gemini AO Working Group and the GNAO Science Team, as well as external instrument teams. The modular design of GNAO will enable a planned upgrade to a Ground Layer AO (GLAO) mode when combined with an Adaptive Secondary Mirror (ASM). By enhancing the natural seeing by an expected factor of two, GLAO will vastly improve Gemini North's observing efficiency for seeing-limited instruments and strengthen its survey capabilities for multi-messenger astronomy.

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.

HD163296 is a Herbig Ae/Be star known to host a protoplanetary disk with a ringed structure. To explain the disk features, previous works proposed the presence of planets embedded into the disk. We have observed HD163296 with the near-infrared (NIR) branch of SPHERE composed by IRDIS and IFS with the aim to put tight constraints on the presence of substellar companions around this star. Despite the low rotation of the field of view during our observation we were able to put upper mass limits of few M_Jup around this object. These limits do not allow to give any definitive conclusion about the planets proposed through the disk characteristics. On the other hand, our results seem to exclude the presence of the only candidate proposed until now using direct imaging in the NIR even if some caution has to be taken considered the different wavelength bands of the two observations.

For extremely large telescopes, adaptive optics will be required to correct the Earth's turbulent atmosphere. The performance of tomographic adaptive optics is strongly dependent on the vertical distribution (profile) of this turbulence. An important way in which this manifests is the tomographic error, arising from imperfect measurement and reconstruction of the turbulent phase at altitude. Conventionally, a small number of reference profiles are used to obtain this error in simulation however these profiles are not constructed to be representative in terms of tomographic error. It is therefore unknown whether these simulations are providing realistic performance estimates. Here, we employ analytical adaptive optics simulation that drastically reduces computation times to compute tomographic error for 10 691 measurements of the turbulence profile gathered by the Stereo-SCIDAR instrument at ESO Paranal. We assess for the first time the impact of the profile on tomographic error in a statistical manner. We find, in agreement with previous work, that the tomographic error is most directly linked with the distribution of turbulence into discrete, stratified layers. Reference profiles are found to provide mostly higher tomographic error than expected, which we attribute to the fact that these profiles are primarily composed of averages of many measurements resulting in unrealistic, continuous distributions of turbulence. We propose that a representative profile should be defined with respect to a particular system, and that as such simulations with a large statistical sample of profiles must be an important step in the design process.

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.

In order to enhance accuracy of astrophysical estimates obtained on Adaptive-optics (AO) images, such as photometry and astrometry, we investigate a new concept to constrain the Point Spread Function (PSF) model called PSF Reconstruction and Identification for Multi-sources characterization Enhancement (PRIME), that handles jointly the science image and the AO control loop data. We present in this paper the concept of PRIME and validate it on Keck II telescope NIRC2 images. We show that by calibrating the PSF model over the scientific image, PSF reconstruction achieves 1\% and 3 mas of accuracy on respectively the Strehl-ratio and the PSF full width at half maximum. We show on NIRC2 binary images that PRIME is sufficiently robust to noise to retain photometry and astrometry precision below 0.005 mag and 100$\mu$as on a $m_H=$ 14 mag object. Finally, we also validate that PRIME performs a PSF calibration on the triple system Gl569BAB which provides a separation of 66.73$\pm 1.02$ and a differential photometry of 0.538$\pm 0.048$, compared to the reference values obtained with the extracted PSF which are 66.76 mas $\pm$ 0.94 and 0.532 mag $\pm$ 0.041.

Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE as inputs. We analyzed the SPHERE images to infer the asteroid's global shape and the morphology of its main craters. We present the volume-equivalent diameter D$_{{\rm eq}}$=214$\pm$5 km, and bulk density $\rho$=2.7$\pm$0.3 g cm$^{-3}$ of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20--40 km in diameter detected at several epochs, which we interpret as impact craters. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. The bulk density of Iris is consistent with that of its meteoritic analog, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris's shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities.

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.

Several telescopes include large Deformable Mirrors (DM) located directly inside the telescope. These adaptive telescopes trigger new constraints for the calibration of the Adaptive Optics (AO) systems as they usually offer no access to an artificial calibration source for the interaction matrix measurement. Moreover, the optical propagation between the DM and the Wave-Front Sensor (WFS) may evolve during the operation, resulting in mis-registrations that highly affect the AO performance and thus the scientific observation. They have to be measured and compensated, for instance by updating the calibration. A new strategy consists of estimating the mis-registrations and injecting them into synthetic models to generate noise-free interaction matrices. This pseudo-synthetic approach is the baseline for the Adaptive Optics Facility working with a Shack-Hartmann WFS and seems particularly suited for the future Extremely Large Telescope as the calibration will have to be regularly updated, for a large numbers of actuators. In this paper, the feasibility of a pseudo synthetic calibration with Pyramid WFS at the Large Binocular Telescope (LBT) is investigated. A synthetic model of the LBT AO systems is developed, and the procedure to adjust the mis-registrations parameters is introduced, extracting them from an experimental interaction matrix. We successfully tested an interaction matrix generated from the model on the real system in high-order AO mode. We recorded a slightly better performance with respect to the experimental one. This work demonstrates that a high accuracy calibration can be obtained using the pseudo synthetic approach with pyramid WFS.

Context. CM-like asteroids (Ch and Cgh classes) are a major population within the broader C-complex, encompassing about 10% of the mass of the main asteroid belt. Their internal structure has been predicted to be homogeneous, based on their compositional similarity as inferred from spectroscopy (Vernazza et al., 2016, AJ 152, 154) and numerical modeling of their early thermal evolution (Bland & Travis, 2017, Sci. Adv. 3, e1602514). Aims. Here we aim to test this hypothesis by deriving the density of the CM-like asteroid (41) Daphne from detailed modeling of its shape and the orbit of its small satellite. Methods. We observed Daphne and its satellite within our imaging survey with the Very Large Telescope extreme adaptive-optics SPHERE/ZIMPOL camera (ID 199.C-0074, PI P. Vernazza) and complemented this data set with earlier Keck/NIRC2 and VLT/NACO observations. We analyzed the dynamics of the satellite with our Genoid meta-heuristic algorithm. Combining our high-angular resolution images with optical lightcurves and stellar occultations, we determine the spin period, orientation, and 3-D shape, using our ADAM shape modeling algorithm. Results. The satellite orbits Daphne on an equatorial, quasi-circular, prograde orbit, like the satellites of many other large main-belt asteroids. The shape model of Daphne reveals several large flat areas that could be large impact craters. The mass determined from this orbit combined with the volume computed from the shape model implies a density for Daphne of 1.77+/-0.26 g/cm3 (3 {\sigma}). This density is consistent with a primordial CM-like homogeneous internal structure with some level of macroporosity (~17%). Conclusions. Based on our analysis of the density of Daphne and 75 other Ch/Cgh-type asteroids gathered from the literature, we conclude that the primordial internal structure of the CM parent bodies was homogeneous.

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.

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

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.

Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche's three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. We derived a bulk density of 3.99\,$\pm$\,0.26\,g$\cdot$cm$^{-3}$ for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites ($\sim$7.8\,g$\cdot$cm$^{-3}$), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density $\sim$4.25\,$\cdot$cm$^{-3}$). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460\,$\pm$\,200}\,m at 150\,km from Psyche (0.2\%\,$\times$\,R$_{Hill}$, the Hill radius) and 800\,$\pm$\,200\,m at 2,000\,km (3\%\,$\times$\,$R_{Hill}$). Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites.

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.

Knowledge of the atmospheric turbulence in the telescope line-of-sight is crucial for wide-field observations assisted by adaptive optics (AO), for which the Point Spread Function (PSF) becomes strongly elongated due to the anisoplanatism effect. This one must be modelled accurately to extrapolate the PSF anywhere across the Field of view (FOV) and improve the science exploitation. However, anisoplanatism is a function of the Cn2(h) profile, that is not directly accessible from single conjugate AO telemetry. One may rely on external profilers, but recent studies have highlighted more than 10% of discrepancies with AO internal measurements, while we aim better than 1% of accuracy for PSF modelling. To tackle this existing limitation, we present the Focal plane profiling (FPP), as a $C_n^2(h)$ profiling method that relies on post-AO focal plane images. We demonstrate such an approach complies with a 1%-level of accuracy on the $C_n^2(h)$ estimation and establish how this accuracy varies regarding the calibration stars magnitudes and positions in the field. We highlight that photometry and astrometry errors due to PSF mis-modelling reaches respectively 1% and 50{\mu}as using FPP on a Keck baseline, with a preliminary calibration using a star of magnitude H=14 at 20". We also validate this concept using Canadas NRC-Herzberg HENOS testbed images in comparing FPP retrieval with alternative $C_n^2(h)$ measurements on HeNOS. The FPP approach allows to profile the $C_n^2(h)$ using the SCAO systems and improve significantly the PSF characterisation. Such a methodology is also ELT-size compliant and will be extrapolated to tomographic systems in a near future.

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

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

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

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.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

When combined with the huge collecting area of the ELT, MOSAIC will be the most effective and flexible Multi-Object Spectrograph (MOS) facility in the world, having both a high multiplex and a multi-Integral Field Unit (Multi-IFU) capability. It will be the fastest way to spectroscopically follow-up the faintest sources, probing the reionisation epoch, as well as evaluating the evolution of the dwarf mass function over most of the age of the Universe. MOSAIC will be world-leading in generating an inventory of both the dark matter (from realistic rotation curves with MOAO fed NIR IFUs) and the cool to warm-hot gas phases in z=3.5 galactic haloes (with visible wavelenth IFUs). Galactic archaeology and the first massive black holes are additional targets for which MOSAIC will also be revolutionary. MOAO and accurate sky subtraction with fibres have now been demonstrated on sky, removing all low Technical Readiness Level (TRL) items from the instrument. A prompt implementation of MOSAIC is feasible, and indeed could increase the robustness and reduce risk on the ELT, since it does not require diffraction limited adaptive optics performance. Science programmes and survey strategies are currently being investigated by the Consortium, which is also hoping to welcome a few new partners in the next two years.

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.

Tomographic wave-front reconstruction is the main computational bottleneck to realize real-time correction for turbulence-induced wave-front aberrations in future laser-assisted tomographic adaptive-optics (AO) systems for ground-based Giant Segmented Mirror Telescopes (GSMT), because of its unprecedented number of degrees of freedom, $N$, i.e. the number of measurements from wave-front sensors (WFS). In this paper, we provide an efficient implementation of the minimum-mean-square error (MMSE) tomographic wave-front reconstruction mainly useful for some classes of AO systems not requiring a multi-conjugation, such as laser-tomographic AO (LTAO), multi-object AO (MOAO) and ground-layer AO (GLAO) systems, but also applicable to multi-conjugate AO (MCAO) systems. This work expands that by R. Conan [ProcSPIE, 9148, 91480R (2014)] to the multi-wave-front, tomographic case using natural and laser guide stars. The new implementation exploits the Toeplitz structure of covariance matrices used in a MMSE reconstructor, which leads to an overall $O(N\log N)$ real-time complexity compared to $O(N^2)$ of the original implementation using straight vector-matrix multiplication. We show that the Toeplitz-based algorithm leads to 60\,nm rms wave-front error improvement for the European Extremely Large Telescope Laser-Tomography AO system over a well-known sparse-based tomographic reconstruction, but the number of iterations required for suitable performance is still beyond what a real-time system can accommodate to keep up with the time-varying turbulence

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.

Adaptive optics (AO) restore the angular resolution of ground-based telescopes, but at the cost of delivering a time- and space-varying point spread function (PSF) with a complex shape. PSF knowledge is crucial for breaking existing limits on the measured accuracy of photometry and astrometry in science observations. In this paper, we concentrate our analyses on anisoplanatism signature only onto PSF: for large-field observations (20") with single- conjugated AO, PSFs are strongly elongated due to anisoplanatism that manifests itself as three different terms for Laser-guide star (LGS) systems: angular, focal and tilt. We propose a generalized model that relies on a point-wise decomposition of the phase and encompasses the non-stationarity of LGS systems. We demonstrate it is more accurate and less computationally demanding than existing models: it agrees with end-to-end physical-optics simulations to within 0.1% of PSF measurables, such as Strehl-ratio, FWHM and fraction of variance unexplained.Secondly, we study off-axis PSF modelling is with respect to $C_n^2(h)$ profile (heights and fractional weights). For 10m class telescope, PSF morphology is estimated at 1% level as long as we model the atmosphere with at least 7 layers whose heights and weights are known respectively with 200m and 10% precision. As a verification test we used the Canada's NRC-Herzberg HeNOS testbed data, featuring four lasers. We highlight capability of retrieving off-axis PSF characteristics within 10% of fraction of variance unexplained, which complies with the expected range from the sensitivity analysis. Our new off-axis PSF modelling method lays the ground-work for testing on-sky in the near future.

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.

We present observations of the Herbig Ae star HD169142 with VLT/SPHERE instruments InfraRed Dual-band Imager and Spectrograph (IRDIS) ($K1K2$ and $H2H3$ bands) and the Integral Field Spectrograph (IFS) ($Y$, $J$ and $H$ bands). We detect several bright blobs at $\sim$180 mas separation from the star, and a faint arc-like structure in the IFS data. Our reference differential imaging (RDI) data analysis also finds a bright ring at the same separation. We show, using a simulation based on polarized light data, that these blobs are actually part of the ring at 180 mas. These results demonstrate that the earlier detections of blobs in the $H$ and $K_S$ bands at these separations in Biller et al. as potential planet/substellar companions are actually tracing a bright ring with a Keplerian motion. Moreover, we detect in the images an additional bright structure at $\sim$93 mas separation and position angle of 355$^{\circ}$, at a location very close to previous detections. It appears point-like in the $YJ$ and $K$ bands but is more extended in the $H$ band. We also marginally detect an inner ring in the RDI data at $\sim$100 mas. Follow-up observations are necessary to confirm the detection and the nature of this source and structure.

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

Context. The high-angular-resolution capability of the new-generation ground-based adaptive-optics camera SPHERE at ESO VLT allows us to assess, for the very first time, the cratering record of medium-sized (D~100-200 km) asteroids from the ground, opening the prospect of a new era of investigation of the asteroid belt's collisional history. Aims. We investigate here the collisional history of asteroid (6) Hebe and challenge the idea that Hebe may be the parent body of ordinary H chondrites, the most common type of meteorites found on Earth (~34% of the falls). Methods. We observed Hebe with SPHERE as part of the science verification of the instrument. Combined with earlier adaptive-optics images and optical light curves, we model the spin and three-dimensional (3D) shape of Hebe and check the consistency of the derived model against available stellar occultations and thermal measurements. Results. Our 3D shape model fits the images with sub-pixel residuals and the light curves to 0.02 mag. The rotation period (7.274 47 h), spin (343 deg,+47 deg), and volume-equivalent diameter (193 +/- 6km) are consistent with previous determinations and thermophysical modeling. Hebe's inferred density is 3.48 +/- 0.64 g.cm-3 , in agreement with an intact interior based on its H-chondrite composition. Using the 3D shape model to derive the volume of the largest depression (likely impact crater), it appears that the latter is significantly smaller than the total volume of close-by S-type H-chondrite-like asteroid families. Conclusions. Our results imply that (6) Hebe is not the most likely source of H chondrites. Over the coming years, our team will collect similar high-precision shape measurements with VLT/SPHERE for ~40 asteroids covering the main compositional classes, thus providing an unprecedented dataset to investigate the origin and collisional evolution of the asteroid belt.

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.

Context. The SAO 206462 (HD 135344B) disk is one of the few known transitional disks showing asymmetric features in scattered light and thermal emission. Near-infrared scattered-light images revealed two bright outer spiral arms and an inner cavity depleted in dust. Giant protoplanets have been proposed to account for the disk morphology. Aims. We aim to search for giant planets responsible for the disk features and, in the case of non-detection, to constrain recent planet predictions using the data detection limits. Methods. We obtained new high-contrast and high-resolution total intensity images of the target spanning the Y to the K bands (0.95-2.3 mic) using the VLT/SPHERE near-infrared camera and integral field spectrometer. Results. The spiral arms and the outer cavity edge are revealed at high resolutions and sensitivities without the need for image post-processing techniques, which introduce photometric biases. We do not detect any close-in companions. For the derivation of the detection limits on putative giant planets embedded in the disk, we show that the knowledge of the disk aspect ratio and viscosity is critical for the estimation of the attenuation of a planet signal by the protoplanetary dust because of the gaps that these putative planets may open. Given assumptions on these parameters, the mass limits can vary from ~2-5 to ~4-7 Jupiter masses at separations beyond the disk spiral arms. The SPHERE detection limits are more stringent than those derived from archival NaCo/L' data and provide new constraints on a few recent predictions of massive planets (4-15 MJ) based on the spiral density wave theory. The SPHERE and ALMA data do not favor the hypotheses on massive giant planets in the outer disk (beyond 0.6). There could still be low-mass planets in the outer disk and/or planets inside the cavity.

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.

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

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

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.

Canary is the Multi-Object Adaptive Optics (MOAO) pathfinder for the future MOAO-assisted Integral-Field Units (IFU) proposed for Extremely Large Telescopes (ELT). The MOAO concept relies on tomographically reconstructing the turbulence using multiple measurements along different lines of sight. Tomography requires the knowledge of the statistical turbulence parameters, commonly recovered from the system telemetry using a dedicated profiling technique. For demonstration purposes with the MOAO pathfinder Canary , this identification is performed thanks to the Learn & Apply (L&A) algorithm, that consists in model- fitting the covariance matrix of WFS measurements dependent on relevant parameters: $C_n^2(h)$ profile, outer scale profile and system mis-registration. We explore an upgrade of this algorithm, the Learn 3 Steps (L3S) approach, that allows one to dissociate the identification of the altitude layers from the ground in order to mitigate the lack of convergence of the required empirical covariance matrices therefore reducing the required length of data time-series for reaching a given accuracy. For nominal observation conditions, the L3S can reach the same level of tomographic error in using five times less data frames than the L&A approach. The L3S technique has been applied over a large amount of Canary data to characterize the turbu- lence above the William Herschel Telescope (WHT). These data have been acquired the 13th, 15th, 16th, 17th and 18th September 2013 and we find 0.67"/8.9m/3.07m/s of total seeing/outer scale/wind-speed, with 0.552"/9.2m/2.89m/s below 1.5 km and 0.263"/10.3m/5.22m/s between 1.5 and 20 km. We have also de- termined the high altitude layers above 20 km, missed by the tomographic reconstruction on Canary , have a median seeing of 0.187" and have occurred 16% of observation time.

In preparation of future Multi-Object Spectrographs (MOS) whose one of the major role is to provide an extensive statistical studies of high redshifted galaxies surveyed, the demonstrator Canary has been designed to tackle technical challenges related to open-loop Adaptive-Optics (AO) control with jointed Natural Guide Star (NGS) and Laser Guide Star (LGS) tomography. We have developed a Point Spread Function (PSF)-Reconstruction algorithm dedicated to MOAO systems using system telemetry to estimate the PSF potentially anywhere in the observed field, a prerequisite to post- process AO-corrected observations in Integral Field Spectroscopy (IFS). In this paper we show how to handle off-axis data to estimate the PSF using atmospheric tomography and compare it to a classical approach that uses on-axis residual phase from a truth sensor observing a natural bright source. We have reconstructed over 450 on-sky Canary PSFs and we get bias/1-${\sigma}$ standard-deviation (std) of 1.3/4.8 on the H-band Strehl ratio (SR) with 92.3% of correlation between reconstructed and sky SR. On the Full Width at Half Maximum (FWHM), we get respectively 2.94 mas, 19.9 mas and 88.3% for the bias, std and correlation. The reference method achieves 0.4/3.5/95% on the SR and 2.71 mas/14.9 mas/92.5% on the FWHM for the bias/std/correlation.

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.

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.

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.

Aims. The formation and properties of star clusters at the edge of H II regions are poorly known, partly due to limitations in angular resolution and sensitivity, which become particularly critical when dealing with extragalactic clusters. In this paper we study the stellar content and star-formation processes in the young N159W region in the Large Magellanic Cloud. Methods. We investigate the star-forming sites in N159W at unprecedented spatial resolution using JHKs-band images obtained with the GeMS/GSAOI instrument on the Gemini South telescope. The typical angular resolution of the images is of 100 mas, with a limiting magnitude in H of 22 mag (90 percent completeness). Photometry from our images is used to identify candidate young stellar objects (YSOs) in N159W. We also determine the H-band luminosity function of the star cluster at the centre of the H II region and use this to estimate its initial mass function (IMF). Results. We estimate an age of 2 + or - 1 Myr for the central cluster, with its IMF described by a power-law with an index of gamma = - 1.05 + or - 0.2 , and with a total estimated mass of 1300 solar mass. We also identify 104 candidate YSOs, which are concentrated in clumps and subclusters of stars, principally at the edges of the H II region. These clusters display signs of recent and active star-formation such as ultra-compact H II regions, and molecular outflows. This suggests that the YSOs are typically younger than the central cluster, pointing to sequential star-formation in N159W, which has probably been influenced by interactions with the expanding H II bubble.

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

Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time control systems of two wide-field AO systems. These are CANARY, which has been operating on-sky since 2010, and DRAGON, which is a laboratory adaptive optics real-time demonstrator instrument. We detail the issues and difficulties that arose, along with the solutions we developed. We also provide recommendations for consideration when developing future wide-field AO systems.

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

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.

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.

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]

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 present HSIM: a dedicated pipeline for simulating observations with the HARMONI integral field spectrograph on the European Extremely Large Telescope. HSIM takes high spectral and spatial resolution input data-cubes, encoding physical descriptions of astrophysical sources, and generates mock observed data-cubes. The simulations incorporate detailed models of the sky, telescope and instrument to produce realistic mock data. Further, we employ a new method of incorporating the strongly wavelength dependent adaptive optics point spread functions. HSIM provides a step beyond traditional exposure time calculators and allows us to both predict the feasibility of a given observing programme with HARMONI, as well as perform instrument design trade-offs. In this paper we concentrate on quantitative measures of the feasibility of planned observations. We give a detailed description of HSIM and present two studies: estimates of point source sensitivities along with simulations of star-forming emission-line galaxies at $z\sim 2-3$. We show that HARMONI will provide exquisite resolved spectroscopy of these objects on sub-kpc scales, probing and deriving properties of individual star-forming regions.

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 investigate the star formation activity in a young star forming cluster embedded at the edge of the RCW 41 HII region. As a complementary goal, we aim at demonstrating the gain provided by Wide-Field Adaptive Optics instruments to study young clusters. We used deep, JHKs images from the newly commissioned Gemini-GeMS/GSAOI instrument, complemented with Spitzer IRAC observations, in order to study the photometric properties of the young stellar cluster. GeMS is an AO instrument, delivering almost diffraction limited images over a field of 2' across. The exquisite angular resolution allows us to reach a limiting magnitude of J = 22 for 98% completeness. The combination of the IRAC photometry with our JHKs catalog is used to build color-color diagrams, and select Young Stellar Objects (YSOs) candidates. We detect the presence of 80 Young Stellar Object (YSO) candidates. Those YSOs are used to infer the cluster age, which is found to be in the range 1 to 5 Myr. We find that 1/3 of the YSOs are in a range between 3 to 5 Myr, while 2/3 of the YSO are < 3 Myr. When looking at the spatial distribution of these two populations, we evidence a potential age gradient across the field, suggesting sequential star formation. We construct the IMF, and show that we can sample the mass distribution well into the brown dwarf regime (down to 0.01 Msun). The logarithmic mass function rises to peak at 0.3 Msun, before turning over and declining into the brown dwarf regime. The total cluster mass derived is estimated to be 78 +/- 18 Msun, while the ratio of brown dwarfs to star derived is 18 p/- 5 %. When comparing with other young clusters, we find that the IMF shape of the young cluster embedded within RCW 41 is consistent with those of Trapezium, IC 348 or Chamaeleon I, except for the IMF peak, which happens to be at higher mass. This characteristic is also seen in clusters like NGC 6611 or even Taurus.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

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.

We propose a new algorithm for an adaptive optics system control law, based on the Linear Quadratic Gaussian approach and a Kalman Filter adaptation with localizations. It allows to handle non-stationary behaviors, to obtain performance close to the optimality defined with the residual phase variance minimization criterion, and to reduce the computational burden with an intrinsically parallel implementation on the Extremely Large Telescopes (ELTs).

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

We propose a new algorithm for an adaptive optics system control law which allows to reduce the computational burden in the case of an Extremely Large Telescope (ELT) and to deal with non-stationary behaviors of the turbulence. This approach, using Ensemble Transform Kalman Filter and localizations by domain decomposition is called the local ETKF: the pupil of the telescope is split up into various local domains and calculations for the update estimate of the turbulent phase on each domain are performed independently. This data assimilation scheme enables parallel computation of markedly less data during this update step. This adapts the Kalman Filter to large scale systems with a non-stationary turbulence model when the explicit storage and manipulation of extremely large covariance matrices are impossible. First simulation results are given in order to assess the theoretical analysis and to demonstrate the potentiality of this new control law for complex adaptive optics systems on ELTs.

Context. Several exoplanet direct imaging instruments will soon be in operation. They use an extreme adaptive optics (XAO) system to correct the atmospheric turbulence and provide a highly-corrected beam to a near-infrared (NIR) coronagraph for starlight suppression. The performance of the coronagraph is however limited by the non-common path aberrations (NCPA) due to the differential wavefront errors existing between the visible XAO sensing path and the NIR science path, leading to residual speckles in the coronagraphic image. Aims. Several approaches have been developed in the past few years to accurately calibrate the NCPA, correct the quasi-static speckles and allow the observation of exoplanets at least 1e6 fainter than their host star. We propose an approach based on the Zernike phase-contrast method for the measurements of the NCPA between the optical path seen by the visible XAO wavefront sensor and that seen by the near-IR coronagraph. Methods. This approach uses a focal plane phase mask of size {\lambda}/D, where {\lambda} and D denote the wavelength and the telescope aperture diameter, respectively, to measure the quasi-static aberrations in the upstream pupil plane by encoding them into intensity variations in the downstream pupil image. We develop a rigorous formalism, leading to highly accurate measurement of the NCPA, in a quasi-linear way during the observation. Results. For a static phase map of standard deviation 44 nm rms at {\lambda} = 1.625 {\mu}m (0.026 {\lambda}), we estimate a possible reduction of the chromatic NCPA by a factor ranging from 3 to 10 in the presence of AO residuals compared with the expected performance of a typical current-generation system. This would allow a reduction of the level of quasi-static speckles in the detected images by a factor 10 to 100 hence, correspondingly improving the capacity to observe exoplanets.

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.

The precise reconstruction of the turbulent volume is a key point in the development of new-generation Adaptive Optics systems. We propose a new Cn2 profilometry method named CO-SLIDAR (COupled Slope and scIntillation Detection And Ranging), that uses correlations of slopes and scintillation indexes recorded on a Shack-Hartmann from two separated stars. CO-SLIDAR leads to an accurate Cn2 retrieval for both low and high altitude layers. Here, we present an end-to-end simulation of the Cn2 profile measurement. Two Shack-Hartmann geometries are considered. The detection noises are taken into account and a method to subtract the bias is proposed. Results are compared to Cn2 profiles obtained from correlations of slopes only or correlations of scintillation indexes only.

The EAGLE instrument is a Multi-Object Adaptive Optics (MOAO) fed, multiple Integral Field Spectrograph (IFS), working in the Near Infra-Red (NIR), on the European Extremely Large Telescope (E-ELT). A Phase A design study was delivered to the European Southern Observatory (ESO) leading to a successful review in October 2009. Since that time there have been a number of developments, which we summarize here. Some of these developments are also described in more detail in other submissions at this meeting. The science case for the instrument, while broad, highlighted in particular: understanding the stellar populations of galaxies in the nearby universe, the observation of the evolution of galaxies during the period of rapid stellar build-up between redshifts of 2-5, and the search for 'first light' in the universe at redshifts beyond 7. In the last 2 years substantial progress has been made in these areas, and we have updated our science case to show that EAGLE is still an essential facility for the E-ELT. This in turn allowed us to revisit the science requirements for the instrument, confirming most of the original decisions, but with one modification. The original location considered for the instrument (a gravity invariant focal station) is no longer in the E-ELT Construction Proposal, and so we have performed some preliminary analyses to show that the instrument can be simply adapted to work at the E-ELT Nasmyth platform. Since the delivery of the Phase A documentation, MOAO has been demonstrated on-sky by the CANARY experiment at the William Herschel Telescope.

EAGLE is the multi-object, spatially-resolved, near-IR spectrograph instrument concept for the E-ELT, relying on a distributed Adaptive Optics, so-called Multi Object Adaptive Optics. This paper presents the results of a phase A study. Using 84x84 actuator deformable mirrors, the performed analysis demonstrates that 6 laser guide stars and up to 5 natural guide stars of magnitude R<17, picked-up in a 7.3' diameter patrol field of view, allow us to obtain an overall performance in terms of Ensquared Energy of 35% in a 75x75 mas^2 spaxel at H band, whatever the target direction in the centred 5' science field for median seeing conditions. The computed sky coverage at galactic latitudes |b|~60 is close to 90%.

We acquired and analyzed adaptive-optics imaging observations of asteroid (2) Pallas from Keck II and the Very Large Telescope taken during four Pallas oppositions between 2003 and 2007, with spatial resolution spanning 32-88 km (image scales 13-20 km/pix). We improve our determination of the size, shape, and pole by a novel method that combines our AO data with 51 visual light-curves spanning 34 years of observations as well as occultation data. The shape model of Pallas derived here reproduces well both the projected shape of Pallas on the sky and light-curve behavior at all the epochs considered. We resolved the pole ambiguity and found the spin-vector coordinates to be within 5 deg. of [long, lat] = [30 deg., -16 deg.] in the ECJ2000.0 reference frame, indicating a high obliquity of ~84 deg., leading to high seasonal contrast. The best triaxial-ellipsoid fit returns radii of a=275 km, b= 258 km, and c= 238 km. From the mass of Pallas determined by gravitational perturbation on other minor bodies [(1.2 +/- 0.3) x 10-10 Solar Masses], we derive a density of 3.4 +/- 0.9 g.cm-3 significantly different from the density of C-type (1) Ceres of 2.2 +/- 0.1 g.cm-3. Considering the spectral similarities of Pallas and Ceres at visible and near-infrared wavelengths, this may point to fundamental differences in the interior composition or structure of these two bodies. We define a planetocentric longitude system for Pallas, following IAU guidelines. We also present the first albedo maps of Pallas covering ~80% of the surface in K-band. These maps reveal features with diameters in the 70-180 km range and an albedo contrast of about 6% wrt the mean surface albedo.

With the development of high contrast imaging techniques and infrared detectors, vast efforts have been devoted during the past decade to detect and characterize lighter, cooler and closer companions to nearby stars, and ultimately image new planetary systems. Complementary to other observing techniques (radial velocity, transit, micro-lensing, pulsar-timing), this approach has opened a new astrophysical window to study the physical properties and the formation mechanisms of brown dwarfs and planets. I here will briefly present the observing challenge, the different observing techniques, strategies and samples of current exoplanet imaging searches that have been selected in the context of the LAOG-Planet Imaging Surveys. I will finally describe the most recent results that led to the discovery of giant planets probably formed like the ones of our solar system, offering exciting and attractive perspectives for the future generation of deep imaging instruments.

One of the main science goal of the future European Extremely Large Telescope will be to understand the mass assembly process in galaxies as a function of cosmic time. To this aim, a multi-object, AO-assisted integral field spectrograph will be required to map the physical and chemical properties of very distant galaxies. In this paper, we examine the ability of such an instrument to obtain spatially resolved spectroscopy of a large sample of massive (0.1<Mstellar<5e11Mo) galaxies at 2<z<6, selected from future large area optical-near IR surveys. We produced a set of about one thousand numerical simulations of 3D observations using reasonable assumptions about the site, telescope, and instrument, and about the physics of distant galaxies. These data-cubes were analysed as real data to produce realistic kinematic measurements of very distant galaxies. We then studied how sensible the scientific goals are to the observational (i.e., site-, telescope-, and instrument-related) and physical (i.e., galaxy-related) parameters. We specifically investigated the impact of AO performance on the science goal. We did not identify any breaking points with respect to the parameters (e.g., the telescope diameter), with the exception of the telescope thermal background, which strongly limits the performance in the highest (z>5) redshift bin. We find that a survey of Ngal galaxies that fulfil the range of science goals can be achieved with a ~90 nights program on the E-ELT, provided a multiplex capability M Ngal/8.

One of the prominent science goal of the ELTs will be to study the physics and mass assembly of galaxies at very high redshifts. Here, we present the galaxy evolution science case for EAGLE, which is a NIR multi-integral field spectrograph for the E-ELT currently under phase A study. We summarize results of simulations conducted to derive high-level requirements. In particular, we show how we have derived the specifications for the ensquared energy that the AO system needs to provide to reach the scientific goals of the instrument. Finally, we present future strategies to conduct galaxy surveys with EAGLE.

We recently reported on the detection of a possible planetary-mass companion to Beta Pictoris at a projected separation of 8 AU from the star, using data taken in November 2003 with NaCo, the adaptive-optics system installed on the Very Large Telescope UT4. Eventhough no second epoch detection was available, there are strong arguments to favor a gravitationally bound companion rather than a background object. If confirmed and located at a physical separation of 8 AU, this young, hot (~1500 K), massive Jovian companion (~8 Mjup) would be the closest planet to its star ever imaged, could be formed via core-accretion, and could explain the main morphological and dynamical properties of the dust disk. Our goal was to return to Beta Pic five years later to obtain a second-epoch observation of the companion or, in case of a non-detection, constrain its orbit. Deep adaptive-optics L'-band direct images of Beta Pic and Ks-band Four-Quadrant-Phase-Mask (4QPM) coronagraphic images were recorded with NaCo in January and February 2009. We also use 4QPM data taken in November 2004. No point-like signal with the brightness of the companion candidate (apparent magnitudes L'=11.2 or Ks ~ 12.5) is detected at projected distances down to 6.5 AU from the star in the 2009 data. As expected, the non-detection does not allow to rule out a background object; however, we show that it is consistent with the orbital motion of a bound companion that got closer to the star since first observed in 2003 and that is just emerging from behind the star at the present epoch. We place strong constraints on the possible orbits of the companion and discuss future observing prospects.

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.

Several Wide Field of view Adaptive Optics (WFAO) concepts like Multi-Conjugate AO (MCAO), Multi-Object AO (MOAO) or Ground-Layer AO (GLAO) are currently studied for the next generation of Extremely Large Telescopes (ELTs). All these concepts will use atmospheric tomography to reconstruct the turbulent phase volume. In this paper, we explore different reconstruction algorithms and their fundamental limitations. We conduct this analysis in the Fourier domain. This approach allows us to derive simple analytical formulations for the different configurations, and brings a comprehensive view of WFAO limitations. We then investigate model and statistical errors and their impact on the phase reconstruction. Finally, we show some examples of different WFAO systems and their expected performance on a 42m telescope case.

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.

We compare coronagraph concepts and investigate their behavior and suitability for planet finder projects with Extremely Large Telescopes (ELTs, 30-42 meters class telescopes). For this task, we analyze the impact of major error sources that occur in a coronagraphic telescope (central obscuration, secondary support, low-order segment aberrations, segment reflectivity variations, pointing errors) for phase, amplitude and interferometric type coronagraphs. This analysis is performed at two different levels of the detection process: under residual phase left uncorrected by an eXtreme Adaptive Optics system (XAO) for a large range of Strehl ratio and after a general and simple model of speckle calibration, assuming common phase aberrations between the XAO and the coronagraph (static phase aberrations of the instrument) and non-common phase aberrations downstream of the coronagraph (differential aberrations provided by the calibration unit). We derive critical parameters that each concept will have to cope with by order of importance. We evidence three coronagraph categories as function of the accessible angular separation and proposed optimal one in each case. Most of the time amplitude concepts appear more favorable and specifically, the Apodized Pupil Lyot Coronagraph gathers the adequate characteristics to be a baseline design for ELTs.

[Abridged] We have developed an end-to-end simulation to specify the science requirements of a MOAO-fed integral field spectrograph on either an 8m or 42m telescope. Our simulations re-scales observations of local galaxies or results from numerical simulations of disk or interacting galaxies. For the current analysis, we limit ourselves to a local disk galaxy which exhibits simple rotation and a simulation of a merger. We have attempted to generalize our results by introducing the simple concepts of "PSF contrast" which is the amount of light polluting adjacent spectra which we find drives the smallest EE at a given spatial scale. The choice of the spatial sampling is driven by the "scale-coupling", i.e., the relationship between the IFU pixel scale and the size of the features that need to be recovered by 3D spectroscopy in order to understand the nature of the galaxy and its substructure. Because the dynamical nature of galaxies are mostly reflected in their large-scale motions, a relatively coarse spatial resolution is enough to distinguish between a rotating disk and a major merger. Although we used a limited number of morpho-kinematic cases, our simulations suggest that, on a 42m telescope, the choice of an IFU pixel scale of 50-75 mas seems to be sufficient. Such a coarse sampling has the benefit of lowering the exposure time to reach a specific signal-to-noise as well as relaxing the performance of the MOAO system. On the other hand, recovering the full 2D-kinematics of z~4 galaxies requires high signal-to-noise and at least an EE of 34% in 150 mas (2 pixels of 75 mas). Finally, we carried out a similar study at z=1.6 with a MOAO-fed spectrograph for an 8m, and find that at least an EE of 30% at 0.25 arcsec spatial sampling is required to understand the nature of disks and mergers.

SPHERE (Spectro Polarimetric High contrast Exoplanet REsearch), the planet finder instrument for the VLT is designed to study relatively bright extrasolar giant planets around young or nearby stars. SPHERE is a set of three instruments fed by the same AO-system, two of them share the same coronagraph. This complex system has been modeled with Fourier Optics to investigate the performance of the whole instrument. In turns, this end-to-end model was useful to analyze the sensitivity to various parameters (WFE, alignment of the coronagraph, differential aberrations) and to put some specifications on the sub-systems. This paper presents some example of sensitivity analysis and some contrast performance of the instruments as a function of the flux for the main observing mode of SPHERE: the Dual Band Imaging (DBI), equivalent to the Spectral Differential Imaging technique.

The Wide-field High-resolution Infrared TElescope (WHITE) will be dedicated in the first years of its life to carrying out a few (well focused in terms of science objectives and time) legacy surveys. WHITE would have an angular resolution of ~0.3'' homogeneous over ~0.7 sq. deg. in the wavelength range 1 - 5 um, which means that we will very efficiently use all the available observational time during night time and day time. Moreover, the deepest observations will be performed by summing up shorter individual frames. We will have a temporal information that can be used to study variable objects. The three key science objectives of WHITE are : 1) A complete survey of the Magellanic Clouds to make a complete census of young stellar objects in the clouds and in the bridge and to study their star formation history and the link with the Milky Way. The interaction of the two clouds with our Galaxy might the closest example of a minor merging event that could be the main driver of galaxy evolution in the last 5 Gyrs. 2) The building of the first sample of dusty supernovae at z<1.2 in the near infrared range (1-5 um) to constrain the equation of state from these obscured objects, study the formation of dust in galaxies and build the first high resolution sample of high redshift galaxies observed in their optical frame 3) A very wide weak lensing survey over that would allow to estimate the equation of state in a way that would favourably compete with space projects.

We propose a new concept of pupil motion sensor for astronomical adaptive optics systems and present experimental results obtained during the first laboratory validation of this concept. Pupil motion is an important issue in the case of extreme adaptive optics, high contrast systems, such as the proposed Planet Finder instruments for the ESO and Gemini 8-meter telescopes. Such high contrast imaging instruments will definitively require pupil stabilization to minimize the effect of quasi-static aberrations. The concept for pupil stabilization we propose uses the flux information from the AO system wave-front sensor to drive in closed loop a pupil tip-tilt mirror located in a focal plane. A laboratory experiment validates this concept and demonstrates its interest for high contrast imaging instrument.

We study the physical characteristics (shape, dimensions, spin axis direction, albedo maps, mineralogy) of the dwarf-planet Ceres based on high-angular resolution near-infrared observations. We analyze adaptive optics J/H/K imaging observations of Ceres performed at Keck II Observatory in September 2002 with an equivalent spatial resolution of ~50 km. The spectral behavior of the main geological features present on Ceres is compared with laboratory samples. Ceres' shape can be described by an oblate spheroid (a = b = 479.7 +/- 2.3 km, c = 444.4 +/- 2.1 km) with EQJ2000.0 spin vector coordinates RA = 288 +/- 5 deg. and DEC = +66 +/- 5 deg. Ceres sidereal period is measured to be 9.0741 +/- 0.0001 h. We image surface features with diameters in the 50-180 km range and an albedo contrast of ~6% with respect to the average Ceres albedo. The spectral behavior of the brightest regions on Ceres is consistent with phyllosilicates and carbonate compounds. Darker isolated regions could be related to the presence of frost.

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.

Surveys provide a wealth of data to the astronomical community that are used well after their completion. In this paper, we propose a project that would take the maximum benefit of Dome C in Antarctica by performing two surveys, in the wavelength range from 1-5 micron, complementary to SNAP space surveys. The first one over 1000 sq. deg. (1 KdF) for 4 years and the second one over 15 sq. deg (SNAP-IR) for the next 4 years at the same time as SNAP 0.35-1.7 microns survey. By using a Ground-Layer Adaptive Optics system, we would be able to recover, at the ice level and over at least half a degree in radius, the 300 mas angular resolution available above the 30-m high turbulent layer. Such a survey, combining a high angular resolution with high sensitivities in the NIR and MIR, should also play the role of a pre-survey for JWST and ALMA.

We report the discovery of companions to 5 nearby late M dwarfs (>M5), LHS1901, LHS4009, LHS6167, LP869-26 and WT460, and we confirm that the recently discovered mid-T brown dwarf companion to SCR1845-6357 is physically bound to that star. These discoveries result from our adaptive optics survey of all M dwarfs within 12 pc. The new companions have spectral types M5 to L1, and orbital separations between 1 and 10 AU. They add significantly to the number of late M dwarfs binaries in the immediate solar neighbourhood, and will improve the multiplicity statistics of late M dwarfs. The expected periods range from 3 to 130 years. Several pairs thus have good potential for accurate mass determination in this poorly sampled mass range.

Adaptive optics (AO) allows one to derive the point spread function (PSF) simultaneously to the science image, which is a major advantage in post-processing tasks such as astrometry/photometry or deconvolution. Based on the algorithm of \citet{veran97}, PSF reconstruction has been developed for four different AO systems so far: PUEO, ALFA, Lick-AO and Altair. A similar effort is undertaken for NAOS/VLT in a collaboration between the group PHASE (Onera and Observatoire de Paris/LESIA) and ESO. In this paper, we first introduce two new algorithms that prevent the use of the so-called "$U\_{ij}$ functions" to: (1) avoid the storage of a large amount of data (for both new algorithms), (2) shorten the PSF reconstruction computation time (for one of the two) and (3) provide an estimation of the PSF variability (for the other one). We then identify and explain issues in the exploitation of real-time Shack-Hartmann (SH) data for PSF reconstruction, emphasising the large impact of thresholding in the accuracy of the phase residual estimation. Finally, we present the data provided by the NAOS real-time computer (RTC) to reconstruct PSF ({\em (1)} the data presently available, {\em (2)} two NAOS software modifications that would provide new data to increase the accuracy of the PSF reconstruction and {\em (3)} the tests of these modifications) and the PSF reconstruction algorithms we are developing for NAOS on that basis.

Context. The knowledge of the point-spread function compensated by adaptive optics is of prime importance in several image restoration techniques such as deconvolution and astrometric/photometric algorithms. Wavefront-related data from the adaptive optics real-time computer can be used to accurately estimate the point-spread function in adaptive optics observations. The only point-spread function reconstruction algorithm implemented on astronomical adaptive optics system makes use of particular functions, named $U\_{ij}$. These $U\_{ij}$ functions are derived from the mirror modes, and their number is proportional to the square number of these mirror modes. Aims. We present here two new algorithms for point-spread function reconstruction that aim at suppressing the use of these $U\_{ij}$ functions to avoid the storage of a large amount of data and to shorten the computation time of this PSF reconstruction. Methods. Both algorithms take advantage of the eigen decomposition of the residual parallel phase covariance matrix. In the first algorithm, the use of a basis in which the latter matrix is diagonal reduces the number of $U\_{ij}$ functions to the number of mirror modes. In the second algorithm, this eigen decomposition is used to compute phase screens that follow the same statistics as the residual parallel phase covariance matrix, and thus suppress the need for these $U\_{ij}$ functions. Results. Our algorithms dramatically reduce the number of $U\_{ij}$ functions to be computed for the point-spread function reconstruction. Adaptive optics simulations show the good accuracy of both algorithms to reconstruct the point-spread function.

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.

In the last few years, new Adaptive Optics [AO] techniques have emerged to answer new astronomical challenges: Ground-Layer AO [GLAO] and Multi-Conjugate AO [MCAO] to access a wider Field of View [FoV], Multi-Object AO [MOAO] for the simultaneous observation of several faint galaxies, eXtreme AO [XAO] for the detection of faint companions. In this paper, we focus our study to one of these applications : high red-shift galaxy observations using MOAO techniques in the framework of Extremely Large Telescopes [ELTs]. We present the high-level specifications of a dedicated instrument. We choose to describe the scientific requirements with the following criteria : 40% of Ensquared Energy [EE] in H band (1.65um) and in an aperture size from 25 to 150 mas. Considering these specifications we investigate different AO solutions thanks to Fourier based simulations. Sky Coverage [SC] is computed for Natural and Laser Guide Stars [NGS, LGS] systems. We show that specifications are met for NGS-based systems at the cost of an extremely low SC. For the LGS approach, the option of low order correction with a faint NGS is discussed. We demonstrate that, this last solution allows the scientific requirements to be met together with a quasi full SC.

We present the results of a deep imaging survey of stars surrounded by planets detected with the radial velocity technique. The purpose is to search for and to characterize long-period stellar and substellar companions. The sample contains a total of 26 stars, among which 6 exhibit additional radial velocity drifts. We used NACO, at the ESO Very Large Telescope, and PUEO-KIR, at the Candian French Hawaiian Telescope, to conduct a near-infrared coronographic survey with adaptive optics of the faint circumstellar environment of the planetary hosts. The domain investigated ranges between 0.1 to 15" (i.e. about 3 to 500 AU, according to the mean distance of the sample). The survey is sensitive to companions within the stellar and the substellar domains, depending on the distance to the central stars and on the star properties. The images of 14 stars do not reveal any companions once the field objects are removed. 8 stars have close potential companions that need to be re-observed within 1-2 years to check for physical companionship. 4 stars are surrounded by faint objects which are confirmed or very probable companions. The companion to HD13445 (Gliese 86) is already known. The HD196885 star is a new close visual binary system with a high probability of being bound. The 2 newly discovered companions, HD1237 B and HD27442 B, share common proper motions with the central stars. Orbital motion is detected for HD1237 B. HD1237 B is likely a low-mass M star, located at 70 AU (projected distance) from the primary. HD27442 B is most probably a white dwarf companion located at about 240 AU (projected distance).

In this paper, we present a first comparison of different Adaptive Optics (AO) concepts to reach a given scientific specification for 3D spectroscopy on Extremely Large Telescope (ELT). We consider that a range of 30%-50% of Ensquarred Energy (EE) in H band (1.65um) and in an aperture size from 25 to 100mas is representative of the scientific requirements. From these preliminary choices, different kinds of AO concepts are investigated : Ground Layer Adaptive Optics (GLAO), Multi-Object AO (MOAO) and Laser Guide Stars AO (LGS). Using Fourier based simulations we study the performance of these AO systems depending on the telescope diameter.

In the context of extrasolar planet direct detection, we evaluated the performance of differential imaging with ground-based telescopes. This study was carried out in the framework of the VLT-Planet Finder project and is further extended to the case of Extremely Large Telescopes. Our analysis is providing critical specifications for future instruments mostly in terms of phase aberrations but also regarding alignments of the instrument optics or offset pointing on the coronagraph. It is found that Planet Finder projects on 8m class telescopes can be successful at detecting Extrasolar Giant Planets providing phase aberrations, alignments and pointing are accurately controlled. The situation is more pessimistic for the detection of terrestrial planets with Extremely Large Telescopes for which phase aberrations must be lowered at a very challenging level.

We analyse the fundamental limitations for the detection of extraterrestrial planets with Extremely Large Telescopes. For this task, a coronagraphic device combined to a very high order wavefront correction system is required but not sufficient to achieve the $10^{-10}$ contrast level needed for detecting an Earth-like planet. The stellar residuals left uncorrected by the wavefront correction system need to be calibrated and subtracted. In this paper, we consider a general model including the dynamic phase aberrations downstream the wavefront correction system, the static phase aberrations of the instrument and some differential aberrations provided by the calibration unit. A rather optimistic case of a filled circular pupil and of a perfect coronagraph is elsewhere assumed. As a result of the analytical study, the limitation mostly comes from the static aberrations. Using numerical simulations we confirm this result and evaluate the requirements in terms of phase aberrations to detect Earth-like planets on Extremely Large Telescopes.

We present the discovery of a tight M8V binary, with a separation of only 1.2 astronomical units, obtained with the PUEO and NACO adaptive optics systems, respectively at the CFHT and VLT telescopes. The estimated period of LP 349-25 is approximately 5 years, and this makes it an excellent candidate for a precise mass measurement.

We report VLT/NACO imaging observations of the stars GSC 08047-00232 and HIP 6856, probable members of the large Tucana-Horologium association. During our previous ADONIS/SHARPII deep imaging survey, a substellar candidate companion was discovered around each star. Based on VLT/NACO astrometric measurements, we find that GSC 08047-00232 and the faint candidate companion near to it share the same proper motion with a significance of 3.1sigma. On the contrary, the candidate companion to HIP 6856 is probably a background object with a significance of 4.3sigma. We also detect a new fainter and closer candidate companion to HIP 6856, but which is likely a background object too with a significance of 4sigma. Recent VLT/NACO spectroscopic measurements of GSC 08047-00232 finally confirm the substellar nature of this young brown dwarf with a derived spectral type M9.5. GSC 08047-00232 B, with an estimated mass of 25MJup and an effective temperature of 2100K, is presently the third substellar companion identified among young, nearby associations.

We report the detection of a possible subsolar mass companion to the triple young system TY CrA using the NACO instrument at the VLT UT4 during its commissioning. Assuming for TY CrA a distance similar to that of the close binary system HD 176386, the photometric spectral type of this fourth stellar component candidate is consistent with an ~M4 star. We discuss the dynamical stability of this possible quadruple system as well as the possible location of dusty particles inside or outside the system.

We present the results of an adaptive optics (AO) imaging survey of the common associations of Tucana and Horologium, carried out at the ESO 3.6m telescope with the ADONIS/SHARPII system. Based on our observations of two dozen probable association members, HIP 1910 and HIP 108422 appear to have low-mass stellar companions, while HIP 6856 and GSC 8047-0232 have possible sub-stellar candidate companions. Astrometric measurements, performed in November 2000 and October 2001, indicate that HIP 1910 B likely is bound to its primary, but are inconclusive in the case of the candidate companion to HIP 6856. Additional observations are needed to confirm the HIP 6856 companionship as well as for HIP 108422 and GSC 8047-0232.

We report adaptive optics (AO) observations of the young and nearby association MBM 12 obtained with the Canada-France-Hawaii Telescope. Our main observational result is the discovery of six new binary systems, LkHa 264, E 0255+2018, RX J0255.4+2005, S18, MBM 12-10, RX J0255.3+1915, and the confirmation of HD 17332, already known as a binary. We also detected a possible quadruple system. It is composed of the close binary LkHa 263 AB (separation of 0.41 ''), of LkH\alpha 262 located 15.25 '' from LkHa 263 A, and of LkHa 263 C, located 4.1 '' from LkH\alpha 263 A. A preliminary study of the binary fraction suggests a binary excess in the MBM 12 association as compared to the field and IC 348. Because of the high binarity rate, previous estimations of spectral types and measurements of IR excesses for several candidate members of MBM 12 have to be revised. LkH\alpha 263 C is a nebulous object that we interpret as a disk oriented almost perfectly edge-on and seen in scattered light. This object has already been reported by Jayawardhana et al. (2002). Scattered light models allow us to estimate some of the structural parameters (i.e. inclination, diameter and to a lesser extent dust mass) of the circumstellar disk. We find an inclination of 89^o and a outer radius for the disk, 165 AU if the distance to MBM 12 is 275 pc. With the present data set, we do not attempt to re-assess the distance to MBM 12. We estimate however that the distance to the candidate member RX J0255.3+1915 is d > 175 pc.

We present new, high angular resolution images of HD 199143 in the Capricornus association, obtained with the adaptive optics system ADONIS+SHARPII at the ESO 3.6m Telescope of La Silla Observatory. HD 199143 and its neighbour star HD 358623 (separation 5'away) have previously been imaged with adaptive optics. For each star, a companion has been detected in the J and K bands at respective separations of 1.1'' and 2.2'' (Jayawardhana & Brandeker 2001). Our new photometry of HD 199143 B suggests that it is a M2 star and that the presence of circumstellar dust proposed by Van den Ancker et al. (2000) is no longer necessary. We show that the 12 microns flux detected by IRAS previously interpreted as an IR excess, can be explained by the presence of the late-type companion.

A large fraction of the present-day stellar mass was formed between z=0.5 and z~3 and our understanding of the formation mechanisms at work at these epochs requires both high spatial and high spectral resolution: one shall simultaneously} obtain images of objects with typical sizes as small as 1-2kpc(~0''.1), while achieving 20-50 km/s (R >= 5000) spectral resolution. The obvious instrumental solution to adopt in order to tackle the science goal is therefore a combination of multi-object 3D spectrograph with multi-conjugate adaptive optics in large fields. A partial, but still competitive correction shall be prefered, over a much wider field of view. This can be done by estimating the turbulent volume from sets of natural guide stars, by optimizing the correction to several and discrete small areas of few arcsec2 selected in a large field (Nasmyth field of 25 arcmin) and by correcting up to the 6th, and eventually, up to the 60th Zernike modes. Simulations on real extragalactic fields, show that for most sources (>80%), the recovered resolution could reach 0".15-0".25 in the J and H bands. Detection of point-like objects is improved by factors from 3 to >10, when compared with an instrument without adaptive correction. The proposed instrument concept, FALCON, is equiped with deployable mini-integral field units (IFUs), achieving spectral resolutions between R=5000 and 20000. Its multiplex capability, combined with high spatial and spectral resolution characteristics, is a natural ground based complement to the next generation of space telescopes.