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.

We present the discovery of a white dwarf companion at 3.6" from GJ3346, a nearby ($\pi\sim$42 mas) K star observed with SPHERE@VLT as part of an open time survey for faint companions to objects with significant proper motion discrepancies ($\Delta\mu$) between Gaia DR1 and Tycho-2. Syrius-like systems like GJ3346AB, which include a main sequence star and a white dwarf, can be difficult to detect because of the intrinsic faintness of the latter. They have, however, been found to be common contaminants for direct imaging searches. White dwarfs have in fact similar brightness to sub-stellar companions in the infrared, while being much brighter in the visible bands like those used by Gaia. Combining our observations with Gaia DR2 and with several additional archival data sets, we were able to fully constrain the physical properties of GJ3346B, such as its effective temperature (11$\times$10$^3\pm$500 K) as well as the cooling age of the system (648$\pm$58 Myrs). This allowed us to better understand the system history and to partially explains the discrepancies previously noted in the age indicators for this objects. Although further investigation is still needed, it seems that GJ3346, which was previously classified as young, is in fact most likely to be older than 4 Gyrs. Finally, given that the mass (0.58$\pm$0.01$M_{\odot}$)} and separation (85 au) of GJ3346B are compatible with the observed $\Delta\mu$, this discovery represents a further confirmation of the potential of this kind of dynamical signatures as selection methods for direct imaging surveys targeting faint, sub-stellar companions.

In the second paper of this series we perfected our method of linking high precision Hubble Space Telescope astrometry to the high-accuracy Gaia DR2 absolute reference system to overcome the limitations of relative astrometry with narrow-field cameras. Our test case here is the Y brown dwarf WISE J163940.83-684738.6, observed at different epochs spread over a 6-yr time baseline with the Infra-Red channel of the Wide Field Camera 3. We derived significantly improved astrometric parameters compared to previous determinations, finding: (mu_RAcosDc,mu_DC,parallax) = (577.21+/-0.24mas/yr,-3108.39+/-0.27mas/yr,210.4+/-1.8mas). In particular, our derived absolute parallax corresponds to a distance of 4.75+/-0.05pc for the faint ultracool dwarf.

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.

Stellar multiplicity is believed to influence planetary formation and evolution, although the precise nature and extent of this role remain ambiguous. We present a study aimed at testing the role of stellar multiplicity in the formation and/or evolution of the most massive, close-in planetary and substellar companions. Using direct imaging observations, as well as the Gaia DR2 catalogue, we searched for wide binary companions to 38 stars hosting massive giant planets or brown dwarfs (M > 7 MJup) on orbits shorter than ~1 AU. We report the discovery of a new component in the WASP-14 system, and present an independent confirmation of a comoving companion to WASP-18. From a robust Bayesian statistical analysis, we derived a binary fraction of 79.0+13.2-14.7% between 20-10,000 AU for our sample, twice as high as for field stars with a 3-{\sigma} significance. This binary frequency was found to be larger than for lower-mass planets on similar orbits, and we observed a marginally higher binary rate for inner companions with periods shorter than 10 days. These results demonstrate that stellar companions greatly influence the formation and/or evolution of these systems, suggesting that the role played by binary companions becomes more important for higher-mass planets, and that this trend may be enhanced for systems with tighter orbits. Our analysis also revealed a peak in binary separation at 250 AU, highlighting a shortfall of close binaries among our sample. This indicates that the mechanisms affecting planet and brown dwarf formation or evolution in binaries must operate from wide separations, although we found that the Kozai-Lidov mechanism is unlikely to be the dominant underlying process. We conclude that binarity plays a crucial role in the existence of very massive short-period giant planets and brown dwarf desert inhabitants, which are almost exclusively observed in multiple systems.

We detect the accreting low-mass companion HD142527B at a separation of 73 mas (11.4 au) from the star. No other companions with mass greater than 10 MJ are visible in the field of view of IFS (\sim 100 au centered on the star) or in the IRDIS field of view (\sim 400 au centered on the star). Measurements from IFS, SAM IFS, and IRDIS suggest an M6 spectral type for HD142527B, with an uncertainty of one spectral subtype, compatible with an object of M=0.11 \pm 0.06 MSun and R=0.15 \pm 0.07 RSun. The determination of the mass remains a challenge using contemporary evolutionary models, as they do not account for the energy input due to accretion from infalling material. We consider that the spectral type of the secondary may also be earlier than the type we derived from IFS spectra. From dynamical considerations, we further constrain the mass to 0.26^{+0.16}_{-0.14} MSun , which is consistent with both our spectroscopic analysis and the values reported in the literature. Following previous methods, the lower and upper dynamical mass values correspond to a spectral type between M2.5 and M5.5 for the companion. By fitting the astrometric points, we find the following orbital parameters: a period of P=35-137 yr; an inclination of i=121-130 deg.; , a value of Omega=124-135 deg for the longitude of node, and an 68% confidence interval of \sim 18 - 57 au for the separation at periapsis. Eccentricity and time at periapsis passage exhibit two groups of values: \sim0.2-0.45 and \sim0.45-0.7 for e, and \sim 2015-2020 and \sim2020-2022 for T_0. While these orbital parameters might at first suggest that HD142527B is not the companion responsible for the outer disk truncation, a previous hydrodynamical analysis of this system showed that they are compatible with a companion that is able to produce the large cavity and other observed features.

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.

One field containing WISE J154151.65-225024.9 was observed by Hubble Space Telescope at three different epochs taken in ~5 yrs. We measured positions of sources in all images and successfully linked these positions to the Gaia DR2 absolute system to derive the astrometric parameters for this faint close-by Y1 brown dwarf. The developed procedure avoids traditional limitations of relative imaging-astrometry with narrow-field cameras, extending Gaia DR2 to fainter magnitudes. We found (mu_RA,mu_d,parallax) = (-902.62+/-0.35mas/yr, -88.26+/-0.35mas/yr,168.4+/-2.2mas), which represent a sizable improvement over recent determinations in the literature. Applying a correction from relative to absolute parallax we found an absolute parallax of 169+/-2 mas, corresponding to a distance of 5.9+/-0.1 pc.

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.

Binary statistics of the latest-type T and Y brown dwarfs are sparse and it is unclear whether the trends seen in the multiplicity properties of their more massive counterparts hold for the very coolest brown dwarfs. We present results from a search for substellar and planetary-mass companions to a sample of 12 ultracool T8$-$Y0 field brown dwarfs with the Hubble Space Telescope/Wide Field Camera 3. We find no evidence for resolved binary companions among our sample down to separations of 0.7$-$2.5 AU. Combining our survey with prior searches, we place some of the first statistically robust constraints to date on the multiplicity properties of the coolest, lowest-mass brown dwarfs in the field. Accounting for observational biases and incompleteness, we derive a binary frequency of $f = 5.5^{+5.2}_{-3.3}$% for T5$-$Y0 brown dwarfs at separations of 1.5$-$1000 AU, for an overall binary fraction of $f_\mathrm{tot} = 8\pm6$%. Modelling the projected separation as a lognormal distribution, we find a peak in separation at $\rho_0 = 2.9^{+0.8}_{-1.4}$ AU with a logarithmic width of $\sigma = 0.21^{+0.14}_{-0.08}$. We infer a mass ratio distribution peaking strongly towards unity, with a power law index of $\gamma = 6.1^{+4.0}_{-2.7}$, reinforcing the significance of the detection of a tighter and higher mass ratio companion population around lower-mass primaries. These results are consistent with prior studies and support the idea of a decreasing binary frequency with spectral type in the Galactic field.

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.

In this paper we present the results of the SPHERE observation of the HD 284149 system, aimed at a more detailed characterisation of both the primary and its brown dwarf companion. We observed HD 284149 in the near-infrared with SPHERE, using the imaging mode (IRDIS+IFS) and the long-slit spectroscopy mode (IRDIS-LSS). The data were reduced using the dedicated SPHERE pipeline, and algorithms such as PCA and TLOCI were applied to reduce the speckle pattern. The IFS images revealed a previously unknown low-mass (~0.16$M_{\odot}$) stellar companion (HD 294149 B) at ~0.1$^{\prime\prime}$, compatible with previously observed radial velocity differences, as well as proper motion differences between Gaia and Tycho-2 measurements. The known brown dwarf companion (HD 284149 b) is clearly visible in the IRDIS images. This allowed us to refine both its photometry and astrometry. The analysis of the medium resolution IRDIS long slit spectra also allowed a refinement of temperature and spectral type estimates. A full reassessment of the age and distance of the system was also performed, leading to more precise values of both mass and semi-major axis. As a result of this study, HD 284149 ABb therefore becomes the latest addition to the (short) list of brown dwarfs on wide circumbinary orbits, providing new evidence to support recent claims that object in such configuration occur with a similar frequency to wide companions to single stars.

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.