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.

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.

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.

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

We present observations of the Herbig Ae star HD169142 with VLT/SPHERE instruments IRDIS (K1K2 and H2H3 bands) and IFS (Y-, J- and H-bands). We detect several bright blobs at ~180 mas separation from the star, and a faint arc-like structure in the IFS data. Our Radial Differential Imaging data analysis also displays 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 H- and K$_s$-bands at these separations in Biller et al. (2014) 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 ~93 mas separation and PA=355{\deg}, 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 ~100 mas. Follow-up observations are necessary to confirm the detection and the nature of this source and structure.

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.

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.

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.

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.

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

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

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

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