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.

Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere.

We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the $JHK$ bands, possibly correlated with the presence of the shadows. Analysis of archival VLTI/PIONIER $H$-band visibilities constrain the orientation of the inner disk to $i = 18.2\deg^{+3.4}_{-4.1}$ and ${\rm PA} = 57.3\deg \pm 5.7\deg$, consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere ($H/r \lesssim 0.2$). Possible mechanisms that may cause asymmetric variations in the optical depth ($\Delta\tau\lesssim1$) through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon.

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 discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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.

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the S-type AGB star chi Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of chi Cyg. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been performed to derive the temperature, density, and velocity structure of the circumstellar envelope. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high-J CO lines to derive the parameters for the circumstellar envelope, we modelled both the ortho- and para-H2O lines. Our modelling results are consistent with the velocity structure expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 +/- 0.2)E-5, much lower than that in O-rich stars. The derived ortho-to-para ratio of 2.1 +/- 0.6 is close to the high-temperature equilibrium limit, consistent with H2O being formed in the photosphere.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

[abridged] Our knowledge of the inner structure of embedded massive young stellar objects is still quite limited. We attempt here to overcome the spatial resolution limitations of conventional thermal infrared imaging. We employed mid-infrared interferometry using the MIDI instrument on the ESO/VLTI facility to investigate M8E-IR, a well-known massive young stellar object suspected of containing a circumstellar disk. Spectrally dispersed visibilities in the 8-13 micron range were obtained at seven interferometric baselines. We resolve the mid-infrared emission of M8E-IR and find typical sizes of the emission regions of the order of 30 milli-arcseconds (~45 AU). Radiative transfer simulations have been performed to interpret the data. The fitting of the spectral energy distribution, in combination with the measured visibilities, does not provide evidence for an extended circumstellar disk with sizes > 100 AU but requires the presence of an extended envelope. The data are not able to constrain the presence of a small-scale disk in addition to an envelope. In either case, the interferometry measurements indicate the existence of a strongly bloated, relatively cool central object, possibly tracing the recent accretion history of M8E-IR. In addition, we present 24.5 micron images that clearly distinguish between M8E-IR and the neighbouring ultracompact HII region and which show the cometary-shaped infrared morphology of the latter source. Our results show that IR interferometry, combined with radiative transfer modelling, can be a viable tool to reveal crucial structure information on embedded massive young stellar objects and to resolve ambiguities arising from fitting the SED.

Based on the SED, Herbig stars have been categorized into two observational groups, reflecting their overall disk structure: group I members have disks with a higher degree of flaring than their group II counterparts. We investigate the 5-35 um Spitzer IRS spectra of a sample of 13 group I sources and 20 group II sources. We focus on the continuum emission to study the underlying disk geometry. We have determined the [30/13.5] and [13.5/7] continuum flux ratios. The 7-um flux excess with respect to the stellar photosphere is measured, as a marker for the strength of the near-IR emission produced by the inner disk. We have compared our data to self-consistent passive-disk model spectra, for which the same quantities were derived. We confirm the literature result that the difference in continuum emission between group I and II sources can largely be explained by a different amount of small dust grains. However, we report a strong correlation between the [30/13.5] and [13.5/7] flux ratios for Meeus group II sources. Moreover, the [30/13.5] flux ratio decreases with increasing 7-um excess for all targets in the sample. To explain these correlations with the models, we need to introduce an artificial scaling factor for the inner disk height. In roughly 50% of the Herbig Ae/Be stars in our sample, the inner disk must be inflated by a factor 2 to 3 beyond what hydrostatic calculations predict. The total disk mass in small dust grains determines the degree of flaring. We conclude, however, that for any given disk mass in small dust grains, the shadowing of the outer (tens of AU) disk is determined by the scale height of the inner disk (1 AU). The inner disk partially obscures the outer disk, reducing the disk surface temperature. Here, for the first time, we prove these effects observationally.

We have analysed near-infrared NTT/SofI observations of the starburst cluster Westerlund 1, which is among the most massive young clusters in the Milky Way. A comparison of colour-magnitude diagrams with theoretical main-sequence and pre-main sequence evolutionary tracks yields improved extinction and distance estimates of A_Ks = 1.13+-0.03 mag and d = 3.55+-0.17 kpc (DM = 12.75+-0.10 mag). The pre-main sequence population is best fit by a Palla & Stahler isochrone for an age of 3.2 Myr, while the main sequence population is in agreement with a cluster age of 3 to 5 Myr. An analysis of the structural parameters of the cluster yields that the half-mass radius of the cluster population increases towards lower mass, indicative of the presence of mass segregation. The cluster is clearly elongated with an eccentricity of 0.20 for stars with masses between 10 and 32 Msun, and 0.15 for stars with masses in the range 3 to 10 Msun. We derive the slope of the stellar mass function for stars with masses between 3.4 and 27 Msun. In an annulus with radii between 0.75 and 1.5 pc from the cluster centre, we obtain a slope of Gamma = -1.3. Closer in, the mass function of Westerlund 1 is shallower with Gamma = -0.6. The extrapolation of the mass function for stars with masses from 0.08 to 120 Msun yields an initial total stellar mass of ~52,000 Msun, and a present-day mass of 20,000 to 45,000 Msun (about 10 times the stellar mass of the Orion Nebula Cluster, and 2 to 4 times the mass of the NGC 3603 young cluster), indicating that Westerlund 1 is the most massive starburst cluster identified to date in the Milky Way.

We report on observations of circumstellar disks around young stars that have been obtained with the MIDI instrument, which is mounted on the VLT Interferometer and operates in the 10 micrometer atmospheric window. The maximum spatial resolution of 5 milli-arcsec corresponds to sub-AU scales at the distance to nearby star formation regions. Thus, we can study the disks on the spatial scales at which important processes occur, such as accretion, dust processing, and planet formation. The main results obtained so far can be summarized as follows: 1. The measured interferometric visibilities are in good qualitative agreement with those predicted by models of circumstellar disks. In particular, a predicted correlation between the strength of the far-infrared excess and the spatial structure of the disk is confirmed by direct measurements; 2. In several objects strong evidence for deviations from circular symmetry is present, indicating that an inclined disk is indeed the dominant component seen in the mid-infrared; 3. The dust properties are not uniform over the disk, but are instead a strong function of distance to the central star. The dust in the innermost disk regions is observed to be more ``processed'' than the dust further out, both in Herbig Ae star disks and in those around T-Tauri stars.

Observations of extrasolar planets using Integral Field Spectroscopy (IFS), if coupled with an extreme Adaptive Optics system and analyzed with a Simultaneous Differential Imaging technique (SDI), are a powerful tool to detect and characterize extrasolar planets directly; they enhance the signal of the planet and, at the same time, reduces the impact of stellar light and consequently important noise sources like speckles. In order to verify the efficiency of such a technique, we developed a simulation code able to test the capabilities of this IFS-SDI technique for different kinds of planets and telescopes, modelling the atmospheric and instrumental noise sources. The first results obtained by the simulations show that many significant extrasolar planet detections are indeed possible using the present 8m-class telescopes within a few hours of exposure time. The procedure adopted to simulate IFS observations is presented here in detail, explaining in particular how we obtain estimates of the speckle noise, Adaptive Optics corrections, specific instrumental features, and how we test the efficiency of the SDI technique to increase the signal-to-noise ratio of the planet detection. The most important results achieved by simulations of various objects, from 1 M_J to brown dwarfs of 30 M_J, for observations with an 8 meter telescope, are then presented and discussed.

We review recent advances in our understanding of the innermost regions of the circumstellar environment around young stars, made possible by the technique of long baseline interferometry at infrared wavelengths. Near-infrared observations directly probe the location of the hottest dust. The characteristic sizes found are much larger than previously thought, and strongly correlate with the luminosity of the central young stars. This relation has motivated in part a new class of models of the inner disk structure. The first mid-infrared observations have probed disk emission over a larger range of scales, and spectrally resolved interferometry has for the first time revealed mineralogy gradients in the disk. These new measurements provide crucial information on the structure and physical properties of young circumstellar disks, as initial conditions for planet formation.

We review the properties of dust in protoplanetary disks around optically visible pre-main sequence stars obtained with a variety of observational techniques, from measurements of scattered light at visual and infrared wavelengths to mid-infrared spectroscopy and millimeter interferometry. A general result is that grains in disks are on average much larger than in the diffuse interstellar medium (ISM). In many disks, there is evidence that a large mass of dust is in grains with millimeter and centimeter sizes, more similar to "sand and pebbles" than to grains. Smaller grains (with micron-sizes) exist closer to the disk surface, which also contains much smaller particles, e.g., polycyclic aromatic hydrocarbons. There is some evidence of a vertical stratification, with smaller grains closer to the surface. Another difference with ISM is the higher fraction of crystalline relative to amorphous silicates found in disk surfaces. There is a large scatter in dust properties among different sources, but no evidence of correlation with the stellar properties, for samples that include objects from intermediate to solar mass stars and brown dwarfs. There is also no apparent correlation with the age of the central object, over a range roughly between 1 and 10 Myr. This suggests a scenario where significant grain processing may occur very early in the disk evolution, possibly when it is accreting matter from the parental molecular core. Further evolution may occur, but not necessarily rapidly, since we have evidence that large amounts of grains, from micron to centimeter size, can survive for periods as long as 10 Myr.

High-quality K-band spectra of point sources, deeply embedded in massive star-forming regions, have revealed a population of 20 young massive stars showing no photospheric absorption lines, but only emission lines. The K-band spectra exhibit one or more features commonly associated with massive Young Stellar Objects surrounded by circumstellar material: a very red color (J-K) = 2, CO bandhead emission, hydrogen emission lines (sometimes doubly peaked), and FeII and/or MgII emission lines. The CO emission comes from a relatively dense (~10^10 cm^(-3)) and hot (T ~ 2000-5000 K) region, sufficiently shielded from the intense UV radiation field of the young massive star. Modeling of the CO-first overtone emission shows that the CO gas is located within 5 AU of the star. The hydrogen emission is produced in an ionized medium exposed to UV radiation. The best geometrical configuration is a dense and neutral circumstellar disk causing the CO bandhead emission, and an ionized upper layer where the hydrogen lines are produced. We argue that the circumstellar disk is likely a remnant of the accretion via a circumstellar disk.

We report the discovery of a symmetric bow shock around the well-known high-mass X-ray binary (HMXB) Vela X-1. Wind bow shocks are a ubiquitous phenomenon around OB-runaway stars, but now such a structure is found around a HMXB. The presence of a bow shock indicates that the system has a high (supersonic) velocity with respect to the interstellar medium. From the symmetry of the bow shock, the direction of motion and, moreover, the origin and age of the system can be derived. Our observation supports Blaauw's scenario for the formation of an OB-runaway star by the supernova explosion of the binary companion.