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

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

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 nearby Herbig Be star HD100546 is known to be a laboratory for the study of protoplanets and their relation with the circumstellar disk that is carved by at least 2 gaps. We observed the HD100546 environment with high contrast imaging exploiting several different observing modes of SPHERE, including datasets with/without coronagraphs, dual band imaging, integral field spectroscopy and polarimetry. The picture emerging from these different data sets is complex. Flux-conservative algorithms images clearly show the disk up to 200au. More aggressive algorithms reveal several rings and warped arms overlapping the main disk. The bright parts of this ring lie at considerable height over the disk mid-plane at about 30au. Our images demonstrate that the brightest wings close to the star in the near side of the disk are a unique structure, corresponding to the outer edge of the intermediate disk at ~40au. Modeling of the scattered light from the disk with a geometrical algorithm reveals that a moderately thin structure can well reproduce the light distribution in the flux-conservative images. We suggest that the gap between 44 and 113 au span between the 1:2 and 3:2 resonance orbits of a massive body located at ~70au that might coincide with the candidate planet HD100546b detected with previous thermal IR observations. In this picture, the two wings can be the near side of a ring formed by disk material brought out of the disk at the 1:2 resonance with the same massive object. While we find no clear evidence confirming detection of the planet candidate HD100546c in our data, we find a diffuse emission close to the expected position of HD100546b. This source can be described as an extremely reddened substellar object surrounded by a dust cloud or its circumplanetary disk. Its astrometry is broadly consistent with a circular orbital motion on the disk plane.

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

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

We present observations with the planet finder SPHERE of a selected sample of the most promising radial velocity (RV) companions for high-contrast imaging. Using a Monte Carlo simulation to explore all the possible inclinations of the orbit of wide RV companions, we identified the systems with companions that could potentially be detected with SPHERE. We found the most favorable RV systems to observe are : HD\,142, GJ\,676, HD\,39091, HIP\,70849, and HD\,30177 and carried out observations of these systems during SPHERE Guaranteed Time Observing (GTO). To reduce the intensity of the starlight and reveal faint companions, we used Principle Component Analysis (PCA) algorithms alongside angular and spectral differential imaging. We injected synthetic planets with known flux to evaluate the self-subtraction caused by our data reduction and to determine the 5$\sigma$ contrast in the J band $vs$ separation for our reduced images. We estimated the upper limit on detectable companion mass around the selected stars from the contrast plot obtained from our data reduction. Although our observations enabled contrasts larger than 15 mag at a few tenths of arcsec from the host stars, we detected no planets. However, we were able to set upper mass limits around the stars using AMES-COND evolutionary models. We can exclude the presence of companions more massive than 25-28 \MJup around these stars, confirming the substellar nature of these RV companions.

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

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

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

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.

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

The transition disk around the T Tauri star T Cha possesses a large gap, making it a prime target for high-resolution imaging in the context of planet formation. We aim to find signs of disk evolutionary processes by studying the disk geometry and the dust grain properties at its surface, and to search for companion candidates. We analyze a set of VLT/SPHERE data at near-infrared and optical wavelengths. We performed polarimetric imaging of T Cha with IRDIS (1.6 $\mu$m) and ZIMPOL (0.5-0.9 $\mu$m), and obtained intensity images from IRDIS dual-band imaging with simultaneous spectro-imaging with IFS (0.9-1.3 $\mu$m). The disk around T Cha is detected in all observing modes and its outer disk is resolved in scattered light with unprecedented angular resolution and signal-to-noise. The images reveal a highly inclined disk with a noticeable east-west brightness asymmetry. The significant amount of non-azimuthal polarization signal in the $U_{\phi}$ images, with a $U_{\phi}$/$Q_{\phi}$ peak-to-peak value of 14%, is in accordance with theoretical studies on multiple scattering. Our optimal axisymmetric radiative transfer model considers two coplanar inner and outer disks, separated by a gap of 0.28" (~30au) in size. We derive a disk inclination of ~69 deg and PA of ~114 deg. In order to self-consistently reproduce the intensity and polarimetric images, the dust grains, responsible for the scattered light, need to be dominated by sizes of around ten microns. A point source is detected at an angular distance of 3.5" from the central star. It is, however, found not to be co-moving. We confirm that the dominant source of emission is forward scattered light from the near edge of the outer disk. Our point source analysis rules out the presence of a companion heavier than ~8.5 $M_{\mathrm{jup}}$ between 0.1" and 0.3". The detection limit decreases to ~2 $M_{\mathrm{jup}}$ for 0.3" to 4.0".

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

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

We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets. We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the Pynpoint package. The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate 'c' might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of 'b'. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to 'b'.